201
|
Lu XR, Duan AQ, Li WQ, Abdel-Shafy H, Rushdi HE, Liang SS, Ma XY, Liang XW, Deng TX. Genome-wide analysis reveals genetic diversity, linkage disequilibrium, and selection for milk production traits in Chinese buffalo breeds. J Dairy Sci 2020; 103:4545-4556. [PMID: 32147265 DOI: 10.3168/jds.2019-17364] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2019] [Accepted: 01/13/2020] [Indexed: 11/19/2022]
Abstract
The water buffalo is an important dual-purpose livestock that is widespread throughout central and southern China. However, there has been no characterization of the population genetics of Chinese buffalo. Using an Axiom buffalo genotyping array (Thermo Fisher Scientific, Wilmington, DE), we analyzed the genetic diversity, linkage disequilibrium pattern, and signature of selection in 176 Chinese buffaloes from 13 breeds. A total of 35,547 SNP passed quality control and were used for further analyses. Population genetic analysis revealed a clear separation between swamp and river types. Ten Chinese indigenous breeds were clustered into the swamp group, the Murrah and Nili-Ravi breeds were clustered into the river group, and the crossbred breed was closer to the river group. Genetic diversity analysis showed that the swamp group had a lower average expected heterozygosity. Linkage disequilibrium decay distance was much shorter in the swamp group compared with the river group, with an average square of correlation coefficient value of 0.2 of approximately 50 kb. Analysis of runs of homozygosity indicated extensive remote and recent inbreeding within swamp and river groups, respectively. Moreover, one genomic region under selection was detected between the river and swamp groups. Our findings contribute to our understanding of the characterization of population genetics in Chinese buffaloes, which in turn may be used in buffalo breeding programs.
Collapse
Affiliation(s)
- X R Lu
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction Technology, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning 530001, China
| | - A Q Duan
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction Technology, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning 530001, China
| | - W Q Li
- BGI Genomics, BGI-Shenzhen, Shenzhen 518083, China
| | - H Abdel-Shafy
- Department of Animal Production, Faculty of Agriculture, Cairo University, 12613 Giza, Egypt
| | - H E Rushdi
- Department of Animal Production, Faculty of Agriculture, Cairo University, 12613 Giza, Egypt
| | - S S Liang
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction Technology, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning 530001, China
| | - X Y Ma
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction Technology, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning 530001, China
| | - X W Liang
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction Technology, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning 530001, China
| | - T X Deng
- Key Laboratory of Buffalo Genetics, Breeding and Reproduction Technology, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning 530001, China.
| |
Collapse
|
202
|
Ghoreishifar SM, Moradi-Shahrbabak H, Fallahi MH, Jalil Sarghale A, Moradi-Shahrbabak M, Abdollahi-Arpanahi R, Khansefid M. Genomic measures of inbreeding coefficients and genome-wide scan for runs of homozygosity islands in Iranian river buffalo, Bubalus bubalis. BMC Genet 2020; 21:16. [PMID: 32041535 PMCID: PMC7011551 DOI: 10.1186/s12863-020-0824-y] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2019] [Accepted: 02/04/2020] [Indexed: 01/06/2023] Open
Abstract
Background Consecutive homozygous fragments of a genome inherited by offspring from a common ancestor are known as runs of homozygosity (ROH). ROH can be used to calculate genomic inbreeding and to identify genomic regions that are potentially under historical selection pressure. The dataset of our study consisted of 254 Azeri (AZ) and 115 Khuzestani (KHZ) river buffalo genotyped for ~ 65,000 SNPs for the following two purposes: 1) to estimate and compare inbreeding calculated using ROH (FROH), excess of homozygosity (FHOM), correlation between uniting gametes (FUNI), and diagonal elements of the genomic relationship matrix (FGRM); 2) to identify frequently occurring ROH (i.e. ROH islands) for our selection signature and gene enrichment studies. Results In this study, 9102 ROH were identified, with an average number of 21.2 ± 13.1 and 33.2 ± 15.9 segments per animal in AZ and KHZ breeds, respectively. On average in AZ, 4.35% (108.8 ± 120.3 Mb), and in KHZ, 5.96% (149.1 ± 107.7 Mb) of the genome was autozygous. The estimated inbreeding values based on FHOM, FUNI and FGRM were higher in AZ than they were in KHZ, which was in contrast to the FROH estimates. We identified 11 ROH islands (four in AZ and seven in KHZ). In the KHZ breed, the genes located in ROH islands were enriched for multiple Gene Ontology (GO) terms (P ≤ 0.05). The genes located in ROH islands were associated with diverse biological functions and traits such as body size and muscle development (BMP2), immune response (CYP27B1), milk production and components (MARS, ADRA1A, and KCTD16), coat colour and pigmentation (PMEL and MYO1A), reproductive traits (INHBC, INHBE, STAT6 and PCNA), and bone development (SUOX). Conclusion The calculated FROH was in line with expected higher inbreeding in KHZ than in AZ because of the smaller effective population size of KHZ. Thus, we find that FROH can be used as a robust estimate of genomic inbreeding. Further, the majority of ROH peaks were overlapped with or in close proximity to the previously reported genomic regions with signatures of selection. This tells us that it is likely that the genes in the ROH islands have been subject to artificial or natural selection.
Collapse
Affiliation(s)
- Seyed Mohammad Ghoreishifar
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Hossein Moradi-Shahrbabak
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran.
| | - Mohammad Hossein Fallahi
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Ali Jalil Sarghale
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Mohammad Moradi-Shahrbabak
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Rostam Abdollahi-Arpanahi
- Departments of Animal and Poultry Science, College of Aburaihan, University of Tehran, Pakdasht, 33916-53755, Iran
| | - Majid Khansefid
- AgriBio Centre for AgriBioscience, Agriculture Victoria, Bundoora, VIC, 3083, Australia
| |
Collapse
|
203
|
Getachew T, Haile A, Mészáros G, Rischkowsky B, Huson H, Gizaw S, Wurzinger M, Mwai A, Sölkner J. Genetic diversity, population structure and runs of homozygosity in Ethiopian short fat-tailed and Awassi sheep breeds using genome-wide 50k SNP markers. Livest Sci 2020. [DOI: 10.1016/j.livsci.2019.103899] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
204
|
Meyermans R, Gorssen W, Buys N, Janssens S. How to study runs of homozygosity using PLINK? A guide for analyzing medium density SNP data in livestock and pet species. BMC Genomics 2020; 21:94. [PMID: 31996125 PMCID: PMC6990544 DOI: 10.1186/s12864-020-6463-x] [Citation(s) in RCA: 111] [Impact Index Per Article: 27.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2019] [Accepted: 01/08/2020] [Indexed: 12/01/2022] Open
Abstract
Background PLINK is probably the most used program for analyzing SNP genotypes and runs of homozygosity (ROH), both in human and in animal populations. The last decade, ROH analyses have become the state-of-the-art method for inbreeding assessment. In PLINK, the --homozyg function is used to perform ROH analyses and relies on several input settings. These settings can have a large impact on the outcome and default values are not always appropriate for medium density SNP array data. Guidelines for a robust and uniform ROH analysis in PLINK using medium density data are lacking, albeit these guidelines are vital for comparing different ROH studies. In this study, 8 populations of different livestock and pet species are used to demonstrate the importance of PLINK input settings. Moreover, the effects of pruning SNPs for low minor allele frequencies and linkage disequilibrium on ROH detection are shown. Results We introduce the genome coverage parameter to appropriately estimate FROH and to check the validity of ROH analyses. The effect of pruning for linkage disequilibrium and low minor allele frequencies on ROH analyses is highly population dependent and such pruning may result in missed ROH. PLINK’s minimal density requirement is crucial for medium density genotypes and if set too low, genome coverage of the ROH analysis is limited. Finally, we provide recommendations for the maximal gap, scanning window length and threshold settings. Conclusions In this study, we present guidelines for an adequate and robust ROH analysis in PLINK on medium density SNP data. Furthermore, we advise to report parameter settings in publications, and to validate them prior to analysis. Moreover, we encourage authors to report genome coverage to reflect the ROH analysis’ validity. Implementing these guidelines will substantially improve the overall quality and uniformity of ROH analyses.
Collapse
Affiliation(s)
- R Meyermans
- Department of Biosystems, Livestock Genetics, KU Leuven, Kasteelpark Arenberg 30 - Box 2472, 3001, Leuven, Belgium
| | - W Gorssen
- Department of Biosystems, Livestock Genetics, KU Leuven, Kasteelpark Arenberg 30 - Box 2472, 3001, Leuven, Belgium
| | - N Buys
- Department of Biosystems, Livestock Genetics, KU Leuven, Kasteelpark Arenberg 30 - Box 2472, 3001, Leuven, Belgium
| | - S Janssens
- Department of Biosystems, Livestock Genetics, KU Leuven, Kasteelpark Arenberg 30 - Box 2472, 3001, Leuven, Belgium.
| |
Collapse
|
205
|
Barbato M, Hailer F, Upadhyay M, Del Corvo M, Colli L, Negrini R, Kim ES, Crooijmans RPMA, Sonstegard T, Ajmone-Marsan P. Adaptive introgression from indicine cattle into white cattle breeds from Central Italy. Sci Rep 2020; 10:1279. [PMID: 31992729 PMCID: PMC6987186 DOI: 10.1038/s41598-020-57880-4] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2019] [Accepted: 12/26/2019] [Indexed: 11/19/2022] Open
Abstract
Cattle domestication occurred at least twice independently and gave rise to the modern taurine and indicine cattle breeds. European cattle diversity is generally dominated by taurine cattle, although elevated levels of indicine ancestry have been recorded in several breeds from southern Europe. Here we use genome-wide high-density SNP genotyping data to investigate the taurine and indicine ancestry in southern European cattle, based on a dataset comprising 508 individuals from 23 cattle breeds of taurine, indicine and mixed ancestry, including three breeds from Central Italy known to exhibit the highest levels of indicine introgression among southern European breeds. Based on local genomic ancestry analyses, we reconstruct taurine and indicine ancestry genome-wide and along chromosomes. We scrutinise local genomic introgression signals and identify genomic regions that have introgressed from indicine into taurine cattle under positive selection, harbouring genes with functions related to body size and feed efficiency. These findings suggest that indicine-derived traits helped enhance Central Italian cattle through adaptive introgression. The identified genes could provide genomic targets for selection for improved cattle performance. Our findings elucidate the key role of adaptive introgression in shaping the phenotypic features of modern cattle, aided by cultural and livestock exchange among historic human societies.
Collapse
Affiliation(s)
- Mario Barbato
- Università Cattolica del Sacro Cuore, Department of Animal Science Food and Nutrition - DIANA, Nutrigenomics and Proteomics Research Centre - PRONUTRIGEN, and Biodiversity and Ancient DNA Research Centre, Università Cattolica del Sacro Cuore, Piacenza, Italy.
| | - Frank Hailer
- School of Biosciences, Cardiff University, Cardiff, Wales, UK
| | - Maulik Upadhyay
- Animal Breeding and Genomics, Wageningen University & Research, Wageningen, The Netherlands.,Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Marcello Del Corvo
- Università Cattolica del Sacro Cuore, Department of Animal Science Food and Nutrition - DIANA, Nutrigenomics and Proteomics Research Centre - PRONUTRIGEN, and Biodiversity and Ancient DNA Research Centre, Università Cattolica del Sacro Cuore, Piacenza, Italy
| | - Licia Colli
- Università Cattolica del Sacro Cuore, Department of Animal Science Food and Nutrition - DIANA, Nutrigenomics and Proteomics Research Centre - PRONUTRIGEN, and Biodiversity and Ancient DNA Research Centre, Università Cattolica del Sacro Cuore, Piacenza, Italy
| | - Riccardo Negrini
- Università Cattolica del Sacro Cuore, Department of Animal Science Food and Nutrition - DIANA, Nutrigenomics and Proteomics Research Centre - PRONUTRIGEN, and Biodiversity and Ancient DNA Research Centre, Università Cattolica del Sacro Cuore, Piacenza, Italy
| | | | | | | | - Paolo Ajmone-Marsan
- Università Cattolica del Sacro Cuore, Department of Animal Science Food and Nutrition - DIANA, Nutrigenomics and Proteomics Research Centre - PRONUTRIGEN, and Biodiversity and Ancient DNA Research Centre, Università Cattolica del Sacro Cuore, Piacenza, Italy.
| |
Collapse
|
206
|
Davenport KM, Hiemke C, McKay SD, Thorne JW, Lewis RM, Taylor T, Murdoch BM. Genetic structure and admixture in sheep from terminal breeds in the United States. Anim Genet 2020; 51:284-291. [PMID: 31970815 PMCID: PMC7065203 DOI: 10.1111/age.12905] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/12/2019] [Indexed: 11/29/2022]
Abstract
Selection for performance in diverse production settings has resulted in variation across sheep breeds worldwide. Although sheep are an important species to the United States, the current genetic relationship among many terminal sire breeds is not well characterized. Suffolk, Hampshire, Shropshire and Oxford (terminal) and Rambouillet (dual purpose) sheep (n = 248) sampled from different flocks were genotyped using the Applied Biosystems Axiom Ovine Genotyping Array (50K), and additional Shropshire sheep (n = 26) using the Illumina Ovine SNP50 BeadChip. Relationships were investigated by calculating observed heterozygosity, inbreeding coefficients, eigenvalues, pairwise Wright’s FST estimates and an identity by state matrix. The mean observed heterozygosity for each breed ranged from 0.30 to 0.35 and was consistent with data reported in other US and Australian sheep. Suffolk from two different regions of the United States (Midwest and West) clustered separately in eigenvalue plots and the rectangular cladogram. Further, divergence was detected between Suffolk from different regions with Wright’s FST estimate. Shropshire animals showed the greatest divergence from other terminal breeds in this study. Admixture between breeds was examined using admixture, and based on cross‐validation estimates, the best fit number of populations (clusters) was K = 6. The greatest admixture was observed within Hampshire, Suffolk, and Shropshire breeds. When plotting eigenvalues, US terminal breeds clustered separately in comparison with sheep from other locations of the world. Understanding the genetic relationships between terminal sire breeds in sheep will inform us about the potential applicability of markers derived in one breed to other breeds based on relatedness.
Collapse
Affiliation(s)
- K M Davenport
- Department of Animal and Veterinary Science, University of Idaho, Moscow, ID, 83844, USA
| | - C Hiemke
- Niman Ranch and Mapleton Mynd Shropshires, Stoughton, MA, 53589, USA
| | - S D McKay
- Department of Animal and Veterinary Sciences, University of Vermont, Burlington, VT, 05405, USA
| | - J W Thorne
- Department of Animal and Veterinary Science, University of Idaho, Moscow, ID, 83844, USA.,Texas A&M AgriLife Extension, San Angelo, TX, 76901, USA
| | - R M Lewis
- Department of Animal Science, University of Nebraska-Lincoln, Lincoln, NE, 68583, USA
| | - T Taylor
- Department of Animal Science, Arlington Research Station, University of Wisconsin-Madison, Arlington, WI, 53911, USA
| | - B M Murdoch
- Department of Animal and Veterinary Science, University of Idaho, Moscow, ID, 83844, USA
| |
Collapse
|
207
|
Li R, Li C, Chen H, Li R, Chong Q, Xiao H, Chen S. Genome-wide scan of selection signatures in Dehong humped cattle for heat tolerance and disease resistance. Anim Genet 2019; 51:292-299. [PMID: 31887783 DOI: 10.1111/age.12896] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/30/2019] [Indexed: 01/11/2023]
Abstract
Dehong humped cattle (DHH) is an indigenous zebu breed from southwestern China that possesses characteristics of heat tolerance and strong disease resistance and adapts well to the local tropical and subtropical climatic conditions. However, information on selection signatures of DHH is scarce. Herein, we compared the genomes of DHH and each of Diqing and Zhaotong cattle breeds using the population differentiation index (FST ), cross-population extended haplotype homozygosity (XP-EHH) and cross-population composite likelihood ratio (XP-CLR) methods to explore the genomic signatures of heat tolerance and disease resistance in DHH. Several pathways and genes carried selection signatures, including thermal sweating (calcium signaling pathway), heat shock (HSF1) and oxidative stress response (PLCB1, PLCB4), coat color (RAB31), feed intake (ATP8A1, SHC3) and reproduction (TP63, MAP3K13, PTPN4, PPP3CC, ADAMTSL1, SS18L1, OSBPL2, TOX, RREB1, GRK2). These identified pathways and genes may contribute to heat tolerance in DHH. Simultaneously, we also identified LIPH, TP63 and CBFA2T3 genes under positive selection that were associated with immunity.
Collapse
Affiliation(s)
- R Li
- School of Life Sciences, Yunnan University, Kunming, Yunnan, 650500, China
| | - C Li
- School of Life Sciences, Yunnan University, Kunming, Yunnan, 650500, China.,National Demonstration Center for Experimental Life Sciences Education, Yunnan University, Kunming, Yunnan, 650500, China
| | - H Chen
- School of Life Sciences, Yunnan University, Kunming, Yunnan, 650500, China
| | - R Li
- School of Life Sciences, Yunnan University, Kunming, Yunnan, 650500, China
| | - Q Chong
- School of Life Sciences, Yunnan University, Kunming, Yunnan, 650500, China
| | - H Xiao
- School of Life Sciences, Yunnan University, Kunming, Yunnan, 650500, China
| | - S Chen
- School of Life Sciences, Yunnan University, Kunming, Yunnan, 650500, China.,National Demonstration Center for Experimental Life Sciences Education, Yunnan University, Kunming, Yunnan, 650500, China
| |
Collapse
|
208
|
Meyermans R, Gorssen W, Wijnrocx K, Lenstra JA, Vellema P, Buys N, Janssens S. Unraveling the genetic diversity of Belgian Milk Sheep using medium-density SNP genotypes. Anim Genet 2019; 51:258-265. [PMID: 31881555 PMCID: PMC7065072 DOI: 10.1111/age.12891] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2019] [Revised: 07/23/2019] [Accepted: 11/13/2019] [Indexed: 01/05/2023]
Abstract
The present study focuses on the Belgian Milk Sheep in Flanders (Belgium) and compares its genetic diversity and relationship with the Flemish Sheep, the Friesian Milk Sheep, the French Lacaune dairy sheep and other Northern European breeds. For this study, 94 Belgian Milk Sheep, 23 Flemish Sheep and 22 Friesian Milk Sheep were genotyped with the OvineSNP50 array. In addition, 29 unregistered animals phenotypically similar to Belgian Milk Sheep were genotyped using the 15K ISGC chip. Both Belgian and Friesian Milk Sheep as well as the East Friesian Sheep were found to be less diverse than the other seven breeds included in this study. Genomic inbreeding coefficients based on runs of homozygosity (ROH) were estimated at 14.5, 12.4 and 10.2% for Belgian Milk Sheep, Flemish Sheep and Friesian Milk Sheep respectively. Out of 29 unregistered Belgian Milk Sheep, 28 mapped in the registered Belgian Milk Sheep population. Ancestry analysis, PCA and FST calculations showed that Belgian Milk Sheep are more related to Friesian Milk Sheep than to Flemish Sheep, which was contrary to the breeders' expectations. Consequently, breeders may prefer to crossbreed Belgian Milk Sheep with Friesian sheep populations (Friesian Milk Sheep or East Friesian Sheep) in order to increase diversity. This research underlines the usefulness of SNP chip genotyping and ROH analyses for monitoring genetic diversity and studying genetic links in small livestock populations, profiting from internationally available genotypes. As assessment of genetic diversity is vital for long-term breed survival, these results will aid flockbooks to preserve genetic diversity.
Collapse
Affiliation(s)
- R Meyermans
- Livestock Genetics, Department of Biosystems, KU Leuven, Kasteelpark Arenberg 30 - 2472, 3001, Leuven, Belgium
| | - W Gorssen
- Livestock Genetics, Department of Biosystems, KU Leuven, Kasteelpark Arenberg 30 - 2472, 3001, Leuven, Belgium
| | - K Wijnrocx
- Livestock Genetics, Department of Biosystems, KU Leuven, Kasteelpark Arenberg 30 - 2472, 3001, Leuven, Belgium
| | - J A Lenstra
- Faculty of Veterinary Medicine, Utrecht University, 3584CM, Utrecht, The Netherlands
| | - P Vellema
- Department of Small Ruminant Health, GD Animal Health, PO Box 9, 7400 AA, Deventer, The Netherlands
| | - N Buys
- Livestock Genetics, Department of Biosystems, KU Leuven, Kasteelpark Arenberg 30 - 2472, 3001, Leuven, Belgium
| | - S Janssens
- Livestock Genetics, Department of Biosystems, KU Leuven, Kasteelpark Arenberg 30 - 2472, 3001, Leuven, Belgium
| |
Collapse
|
209
|
Gebreselassie G, Berihulay H, Jiang L, Ma Y. Review on Genomic Regions and Candidate Genes Associated with Economically Important Production and Reproduction Traits in Sheep ( Ovies aries). Animals (Basel) 2019; 10:E33. [PMID: 31877963 PMCID: PMC7022721 DOI: 10.3390/ani10010033] [Citation(s) in RCA: 44] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2019] [Revised: 11/24/2019] [Accepted: 12/20/2019] [Indexed: 12/15/2022] Open
Abstract
Sheep (Ovis aries) is one of the most economically, culturally, and socially important domestic animals. They are reared primarily for meat, milk, wool, and fur production. Sheep were reared using natural selection for a long period of time to offer these traits. In fact, this production system has been slowing the productivity and production potential of the sheep. To improve production efficiency and productivity of this animal through genetic improvement technologies, understanding the genetic background of traits such as body growth, weight, carcass quality, fat percent, fertility, milk yield, wool quality, horn type, and coat color is essential. With the development and utilization of animal genotyping technologies and gene identification methods, many functional genes and genetic variants associated with economically important phenotypic traits have been identified and annotated. This is useful and presented an opportunity to increase the pace of animal genetic gain. Quantitative trait loci and genome wide association study have been playing an important role in identifying candidate genes and animal characterization. This review provides comprehensive information on the identified genomic regions and candidate genes associated with production and reproduction traits, and gene function in sheep.
Collapse
Affiliation(s)
- Gebremedhin Gebreselassie
- Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China; (G.G.); (H.B.); (L.J.)
- National Germplasm Center of Domestic Animal Resources, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China
- Department of Agricultural Biotechnology, Biotechnology Center, Ethiopian Biotechnology Institute, Ministry of Innovation and Technology, Addis Ababa 1000, Ethiopia
| | - Haile Berihulay
- Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China; (G.G.); (H.B.); (L.J.)
- National Germplasm Center of Domestic Animal Resources, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China
| | - Lin Jiang
- Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China; (G.G.); (H.B.); (L.J.)
- National Germplasm Center of Domestic Animal Resources, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China
| | - Yuehui Ma
- Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China; (G.G.); (H.B.); (L.J.)
- National Germplasm Center of Domestic Animal Resources, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China
| |
Collapse
|
210
|
Baazaoui I, McEwan J, Anderson R, Brauning R, McCulloch A, Van Stijn T, Bedhiaf-Romdhani S. GBS Data Identify Pigmentation-Specific Genes of Potential Role in Skin-Photosensitization in Two Tunisian Sheep Breeds. Animals (Basel) 2019; 10:ani10010005. [PMID: 31861491 PMCID: PMC7022847 DOI: 10.3390/ani10010005] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2019] [Revised: 12/09/2019] [Accepted: 12/11/2019] [Indexed: 02/06/2023] Open
Abstract
The Tunisian Noire de Thibar sheep breed is a composite breed, recently selected to create animals that are uniformly black in order to avoid skin photosensitization after the ingestion of toxic "hypericum perforatum" weeds, which causes a major economic loss to sheep farmers. We assessed genetic differentiation and estimated marker FST using genotyping-by-sequencing (GBS) data in black (Noire de Thibar) and related white-coated (Queue fine de l'ouest) sheep breeds to identify signals of artificial selection. The results revealed the selection signatures within candidate genes related to coat color, which are assumed to be indirectly involved in the mechanism of photosensitization in sheep. The identified genes could provide important information for molecular breeding.
Collapse
Affiliation(s)
- Imen Baazaoui
- Faculty of Science of Bizerte, University of Carthage, Carthage 1054, Tunisia
| | - John McEwan
- AgResearch Ltd., Invermay Agricultural Centre; Mosgiel 9092, New Zealand
| | - Rayna Anderson
- AgResearch Ltd., Invermay Agricultural Centre; Mosgiel 9092, New Zealand
| | - Rudiger Brauning
- AgResearch Ltd., Invermay Agricultural Centre; Mosgiel 9092, New Zealand
| | - Alan McCulloch
- AgResearch Ltd., Invermay Agricultural Centre; Mosgiel 9092, New Zealand
| | - Tracey Van Stijn
- AgResearch Ltd., Invermay Agricultural Centre; Mosgiel 9092, New Zealand
| | - Sonia Bedhiaf-Romdhani
- National Agricultural Research Institute of Tunisia, Laboratory of Animal and forage Production, University of Carthage, Ariana 1004, Tunisia
- Correspondence: ; Tel.: +216-25-113-344
| |
Collapse
|
211
|
Deniskova T, Dotsev A, Lushihina E, Shakhin A, Kunz E, Medugorac I, Reyer H, Wimmers K, Khayatzadeh N, Sölkner J, Sermyagin A, Zhunushev A, Brem G, Zinovieva N. Population Structure and Genetic Diversity of Sheep Breeds in the Kyrgyzstan. Front Genet 2019; 10:1311. [PMID: 31921318 PMCID: PMC6922024 DOI: 10.3389/fgene.2019.01311] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2019] [Accepted: 11/28/2019] [Indexed: 11/29/2022] Open
Abstract
Sheep are a main livestock species of Kyrgyzstan, a Central Asian country with predominating mountain terrain. The current gene pool of local sheep resources has been forming under diverse climate conditions from the era of the trading caravans of the Great Silk Road, through the Soviet period of large-scale livestock improvements, which was followed by the deep crisis at the end of the 20th century, up to now. However, not much is known about the genetic background and variability of the local sheep populations. Therefore, our aims were to provide a characterization of the population structure and genetic relations within the Kyrgyz sheep breeds and to study their genetic connections with the global sheep breeds using SNP analysis. Samples of the Alai (n = 31), Gissar (n = 30), Kyrgyz coarse wool (n = 13), Aykol (n = 31), and Tien-Shan (n = 24) breeds were genotyped with the OvineSNP50 BeadChip or the Ovine Infinium HD BeadChip (Illumina Inc., USA). The measure of inbreeding based on runs of homozygosity showed a minimum value in the Aykol breed (FROH = 0.034), while the maximum was found in the Alai breed (FROH = 0.071). Short ROH segments (ROH ≤ 4 Mb) were predominant in all breeds. Long ROH segments (ROH > 16 Mb) were absent in the Gissar breed. The Gissar and Aykol breeds had the highest values of the effective population sizes estimated for five generations ago (Ne5 = 660 and 563), whereas the Alai and Kyrgyz coarse wool displayed lower values (Ne5 = 176 and 128, respectively). The synthetic origin of the Aykol breed was clearly evidenced by all analyses applied. Based on the network and admixture analyses of the Kyrgyz and global sheep breeds, the Tien-Shan and the Russian semi-fine wool breeds demonstrated a common ancestry that most likely is due to a contribution of the Lincoln breed. The Gissar, Aykol, and Kyrgyz coarse wool breeds showed a genetic background predominating in sheep populations from Iran and China whereas the Alai demonstrated the different ancestry type. The revealed admixture patterns probably resulted from the exchange and trade during the era of the Great Silk Road, which partly overlapped with historical and archeological findings.
Collapse
Affiliation(s)
- Tatiana Deniskova
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, Russia
| | - Arsen Dotsev
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, Russia
| | - Eugenia Lushihina
- Institute of Biotechnology, National Academy of Science of Kyrgyz Republic, Bishkek, Kyrgyzstan
| | - Alexey Shakhin
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, Russia
| | - Elisabeth Kunz
- Population Genomics Group, Department of Veterinary Sciences, Ludwig Maximilians University of Munich, Munich, Germany
| | - Ivica Medugorac
- Population Genomics Group, Department of Veterinary Sciences, Ludwig Maximilians University of Munich, Munich, Germany
| | - Henry Reyer
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology, Dummerstorf, Germany
| | - Klaus Wimmers
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology, Dummerstorf, Germany
| | - Negar Khayatzadeh
- Division of Livestock Sciences, University of Natural Resources and Life Sciences Vienna, Vienna, Austria
| | - Johann Sölkner
- Division of Livestock Sciences, University of Natural Resources and Life Sciences Vienna, Vienna, Austria
| | | | - Asankadyr Zhunushev
- Institute of Biotechnology, National Academy of Science of Kyrgyz Republic, Bishkek, Kyrgyzstan
| | - Gottfried Brem
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, Russia.,Institute of Animal Breeding and Genetics, University of Veterinary Medicine, Vienna, Austria
| | - Natalia Zinovieva
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, Russia
| |
Collapse
|
212
|
Eck K, Kunz E, Mendel C, Lühken G, Medugorac I. Morphometric measurements in lambs as a basis for future mapping studies. Small Rumin Res 2019. [DOI: 10.1016/j.smallrumres.2019.04.007] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
213
|
Yaro M, Munyard KA, Morgan E, Allcock RJN, Stear MJ, Groth DM. Analysis of pooled genome sequences from Djallonke and Sahelian sheep of Ghana reveals co-localisation of regions of reduced heterozygosity with candidate genes for disease resistance and adaptation to a tropical environment. BMC Genomics 2019; 20:816. [PMID: 31699027 PMCID: PMC6836352 DOI: 10.1186/s12864-019-6198-8] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2018] [Accepted: 10/16/2019] [Indexed: 12/17/2022] Open
Abstract
BACKGROUND The Djallonke sheep is well adapted to harsh environmental conditions, and is relatively resistant to Haemonchosis and resilient to animal trypanosomiasis. The larger Sahelian sheep, which cohabit the same region, is less well adapted to these disease challenges. Haemonchosis and Trypanosomiasis collectively cost the worldwide animal industry billions of dollars in production losses annually. RESULTS Here, we separately sequenced and then pooled according to breed the genomes from five unrelated individuals from each of the Djallonke and Sahelian sheep breeds (sourced from Ghana), at greater than 22-fold combined coverage for each breed. A total of approximately 404 million (97%) and 343 million (97%) sequence reads from the Djallonke and Sahelian breeds respectively, were successfully mapped to the sheep reference genome Oar v3.1. We identified approximately 11.1 million and 10.9 million single nucleotide polymorphisms (SNPs) in the Djallonke and Sahelian breeds, with approximately 15 and 16% respectively of these not previously reported in sheep. Multiple regions of reduced heterozygosity were also found; 70 co-localised within genomic regions harbouring genes that mediate disease resistance, immune response and adaptation in sheep or cattle. Thirty- three of the regions of reduced heterozygosity co-localised with previously reported genes for resistance to haemonchosis and trypanosomiasis. CONCLUSIONS Our analyses suggest that these regions of reduced heterozygosity may be signatures of selection for these economically important diseases.
Collapse
Affiliation(s)
- M. Yaro
- School of Biomedical Sciences, CHIRI Biosciences Research Precinct, Faculty of Health Sciences, Curtin University, GPO Box U1987, Perth, WA 6845 Australia
| | - K. A. Munyard
- School of Biomedical Sciences, CHIRI Biosciences Research Precinct, Faculty of Health Sciences, Curtin University, GPO Box U1987, Perth, WA 6845 Australia
| | - E. Morgan
- School of Biomedical Sciences, CHIRI Biosciences Research Precinct, Faculty of Health Sciences, Curtin University, GPO Box U1987, Perth, WA 6845 Australia
| | - R. J. N. Allcock
- The University of Western Australia, 35 Stirling Highway, Crawley, Perth, WA Australia
- Pathwest Laboratory Medicine WA, QEII Medical Centre, Monash Avenue, Nedlands, 6009 Australia
| | - M. J. Stear
- Agribio centre for Agribioscience, La Trobe University, Melbourne, Australia
- Institute of Biodiversity, Animal Health and Comparative Medicine University of Glasgow, Bearsden Road, Glasgow, G61 1QH UK
| | - D. M. Groth
- School of Biomedical Sciences, CHIRI Biosciences Research Precinct, Faculty of Health Sciences, Curtin University, GPO Box U1987, Perth, WA 6845 Australia
| |
Collapse
|
214
|
Sparks AM, Watt K, Sinclair R, Pilkington JG, Pemberton JM, McNeilly TN, Nussey DH, Johnston SE. The genetic architecture of helminth-specific immune responses in a wild population of Soay sheep (Ovis aries). PLoS Genet 2019; 15:e1008461. [PMID: 31697674 PMCID: PMC6863570 DOI: 10.1371/journal.pgen.1008461] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2019] [Revised: 11/19/2019] [Accepted: 10/04/2019] [Indexed: 12/22/2022] Open
Abstract
Much of our knowledge of the drivers of immune variation, and how these responses vary over time, comes from humans, domesticated livestock or laboratory organisms. While the genetic basis of variation in immune responses have been investigated in these systems, there is a poor understanding of how genetic variation influences immunity in natural, untreated populations living in complex environments. Here, we examine the genetic architecture of variation in immune traits in the Soay sheep of St Kilda, an unmanaged population of sheep infected with strongyle gastrointestinal nematodes. We assayed IgA, IgE and IgG antibodies against the prevalent nematode Teladorsagia circumcincta in the blood plasma of > 3,000 sheep collected over 26 years. Antibody levels were significantly heritable (h2 = 0.21 to 0.57) and highly stable over an individual’s lifespan. IgA levels were strongly associated with a region on chromosome 24 explaining 21.1% and 24.5% of heritable variation in lambs and adults, respectively. This region was adjacent to two candidate loci, Class II Major Histocompatibility Complex Transactivator (CIITA) and C-Type Lectin Domain Containing 16A (CLEC16A). Lamb IgA levels were also associated with the immunoglobulin heavy constant loci (IGH) complex, and adult IgE levels and lamb IgA and IgG levels were associated with the major histocompatibility complex (MHC). This study provides evidence of high heritability of a complex immunological trait under natural conditions and provides the first evidence from a genome-wide study that large effect genes located outside the MHC region exist for immune traits in the wild. Understanding how immune responses vary in natural populations can give an insight into how infection affects the ability of hosts and parasites to survive and reproduce, and how this drives evolutionary and ecological dynamics. Yet, very little is known about how immune responses vary over an individual’s lifetime and how genes contribute to this variation under natural conditions. Our study investigates the genetic architecture of variation in three antibody types, IgA, IgE and IgG in a wild population of Soay sheep on the St Kilda archipelago in North-West Scotland. Using data collected over 26 years, we show that antibody levels have a heritable basis in lambs and adults and are stable over an individual’s lifetime. We also identify several genomic regions with large effects on immune responses. Our study offers the first insights into the genetic control of immunity in a wild population, which is essential to understand how immune profiles vary in challenging natural conditions and how natural selection maintains genetic variation in complex immune traits.
Collapse
Affiliation(s)
- Alexandra M. Sparks
- Institutes of Evolutionary Biology and Immunology and Infection Research, School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom
- Faculty of Biological Sciences, School of Biology, University of Leeds, Leeds, United Kingdom
- * E-mail:
| | - Kathryn Watt
- Institutes of Evolutionary Biology and Immunology and Infection Research, School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom
| | - Rona Sinclair
- Institutes of Evolutionary Biology and Immunology and Infection Research, School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom
| | - Jill G. Pilkington
- Institutes of Evolutionary Biology and Immunology and Infection Research, School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom
| | - Josephine M. Pemberton
- Institutes of Evolutionary Biology and Immunology and Infection Research, School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom
| | - Tom N. McNeilly
- Moredun Research Institute, Pentlands Science Park, Bush Loan, Midlothian, United Kingdom
| | - Daniel H. Nussey
- Institutes of Evolutionary Biology and Immunology and Infection Research, School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom
| | - Susan E. Johnston
- Institutes of Evolutionary Biology and Immunology and Infection Research, School of Biological Sciences, University of Edinburgh, Edinburgh, United Kingdom
| |
Collapse
|
215
|
Benjelloun B, Boyer F, Streeter I, Zamani W, Engelen S, Alberti A, Alberto FJ, BenBati M, Ibnelbachyr M, Chentouf M, Bechchari A, Rezaei HR, Naderi S, Stella A, Chikhi A, Clarke L, Kijas J, Flicek P, Taberlet P, Pompanon F. An evaluation of sequencing coverage and genotyping strategies to assess neutral and adaptive diversity. Mol Ecol Resour 2019; 19:1497-1515. [PMID: 31359622 PMCID: PMC7115901 DOI: 10.1111/1755-0998.13070] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2018] [Revised: 06/30/2019] [Accepted: 07/08/2019] [Indexed: 12/12/2022]
Abstract
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.
Collapse
Affiliation(s)
- Badr Benjelloun
- Univ. Grenoble-Alpes, Univ. Savoie Mont Blanc, CNRS, LECA, F-38000 Grenoble, France
- National Institute of Agronomic Research (INRA Maroc), Regional Centre of Agronomic Research, 23000 Beni-Mellal, Morocco
| | - Frédéric Boyer
- Univ. Grenoble-Alpes, Univ. Savoie Mont Blanc, CNRS, LECA, F-38000 Grenoble, France
| | - Ian Streeter
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SD UK
| | - Wahid Zamani
- Univ. Grenoble-Alpes, Univ. Savoie Mont Blanc, CNRS, LECA, F-38000 Grenoble, France
- Department of Environmental Sciences, Faculty of Natural Resources and Marine Sciences, Tarbiat Modares University, 46417-76489 Noor, Mazandaran, Iran
| | - Stefan Engelen
- CEA - Institut de biologie François-Jacob, Genoscope, 2 Rue Gaston Cremieux 91057 Evry Cedex, France
| | - Adriana Alberti
- CEA - Institut de biologie François-Jacob, Genoscope, 2 Rue Gaston Cremieux 91057 Evry Cedex, France
| | - Florian J. Alberto
- Univ. Grenoble-Alpes, Univ. Savoie Mont Blanc, CNRS, LECA, F-38000 Grenoble, France
| | - Mohamed BenBati
- National Institute of Agronomic Research (INRA Maroc), Regional Centre of Agronomic Research, 23000 Beni-Mellal, Morocco
| | - Mustapha Ibnelbachyr
- National Institute of Agronomic Research (INRA Maroc), CRRA Errachidia, 52000 Errachidia, Morocco
| | - Mouad Chentouf
- National Institute of Agronomic Research (INRA Maroc), CRRA Tangier, 90010 Tangier, Morocco
| | - Abdelmajid Bechchari
- National Institute of Agronomic Research (INRA Maroc), CRRA Oujda, 60000 Oujda, Morocco
| | - Hamid R. Rezaei
- Department of Environmental Sci, Gorgan University of Agricultural Sciences & Natural Resources, 41996-13776 Gorgan, Iran
| | - Saeid Naderi
- Environmental Sciences Department, Natural Resources Faculty, University of Guilan, 49138-15749 Guilan, Iran
| | - Alessandra Stella
- PTP Science Park, Bioinformatics Unit, Via Einstein-Loc. Cascina Codazza, 26900 Lodi, Italy
| | - Abdelkader Chikhi
- National Institute of Agronomic Research (INRA Maroc), CRRA Errachidia, 52000 Errachidia, Morocco
| | - Laura Clarke
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SD UK
| | - James Kijas
- Commonwealth Scientific and Industrial Research Organisation Animal Food and Health Sciences, St Lucia, QLD 4067, Australia
| | - Paul Flicek
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SD UK
| | - Pierre Taberlet
- Univ. Grenoble-Alpes, Univ. Savoie Mont Blanc, CNRS, LECA, F-38000 Grenoble, France
| | - François Pompanon
- Univ. Grenoble-Alpes, Univ. Savoie Mont Blanc, CNRS, LECA, F-38000 Grenoble, France
| |
Collapse
|
216
|
Megdiche S, Mastrangelo S, Ben Hamouda M, Lenstra JA, Ciani E. A Combined Multi-Cohort Approach Reveals Novel and Known Genome-Wide Selection Signatures for Wool Traits in Merino and Merino-Derived Sheep Breeds. Front Genet 2019; 10:1025. [PMID: 31708969 PMCID: PMC6824410 DOI: 10.3389/fgene.2019.01025] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2019] [Accepted: 09/24/2019] [Indexed: 12/24/2022] Open
Abstract
Merino sheep represents a valuable genetic resource worldwide. In this study, we investigated selection signatures in Merino (and Merino-derived) sheep breeds using genome-wide SNP data and two different approaches: a classical FST-outlier method and an approach based on the analysis of local ancestry in admixed populations. In order to capture the most reliable signals, we adopted a combined, multi-cohort approach. In particular, scenarios involving four Merino breeds (Spanish Merino, Australian Merino, Chinese Merino, and Sopravissana) were tested via the local ancestry approach, while nine pair-wise breed comparisons contrasting the above breeds, as well as the Gentile di Puglia breed, with non-Merino breeds from the same geographic area were tested via the FST-outlier method. Signals observed using both methods were compared with genome-wide patterns of distribution of runs of homozygosity (ROH) islands. Novel and known selection signatures were detected. The most reliable signals were observed on OAR 3 (MSRB3 and LEMD3), OAR10 (FRY and RXFP2), OAR 13 (RALY), OAR17 (FAM101A), and OAR18 (NFKBIA, SEC23A, and PAX9). All the above overlapped with known QTLs for wool traits, and evidences from the literature of their involvement in skin/hair/wool biology, as well as gene network analysis, further corroborated these results. The signal on OAR10 also contains well known evidence for association with horn morphology and polledness. More elusive biological evidences of association with the Merino phenotype were observed for a number of other genes, notably LOC101120019 and TMEM132B (OAR17), LOC105609948 (OAR3), LOC101110773 (OAR10), and EIF2S2 (OAR17). Taken together, the above results further contribute to decipher the genetic basis underlying the Merino phenotype.
Collapse
Affiliation(s)
- Sami Megdiche
- Départment des Ressources Animales, Agroalimentaire et Développement Rural, Institut Supérieur Agronomique de Chott-Mariem, Université de Sousse, Sousse, Tunisia.,Dipartimento di Bioscienze, Biotecnologie e Biofarmaceutica, University of Bari "Aldo Moro," Bari, Italy
| | - Salvatore Mastrangelo
- Dipartimento di Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | | | - Johannes A Lenstra
- Faculty of Veterinary Medicine, Utrecht University, Utrecht, Netherlands
| | - Elena Ciani
- Dipartimento di Bioscienze, Biotecnologie e Biofarmaceutica, University of Bari "Aldo Moro," Bari, Italy
| |
Collapse
|
217
|
Population structure and breed composition prediction in a multi-breed sheep population using genome-wide single nucleotide polymorphism genotypes. Animal 2019; 14:464-474. [PMID: 31610818 DOI: 10.1017/s1751731119002398] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
Knowledge of population structure and breed composition of a population can be advantageous for a number of reasons; these include designing optimal (cross)breeding strategies in order to maximise non-additive genetic effects, maintaining flockbook integrity by authenticating animals being registered and as a quality control measure in the genotyping process. The objectives of the present study were to 1) describe the population structure of 24 sheep breeds, 2) quantify the breed composition of both flockbook-recorded and crossbred animals using single nucleotide polymorphism BLUP (SNP-BLUP), and 3) quantify the accuracy of breed composition prediction from low-density genotype panels containing between 2000 and 6000 SNPs. In total, 9334 autosomal SNPs on 11 144 flockbook-recorded animals and 1172 crossbred animals were used. The population structure of all breeds was characterised by principal component analysis (PCA) as well as the pairwise breed fixation index (Fst). The total number of animals, all of which were purebred, included in the calibration population for SNP-BLUP was 2579 with the number of animals per breed ranging from 9 to 500. The remaining 9559 flockbook-recorded animals, composite breeds and crossbred animals represented the test population; three breeds were excluded from breed composition prediction. The breed composition predicted using SNP-BLUP with 9334 SNPs was considered the gold standard prediction. The pairwise breed Fst ranged from 0.040 (between the Irish Blackface and Scottish Blackface) to 0.282 (between the Border Leicester and Suffolk). Principal component analysis revealed that the Suffolk from Ireland and the Suffolk from New Zealand formed distinct, non-overlapping clusters. In contrast, the Texel from Ireland and that from New Zealand formed integrated, overlapping clusters. Composite animals such as the Belclare clustered close to its founder breeds (i.e., Finn, Galway, Lleyn and Texel). When all 9334 SNPs were used to predict breed composition, an animal that had a majority breed proportion predicted to be ≥0.90 was defined as purebred for the present study. As the panel density decreased, the predicted breed proportion threshold, used to identify animals as purebred, also decreased (≥0.85 with 6000 SNPs to ≥0.60 with 2000 SNPs). In all, results from the study suggest that breed composition for purebred and crossbred animals can be determined with SNP-BLUP using ≥5000 SNPs.
Collapse
|
218
|
Sim Z, Coltman DW. Heritability of Horn Size in Thinhorn Sheep. Front Genet 2019; 10:959. [PMID: 31681413 PMCID: PMC6797622 DOI: 10.3389/fgene.2019.00959] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2019] [Accepted: 09/09/2019] [Indexed: 12/31/2022] Open
Abstract
Understanding the genetic basis of fitness-related trait variation has long been of great interest to evolutionary biologists. Secondary sexual characteristics, such as horns in bovids, are particularly intriguing since they can be potentially affected by both natural and sexual selection. Until recently, however, the study of fitness-related quantitative trait variation in wild species has been hampered by a lack of genomic resources, pedigree, and/or phenotype data. Recent innovations in genomic technologies have enabled wildlife researchers to perform marker-based relatedness estimation and acquire adequate loci density, enabling both the “top-down” approach of quantitative genetics and the “bottom-up” approach of association studies to describe the genetic basis of fitness-related traits. Here we combine a cross species application of the OvineHD BeadChip and horn measurements (horn length, base circumference, and volume) from harvested thinhorn sheep to examine the heritability and to perform a genome-wide single-nucleotide polymorphism association study of horn size in the species. Thinhorn sheep are mountain ungulates that reside in the mountainous regions of northwestern North America. Thinhorn sheep males grow massive horns that determine the social rank and mating success. We found horn length, base circumference, and volume to be moderately heritable and two loci to be suggestively associated with horn length.
Collapse
Affiliation(s)
- Zijian Sim
- Department of Biological Sciences, University of Alberta, Edmonton, AB, Canada.,Fish and Wildlife Forensic Unit, Alberta Fish and Wildlife Enforcement Branch, Government of Alberta, Edmonton, AB, Canada
| | - David W Coltman
- Department of Biological Sciences, University of Alberta, Edmonton, AB, Canada
| |
Collapse
|
219
|
McHugo GP, Browett S, Randhawa IAS, Howard DJ, Mullen MP, Richardson IW, Park SDE, Magee DA, Scraggs E, Dover MJ, Correia CN, Hanrahan JP, MacHugh DE. A Population Genomics Analysis of the Native Irish Galway Sheep Breed. Front Genet 2019; 10:927. [PMID: 31649720 PMCID: PMC6792165 DOI: 10.3389/fgene.2019.00927] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2019] [Accepted: 09/05/2019] [Indexed: 12/24/2022] Open
Abstract
The Galway sheep population is the only native Irish sheep breed and this livestock genetic resource is currently categorised as 'at-risk'. In the present study, comparative population genomics analyses of Galway sheep and other sheep populations of European origin were used to investigate the microevolution and recent genetic history of the breed. These analyses support the hypothesis that British Leicester sheep were used in the formation of the Galway. When compared to conventional and endangered breeds, the Galway breed was intermediate in effective population size, genomic inbreeding and runs of homozygosity. This indicates that, although the Galway breed is declining, it is still relatively genetically diverse and that conservation and management plans informed by genomic information may aid its recovery. The Galway breed also exhibited distinct genomic signatures of artificial or natural selection when compared to other breeds, which highlighted candidate genes that may be involved in production and health traits.
Collapse
Affiliation(s)
- Gillian P McHugo
- Animal Genomics Laboratory, UCD School of Agriculture and Food Science, University College Dublin, Dublin, Ireland
| | - Sam Browett
- Animal Genomics Laboratory, UCD School of Agriculture and Food Science, University College Dublin, Dublin, Ireland
| | - Imtiaz A S Randhawa
- Sydney School of Veterinary Science, University of Sydney, Camden, NSW, Australia
| | - Dawn J Howard
- Animal and Grassland Research and Innovation Centre, Athenry, Ireland
| | - Michael P Mullen
- Animal and Grassland Research and Innovation Centre, Athenry, Ireland
| | | | | | - David A Magee
- Animal Genomics Laboratory, UCD School of Agriculture and Food Science, University College Dublin, Dublin, Ireland
| | - Erik Scraggs
- Animal Genomics Laboratory, UCD School of Agriculture and Food Science, University College Dublin, Dublin, Ireland
| | - Michael J Dover
- Animal Genomics Laboratory, UCD School of Agriculture and Food Science, University College Dublin, Dublin, Ireland
| | - Carolina N Correia
- Animal Genomics Laboratory, UCD School of Agriculture and Food Science, University College Dublin, Dublin, Ireland
| | - James P Hanrahan
- Animal and Grassland Research and Innovation Centre, Athenry, Ireland
| | - David E MacHugh
- Animal Genomics Laboratory, UCD School of Agriculture and Food Science, University College Dublin, Dublin, Ireland.,UCD Conway Institute of Biomolecular and Biomedical Research, University College Dublin, Dublin, Ireland
| |
Collapse
|
220
|
Smith CCR, Flaxman SM. Leveraging whole genome sequencing data for demographic inference with approximate Bayesian computation. Mol Ecol Resour 2019; 20:125-139. [DOI: 10.1111/1755-0998.13092] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2019] [Revised: 08/30/2019] [Accepted: 09/06/2019] [Indexed: 01/16/2023]
Affiliation(s)
- Chris C. R. Smith
- Department of Ecology and Evolutionary Biology University of Colorado Boulder CO USA
| | - Samuel M. Flaxman
- Department of Ecology and Evolutionary Biology University of Colorado Boulder CO USA
| |
Collapse
|
221
|
Salavati M, Bush SJ, Palma-Vera S, McCulloch MEB, Hume DA, Clark EL. Elimination of Reference Mapping Bias Reveals Robust Immune Related Allele-Specific Expression in Crossbred Sheep. Front Genet 2019; 10:863. [PMID: 31608110 PMCID: PMC6761296 DOI: 10.3389/fgene.2019.00863] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2019] [Accepted: 08/19/2019] [Indexed: 12/13/2022] Open
Abstract
Pervasive allelic variation at both gene and single nucleotide level (SNV) between individuals is commonly associated with complex traits in humans and animals. Allele-specific expression (ASE) analysis, using RNA-Seq, can provide a detailed annotation of allelic imbalance and infer the existence of cis-acting transcriptional regulation. However, variant detection in RNA-Seq data is compromised by biased mapping of reads to the reference DNA sequence. In this manuscript, we describe an unbiased standardized computational pipeline for allele-specific expression analysis using RNA-Seq data, which we have adapted and developed using tools available under open license. The analysis pipeline we present is designed to minimize reference bias while providing accurate profiling of allele-specific expression across tissues and cell types. Using this methodology, we were able to profile pervasive allelic imbalance across tissues and cell types, at both the gene and SNV level, in Texel×Scottish Blackface sheep, using the sheep gene expression atlas data set. ASE profiles were pervasive in each sheep and across all tissue types investigated. However, ASE profiles shared across tissues were limited, and instead, they tended to be highly tissue-specific. These tissue-specific ASE profiles may underlie the expression of economically important traits and could be utilized as weighted SNVs, for example, to improve the accuracy of genomic selection in breeding programs for sheep. An additional benefit of the pipeline is that it does not require parental genotypes and can therefore be applied to other RNA-Seq data sets for livestock, including those available on the Functional Annotation of Animal Genomes (FAANG) data portal. This study is the first global characterization of moderate to extreme ASE in tissues and cell types from sheep. We have applied a robust methodology for ASE profiling to provide both a novel analysis of the multi-dimensional sheep gene expression atlas data set and a foundation for identifying the regulatory and expressed elements of the genome that are driving complex traits in livestock.
Collapse
Affiliation(s)
- Mazdak Salavati
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, Edinburgh, United Kingdom
| | - Stephen J. Bush
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, Edinburgh, United Kingdom
| | - Sergio Palma-Vera
- Leibniz Institute for Farm Animal Biology (FBN), Institute for Reproductive Biology, Dummerstorf, Germany
| | - Mary E. B. McCulloch
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, Edinburgh, United Kingdom
| | - David A. Hume
- Mater Research Institute-University of Queensland, Translational Research Institute, Woolloongabba, QLD, Australia
| | - Emily L. Clark
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, Edinburgh, United Kingdom
| |
Collapse
|
222
|
Ablondi M, Viklund Å, Lindgren G, Eriksson S, Mikko S. Signatures of selection in the genome of Swedish warmblood horses selected for sport performance. BMC Genomics 2019; 20:717. [PMID: 31533613 PMCID: PMC6751828 DOI: 10.1186/s12864-019-6079-1] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2019] [Accepted: 09/04/2019] [Indexed: 01/09/2023] Open
Abstract
Background A growing demand for improved physical skills and mental attitude in modern sport horses has led to strong selection for performance in many warmblood studbooks. The aim of this study was to detect genomic regions with low diversity, and therefore potentially under selection, in Swedish Warmblood horses (SWB) by analysing high-density SNP data. To investigate if such signatures could be the result of selection for equestrian sport performance, we compared our SWB SNP data with those from Exmoor ponies, a horse breed not selected for sport performance traits. Results The genomic scan for homozygous regions identified long runs of homozygosity (ROH) shared by more than 85% of the genotyped SWB individuals. Such ROH were located on ECA4, ECA6, ECA7, ECA10 and ECA17. Long ROH were instead distributed evenly across the genome of Exmoor ponies in 77% of the chromosomes. Two population differentiation tests (FST and XP-EHH) revealed signatures of selection on ECA1, ECA4, and ECA6 in SWB horses. Conclusions Genes related to behaviour, physical abilities and fertility, appear to be targets of selection in the SWB breed. This study provides a genome-wide map of selection signatures in SWB horses, and ground for further functional studies to unravel the biological mechanisms behind complex traits in horses.
Collapse
Affiliation(s)
- Michela Ablondi
- Dept. of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, PO Box 7023, S-750 07, Uppsala, Sweden.,Department of Veterinary Science, Università degli Studi di Parma, 43126, Parma, Italy
| | - Åsa Viklund
- Dept. of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, PO Box 7023, S-750 07, Uppsala, Sweden
| | - Gabriella Lindgren
- Dept. of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, PO Box 7023, S-750 07, Uppsala, Sweden.,Livestock Genetics, Department of Biosystems, Leuven, KU, Belgium
| | - Susanne Eriksson
- Dept. of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, PO Box 7023, S-750 07, Uppsala, Sweden
| | - Sofia Mikko
- Dept. of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, PO Box 7023, S-750 07, Uppsala, Sweden.
| |
Collapse
|
223
|
Rochus CM, Westberg Sunesson K, Jonas E, Mikko S, Johansson AM. Mutations in ASIP and MC1R: dominant black and recessive black alleles segregate in native Swedish sheep populations. Anim Genet 2019; 50:712-717. [PMID: 31475378 DOI: 10.1111/age.12837] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/08/2019] [Indexed: 01/03/2023]
Abstract
By studying genes associated with coat colour, we can understand the role of these genes in pigmentation but also gain insight into selection history. North European short-tailed sheep, including Swedish breeds, have variation in their coat colour, making them good models to expand current knowledge of mutations associated with coat colour in sheep. We studied ASIP and MC1R, two genes with known roles in pigmentation, and their association with black coat colour. We did this by sequencing the coding regions of ASIP in 149 animals and MC1R in 129 animals from seven native Swedish sheep breeds in individuals with black, white or grey fleece. Previously known mutations in ASIP [recessive black allele: g.100_105del (D5 ) and/or g.5172T>A] were associated with black coat colour in Klövsjö and Roslag sheep breeds and mutations in both ASIP and MC1R (dominant black allele: c.218T>A and/or c.361G>A) were associated with black coat colour in Swedish Finewool. In Gotland, Gute, Värmland and Helsinge sheep breeds, coat colour inheritance was more complex: only 11 of 16 individuals with black fleece had genotypes that could explain their black colour. These breeds have grey individuals in their populations, and grey is believed to be a result of mutations and allelic copy number variation within the ASIP duplication, which could be a possible explanation for the lack of a clear inheritance pattern in these breeds. Finally, we found a novel missense mutation in MC1R (c.452G>A) in Gotland, Gute and Värmland sheep and evidence of a duplication of MC1R in Gotland sheep.
Collapse
Affiliation(s)
- C M Rochus
- Department of Animal Breeding and Genetics, Faculty of Veterinary Medicine and Animal Science, Swedish University of Agricultural Sciences, Box 7923, SE-75007, Uppsala, Sweden.,UFR Génétique, Élevage et Reproduction, Sciences de la Vie et Santé, AgroParisTech, Université Paris Saclay, 16 rue Claude Bernard, F-75231, Paris Cedex 05, France.,Génétique Physiologie Systèmes d'Elevage, Animal Genetics Division, INRA, 24 chemin de Borde-Rouge-Auzeville Tolosane, F-31326 Castanet-Tolosan, France
| | - K Westberg Sunesson
- Department of Animal Breeding and Genetics, Faculty of Veterinary Medicine and Animal Science, Swedish University of Agricultural Sciences, Box 7923, SE-75007, Uppsala, Sweden
| | - E Jonas
- Department of Animal Breeding and Genetics, Faculty of Veterinary Medicine and Animal Science, Swedish University of Agricultural Sciences, Box 7923, SE-75007, Uppsala, Sweden
| | - S Mikko
- Department of Animal Breeding and Genetics, Faculty of Veterinary Medicine and Animal Science, Swedish University of Agricultural Sciences, Box 7923, SE-75007, Uppsala, Sweden
| | - A M Johansson
- Department of Animal Breeding and Genetics, Faculty of Veterinary Medicine and Animal Science, Swedish University of Agricultural Sciences, Box 7923, SE-75007, Uppsala, Sweden
| |
Collapse
|
224
|
Beeson SK, Mickelson JR, McCue ME. Exploration of fine-scale recombination rate variation in the domestic horse. Genome Res 2019; 29:1744-1752. [PMID: 31434677 PMCID: PMC6771410 DOI: 10.1101/gr.243311.118] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2018] [Accepted: 08/15/2019] [Indexed: 01/17/2023]
Abstract
Total genetic map length and local recombination landscapes typically vary within and across populations. As a first step to understanding the recombination landscape in the domestic horse, we calculated population recombination rates and identified likely recombination hotspots using approximately 1.8 million SNP genotypes for 485 horses from 32 distinct breeds. The resulting breed-averaged recombination map spans 2.36 Gb and accounts for 2939.07 cM. Recombination hotspots occur once per 23.8 Mb on average and account for ∼9% of the physical map length. Regions with elevated recombination rates in the entire cohort were enriched for genes in pathways involving interaction with the environment: immune system processes (specifically, MHC class I and class II genes), responses to stimuli, and serotonin receptor pathways. We found significant correlations between differences in local recombination rates and population differentiation quantified by F ST Analysis of breed-specific maps revealed thousands of hotspot regions unique to particular breeds, as well as unique "coldspots," regions where a particular breed showed below-average recombination, whereas all other breeds had evidence of a hotspot. Finally, we identified relative enrichment (P = 5.88 × 10-27) for the in silico-predicted recognition motif for equine PR/SET domain 9 (PRDM9) in recombination hotspots. These results indicate that selective pressures and PRDM9 function contribute to variation in recombination rates across the domestic horse genome.
Collapse
Affiliation(s)
- Samantha K Beeson
- Veterinary Population Medicine Department, University of Minnesota, St. Paul, Minnesota 55108, USA
| | - James R Mickelson
- Veterinary and Biomedical Sciences Department, University of Minnesota, St. Paul, Minnesota 55108, USA
| | - Molly E McCue
- Veterinary Population Medicine Department, University of Minnesota, St. Paul, Minnesota 55108, USA
| |
Collapse
|
225
|
Luigi-Sierra MG, Cardoso TF, Martínez A, Pons A, Bermejo LA, Jordana J, Delgado JV, Adán S, Ugarte E, Arranz JJ, Casellas J, Amills M. Low genome-wide homozygosity in 11 Spanish ovine breeds. Anim Genet 2019; 50:501-511. [PMID: 31393638 DOI: 10.1111/age.12832] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/04/2019] [Indexed: 12/16/2022]
Abstract
The population of Spanish sheep has decreased from 24 to 15 million heads in the last 75 years due to multiple social and economic factors. Such a demographic reduction might have caused an increase in homozygosity and inbreeding, thus limiting the viability of local breeds with excellent adaptations to harsh ecosystems. The main goal of our study was to investigate the homozygosity patterns of 11 Spanish ovine breeds and to elucidate the relationship of these Spanish breeds with reference populations from Europe, Africa and the Near East. By using Ovine SNP50 BeadChip data retrieved from previous publications, we have found that the majority of studied Spanish ovine breeds have close genetic relatedness with other European populations; the one exception is the Canaria de Pelo breed, which is similar to North African breeds. Our analysis has also demonstrated that, with few exceptions, the genomes of Spanish sheep harbor fewer than 50 runs of homozygosity (ROH) with a total length of less than 350 Mb. Moreover, the frequencies of very long ROH (>30 Mb) are very low, and the inbreeding coefficients (FROH ) are generally small (FROH < 0.10), ranging from 0.008 (Rasa Aragonesa) to 0.086 (Canaria de Pelo). The low levels of homozygosity observed in the 11 Spanish sheep under analysis might be due to their extensive management and the high number of small to medium farms.
Collapse
Affiliation(s)
- M G Luigi-Sierra
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus Universitat Autònoma de Barcelona, Bellaterra, 08193, Spain
| | - T F Cardoso
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus Universitat Autònoma de Barcelona, Bellaterra, 08193, Spain.,CAPES Foundation, Ministry of Education of Brazil, Brasilia D.F, 70.040-020, Brazil
| | - A Martínez
- Departamento de Genética, Universidad de Córdoba, Córdoba, 14071, Spain
| | - A Pons
- Unitat de Races Autòctones, Servei de Millora Agrària i Pesquera (SEMILLA), Son Ferriol, 07198, Spain
| | - L A Bermejo
- Departamento de Ingeniería, Producción y Economía Agrarias, Universidad de La Laguna, La Laguna, Tenerife, 38071, Spain
| | - J Jordana
- Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, Bellaterra, 08193, Spain
| | - J V Delgado
- Departamento de Genética, Universidad de Córdoba, Córdoba, 14071, Spain
| | - S Adán
- Pazo de Fontefiz, Federación de Razas Autóctonas de Galicia (BOAGA), 32152, Coles, Ourense, Spain
| | - E Ugarte
- Neiker-Tecnalia, Campus Agroalimentario de Arkaute, apdo 46 E-01080, Vitoria-Gazteiz (Araba), Spain
| | - J J Arranz
- Departamento de Producción Animal, Universidad de León, León, 24071, Spain
| | - J Casellas
- Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, Bellaterra, 08193, Spain
| | - M Amills
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus Universitat Autònoma de Barcelona, Bellaterra, 08193, Spain.,Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, Bellaterra, 08193, Spain
| |
Collapse
|
226
|
Yuan Z, Li W, Li F, Yue X. Selection signature analysis reveals genes underlying sheep milking performance. Arch Anim Breed 2019; 62:501-508. [PMID: 31807661 PMCID: PMC6859915 DOI: 10.5194/aab-62-501-2019] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2019] [Accepted: 07/15/2019] [Indexed: 11/30/2022] Open
Abstract
Sheep milk is the most important feed resource for newborn lambs and an important food resource for humans. Sheep milk production and ingredients are influenced by genetic and environmental factors. In this study, we implemented selection signature analysis using Illumina Ovine SNP50 BeadChip data of 78 meat Lacaune and 103 milk Lacaune sheep, which have similar genetic backgrounds, from the Sheep HapMap project to identify candidate genes related to ovine milk traits. Since different methods can detect different variation types and complement each other, we used a haplotype-based method (hapFLK) to implement selection signature analysis. The results revealed six selection signature regions showing signs of being selected ( P < 0.001 ): chromosomes 1, 2, 3, 6, 13 and 18. In addition, 38 quantitative trait loci (QTLs) related to sheep milk performance were identified in selection signature regions, which contain 334 candidate genes. Of those, SUCNR1 (succinate receptor 1) and PPARGC1A (PPARG coactivator 1 alpha) may be the most significant genes that affect sheep milking performance, which supply a significant indication for future studies to investigate candidate genes that play an important role in milk production and quality.
Collapse
Affiliation(s)
- Zehu Yuan
- State Key Laboratory of Grassland Agro-Ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural
Affairs; College of Pastoral Agriculture Science and Technology, Lanzhou
University, Lanzhou, 730020, P. R. China
| | - Wanhong Li
- State Key Laboratory of Grassland Agro-Ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural
Affairs; College of Pastoral Agriculture Science and Technology, Lanzhou
University, Lanzhou, 730020, P. R. China
| | - Fadi Li
- State Key Laboratory of Grassland Agro-Ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural
Affairs; College of Pastoral Agriculture Science and Technology, Lanzhou
University, Lanzhou, 730020, P. R. China
- Engineering Laboratory of Sheep Breeding and Reproduction
Biotechnology in Gansu Province, Minqin, 733300, P. R. China
| | - Xiangpeng Yue
- State Key Laboratory of Grassland Agro-Ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural
Affairs; College of Pastoral Agriculture Science and Technology, Lanzhou
University, Lanzhou, 730020, P. R. China
| |
Collapse
|
227
|
Ben Jemaa S, Kdidi S, Gdura AM, Dayhum AS, Eldaghayes IM, Boussaha M, Rebours E, Yahyaoui MH. Inferring the population structure of the Maghreb sheep breeds using a medium-density SNP chip. Anim Genet 2019; 50:526-533. [PMID: 31385357 DOI: 10.1111/age.12831] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/04/2019] [Indexed: 01/26/2023]
Abstract
North Africa has a great diversity of indigenous sheep breeds whose origin is linked to its environmental characteristics and to certain historical events that took place in the region. To date, few genome-wide studies have been conducted to investigate the population structure of North African indigenous sheep. The objective of the present study was to provide a detailed assessment of the genetic structure and admixture patterns of six Maghreb sheep populations using the Illumina 50K Ovine BeadChip and comparisons with 22 global populations of sheep and mouflon. Regardless of the method of analysis used, patterns of multiple hybridization events were observed within all North African populations, leading to a heterogeneous genetic architecture that varies according to the breed. The Barbarine population showed the lowest genetic heterogeneity and major southwest Asian ancestry, providing additional support to the Asian origin of the North African fat-tailed sheep. All other breeds presented substantial Merino introgression ranging from 15% for D'man to 31% for Black Thibar. We highlighted several signals of ancestral introgression between North African and southern European sheep. In addition, we identified two opposite gradients of ancestry, southwest Asian and central European, occurring between North Africa and central Europe. Our results provide further evidence of the weak global population structure of sheep resulting from high levels of gene flow among breeds occurring worldwide. At the regional level, signs of recent admixture among North African populations, resulting in a change of the original genomic architecture of minority breeds, were also detected.
Collapse
Affiliation(s)
- S Ben Jemaa
- Laboratoire des Productions Animales et Fourragères, Institut National de la Recherche Agronomique de Tunisie, Université de Carthage, Ariana, Tunisia
| | - S Kdidi
- Livestock and Wildlife Laboratory, Arid Lands Institute, Route Djorf Km 22, 4119, Medenine, Tunisia
| | - A M Gdura
- Animal Improvement Project, Ministry of Agriculture, Tripoli, Libya
| | - A S Dayhum
- Faculty of Veterinary Medicine, University of Tripoli, 13662, Tripoli, Libya
| | - I M Eldaghayes
- Faculty of Veterinary Medicine, University of Tripoli, 13662, Tripoli, Libya
| | - M Boussaha
- GABI, INRA, AgroParisTech, Paris Saclay University, 78350, Jouy-en-Josas, France
| | - E Rebours
- GABI, INRA, AgroParisTech, Paris Saclay University, 78350, Jouy-en-Josas, France
| | - M H Yahyaoui
- Livestock and Wildlife Laboratory, Arid Lands Institute, Route Djorf Km 22, 4119, Medenine, Tunisia
| |
Collapse
|
228
|
Cesarani A, Sechi T, Gaspa G, Usai MG, Sorbolini S, Macciotta NPP, Carta A. Investigation of genetic diversity and selection signatures between Sarda and Sardinian Ancestral black, two related sheep breeds with evident morphological differences. Small Rumin Res 2019. [DOI: 10.1016/j.smallrumres.2019.06.014] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
|
229
|
Signer-Hasler H, Burren A, Ammann P, Drögemüller C, Flury C. Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds. Anim Genet 2019; 50:512-525. [PMID: 31365135 DOI: 10.1111/age.12828] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/13/2019] [Indexed: 01/01/2023]
Abstract
A dataset consisting of 787 animals with high-density SNP chip genotypes (346 774 SNPs) and 939 animals with medium-density SNP chip genotypes (33 828 SNPs) from eight indigenous Swiss sheep breeds was analyzed to characterize population structure, quantify genomic inbreeding based on runs of homozygosity and identify selection signatures. In concordance with the recent known history of these breeds, the highest genetic diversity was observed in Engadine Red sheep and the lowest in Valais Blacknose sheep. Correlation between FPED and FROH was around 0.50 and thereby lower than that found in similar studies in cattle. Mean FROH estimates from medium-density data and HD data were highly correlated (0.95). Signatures of selection and candidate gene analysis revealed that the most prominent signatures of selection were found in the proximity of genes associated with body size (NCAPG, LCORL, LAP3, SPP1, PLAG1, ALOX12, TP53), litter size (SPP1), milk production (ABCG2, SPP1), coat color (KIT, ASIP, TBX3) and horn status (RXFP2). For the Valais Blacknose sheep, the private signatures in proximity of genes/QTL influencing body size, coat color and fatty acid composition were confirmed based on runs of homozygosity analysis. These private signatures underline the genetic uniqueness of the Valais Blacknose sheep breed. In conclusion, we identified differences in the genetic make-up of Swiss sheep breeds, and we present relevant candidate genes responsible for breed differentiation in locally adapted breeds.
Collapse
Affiliation(s)
- H Signer-Hasler
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Länggasse 85, 3052 , Zollikofen, Switzerland
| | - A Burren
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Länggasse 85, 3052 , Zollikofen, Switzerland
| | - P Ammann
- ProSpecieRara, Unter Brüglingen 6, 4052, Basel, Switzerland
| | - C Drögemüller
- Vetsuisse Faculty, Institute of Genetics, University of Bern, Bremgartenstrasse 109a, 3001, Bern, Switzerland
| | - C Flury
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Länggasse 85, 3052 , Zollikofen, Switzerland
| |
Collapse
|
230
|
Attwood GT, Wakelin SA, Leahy SC, Rowe S, Clarke S, Chapman DF, Muirhead R, Jacobs JME. Applications of the Soil, Plant and Rumen Microbiomes in Pastoral Agriculture. Front Nutr 2019; 6:107. [PMID: 31380386 PMCID: PMC6646666 DOI: 10.3389/fnut.2019.00107] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2019] [Accepted: 06/27/2019] [Indexed: 12/14/2022] Open
Abstract
The production of dairy, meat, and fiber by ruminant animals relies on the biological processes occurring in soils, forage plants, and the animals' rumens. Each of these components has an associated microbiome, and these have traditionally been viewed as distinct ecosystems. However, these microbiomes operate under similar ecological principles and are connected via water, energy flows, and the carbon and nitrogen nutrient cycles. Here, we summarize the microbiome research that has been done in each of these three environments (soils, forage plants, animals' rumen) and investigate what additional benefits may be possible through understanding the interactions between the various microbiomes. The challenge for future research is to enhance microbiome function by appropriate matching of plant and animal genotypes with the environment to improve the output and environmental sustainability of pastoral agriculture.
Collapse
Affiliation(s)
| | | | - Sinead C Leahy
- Animal Science, AgResearch, Palmerston North, New Zealand
| | - Suzanne Rowe
- Animal Science, AgResearch, Invermay, New Zealand
| | | | | | | | | |
Collapse
|
231
|
Schönherz AA, Szekeres BD, Nielsen VH, Guldbrandtsen B. Population structure and genetic characterization of two native Danish sheep breeds. ACTA AGR SCAND A-AN 2019. [DOI: 10.1080/09064702.2019.1639804] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- A. A. Schönherz
- Department of Molecular Biology and Genetics, Aarhus University, Tjele, Denmark
| | - B. D. Szekeres
- Department of Molecular Biology and Genetics, Aarhus University, Tjele, Denmark
| | - V. H. Nielsen
- Danish Centre for Food and Agriculture, Aarhus University, Tjele, Denmark
| | - B. Guldbrandtsen
- Department of Molecular Biology and Genetics, Aarhus University, Tjele, Denmark
| |
Collapse
|
232
|
Al Kalaldeh M, Gibson J, Lee SH, Gondro C, van der Werf JHJ. Detection of genomic regions underlying resistance to gastrointestinal parasites in Australian sheep. Genet Sel Evol 2019; 51:37. [PMID: 31269896 PMCID: PMC6609385 DOI: 10.1186/s12711-019-0479-1] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2018] [Accepted: 06/19/2019] [Indexed: 12/13/2022] Open
Abstract
BACKGROUND This study aimed at identifying genomic regions that underlie genetic variation of worm egg count, as an indicator trait for parasite resistance in a large population of Australian sheep, which was genotyped with the high-density 600 K Ovine single nucleotide polymorphism array. This study included 7539 sheep from different locations across Australia that underwent a field challenge with mixed gastrointestinal parasite species. Faecal samples were collected and worm egg counts for three strongyle species, i.e. Teladorsagia circumcincta, Haemonchus contortus and Trichostrongylus colubriformis were determined. Data were analysed using genome-wide association studies (GWAS) and regional heritability mapping (RHM). RESULTS Both RHM and GWAS detected a region on Ovis aries (OAR) chromosome 2 that was highly significantly associated with parasite resistance at a genome-wise false discovery rate of 5%. RHM revealed additional significant regions on OAR6, 18, and 24. Pathway analysis revealed 13 genes within these significant regions (SH3RF1, HERC2, MAP3K, CYFIP1, PTPN1, BIN1, HERC3, HERC5, HERC6, IBSP, SPP1, ISG20, and DET1), which have various roles in innate and acquired immune response mechanisms, as well as cytokine signalling. Other genes involved in haemostasis regulation and mucosal defence were also detected, which are important for protection of sheep against invading parasites. CONCLUSIONS This study identified significant genomic regions on OAR2, 6, 18, and 24 that are associated with parasite resistance in sheep. RHM was more powerful in detecting regions that affect parasite resistance than GWAS. Our results support the hypothesis that parasite resistance is a complex trait and is determined by a large number of genes with small effects, rather than by a few major genes with large effects.
Collapse
Affiliation(s)
- Mohammad Al Kalaldeh
- Cooperative Research Centre for Sheep Industry Innovation, Armidale, NSW, 2351, Australia. .,School of Environmental and Rural Science, University of New England, Armidale, NSW, 2351, Australia.
| | - John Gibson
- Cooperative Research Centre for Sheep Industry Innovation, Armidale, NSW, 2351, Australia.,School of Environmental and Rural Science, University of New England, Armidale, NSW, 2351, Australia
| | - Sang Hong Lee
- Australian Centre for Precision Health, University of South Australia Cancer Research Institute, University of South Australia, Adelaide, SA, 5000, Australia
| | - Cedric Gondro
- School of Environmental and Rural Science, University of New England, Armidale, NSW, 2351, Australia.,College of Agriculture and Natural Resources, Michigan State University, East Lansing, MI, 48824, USA
| | - Julius H J van der Werf
- Cooperative Research Centre for Sheep Industry Innovation, Armidale, NSW, 2351, Australia.,School of Environmental and Rural Science, University of New England, Armidale, NSW, 2351, Australia
| |
Collapse
|
233
|
Hu XJ, Yang J, Xie XL, Lv FH, Cao YH, Li WR, Liu MJ, Wang YT, Li JQ, Liu YG, Ren YL, Shen ZQ, Wang F, Hehua EE, Han JL, Li MH. The Genome Landscape of Tibetan Sheep Reveals Adaptive Introgression from Argali and the History of Early Human Settlements on the Qinghai-Tibetan Plateau. Mol Biol Evol 2019; 36:283-303. [PMID: 30445533 PMCID: PMC6367989 DOI: 10.1093/molbev/msy208] [Citation(s) in RCA: 65] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Tibetan sheep are the most common and widespread domesticated animals on the Qinghai-Tibetan Plateau (QTP) and have played an essential role in the permanent human occupation of this high-altitude region. However, the precise timing, route, and process of sheep pastoralism in the QTP region remain poorly established, and little is known about the underlying genomic changes that occurred during the process. Here, we investigate the genomic variation in Tibetan sheep using whole-genome sequences, single nucleotide polymorphism arrays, mitochondrial DNA, and Y-chromosomal variants in 986 samples throughout their distribution range. We detect strong signatures of selection in genes involved in the hypoxia and ultraviolet signaling pathways (e.g., HIF-1 pathway and HBB and MITF genes) and in genes associated with morphological traits such as horn size and shape (e.g., RXFP2). We identify clear signals of argali (Ovis ammon) introgression into sympatric Tibetan sheep, covering 5.23-5.79% of their genomes. The introgressed genomic regions are enriched in genes related to oxygen transportation system, sensory perception, and morphological phenotypes, in particular the genes HBB and RXFP2 with strong signs of adaptive introgression. The spatial distribution of genomic diversity and demographic reconstruction of the history of Tibetan sheep show a stepwise pattern of colonization with their initial spread onto the QTP from its northeastern part ∼3,100 years ago, followed by further southwest expansion to the central QTP ∼1,300 years ago. Together with archeological evidence, the date and route reveal the history of human expansions on the QTP by the Tang-Bo Ancient Road during the late Holocene. Our findings contribute to a depth understanding of early pastoralism and the local adaptation of Tibetan sheep as well as the late-Holocene human occupation of the QTP.
Collapse
Affiliation(s)
- Xiao-Ju Hu
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences (CAS), Beijing, China.,College of Life Sciences, University of Chinese Academy of Sciences (UCAS), Beijing, China
| | - Ji Yang
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences (CAS), Beijing, China
| | - Xing-Long Xie
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences (CAS), Beijing, China.,College of Life Sciences, University of Chinese Academy of Sciences (UCAS), Beijing, China
| | - Feng-Hua Lv
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences (CAS), Beijing, China
| | - Yin-Hong Cao
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences (CAS), Beijing, China.,College of Life Sciences, University of Chinese Academy of Sciences (UCAS), Beijing, China
| | - Wen-Rong Li
- Animal Biotechnological Research Center, Xinjiang Academy of Animal Science, Urumqi, China
| | - Ming-Jun Liu
- Animal Biotechnological Research Center, Xinjiang Academy of Animal Science, Urumqi, China
| | - Yu-Tao Wang
- College of Life and Geographic Sciences, Kashi University, Kashi, China
| | - Jin-Quan Li
- College of Animal Science, Inner Mongolia Agricultural University, Hohhot, China
| | - Yong-Gang Liu
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China
| | - Yan-Lin Ren
- Shandong Binzhou Academy of Animal Science and Veterinary Medicine, Binzhou, China
| | - Zhi-Qiang Shen
- Shandong Binzhou Academy of Animal Science and Veterinary Medicine, Binzhou, China
| | - Feng Wang
- Institute of Sheep and Goat Science, Nanjing Agricultural University, Nanjing, China
| | - EEr Hehua
- Grass-Feeding Livestock Engineering Technology Research Center, Ningxia Academy of Agriculture and Forestry Sciences, Yinchuan, China
| | - Jian-Lin Han
- CAAS-ILRI Joint Laboratory on Livestock and Forage Genetic Resources, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, China.,Livestock Genetics Program, International Livestock Research Institute (ILRI), Nairobi, Kenya
| | - Meng-Hua Li
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences (CAS), Beijing, China.,College of Life Sciences, University of Chinese Academy of Sciences (UCAS), Beijing, China
| |
Collapse
|
234
|
Mastrangelo S, Ben Jemaa S, Sottile G, Casu S, Portolano B, Ciani E, Pilla F. Combined approaches to identify genomic regions involved in phenotypic differentiation between low divergent breeds: Application in Sardinian sheep populations. J Anim Breed Genet 2019; 136:526-534. [PMID: 31206848 DOI: 10.1111/jbg.12422] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2019] [Revised: 05/17/2019] [Accepted: 05/23/2019] [Indexed: 12/17/2022]
Abstract
Selective breeding has led to modifications in the genome of many livestock breeds. In this study, we identified the genomic regions that may explain some of the phenotypic differences between two closely related breeds from Sardinia. A total of 44 animals, 20 Sardinian Ancestral Black (SAB) and 24 Sardinian White (SW), were genotyped using the Illumina Ovine 50K array. A total of 68, 38 and 15 significant markers were identified using the case-control genome-wide association study (GWAS), the Bayesian population differentiation analysis (FST ) and the Rsb metric, respectively. Comparisons among the approaches revealed a total of 22 overlapping markers between GWAS and FST and one marker between GWAS and Rsb. Three markers detected by Rsb were also located near (<2 Mb) to highly significant regions identified by GWAS and FST analyses. Moreover, one candidate marker identified by GWAS and FST approaches was located in a run of homozygosity island that was shared by both breeds. We identified several genes involved in many phenotypic differences (such as stature and growth, reproduction, ear size, coat colour, behaviour) between the two analysed breeds. This study shows that combining several genome-wide approaches could improve discovery of regions involved in the variability of breeding traits and responsible for the phenotypic diversity even between closely related breeds. Overall, the combination of such genome-wide methods can be extended to other livestock breeds that share between them a similar genetic background, to understand the process that shapes the patterns of genetic variability between closely related populations.
Collapse
Affiliation(s)
- Salvatore Mastrangelo
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Slim Ben Jemaa
- Laboratoire des Productions Animales et Fourragères, Institut National de la Recherche Agronomique de Tunisie, Université de Carthage, Ariana, Tunisia
| | - Gianluca Sottile
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Sara Casu
- Unità di Ricerca di Genetica e Biotecnologie, Agris Sardegna, Sassari, Italy
| | - Baldassare Portolano
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Elena Ciani
- Dipartimento di Bioscienze Biotecnologie e Biofarmaceutica, University of Bari, Bari, Italy
| | - Fabio Pilla
- Dipartimento di Agricoltura, Ambiente e Alimenti, University of Molise, Campobasso, Italy
| |
Collapse
|
235
|
Mastrangelo S, Bahbahani H, Moioli B, Ahbara A, Al Abri M, Almathen F, da Silva A, Belabdi I, Portolano B, Mwacharo JM, Hanotte O, Pilla F, Ciani E. Novel and known signals of selection for fat deposition in domestic sheep breeds from Africa and Eurasia. PLoS One 2019; 14:e0209632. [PMID: 31199810 PMCID: PMC6568386 DOI: 10.1371/journal.pone.0209632] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2018] [Accepted: 05/28/2019] [Indexed: 01/24/2023] Open
Abstract
Genomic regions subjected to selection frequently show signatures such as within-population reduced nucleotide diversity and outlier values of differentiation among differentially selected populations. In this study, we analyzed 50K SNP genotype data of 373 animals belonging to 23 sheep breeds of different geographic origins using the Rsb (extended haplotype homozygosity) and FST statistical approaches, to identify loci associated with the fat-tail phenotype. We also checked if these putative selection signatures overlapped with regions of high-homozygosity (ROH). The analyses identified novel signals and confirmed the presence of selection signature in genomic regions that harbor candidate genes known to affect fat deposition. Several genomic regions that frequently appeared in ROH were also identified within each breed, but only two ROH islands overlapped with the putative selection signatures. The results reported herein provide the most complete genome-wide study of selection signatures for fat-tail in African and Eurasian sheep breeds; they also contribute insights into the genetic basis for the fat tail phenotype in sheep, and confirm the great complexity of the mechanisms that underlie quantitative traits, such as the fat-tail.
Collapse
Affiliation(s)
- Salvatore Mastrangelo
- Dipartimento di Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Hussain Bahbahani
- Department of Biological Sciences, Faculty of Science, Kuwait University, Safat, Kuwait
| | - Bianca Moioli
- Consiglio per la ricerca in agricoltura e l’analisi dell’economia agraria (CREA), Monterotondo, Italy
| | - Abulgasim Ahbara
- School of Life Sciences, University of Nottingham, University Park, Nottingham, United Kingdom
- Department of Zoology, Faculty of Sciences, Misurata University, Misurata, Libya
| | - Mohammed Al Abri
- Department of Animal and Veterinary Sciences, Sultan Qaboos University, Oman
| | - Faisal Almathen
- Department of Public Health and Animal Welfare, College of Veterinary Medicine, King Faisal University, Alhufuf, Al-Ahsa, Saudi Arabia
| | - Anne da Silva
- Université de Limoges, INRA, PEREINE EA7500, USC1061 GAMAA, Limoges, France
| | - Ibrahim Belabdi
- Science Veterinary Institute, University of Blida, Blida, Algeria
- Laboratory of Biotechnology related to Animal Reproduction (LBRA), University of Blida, Blida, Algeria
| | - Baldassare Portolano
- Dipartimento di Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Joram M. Mwacharo
- Small Ruminant Genomics, International Center for Agricultural Research in the Dry Areas (ICARDA), Addis Ababa, Ethiopia
| | - Olivier Hanotte
- School of Life Sciences, University of Nottingham, University Park, Nottingham, United Kingdom
| | - Fabio Pilla
- Dipartimento Agricoltura, Ambiente e Alimenti, Università degli Studi del Molise, Campobasso, Italy
| | - Elena Ciani
- Dipartimento di Bioscienze, Biotecnologie e Biofarmaceutica, Università degli Studi di Bari Aldo Moro, Bari, Italy
- * E-mail:
| |
Collapse
|
236
|
Genome-wide scan reveals genetic divergence and diverse adaptive selection in Chinese local cattle. BMC Genomics 2019; 20:494. [PMID: 31200634 PMCID: PMC6570941 DOI: 10.1186/s12864-019-5822-y] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2018] [Accepted: 05/21/2019] [Indexed: 01/18/2023] Open
Abstract
Background Understanding the population structure and genetic bases of well-adapted cattle breeds to local environments is one of the most essential tasks to develop appropriate genetic improvement programs. Results We performed a comprehensive study to investigate the population structure, divergence and selection signatures at genome-wide level in diverse Chinese local cattle using Bovine HD SNPs array, including two breeds from North China, one breed from Northwest China, three breeds from Southwest China and two breeds from South China. Population genetic analyses revealed the genetic structures of these populations were mostly related to the geographic locations. Notably, we detected 294 and 1263 candidate regions under selection using the di and iHS approaches, respectively. A series of group-specific and breed-specific candidate genes were identified, which are involved in immune response, sexual maturation, stature related, birth and bone weight, embryonic development, coat colors and adaptation. Furthermore, haplotype diversity and network pattern for candidate genes, including LPGAT1, LCORL, PPP1R8, RXFP2 and FANCA, suggest that these genes have been under differential selection pressures in various environmental conditions. Conclusions Our results shed insights into diverse selection during breed formation in Chinese local cattle. These findings may promote the application of genome-assisted breeding for well-adapted local breeds with economic and ecological importance. Electronic supplementary material The online version of this article (10.1186/s12864-019-5822-y) contains supplementary material, which is available to authorized users.
Collapse
|
237
|
Genetic homogenization of indigenous sheep breeds in Northwest Africa. Sci Rep 2019; 9:7920. [PMID: 31138837 PMCID: PMC6538629 DOI: 10.1038/s41598-019-44137-y] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2018] [Accepted: 05/07/2019] [Indexed: 11/17/2022] Open
Abstract
Northwest-African sheep represent an ideal case-study for assessing the potential impact of genetic homogenization as a threat to the future of traditional breeds that are adapted to local conditions. We studied ten Algerian and Moroccan breeds of sheep, including three transboundary breeds, distributed over a large part of the Maghreb region, which represents a geographically and historically coherent unit. Our analysis of the dataset that involved carrying out Genome-wide SNP genotyping, revealed a high level of homogenization (ADMIXTURE, NetView, fineSTRUCTURE and IBD segments analyses), in such a way that some breeds from different origins appeared genetically undistinguished: by grouping the eight most admixed populations, we obtained a mean global FST value of 0.0024. The sPCA analysis revealed that the major part of Morocco and the Northern part of Algeria were affected by the phenomenon, including most of the breeds considered. Unsupervised cross-breeding with the popular Ouled-Djellal breed was identified as a proximate cause of this homogenization. The issue of transboundary breeds was investigated, and the Hamra breed in particular was examined via ROH fragments analysis. Genetic diversity was considered in the light of historical archives and anthropological works. All of these elements taken together suggest that homogenization as a factor affecting the Maghrebin sheep stock, has been particularly significant over the last few decades, although this process probably started much earlier. In particular, we have identified the policies set by the French administration during the colonial period of the region’s history as a causal factor that probably contributed significantly to this process. The genetic homogenization that we have observed calls into question the integrity of the farm animal genomic resources represented by these local breeds, whose conservation is of critical importance to the future of the livestock sector.
Collapse
|
238
|
Yurchenko AA, Deniskova TE, Yudin NS, Dotsev AV, Khamiruev TN, Selionova MI, Egorov SV, Reyer H, Wimmers K, Brem G, Zinovieva NA, Larkin DM. High-density genotyping reveals signatures of selection related to acclimation and economically important traits in 15 local sheep breeds from Russia. BMC Genomics 2019; 20:294. [PMID: 32039702 PMCID: PMC7227232 DOI: 10.1186/s12864-019-5537-0] [Citation(s) in RCA: 46] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Background Domestication and centuries of selective breeding have changed genomes of sheep breeds to respond to environmental challenges and human needs. The genomes of local breeds, therefore, are valuable sources of genomic variants to be used to understand mechanisms of response to adaptation and artificial selection. As a step toward this we performed a high-density genotyping and comprehensive scans for signatures of selection in the genomes from 15 local sheep breeds reared across Russia. Results Results demonstrated that the genomes of Russian sheep breeds contain multiple regions under putative selection. More than 50% of these regions matched with intervals identified in previous scans for selective sweeps in sheep genomes. These regions contain well-known candidate genes related to morphology, adaptation, and domestication (e.g., KITLG, KIT, MITF, and MC1R), wool quality and quantity (e.g., DSG@, DSC@, and KRT@), growth and feed intake (e.g., HOXA@, HOXC@, LCORL, NCAPG, LAP3, and CCSER1), reproduction (e.g., CMTM6, HTRA1, GNAQ, UBQLN1, and IFT88), and milk-related traits (e.g., ABCG2, SPP1, ACSS1, and ACSS2). In addition, multiple genes that are putatively related to environmental adaptations were top-ranked in selected intervals (e.g., EGFR, HSPH1, NMUR1, EDNRB, PRL, TSHR, and ADAMTS5). Moreover, we observed that multiple key genes involved in human hereditary sensory and autonomic neuropathies, and genetic disorders accompanied with an inability to feel pain and environmental temperatures, were top-ranked in multiple or individual sheep breeds from Russia pointing to a possible mechanism of adaptation to harsh climatic conditions. Conclusions Our work represents the first comprehensive scan for signatures of selection in genomes of local sheep breeds from the Russian Federation of both European and Asian origins. We confirmed that the genomes of Russian sheep contain previously identified signatures of selection, demonstrating the robustness of our integrative approach. Multiple novel signatures of selection were found near genes which could be related to adaptation to the harsh environments of Russia. Our study forms a basis for future work on using Russian sheep genomes to spot specific genetic variants or haplotypes to be used in efforts on developing next-generation highly productive breeds, better suited to diverse Eurasian environments. Electronic supplementary material The online version of this article (10.1186/s12864-019-5537-0) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Andrey A Yurchenko
- The Federal Research Center Institute of Cytology and Genetics, The Siberian Branch of the Russian Academy of Sciences (ICG SB RAS), Novosibirsk, Russia
| | - Tatiana E Deniskova
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, 142132, Russia
| | - Nikolay S Yudin
- The Federal Research Center Institute of Cytology and Genetics, The Siberian Branch of the Russian Academy of Sciences (ICG SB RAS), Novosibirsk, Russia.,Novosibirsk State University, Novosibirsk, 630090, Russia
| | - Arsen V Dotsev
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, 142132, Russia
| | - Timur N Khamiruev
- Research Institute of Veterinary Medicine of Eastern Siberia, The Branch of the Siberian Federal Scientific Center for Agrobiotechnologies of the Russian Academy of Sciences, Chita, Russia
| | - Marina I Selionova
- All-Russian Research Institute of Sheep and Goat Breeding - branch of the Federal State Budgetary Scientific Institution North Caucasian Agrarian Center, Stavropol, 355017, Russia
| | - Sergey V Egorov
- Siberian Research Institute of Animal Husbandry, Krasnoobsk, Russia
| | - Henry Reyer
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
| | - Klaus Wimmers
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
| | - Gottfried Brem
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, 142132, Russia.,Institute of Animal Breeding and Genetics, University of Veterinary Medicine, Vienna, Austria
| | - Natalia A Zinovieva
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, 142132, Russia.
| | - Denis M Larkin
- The Federal Research Center Institute of Cytology and Genetics, The Siberian Branch of the Russian Academy of Sciences (ICG SB RAS), Novosibirsk, Russia. .,Royal Veterinary College, University of London, London, UK.
| |
Collapse
|
239
|
Selection signatures in four German warmblood horse breeds: Tracing breeding history in the modern sport horse. PLoS One 2019; 14:e0215913. [PMID: 31022261 PMCID: PMC6483353 DOI: 10.1371/journal.pone.0215913] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2018] [Accepted: 04/10/2019] [Indexed: 12/04/2022] Open
Abstract
The study of selection signatures helps to find genomic regions that have been under selective pressure and might host genes or variants that modulate important phenotypes. Such knowledge improves our understanding of how breeding programmes have shaped the genomes of livestock. In this study, 942 stallions were included from four, exemplarily chosen, German warmblood breeds with divergent historical and recent selection focus and different crossbreeding policies: Trakehner (N = 44), Holsteiner (N = 358), Hanoverian (N = 319) and Oldenburger (N = 221). Those breeds are nowadays bred for athletic performance and aptitude for show-jumping, dressage or eventing, with a particular focus of Holsteiner on the first discipline. Blood samples were collected during the health exams of the stallion preselections before licensing and were genotyped with the Illumina EquineSNP50 BeadChip. Autosomal markers were used for a multi-method search for signals of positive selection. Analyses within and across breeds were conducted by using the integrated Haplotype Score (iHS), cross-population Extended Haplotype Homozygosity (xpEHH) and Runs of Homozygosity (ROH). Oldenburger and Hanoverian showed very similar iHS signatures, but breed specificities were detected on multiple chromosomes with the xpEHH. The Trakehner clustered as a distinct group in a principal component analysis and also showed the highest number of ROHs, which reflects their historical bottleneck. Beside breed specific differences, we found shared selection signals in an across breed iHS analysis on chromosomes 1, 4 and 7. After investigation of these iHS signals and shared ROH for potential functional candidate genes and affected pathways including enrichment analyses, we suggest that genes affecting muscle functionality (TPM1, TMOD2-3, MYO5A, MYO5C), energy metabolism and growth (AEBP1, RALGAPA2, IGFBP1, IGFBP3-4), embryonic development (HOXB-complex) and fertility (THEGL, ZPBP1-2, TEX14, ZP1, SUN3 and CFAP61) have been targeted by selection in all breeds. Our findings also indicate selection pressure on KITLG, which is well-documented for influencing pigmentation.
Collapse
|
240
|
Chhotaray S, Panigrahi M, Pal D, Ahmad SF, Bhanuprakash V, Kumar H, Parida S, Bhushan B, Gaur GK, Mishra BP, Singh RK. Genome-wide estimation of inbreeding coefficient, effective population size and haplotype blocks in Vrindavani crossbred cattle strain of India. BIOL RHYTHM RES 2019. [DOI: 10.1080/09291016.2019.1600266] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Supriya Chhotaray
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, India
| | - Manjit Panigrahi
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, India
| | - Dhan Pal
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, India
| | - Sheikh Firdous Ahmad
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, India
| | - V. Bhanuprakash
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, India
| | - Harshit Kumar
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, India
| | - Subhashree Parida
- Division of Pharmacology & Toxicology, Indian Veterinary Research Institute, Bareilly, India
| | - Bharat Bhushan
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, India
| | - G. K. Gaur
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, India
| | - B. P. Mishra
- Division of Animal Biotechnology, Indian Veterinary Research Institute, Bareilly, India
| | - R. K. Singh
- Division of Animal Biotechnology, Indian Veterinary Research Institute, Bareilly, India
| |
Collapse
|
241
|
Wang W, Zhang X, Zhou X, Zhang Y, La Y, Zhang Y, Li C, Zhao Y, Li F, Liu B, Jiang Z. Deep Genome Resequencing Reveals Artificial and Natural Selection for Visual Deterioration, Plateau Adaptability and High Prolificacy in Chinese Domestic Sheep. Front Genet 2019; 10:300. [PMID: 31001329 PMCID: PMC6454055 DOI: 10.3389/fgene.2019.00300] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2018] [Accepted: 03/19/2019] [Indexed: 01/02/2023] Open
Abstract
Sheep were one of the earliest domesticated animals. Both artificial and natural selection during domestication has resulted in remarkable changes in behavioral, physiological, and morphological phenotypes; however, the genetic mechanisms underpinning these changes remain unclear, particularly for indigenous Chinese sheep. In the present study, we performed pooled whole-genome resequencing of 338 sheep from five breeds representative of indigenous Chinese breeds and compared them to the wild ancestors of domestic sheep (Asian mouflon, Ovis orientalis) for detection of genome-wide selective sweeps. Comparative genomic analysis between domestic sheep and Asian mouflon showed that selected regions were enriched for genes involved in bone morphogenesis, growth regulation, and embryonic and neural development in domestic sheep. Moreover, we identified several vision-associated genes with funtional mutations, such as PDE6B (c.G2994C/p.A982P and c.C2284A/p.L762M mutations), PANK2, and FOXC1/GMSD in all five Chinese native breeds. Breed-specific selected regions were determined including genes such as CYP17 for hypoxia adaptability in Tibetan sheep and DNAJB5 for heat tolerance in Duolang sheep. Our findings provide insights into the genetic mechanisms underlying important phenotypic changes that have occurred during sheep domestication and subsequent selection.
Collapse
Affiliation(s)
- Weimin Wang
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, China
| | - Xiaoxue Zhang
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, China
| | - Xiang Zhou
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan, China.,Department of Animal Sciences, Washington State University, Pullman, WA, United States
| | - Yangzi Zhang
- Department of Animal Sciences, Washington State University, Pullman, WA, United States
| | - Yongfu La
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, China
| | - Yu Zhang
- Department of Animal Sciences, Washington State University, Pullman, WA, United States
| | - Chong Li
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, China
| | - Youzhang Zhao
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, China
| | - Fadi Li
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, China.,The State Key Laboratory of Grassland Agro-Ecosystems, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, China.,Engineering Laboratory of Sheep Breeding and Reproduction Biotechnology in Gansu Province, Minqin, China
| | - Bang Liu
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan, China
| | - Zhihua Jiang
- Department of Animal Sciences, Washington State University, Pullman, WA, United States
| |
Collapse
|
242
|
Ghoreishifar SM, Moradi-Shahrbabak H, Parna N, Davoudi P, Khansefid M. Linkage disequilibrium and within-breed genetic diversity in Iranian Zandi sheep. Arch Anim Breed 2019; 62:143-151. [PMID: 31807624 PMCID: PMC6852851 DOI: 10.5194/aab-62-143-2019] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2018] [Accepted: 03/08/2019] [Indexed: 12/23/2022] Open
Abstract
This research aimed to measure the extent of linkage disequilibrium (LD),
effective population size (Ne), and runs of homozygosity (ROHs)
in one of the major Iranian sheep breeds (Zandi) using 96 samples genotyped
with Illumina Ovine SNP50 BeadChip. The amount of LD (r2) for
single-nucleotide polymorphism (SNP) pairs in short distances (10–20 kb)
was 0.21±0.25 but rapidly decreased to 0.10±0.16 by increasing the
distance between SNP pairs (40–60 kb). The Ne of Zandi sheep in
past (approximately 3500 generations ago) and recent (five generations ago)
populations was estimated to be 6475 and 122, respectively. The ROH-based
inbreeding was 0.023. We found 558 ROH regions, of which 37 % were
relatively long (>10 Mb). Compared with the rate of LD
reduction in other species (e.g., cattle and pigs), in Zandi, it was reduced
more rapidly by increasing the distance between SNP pairs. According to the
LD pattern and high genetic diversity of Zandi sheep, we need to use an SNP
panel with a higher density than Illumina Ovine SNP50 BeadChip for genomic
selection and genome-wide association studies in this breed.
Collapse
Affiliation(s)
- Seyed Mohammad Ghoreishifar
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Hossein Moradi-Shahrbabak
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Nahid Parna
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Pourya Davoudi
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Majid Khansefid
- Agriculture Victoria, AgriBio Centre for AgriBioscience, Bundoora, VIC 3083, Australia
| |
Collapse
|
243
|
Sejian V, Bagath M, Krishnan G, Rashamol V, Pragna P, Devaraj C, Bhatta R. Genes for resilience to heat stress in small ruminants: A review. Small Rumin Res 2019. [DOI: 10.1016/j.smallrumres.2019.02.009] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
|
244
|
Manzari Z, Mehrabani-Yeganeh H, Nejati-Javaremi A, Moradi MH, Gholizadeh M. Detecting selection signatures in three Iranian sheep breeds. Anim Genet 2019; 50:298-302. [PMID: 30883840 DOI: 10.1111/age.12772] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/14/2018] [Indexed: 12/26/2022]
Abstract
The objective of genome mapping is to achieve valuable insight into the connection between gene variants (genotype) and observed traits (phenotype). Part of that objective is to understand the selective forces that have operated on a population. Finding links between genotype-phenotype changes makes it possible to identify selective sweeps by patterns of genetic variation and linkage disequilibrium. Based on Illumina 50KSNP chip data, two approaches, XP-EHH (cross-population extend haplotype homozygosity) and FST (fixation index), were carried out in this research to identify selective sweeps in the genome of three Iranian local sheep breeds: Baluchi (n = 86), Lori-Bakhtiari (n = 45) and Zel (n = 45). Using both methods, 93 candidate genomic regions were identified as harboring putative selective sweeps. Bioinformatics analysis of the genomic regions showed that signatures of selection related to multiple candidate genes, such as HOXB9, HOXB13, ACAN, NPR2, TRIL, AOX1, CSF2, GHR, TNS2, SPAG8, HINT2, ALS2, AAAS, RARG, SYCP2, CAV1, PPP1R3D, PLA2G7, TTLL7 and C20orf10, that play a role in skeletal system and tail, sugar and energy metabolisms, growth, reproduction, immune and nervous system traits. Our findings indicated diverse genomic selection during the domestication of Iranian sheep breeds.
Collapse
Affiliation(s)
- Z Manzari
- Department of Animal Science, College of Agriculture and Natural Resources, University of Tehran, Karaj, 3158711167-4111, Iran
| | - H Mehrabani-Yeganeh
- Department of Animal Science, College of Agriculture and Natural Resources, University of Tehran, Karaj, 3158711167-4111, Iran
| | - A Nejati-Javaremi
- Department of Animal Science, College of Agriculture and Natural Resources, University of Tehran, Karaj, 3158711167-4111, Iran
| | - M H Moradi
- Department of Animal Science, Faculty of Agriculture and Natural Resources, Arak University, Arak, 38156-8-8349, Iran
| | - M Gholizadeh
- Department of Animal Science, Faculty of Animal and Aquatic Science, Sari Agricultural Sciences and Natural Resources University, Sari, 4818168984, Iran
| |
Collapse
|
245
|
Berihulay H, Li Y, Liu X, Gebreselassie G, Islam R, Liu W, Jiang L, Ma Y. Genetic diversity and population structure in multiple Chinese goat populations using a SNP panel. Anim Genet 2019; 50:242-249. [PMID: 30883837 DOI: 10.1111/age.12776] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 01/16/2019] [Indexed: 11/30/2022]
Abstract
Information about genetic diversity and population structure among goat breeds is essential for genetic improvement, understanding of environmental adaptation as well as utilization and conservation of goat breeds. Here, we measured genetic diversity and population structure in multiple Chinese goat populations, namely, Nanjiang, Qinggeda, Arbas Cashmere, Jining Grey, Luoping Yellow and Guangfeng goats. A total of 193 individuals were genotyped for about 47 401 autosomal single nucleotide polymorphisms (SNPs). We found a high proportion of informative SNPs, ranging from 69.5% in the Luoping Yellow to 93.9% in the Jining Grey goat breeds with an average mean of 84.7%. Diversity, as measured by expected heterozygosity, ranged from 0.371 in Luoping Yellow to 0.405 in Jining Grey goat populations. The average estimated pair-wise genetic differentiation (FST ) among the populations was 8.6%, ranging from 0.2% to 16% and indicating low to moderate genetic differentiation. Principal component analysis, genetic structure and phylogenetic tree analysis revealed a clustering of six Chinese goat populations according to geographic distribution. The results from this study can contribute valuable genetic information and can properly assist with within-breed diversity, which provides a good opportunity for sustainable utilization of and maintenance of genetic resource improvements in the Chinese goat populations.
Collapse
Affiliation(s)
- H Berihulay
- Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China.,The Key Laboratory for Farm Animal Genetic Resources and Utilization of Ministry of Agriculture of China, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China
| | - Y Li
- Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China.,The Key Laboratory for Farm Animal Genetic Resources and Utilization of Ministry of Agriculture of China, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China
| | - X Liu
- Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China.,The Key Laboratory for Farm Animal Genetic Resources and Utilization of Ministry of Agriculture of China, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China
| | - G Gebreselassie
- Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China.,The Key Laboratory for Farm Animal Genetic Resources and Utilization of Ministry of Agriculture of China, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China
| | - R Islam
- Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China.,The Key Laboratory for Farm Animal Genetic Resources and Utilization of Ministry of Agriculture of China, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China
| | - W Liu
- Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China.,The Key Laboratory for Farm Animal Genetic Resources and Utilization of Ministry of Agriculture of China, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China
| | - L Jiang
- Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China.,The Key Laboratory for Farm Animal Genetic Resources and Utilization of Ministry of Agriculture of China, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China
| | - Y Ma
- Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China.,The Key Laboratory for Farm Animal Genetic Resources and Utilization of Ministry of Agriculture of China, Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China
| |
Collapse
|
246
|
Bertrand AR, Kadri NK, Flori L, Gautier M, Druet T. RZooRoH: An R package to characterize individual genomic autozygosity and identify homozygous‐by‐descent segments. Methods Ecol Evol 2019. [DOI: 10.1111/2041-210x.13167] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Amandine R. Bertrand
- Unit of Animal Genomics, GIGA‐R & Faculty of Veterinary MedicineUniversity of Liège Liège Belgium
| | - Naveen K. Kadri
- Unit of Animal Genomics, GIGA‐R & Faculty of Veterinary MedicineUniversity of Liège Liège Belgium
| | - Laurence Flori
- SELMET, INRA, CIRAD, Montpellier SupagroUniversity of Montpellier Montpellier France
| | - Mathieu Gautier
- INRA, UMR CBGP (INRA – IRD – Cirad – Montpellier SupAgro) Montferrier‐sur‐Lez France
| | - Tom Druet
- Unit of Animal Genomics, GIGA‐R & Faculty of Veterinary MedicineUniversity of Liège Liège Belgium
| |
Collapse
|
247
|
Rodríguez-Ramilo ST, Elsen JM, Legarra A. Inbreeding and effective population size in French dairy sheep: Comparison between genomic and pedigree estimates. J Dairy Sci 2019; 102:4227-4237. [PMID: 30827541 DOI: 10.3168/jds.2018-15405] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2018] [Accepted: 12/23/2018] [Indexed: 01/11/2023]
Abstract
Before availability of dense SNP data, genetic diversity was characterized and managed with pedigree-based information. Besides this classical approach, 2 methodologies have been proposed in recent years to characterize and manage diversity from dense SNP data: the SNP-by-SNP approach and the alternative based on runs of homozygosity (ROH). The establishment of criteria to identify ROH is a current constraint in the literature dealing with ROH. The objective of this study was, using a medium-density SNP chip, to quantify by 3 methods (pedigree, SNP-by-SNP, and ROH) the genetic diversity on 5 selected French dairy sheep subpopulations and breeds and to assess the effect of the definition of ROH on these estimates. The data set available included individuals from the breeds Basco-Béarnaise, Manech Tête Noire, Manech Tête Rousse, and 2 subpopulations of Lacaune: Lacaune Confederation and Lacaune Ovitest. Animals were genotyped with the Illumina OvineSNP50 BeadChip (Illumina Inc., San Diego, CA). After filtering, the genomic data included 38,287 autosomal SNP and 8,700 individuals, which comprised 72,803 animals in the pedigree. The results indicated that no significant differences were observed in effective population size estimates obtained from pedigree or genomic (SNP-by-SNP or ROH) information. In general, estimates of effective population size were above 200 in Lacaune Confederation and Lacaune Ovitest subpopulations and below 200 in Basco-Béarnaise, Manech Tête Noire, and Manech Tête Rousse breeds. The minimum length that constituted a ROH, the minimum number of SNP that constituted a ROH, as well as the minimum density and the maximum distance allowed between 2 homozygous SNP are ROH-defining factors with important implications in the estimation of the rate of inbreeding. The ROH-based rates of inbreeding in concordance with those obtained from pedigree information require a specific set of values. This particular set of values is different from that identified to obtain ROH-based rates of inbreeding similar to those obtained on a SNP-by-SNP basis. Factors to define ROH do not change the results much unless extreme values are considered, although further research on ROH-based inbreeding is still required.
Collapse
Affiliation(s)
| | - J M Elsen
- INRA, UMR 1388 GenPhySE, 31326 Castanet Tolosan, France
| | - A Legarra
- INRA, UMR 1388 GenPhySE, 31326 Castanet Tolosan, France
| |
Collapse
|
248
|
New world goat populations are a genetically diverse reservoir for future use. Sci Rep 2019; 9:1476. [PMID: 30728441 PMCID: PMC6365549 DOI: 10.1038/s41598-019-38812-3] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2018] [Accepted: 10/30/2018] [Indexed: 01/02/2023] Open
Abstract
Western hemisphere goats have European, African and Central Asian origins, and some local or rare breeds are reported to be adapted to their environments and economically important. By-in-large these genetic resources have not been quantified. Using 50 K SNP genotypes of 244 animals from 12 goat populations in United States, Costa Rica, Brazil and Argentina, we evaluated the genetic diversity, population structure and selective sweeps documenting goat migration to the "New World". Our findings suggest the concept of breed, particularly among "locally adapted" breeds, is not a meaningful way to characterize goat populations. The USA Spanish goats were found to be an important genetic reservoir, sharing genomic composition with the wild ancestor and with specialized breeds (e.g. Angora, Lamancha and Saanen). Results suggest goats in the Americas have substantial genetic diversity to use in selection and promote environmental adaptation or product driven specialization. These findings highlight the importance of maintaining goat conservation programs and suggest an awaiting reservoir of genetic diversity for breeding and research while simultaneously discarding concerns about breed designations.
Collapse
|
249
|
Bolormaa S, Chamberlain AJ, Khansefid M, Stothard P, Swan AA, Mason B, Prowse-Wilkins CP, Duijvesteijn N, Moghaddar N, van der Werf JH, Daetwyler HD, MacLeod IM. Accuracy of imputation to whole-genome sequence in sheep. Genet Sel Evol 2019; 51:1. [PMID: 30654735 PMCID: PMC6337865 DOI: 10.1186/s12711-018-0443-5] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2018] [Accepted: 12/18/2018] [Indexed: 12/12/2022] Open
Abstract
Background The use of whole-genome sequence (WGS) data for genomic prediction and association studies is highly desirable because the causal mutations should be present in the data. The sequencing of 935 sheep from a range of breeds provides the opportunity to impute sheep genotyped with single nucleotide polymorphism (SNP) arrays to WGS. This study evaluated the accuracy of imputation from SNP genotypes to WGS using this reference population of 935 sequenced sheep. Results The accuracy of imputation from the Ovine Infinium® HD BeadChip SNP (~ 500 k) to WGS was assessed for three target breeds: Merino, Poll Dorset and F1 Border Leicester × Merino. Imputation accuracy was highest for the Poll Dorset breed, although there were more Merino individuals in the sequenced reference population than Poll Dorset individuals. In addition, empirical imputation accuracies were higher (by up to 1.7%) when using larger multi-breed reference populations compared to using a smaller single-breed reference population. The mean accuracy of imputation across target breeds using the Minimac3 or the FImpute software was 0.94. The empirical imputation accuracy varied considerably across the genome; six chromosomes carried regions of one or more Mb with a mean imputation accuracy of < 0.7. Imputation accuracy in five variant annotation classes ranged from 0.87 (missense) up to 0.94 (intronic variants), where lower accuracy corresponded to higher proportions of rare alleles. The imputation quality statistic reported from Minimac3 (R2) had a clear positive relationship with the empirical imputation accuracy. Therefore, by first discarding imputed variants with an R2 below 0.4, the mean empirical accuracy across target breeds increased to 0.97. Although accuracy of genomic prediction was less affected by filtering on R2 in a multi-breed population of sheep with imputed WGS, the genomic heritability clearly tended to be lower when using variants with an R2 ≤ 0.4. Conclusions The mean imputation accuracy was high for all target breeds and was increased by combining smaller breed sets into a multi-breed reference. We found that the Minimac3 software imputation quality statistic (R2) was a useful indicator of empirical imputation accuracy, enabling removal of very poorly imputed variants before downstream analyses. Electronic supplementary material The online version of this article (10.1186/s12711-018-0443-5) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Sunduimijid Bolormaa
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, 5 Ring Rd, Bundoora, VIC, 3083, Australia. .,Cooperative Research Centre for Sheep Industry Innovation, Armidale, NSW, 2351, Australia.
| | - Amanda J Chamberlain
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, 5 Ring Rd, Bundoora, VIC, 3083, Australia
| | - Majid Khansefid
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, 5 Ring Rd, Bundoora, VIC, 3083, Australia.,Cooperative Research Centre for Sheep Industry Innovation, Armidale, NSW, 2351, Australia
| | - Paul Stothard
- Faculty of Agricultural, Life and Environmental Sciences, University of Alberta, Edmonton, AB, T6G 2R3, Canada
| | - Andrew A Swan
- Cooperative Research Centre for Sheep Industry Innovation, Armidale, NSW, 2351, Australia.,Animal Genetics and Breeding Unit, University of New England, Armidale, NSW, 2351, Australia
| | - Brett Mason
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, 5 Ring Rd, Bundoora, VIC, 3083, Australia
| | - Claire P Prowse-Wilkins
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, 5 Ring Rd, Bundoora, VIC, 3083, Australia
| | - Naomi Duijvesteijn
- Cooperative Research Centre for Sheep Industry Innovation, Armidale, NSW, 2351, Australia.,School of Environmental and Rural Science, University of New England, Armidale, NSW, 2351, Australia
| | - Nasir Moghaddar
- Cooperative Research Centre for Sheep Industry Innovation, Armidale, NSW, 2351, Australia.,School of Environmental and Rural Science, University of New England, Armidale, NSW, 2351, Australia
| | - Julius H van der Werf
- Cooperative Research Centre for Sheep Industry Innovation, Armidale, NSW, 2351, Australia.,School of Environmental and Rural Science, University of New England, Armidale, NSW, 2351, Australia
| | - Hans D Daetwyler
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, 5 Ring Rd, Bundoora, VIC, 3083, Australia.,Cooperative Research Centre for Sheep Industry Innovation, Armidale, NSW, 2351, Australia.,School of Applied Systems Biology, La Trobe University, Bundoora, VIC, 3086, Australia
| | - Iona M MacLeod
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, 5 Ring Rd, Bundoora, VIC, 3083, Australia.,Cooperative Research Centre for Sheep Industry Innovation, Armidale, NSW, 2351, Australia
| |
Collapse
|
250
|
Li R, Li C, Chen H, Liu X, Xiao H, Chen S. Genomic diversity and admixture patterns among six Chinese indigenous cattle breeds in Yunnan. ASIAN-AUSTRALASIAN JOURNAL OF ANIMAL SCIENCES 2019; 32:1069-1076. [PMID: 30744361 PMCID: PMC6599958 DOI: 10.5713/ajas.18.0605] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/12/2018] [Accepted: 12/08/2018] [Indexed: 11/27/2022]
Abstract
Objective Yunnan is not only a frontier zone that connects China with South and Southeast Asia, but also represents an admixture zone between taurine (Bos taurus) and zebu (Bos indicus) cattle. The purpose of this study is to understand the level of genomic diversity and the extent of admixture in each Yunnan native cattle breed. Methods All 120 individuals were genotyped using Illumina BovineHD BeadChip (777,962 single nucleotide polymorphisms [SNPs]). Quality control and genomic diversity indexes were calculated using PLINK software. The principal component analysis (PCA) was assessed using SMARTPCA program implemented in EIGENSOFT software. The ADMIXTURE software was used to reveal admixture patterns among breeds. Results A total of 604,630 SNPs was obtained after quality control procedures. Among six breeds, the highest level of mean heterozygosity was found in Zhaotong cattle from Northeastern Yunnan, whereas the lowest level of heterozygosity was detected in Dehong humped cattle from Western Yunnan. The PCA based on a pruned dataset of 233,788 SNPs clearly separated Dehong humped cattle (supposed to be a pure zebu breed) from other five breeds. The admixture analysis further revealed two clusters (K = 2 with the lowest cross validation error), corresponding to taurine and zebu cattle lineages. All six breeds except for Dehong humped cattle showed different degrees of admixture between taurine and zebu cattle. As expected, Dehong humped cattle showed no signature of taurine cattle influence. Conclusion Overall, considerable genomic diversity was found in six Yunnan native cattle breeds except for Dehong humped cattle from Western Yunnan. Dehong humped cattle is a pure zebu breed, while other five breeds had admixed origins with different extents of admixture between taurine and zebu cattle. Such admixture by crossbreeding between zebu and taurine cattle facilitated the spread of zebu cattle from tropical and subtropical regions to other highland regions in Yunnan.
Collapse
Affiliation(s)
- Rong Li
- School of Life Sciences, Yunnan University, Kunming, Yunnan 650500, China
| | - Chunqing Li
- School of Life Sciences, Yunnan University, Kunming, Yunnan 650500, China.,National Demonstration Center for Experimental Life Sciences Education, Yunnan University, Kunming Yunnan 650500, China
| | - Hongyu Chen
- School of Life Sciences, Yunnan University, Kunming, Yunnan 650500, China
| | - Xuehong Liu
- Faculty of Animal Science and Technology, Yunnan Agricultural University, Kunming, Yunnan 650221, China
| | - Heng Xiao
- School of Life Sciences, Yunnan University, Kunming, Yunnan 650500, China
| | - Shanyuan Chen
- School of Life Sciences, Yunnan University, Kunming, Yunnan 650500, China.,National Demonstration Center for Experimental Life Sciences Education, Yunnan University, Kunming Yunnan 650500, China
| |
Collapse
|