1
|
Le SV, de Las Heras-Saldana S, Alexandri P, Olmo L, Walkden-Brown SW, van der Werf JHJ. Genetic diversity, population structure and origin of the native goats in Central Laos. J Anim Breed Genet 2024; 141:531-549. [PMID: 38520124 DOI: 10.1111/jbg.12862] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2023] [Revised: 02/14/2024] [Accepted: 03/10/2024] [Indexed: 03/25/2024]
Abstract
Maintaining genetic diversity and variation in livestock populations is critical for natural and artificial selection promoting genetic improvement while avoiding problems due to inbreeding. In Laos, there are concerns that there has been a decline in genetic diversity and a rise in inbreeding among native goats in their village-based smallholder system. In this study, we investigated the genetic diversity of Lao native goats in Phin, Songkhone and Sepon districts in Central Laos for the first time using Illumina's Goat SNP50 BeadChip. We also explored the genetic relationships between Lao goats with 163 global goat populations from 36 countries. Our results revealled a close genetic relationship between Lao native goats and Chinese, Mongolian and Pakistani goats, sharing ancestries with Guangfen, Jining Grey and Luoping Yellow breeds (China) and Teddi goats (Pakistan). The observed (Ho) and expected (He) heterozygosity were 0.292 and 0.303 (Laos), 0.288 and 0.288 (Sepon), 0.299 and 0.308 (Phin) and 0.289 and 0.305 (Songkhone), respectively. There was low to moderate genetic differentiation (FST: 0.011-0.043) and negligible inbreeding coefficients (FIS: -0.001 to 0.052) between goat districts. The runs of homozygosity (ROH) had an average length of 5.92-6.85 Mb, with short ROH segments (1-5 Mb length) being the most prevalent (66.34%). Longer ROH segments (20-40 and >40 Mb length categories) were less common, comprising only 4.81% and 1.01%, respectively. Lao goats exhibit moderate genetic diversity, low-inbreeding levels and adequate effective population size. Some genetic distinctions between Lao goats may be explained by geographic and cultural features.
Collapse
Affiliation(s)
- Sang V Le
- School of Rural & Environmental Science, University of New England, Armidale, New South Wales, Australia
| | - Sara de Las Heras-Saldana
- AGBU, a Joint Venture of the NSW Department of Primary Industries and the University of New England, Armidale, New South Wales, Australia
| | - Panoraia Alexandri
- AGBU, a Joint Venture of the NSW Department of Primary Industries and the University of New England, Armidale, New South Wales, Australia
| | - Luisa Olmo
- School of Rural & Environmental Science, University of New England, Armidale, New South Wales, Australia
| | - Stephen W Walkden-Brown
- School of Rural & Environmental Science, University of New England, Armidale, New South Wales, Australia
| | - Julius H J van der Werf
- School of Rural & Environmental Science, University of New England, Armidale, New South Wales, Australia
| |
Collapse
|
2
|
Petretto E, Dettori ML, Luigi-Sierra MG, Noce A, Pazzola M, Vacca GM, Molina A, Martínez A, Goyache F, Carolan S, Amills M. Investigating the footprint of post-domestication dispersal on the diversity of modern European, African and Asian goats. Genet Sel Evol 2024; 56:55. [PMID: 39068382 PMCID: PMC11282621 DOI: 10.1186/s12711-024-00923-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2023] [Accepted: 07/05/2024] [Indexed: 07/30/2024] Open
Abstract
BACKGROUND Goats were domesticated in the Fertile Crescent about 10,000 years before present (YBP) and subsequently spread across Eurasia and Africa. This dispersal is expected to generate a gradient of declining genetic diversity with increasing distance from the areas of early livestock management. Previous studies have reported the existence of such genetic cline in European goat populations, but they were based on a limited number of microsatellite markers. Here, we have analyzed data generated by the AdaptMap project and other studies. More specifically, we have used the geographic coordinates and estimates of the observed (Ho) and expected (He) heterozygosities of 1077 European, 1187 African and 617 Asian goats belonging to 38, 43 and 22 different breeds, respectively, to find out whether genetic diversity and distance to Ganj Dareh, a Neolithic settlement in western Iran for which evidence of an early management of domestic goats has been obtained, are significantly correlated. RESULTS Principal component and ADMIXTURE analyses revealed an incomplete regional differentiation of European breeds, but two genetic clusters representing Northern Europe and the British-Irish Isles were remarkably differentiated from the remaining European populations. In African breeds, we observed five main clusters: (1) North Africa, (2) West Africa, (3) East Africa, (4) South Africa, and (5) Madagascar. Regarding Asian breeds, three well differentiated West Asian, South Asian and East Asian groups were observed. For European and Asian goats, no strong evidence of significant correlations between Ho and He and distance to Ganj Dareh was found. In contrast, in African breeds we detected a significant gradient of diversity, which decreased with distance to Ganj Dareh. CONCLUSIONS The detection of a genetic cline associated with distance to the Ganj Dareh in African but not in European or Asian goat breeds might reflect differences in the post-domestication dispersal process and subsequent migratory movements associated with the management of caprine populations from these three continents.
Collapse
Affiliation(s)
- Elena Petretto
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
- Department of Veterinary Medicine, University of Sassari, 07100, Sassari, Italy
- Departament de Ciència Animal i dels Aliments, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Maria Luisa Dettori
- Department of Veterinary Medicine, University of Sassari, 07100, Sassari, Italy
| | - María Gracia Luigi-Sierra
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
- Departament de Ciència Animal i dels Aliments, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Antonia Noce
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Michele Pazzola
- Department of Veterinary Medicine, University of Sassari, 07100, Sassari, Italy
| | | | - Antonio Molina
- Department of Genetics, University of Cordoba, 14071, Córdoba, Spain
| | - Amparo Martínez
- Department of Genetics, University of Cordoba, 14071, Córdoba, Spain
| | - Félix Goyache
- Área de Genética y Reproducción Animal, SERIDA-Deva, Camino de Rioseco 1225, 33394, Gijón, Spain
| | | | - Marcel Amills
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain.
- Departament de Ciència Animal i dels Aliments, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain.
| |
Collapse
|
3
|
Gmel AI, Mikko S, Ricard A, Velie BD, Gerber V, Hamilton NA, Neuditschko M. Using high-density SNP data to unravel the origin of the Franches-Montagnes horse breed. Genet Sel Evol 2024; 56:53. [PMID: 38987703 PMCID: PMC11238448 DOI: 10.1186/s12711-024-00922-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Accepted: 07/03/2024] [Indexed: 07/12/2024] Open
Abstract
BACKGROUND The Franches-Montagnes (FM) is the last native horse breed of Switzerland, established at the end of the 19th century by cross-breeding local mares with Anglo-Norman stallions. We collected high-density SNP genotype data (Axiom™ 670 K Equine genotyping array) from 522 FM horses, including 44 old-type horses (OF), 514 European Warmblood horses (WB) from Sweden and Switzerland (including a stallion used for cross-breeding in 1990), 136 purebred Arabians (AR), 32 Shagya Arabians (SA), and 64 Thoroughbred (TB) horses, as introgressed WB stallions showed TB origin in their pedigrees. The aim of the study was to ascertain fine-scale population structures of the FM breed, including estimation of individual admixture levels and genomic inbreeding (FROH) by means of Runs of Homozygosity. RESULTS To assess fine-scale population structures within the FM breed, we applied a three-step approach, which combined admixture, genetic contribution, and FROH of individuals into a high-resolution network visualization. Based on this approach, we were able to demonstrate that population substructures, as detected by model-based clustering, can be either associated with a different genetic origin or with the progeny of most influential sires. Within the FM breed, admixed horses explained most of the genetic variance of the current breeding population, while OF horses only accounted for a small proportion of the variance. Furthermore, we illustrated that FM horses showed high TB admixture levels and we identified inconsistencies in the origin of FM horses descending from the Arabian stallion Doktryner. With the exception of WB, FM horses were less inbred compared to the other breeds. However, the relatively few but long ROH segments suggested diversity loss in both FM subpopulations. Genes located in FM- and OF-specific ROH islands had known functions involved in conformation and behaviour, two traits that are highly valued by breeders. CONCLUSIONS The FM remains the last native Swiss breed, clearly distinguishable from other historically introgressed breeds, but it suffered bottlenecks due to intensive selection of stallions, restrictive mating choices based on arbitrary definitions of pure breeding, and selection of rare coat colours. To preserve the genetic diversity of FM horses, future conservation managements strategies should involve a well-balanced selection of stallions (e.g., by integrating OF stallions in the FM breeding population) and avoid selection for rare coat colours.
Collapse
Affiliation(s)
- Annik Imogen Gmel
- Animal GenoPhenomics, Agroscope, Route de la Tioleyre 4, 1725, Posieux, Switzerland
- Equine Department, Vetsuisse Faculty, University of Zurich, Winterthurerstrasse 260, 8053, Zurich, Switzerland
| | - Sofia Mikko
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Box 7023, 750 07, Uppsala, Sweden
| | - Anne Ricard
- Institut National de la Recherche Agronomique, Domaine de Vilvert, 78350, Jouy-en-Josas, France
| | - Brandon D Velie
- Equine Genetics and Genomics Group, School of Life and Environmental Sciences, University of Sydney, RMC Gunn B19-603, Sydney, NSW, 2006, Australia
| | - Vinzenz Gerber
- Institut Suisse de Médecine Equine ISME, Vetsuisse Faculty, University of Bern, Länggassstrasse 124, 3012, Bern, Switzerland
| | - Natasha Anne Hamilton
- Sydney School of Veterinary Science, University of Sydney, Sydney, NSW, 2006, Australia
| | - Markus Neuditschko
- Animal GenoPhenomics, Agroscope, Route de la Tioleyre 4, 1725, Posieux, Switzerland.
| |
Collapse
|
4
|
Vostry L, Vostra-Vydrova H, Moravcikova N, Kasarda R, Margetin M, Rychtarova J, Drzaic I, Shihabi M, Cubric-Curik V, Sölkner J, Curik I. Genomic analysis of conservation status, population structure and admixture in local Czech and Slovak dairy goat breeds. J Dairy Sci 2024:S0022-0302(24)00937-8. [PMID: 38908686 DOI: 10.3168/jds.2023-24607] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2023] [Accepted: 05/30/2024] [Indexed: 06/24/2024]
Abstract
While dairy goat production, characterized by traditional production on small farms, is an important source of income in the Czech Republic and Slovakia, locally adapted breeds have not been fully consolidated over the last 100 years due to large fluctuations in population size and inconsistent breeding programs that allowed for different crossbreeding strategies. Our main objective in this study was therefore to assess the conservation status of 4 Czech (Alpine Goat, White Shorthair, Brown Shorthair and Czech Landrace) and one Slovak (Slovak White Shorthair) local goat breeds, to analyze their population structure and admixture, and to estimate their relatedness to several neighboring breeds. Our analyses included 142 goats belonging to 5 local breeds genotyped with the Illumina 50K BeadChip and 618 previously genotyped animals representing 15 goat breeds from Austria and Switzerland (all analyses based on 46,862 autosomal SNPs and 760 animals). In general, the conservation status of the Czech and Slovak local goat breeds was satisfactory, with the exception of the Brown Shorthair goat, as the analyzed parameters (heterozygosity, haplotype richness, ROH-based inbreeding and effective population size) were mostly above the median of 20 breeds. However, for all 5 Czech and Slovakian breeds, an examination of historical effective population size indicated a substantial decline about 8 to 22 generations ago. In addition, our study revealed that the Czech and Slovakian breeds are not fully consolidated; for instance, White Shorthair and Brown Shorthair were not clearly distinguishable. Considerable admixture, especially in Czech Landrace (effective number of parental clusters equal to 4.2), and low but numerous migration rates from other Austrian and Swiss breeds were found. These results provide valuable insights for future breeding programs and genetic diversity management of local Czech and Slovak goat breeds.
Collapse
Affiliation(s)
- Lubos Vostry
- Czech University of Life Science Prague, Kamycka 129, 16500 Prague, Czech Republic.
| | - Hana Vostra-Vydrova
- Czech University of Life Science Prague, Kamycka 129, 16500 Prague, Czech Republic
| | - Nina Moravcikova
- Slovak University of Agriculture in Nitra, Tr. A. Hlinku 2, 94976 Nitra, Slovak Republic
| | - Radovan Kasarda
- Slovak University of Agriculture in Nitra, Tr. A. Hlinku 2, 94976 Nitra, Slovak Republic
| | - Milan Margetin
- Slovak University of Agriculture in Nitra, Tr. A. Hlinku 2, 94976 Nitra, Slovak Republic
| | - Jana Rychtarova
- Institute of Animal Science, Přátelství 815, 104 00 Prague, Czech Republic
| | - Ivana Drzaic
- University of Zagreb, Faculty of Agriculture, Svetošimunska cesta 25, 10000 Zagreb, Croatia
| | - Mario Shihabi
- University of Zagreb, Faculty of Agriculture, Svetošimunska cesta 25, 10000 Zagreb, Croatia
| | - Vlatka Cubric-Curik
- University of Zagreb, Faculty of Agriculture, Svetošimunska cesta 25, 10000 Zagreb, Croatia
| | - Johan Sölkner
- University of Natural Resources & Life Sciences Vienna, Gregor-Mendel-Strasse 33, 1180 Vienna, Austria
| | - Ino Curik
- University of Zagreb, Faculty of Agriculture, Svetošimunska cesta 25, 10000 Zagreb, Croatia; Institute of Animal Sciences, Hungarian University of Agriculture and Life Sciences (MATE), Guba Sándor u. 40, 7400 Kaposvár, Hungary.
| |
Collapse
|
5
|
Pogorevc N, Dotsev A, Upadhyay M, Sandoval-Castellanos E, Hannemann E, Simčič M, Antoniou A, Papachristou D, Koutsouli P, Rahmatalla S, Brockmann G, Sölkner J, Burger P, Lymberakis P, Poulakakis N, Bizelis I, Zinovieva N, Horvat S, Medugorac I. Whole-genome SNP genotyping unveils ancestral and recent introgression in wild and domestic goats. Mol Ecol 2024; 33:e17190. [PMID: 37909668 DOI: 10.1111/mec.17190] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Revised: 10/04/2023] [Accepted: 10/16/2023] [Indexed: 11/03/2023]
Abstract
After the domestication of goats around 10,000 years before the present (BP), humans transported goats far beyond the range of their wild ancestor, the bezoar goat. This brought domestic goats into contact with many wild goat species such as ibex and markhor, enabling introgression between domestic and wild goats. To investigate this, while shedding light on the taxonomic status of wild and domestic goats, we analysed genome-wide SNP data of 613 specimens from 14 taxonomic units, including Capra hircus, C. pyrenaica, C. ibex (from Switzerland, Austria, Germany and Slovenia), C. aegagrus aegagrus, C. a. cretica, C. h. dorcas, C. caucasica caucasica, C. c. severtzovi, C. c. cylindricornis, C. falconeri, C. sibirica sibirica, C. s. alaiana and C. nubiana, as well as Oreamnos americanus (mountain goat) as an outgroup. To trace gene flow between domestic and wild goats, we integrated genotype data of local goat breeds from the Alps as well as from countries such as Spain, Greece, Türkiye, Egypt, Sudan, Iran, Russia (Caucasus and Altai) and Pakistan. Our phylogenetic analyses displayed a clear separation between bezoar-type and ibex-type clades with wild goats from the Greek islands of Crete and Youra clustered within domestic goats, confirming their feral origin. Our analyses also revealed gene flow between the lineages of Caucasian tur and domestic goats that most likely occurred before or during early domestication. Within the clade of domestic goats, analyses inferred gene flow between African and Iberian goats. The detected events of introgression were consistent with previous reports and offered interesting insights into the historical relationships among domestic and wild goats.
Collapse
Affiliation(s)
- Neža Pogorevc
- Department of Animal Science, Biotechnical Faculty, University of Ljubljana, Ljubljana, Slovenia
- Population Genomics Group, Department of Veterinary Sciences, Faculty of Veterinary Medicine, Ludwig-Maximilians-University of Munich, Martinsried/Planegg, Germany
| | - Arsen Dotsev
- L.K. Ernst Federal Research Center for Animal Husbandry, Podolsk, Russia
| | - Maulik Upadhyay
- Population Genomics Group, Department of Veterinary Sciences, Faculty of Veterinary Medicine, Ludwig-Maximilians-University of Munich, Martinsried/Planegg, Germany
| | - Edson Sandoval-Castellanos
- Population Genomics Group, Department of Veterinary Sciences, Faculty of Veterinary Medicine, Ludwig-Maximilians-University of Munich, Martinsried/Planegg, Germany
| | - Elisabeth Hannemann
- Population Genomics Group, Department of Veterinary Sciences, Faculty of Veterinary Medicine, Ludwig-Maximilians-University of Munich, Martinsried/Planegg, Germany
| | - Mojca Simčič
- Department of Animal Science, Biotechnical Faculty, University of Ljubljana, Ljubljana, Slovenia
| | - Aglaia Antoniou
- Institute of Marine Biology, Biotechnology and Aquaculture (IMBBC), Hellenic Centre for Marine Research (HCMR), Crete, Greece
| | - Dimitris Papachristou
- Laboratory of Animal Breeding and Husbandry, Department of Animal Science, School of Animal Biosciences, Agricultural University of Athens, Athens, Greece
| | - Panagiota Koutsouli
- Laboratory of Animal Breeding and Husbandry, Department of Animal Science, School of Animal Biosciences, Agricultural University of Athens, Athens, Greece
| | - Siham Rahmatalla
- Albrecht Daniel Thaer-Institute for Agricultural and Horticultural Sciences, Animal Breeding and Molecular Genetics, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Gudrun Brockmann
- Albrecht Daniel Thaer-Institute for Agricultural and Horticultural Sciences, Animal Breeding and Molecular Genetics, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Johann Sölkner
- Division of Livestock Sciences, Department of Sustainable Agricultural Systems, University of Natural Resources and Life Sciences, Vienna, Austria
| | - Pamela Burger
- Department of Interdisciplinary Life Sciences, Research Institute of Wildlife Ecology, Vetmeduni Vienna, Vienna, Austria
| | - Petros Lymberakis
- Natural History Museum of Crete, School of Sciences and Engineering, University of Crete, Irakleio, Greece
| | - Nikos Poulakakis
- Natural History Museum of Crete, School of Sciences and Engineering, University of Crete, Irakleio, Greece
- Biology Department, School of Sciences and Engineering, University of Crete, Irakleio, Greece
- Institute of Molecular Biology and Biotechnology (IMBB), Foundation for Research and Technology - Hellas (FORTH), Irakleio, Greece
| | - Iosif Bizelis
- Laboratory of Animal Breeding and Husbandry, Department of Animal Science, School of Animal Biosciences, Agricultural University of Athens, Athens, Greece
| | - Natalia Zinovieva
- L.K. Ernst Federal Research Center for Animal Husbandry, Podolsk, Russia
| | - Simon Horvat
- Department of Animal Science, Biotechnical Faculty, University of Ljubljana, Ljubljana, Slovenia
| | - Ivica Medugorac
- Population Genomics Group, Department of Veterinary Sciences, Faculty of Veterinary Medicine, Ludwig-Maximilians-University of Munich, Martinsried/Planegg, Germany
| |
Collapse
|
6
|
Abdoli R, Mirhoseini SZ, Ghavi Hossein-Zadeh N, Zamani P, Moradi MH, Ferdosi MH, Sargolzaei M, Gondro C. Runs of homozygosity and cross-generational inbreeding of Iranian fat-tailed sheep. Heredity (Edinb) 2023; 130:358-367. [PMID: 37016136 PMCID: PMC10238534 DOI: 10.1038/s41437-023-00611-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2021] [Revised: 03/09/2023] [Accepted: 03/13/2023] [Indexed: 04/06/2023] Open
Abstract
The Lori-Bakhtiari fat-tailed sheep is one of the most important heavyweight native breeds of Iran. The breed is robust and well-adapted to semi-arid regions and an important resource for smallholder farms. An established nucleus-based breeding scheme is used to improve their production traits but there is an indication of inbreeding depression and loss of genetic diversity due to selection. Here, we estimated the inbreeding levels and the distribution of runs of homozygosity (ROH) islands in 122 multi-generational female Lori-Bakhtiari from different half-sib families selected from a breeding station that were genotyped on the 50k array. A total of 2404 ROH islands were identified. On average, there were 19.70 ± 1.4 ROH per individual ranging between 6 and 41. The mean length of the ROH was 4.1 ± 0.14 Mb. There were 1999 short ROH of length 1-6 Mb and another 300 in the 6-12 Mb range. Additionally long ROH indicative of inbreeding were found in the ranges of 12-24 Mb (95) and 24-48 Mb (10). The average inbreeding coefficient (FROH) was 0.031 ± 0.003 with estimates varying from 0.006 to 0.083. Across generations, FROH increased from 0.019 ± 0.012 to 0.036 ± 0.007. Signatures of selection were identified on chromosomes 2, 6, and 10, encompassing 55 genes and 23 QTL associated with production traits. Inbreeding coefficients are currently within acceptable levels but across generations, inbreeding is increasing due to selection. The breeding program needs to actively monitor future inbreeding rates and ensure that the breed maintains or improves on its current levels of environmental adaptation.
Collapse
Affiliation(s)
- Ramin Abdoli
- Iran Silk Research Center, Agricultural Research, Education and Extension Organization (AREEO), Gilan, Iran.
| | - Seyed Ziaeddin Mirhoseini
- Department of Animal Science, Faculty of Agricultural Sciences, University of Guilan, Rasht, 41635-1314, Iran.
| | - Navid Ghavi Hossein-Zadeh
- Department of Animal Science, Faculty of Agricultural Sciences, University of Guilan, Rasht, 41635-1314, Iran
| | - Pouya Zamani
- Department of Animal Science, Faculty of Agriculture, Bu-Ali Sina University, Hamedan, 65178-33131, Iran
| | - Mohammad Hossein Moradi
- Department of Animal Science, Faculty of Agriculture and Natural Resources, Arak University, Arāk, 38156-8-8349, Iran
| | - Mohammad Hossein Ferdosi
- AGBU, a joint venture of NSW Department of Primary Industries and University of New England, University of New England, Armidale, NSW, 2351, Australia
| | - Mehdi Sargolzaei
- Department of Pathobiology, Ontario Veterinary College, University of Guelph, Guelph, Canada
| | - Cedric Gondro
- Department of Animal Science, College of Agriculture and Natural Resources, Michigan State University, East Lansing, MI, 48824, USA
| |
Collapse
|
7
|
Sun X, Niu Q, Jiang J, Wang G, Zhou P, Li J, Chen C, Liu L, Xu L, Ren H. Identifying Candidate Genes for Litter Size and Three Morphological Traits in Youzhou Dark Goats Based on Genome-Wide SNP Markers. Genes (Basel) 2023; 14:1183. [PMID: 37372363 DOI: 10.3390/genes14061183] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2023] [Revised: 05/25/2023] [Accepted: 05/25/2023] [Indexed: 06/29/2023] Open
Abstract
This study aimed to reveal the potential genetic basis for litter size, coat colour, black middorsal stripe and skin colour by combining genome-wide association analysis (GWAS) and selection signature analysis and ROH detection within the Youzhou dark (YZD) goat population (n = 206) using the Illumina GoatSNP54 BeadChip. In the GWAS, we identified one SNP (snp54094-scaffold824-899720) on chromosome 11 for litter size, two SNPs on chromosome 26 (snp11508-scaffold142-1990450, SORCS3) and chromosome 12 (snp55048-scaffold842-324525, LOC102187779) for coat colour and one SNP on chromosome 18 (snp56013-scaffold873-22716, TCF25) for the black middorsal stripe. In contrast, no SNPs were identified for skin colour. In selection signature analysis, 295 significant iHS genomic regions with a mean |iHS| score > 2.66, containing selection signatures encompassing 232 candidate genes were detected. In particular, 43 GO terms and one KEGG pathway were significantly enriched in the selected genes, which may contribute to the excellent environmental adaptability and characteristic trait formation during the domestication of YZD goats. In ROH detection, we identified 4446 ROH segments and 282 consensus ROH regions, among which nine common genes overlapped with those detected using the iHS method. Some known candidate genes for economic traits such as reproduction (TSHR, ANGPT4, CENPF, PIBF1, DACH1, DIS3, CHST1, COL4A1, PRKD1 and DNMT3B) and development and growth (TNPO2, IFT80, UCP2, UCP3, GHRHR, SIM1, CCM2L, CTNNA3 and CTNNA1) were revealed by iHS and ROH detection. Overall, this study is limited by the small population size, which affects the results of GWAS to a certain extent. Nevertheless, our findings could provide the first overview of the genetic mechanism underlying these important traits and provide novel insights into the future conservation and utilisation of Chinese goat germplasm resources.
Collapse
Affiliation(s)
- Xiaoyan Sun
- Chongqing Academy of Animal Sciences, Rongchang 402460, China
| | - Qunhao Niu
- Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China
| | - Jing Jiang
- Chongqing Academy of Animal Sciences, Rongchang 402460, China
| | - Gaofu Wang
- Chongqing Academy of Animal Sciences, Rongchang 402460, China
| | - Peng Zhou
- Chongqing Academy of Animal Sciences, Rongchang 402460, China
| | - Jie Li
- Chongqing Academy of Animal Sciences, Rongchang 402460, China
| | - Cancan Chen
- Chongqing Academy of Animal Sciences, Rongchang 402460, China
| | - Liangjia Liu
- Chongqing Academy of Animal Sciences, Rongchang 402460, China
| | - Lingyang Xu
- Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China
| | - Hangxing Ren
- Chongqing Academy of Animal Sciences, Rongchang 402460, China
| |
Collapse
|
8
|
Zhong T, Wang X, Huang C, Yang L, Zhao Q, Chen X, Freitas-de-Melo A, Zhan S, Wang L, Dai D, Cao J, Guo J, Li L, Zhang H, Niu L. A genome-wide perspective on the diversity and selection signatures in indigenous goats using 53 K single nucleotide polymorphism array. Animal 2023; 17:100706. [PMID: 36758301 DOI: 10.1016/j.animal.2023.100706] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2022] [Revised: 01/03/2023] [Accepted: 01/06/2023] [Indexed: 01/15/2023] Open
Abstract
Tibetan goats, Taihang goats, Jining grey goats, and Meigu goats are the representative indigenous goats in China, found in Qinghai-Tibet Plateau, Western pastoral area, Northern and Southern agricultural regions. Very few studies have conducted a comprehensive analysis of the genomic diversity and selection of these breeds. We genotyped 96 unrelated individuals, using goat 53 K Illumina BeadChip array, of the following goat breeds: Tibetan (TG), Taihang (THG), Jining grey (JGG), and Meigu (MGG). A total of 45 951 single nucleotide polymorphisms were filtered to estimate the genetic diversity and selection signatures. All breeds had a high proportion (over 95%) of polymorphic loci. The observed and excepted heterozygosity ranged from 0.338 (MGG) to 0.402 (JGG) and 0.339 (MGG) to 0.395 (JGG), respectively. Clustering analysis displayed a genetically distinct lineage for each breed, and their Fst were greater than 0.25, indicating that they had a higher genetic differentiation between groups. Furthermore, effective population size reduced in all four populations, indicating a loss of genetic diversity. In addition, runs of homozygosity were mainly distributed in 5-10 Mb. Lastly, we identified signature genes, which were closely related to high-altitude adaptation (ADIRF) and prolificity (CNTROB, SMC3, and PTEN). This study provides a valuable resource for future studies on genome-wide perspectives on the diversity and selection signatures of Chinese indigenous goats.
Collapse
Affiliation(s)
- Tao Zhong
- Key Laboratory of Livestock and Poultry Multi-omics, Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
| | - Xinlu Wang
- Key Laboratory of Livestock and Poultry Multi-omics, Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
| | - Chunhua Huang
- Key Laboratory of Livestock and Poultry Multi-omics, Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
| | - Liu Yang
- Key Laboratory of Livestock and Poultry Multi-omics, Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
| | - Qianjun Zhao
- Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China
| | - Xiaoyong Chen
- College of Animal Science and Technology, Hebei Agricultural University, Baoding 071001, China
| | - Aline Freitas-de-Melo
- Departamento de Biociencias Veterinarias, Facultad de Veterinaria, Universidad de la República, Montevideo, Uruguay
| | - Siyuan Zhan
- Key Laboratory of Livestock and Poultry Multi-omics, Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
| | - Linjie Wang
- Key Laboratory of Livestock and Poultry Multi-omics, Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
| | - Dinghui Dai
- Key Laboratory of Livestock and Poultry Multi-omics, Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
| | - Jiaxue Cao
- Key Laboratory of Livestock and Poultry Multi-omics, Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
| | - Jiazhong Guo
- Key Laboratory of Livestock and Poultry Multi-omics, Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
| | - Li Li
- Key Laboratory of Livestock and Poultry Multi-omics, Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
| | - Hongping Zhang
- Key Laboratory of Livestock and Poultry Multi-omics, Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
| | - Lili Niu
- Key Laboratory of Livestock and Poultry Multi-omics, Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China.
| |
Collapse
|
9
|
Gmel AI, Guichard M, Dainat B, Williams GR, Eynard S, Vignal A, Servin B, Neuditschko M. Identification of runs of homozygosity in Western honey bees ( Apis mellifera) using whole-genome sequencing data. Ecol Evol 2023; 13:e9723. [PMID: 36694553 PMCID: PMC9843643 DOI: 10.1002/ece3.9723] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2021] [Revised: 12/15/2022] [Accepted: 12/19/2022] [Indexed: 01/19/2023] Open
Abstract
Runs of homozygosity (ROH) are continuous homozygous segments that arise through the transmission of haplotypes that are identical by descent. The length and distribution of ROH segments provide insights into the genetic diversity of populations and can be associated with selection signatures. Here, we analyzed reconstructed whole-genome queen genotypes, from a pool-seq data experiment including 265 Western honeybee colonies from Apis mellifera mellifera and Apis mellifera carnica. Integrating individual ROH patterns and admixture levels in a dynamic population network visualization allowed us to ascertain major differences between the two subspecies. Within A. m. mellifera, we identified well-defined substructures according to the genetic origin of the queens. Despite the current applied conservation efforts, we pinpointed 79 admixed queens. Genomic inbreeding (F ROH) strongly varied within and between the identified subpopulations. Conserved A. m. mellifera from Switzerland had the highest mean F ROH (3.39%), while queens originating from a conservation area in France, which were also highly admixed, showed significantly lower F ROH (0.45%). The majority of A. m. carnica queens were also highly admixed, except 12 purebred queens with a mean F ROH of 2.33%. Within the breed-specific ROH islands, we identified 14 coding genes for A. m. mellifera and five for A. m. carnica, respectively. Local adaption of A. m. mellifera could be suggested by the identification of genes involved in the response to ultraviolet light (Crh-BP, Uvop) and body size (Hex70a, Hex70b), while the A. m. carnica specific genes Cpr3 and Cpr4 are most likely associated with the lighter striping pattern, a morphological phenotype expected in this subspecies. We demonstrated that queen genotypes derived from pooled workers are useful tool to unravel the population dynamics in A. mellifera and provide fundamental information to conserve native honey bees.
Collapse
Affiliation(s)
- Annik Imogen Gmel
- Animal GenoPhenomics, Animal Production Systems and Animal HealthAgroscopePosieuxSwitzerland
| | - Matthieu Guichard
- Animal GenoPhenomics, Animal Production Systems and Animal HealthAgroscopePosieuxSwitzerland
- Swiss Bee Research CentreAgroscopeLiebefeldSwitzerland
| | | | | | - Sonia Eynard
- GenPhySEINRAE, INPT, INPENVTUniversité de ToulouseCastanet‐TolosanFrance
- UMT PrADEProtection des Abeilles Dans L'EnvironnementAvignonFrance
| | - Alain Vignal
- GenPhySEINRAE, INPT, INPENVTUniversité de ToulouseCastanet‐TolosanFrance
- UMT PrADEProtection des Abeilles Dans L'EnvironnementAvignonFrance
| | - Bertrand Servin
- GenPhySEINRAE, INPT, INPENVTUniversité de ToulouseCastanet‐TolosanFrance
- UMT PrADEProtection des Abeilles Dans L'EnvironnementAvignonFrance
| | | | - Markus Neuditschko
- Animal GenoPhenomics, Animal Production Systems and Animal HealthAgroscopePosieuxSwitzerland
| |
Collapse
|
10
|
Magoro AM, Mtileni B, Hadebe K, Zwane A. Assessment of Genetic Diversity and Conservation in South African Indigenous Goat Ecotypes: A Review. Animals (Basel) 2022; 12:ani12233353. [PMID: 36496874 PMCID: PMC9735466 DOI: 10.3390/ani12233353] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2022] [Revised: 11/21/2022] [Accepted: 11/25/2022] [Indexed: 12/04/2022] Open
Abstract
Goats were amongst the first livestock to be domesticated more than 10,000 years ago for their meat, milk, skin, and fiber. They were introduced to Southern Africa by migrating nations from Central Africa to the south. Due to local adaptation to the different agro-ecological zones and selection, indigenous goats are identified as ecotypes within the indigenous veld goat breed. Their ability to thrive in a resource-limited production system and in challenging environmental conditions makes them valuable animal resources for small-scale and emerging farmers. They play important roles in household agriculture and cultural activities as well as in poverty alleviation. Studies have described the phenotypic and genetic variations in indigenous goats, targeting the major goat-producing regions and the breeds of South Africa. In turn, information is restricted to certain breeds and regions, and the experimental design is often not adequate to inform the conservation status and priorities in changing environments. Advances in genomics technologies have availed more opportunities for the assessment of the biodiversity, demographic histories, and detection regions associated with local adaptation. These technologies are essential for breeding and conservation strategies for sustainable production for food security. This review focuses on the status of indigenous goats in South Africa and the application of genomics technologies for characterization, with emphasis on prioritization for conservation and sustainable utilization.
Collapse
Affiliation(s)
- Aletta Matshidiso Magoro
- Department of Animal Sciences, Tshwane University of Technology, Pretoria 0001, South Africa
- Agricultural Research Council, Animal Production, Irene 0062, South Africa
- Correspondence:
| | - Bohani Mtileni
- Department of Animal Sciences, Tshwane University of Technology, Pretoria 0001, South Africa
| | - Khanyisile Hadebe
- Agricultural Research Council, Biotechnology Platform, Onderstepoort 0110, South Africa
| | - Avhashoni Zwane
- Agricultural Research Council, Animal Production, Irene 0062, South Africa
| |
Collapse
|
11
|
Chen Y, Li R, Sun J, Li C, Xiao H, Chen S. Genome-Wide Population Structure and Selection Signatures of Yunling Goat Based on RAD-seq. Animals (Basel) 2022; 12:ani12182401. [PMID: 36139261 PMCID: PMC9495202 DOI: 10.3390/ani12182401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2022] [Revised: 09/07/2022] [Accepted: 09/10/2022] [Indexed: 12/04/2022] Open
Abstract
Simple Summary Goats are important domestic animals that provide meat, milk, fur, and other products for humans. The demand for these products has increased in recent years. Disease resistance among goat breeds is different, but the genetic basis of the differences in resistance to diseases is still unclear and needs to be further studied. In this study, many genes and pathways related to immunity and diseases were identified to be under positive selection between Yunling and Nubian goats using RAD-seq technology. This study on the selection signatures of Yunling goats provides the scientific basis and technical support for the breeding of domestic goats for disease resistance, which has important social and economic significance. Abstract Animal diseases impose a huge burden on the countries where diseases are endemic. Conventional control strategies of vaccines and veterinary drugs are to control diseases from a pharmaceutical perspective. Another alternative approach is using pre-existing genetic disease resistance or tolerance. We know that the Yunling goat is an excellent local breed from Yunnan, southwestern China, which has characteristics of strong disease resistance and remarkable adaptability. However, genetic information about the selection signatures of Yunling goats is limited. We reasoned that the genes underlying the observed difference in disease resistance might be identified by investigating selection signatures between two different goat breeds. Herein, we selected the Nubian goat as the reference group to perform the population structure and selection signature analysis by using RAD-seq technology. The results showed that two goat breeds were divided into two clusters, but there also existed gene flow. We used Fst (F-statistics) and π (pi/θπ) methods to carry out selection signature analysis. Eight selected regions and 91 candidate genes were identified, in which some genes such as DOK2, TIMM17A, MAVS, and DOCK8 related to disease and immunity and some genes such as SPEFI, CDC25B, and MIR103 were associated with reproduction. Four GO (Gene Ontology) terms (GO:0010591, GO:001601, GO:0038023, and GO:0017166) were associated with cell migration, signal transduction, and immune responses. The KEGG (Kyoto Encyclopedia of Genes and Genomes) signaling pathways were mainly associated with immune responses, inflammatory responses, and stress reactions. This study preliminarily revealed the genetic basis of strong disease resistance and adaptability of Yunling goats. It provides a theoretical basis for the subsequent genetic breeding of disease resistance of goats.
Collapse
Affiliation(s)
- Yuming Chen
- School of Ecology and Environmental Science, Yunnan University, Kunming 650500, China; (Y.C.); (R.L.); (C.L.); (H.X.)
- School of Life Sciences, Yunnan University, Kunming 650500, China;
| | - Rong Li
- School of Ecology and Environmental Science, Yunnan University, Kunming 650500, China; (Y.C.); (R.L.); (C.L.); (H.X.)
- College of Life Science, Yunnan Normal University, Kunming 650500, China
| | - Jianshu Sun
- School of Life Sciences, Yunnan University, Kunming 650500, China;
| | - Chunqing Li
- School of Ecology and Environmental Science, Yunnan University, Kunming 650500, China; (Y.C.); (R.L.); (C.L.); (H.X.)
| | - Heng Xiao
- School of Ecology and Environmental Science, Yunnan University, Kunming 650500, China; (Y.C.); (R.L.); (C.L.); (H.X.)
| | - Shanyuan Chen
- School of Ecology and Environmental Science, Yunnan University, Kunming 650500, China; (Y.C.); (R.L.); (C.L.); (H.X.)
- Correspondence: ; Tel.: +86-18687122260
| |
Collapse
|
12
|
Deniskova TE, Dotsev AV, Selionova MI, Upadhyay M, Medugorac I, Zinovieva NA. Characteristics of Dagestan Local Goat Subpopulations (Capra hircus) Based on the Analysis of the Complete Mitogenome Polymorphism. RUSS J GENET+ 2022. [DOI: 10.1134/s1022795422050040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
13
|
Capturing Genetic Diversity and Selection Signatures of the Endangered Kosovar Balusha Sheep Breed. Genes (Basel) 2022; 13:genes13050866. [PMID: 35627251 PMCID: PMC9140571 DOI: 10.3390/genes13050866] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Revised: 04/23/2022] [Accepted: 05/09/2022] [Indexed: 11/17/2022] Open
Abstract
There is a growing concern about the loss of animal genetic resources. The aim of this study was to analyze the genetic diversity and potential peculiarity of the endangered Kosovar sheep breed Balusha. For this purpose, a dataset consisting of medium-density SNP chip genotypes (39,879 SNPs) from 45 Balusha sheep was generated and compared with SNP chip genotypes from 29 individuals of a second Kosovar breed, Bardhoka. Publicly available SNP genotypes from 39 individuals of the relatively closely located sheep breeds Istrian Pramenka and Ruda were additionally included in the analyses. Analysis of heterozygosity, allelic richness and effective population size was used to assess the genetic diversity. Inbreeding was evaluated using two different methods (FIS, FROH). The standardized FST (di) and cross-population extended haplotype homozygosity (XPEHH) methods were used to detect signatures of selection. We observed the lowest heterozygosity (HO = 0.351) and effective population size (Ne5 = 25, Ne50 = 228) for the Balusha breed. The mean allelic richness levels (1.780–1.876) across all analyzed breeds were similar and also comparable with those in worldwide breeds. FROH estimates (0.023–0.077) were highest for the Balusha population, although evidence of decreased inbreeding was observed in FIS results for the Balusha breed. Two Gene Ontology (GO) TERMs were strongly enriched for Balusha, and involved genes belonging to the melanogenesis and T cell receptor signaling pathways, respectively. This could result from selection for the special coat color pattern of Balusha (black head) and resistance to certain infectious diseases. The analyzed diversity parameters highlight the urgency to preserve the local Kosovar Balusha sheep as it is clearly distinguished from other sheep of Southeastern Europe, has the lowest diversity level and may harbor valuable genetic variants, e.g., for resistance to infectious diseases.
Collapse
|
14
|
Pakpahan S, Widayanti R, Artama WT. Selection signatures and formation of the Samosir goat breed through the cultures of the Batak Toba Tribe in Samosir Island, Indonesia. Vet World 2022; 15:1044-1050. [PMID: 35698517 PMCID: PMC9178575 DOI: 10.14202/vetworld.2022.1044-1050] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2021] [Accepted: 03/11/2022] [Indexed: 11/16/2022] Open
Abstract
Background and Aim: The Samosir goat has a high cultural value and is a source of germplasm in Indonesia. This study aimed to reveal the history and selection signatures of the Samosir goat. Materials and Methods: A total of 25 goats were divided into seven subpopulations of Indonesian goat breeds. Deoxyribonucleic acid (DNA) from blood samples was isolated with the use of the gSYNC™ DNA Mini Kit (Geneaid, Taipei, Taiwan). Cytb gene amplification was performed by the polymerase chain reaction (PCR) method, and the PCR products were sequenced. A phylogenetic tree was constructed by the neighbor-joining method using MEGA 11 software. A questionnaire was used to collect information related to the history and breeding practices of the Samosir goat on Samosir Island. Results: Samosir goats are divided into four groups based on their coat color: Completely white, white with brown spots, white with black spots, and white with brown and black spots. The body form of the Samosir goat is similar to that of the Kacang goat. The space below a traditional Toba Batak house is used as a goat pen. The genetic difference between the Samosir goat and the Kacang goat based on the Cytb gene was approximately 0.1%. Conclusion: Phylogenetic analysis between Samosir goats and other indigenous Indonesian goats revealed that Samosir goats form a single clade, with a very close genetic distance from other local goats, such as the Kacang goat. The Toba Batak culture on Samosir Island has significantly influenced the selection and formation of the Samosir goat breed.
Collapse
Affiliation(s)
- Suhendra Pakpahan
- Museum Zoologicum Bogoriense, Research Center for Applied Zoology, National Research and Innovation Agency (BRIN), Cibinong, West Java, Indonesia
| | - Rini Widayanti
- Department of Biochemistry and Molecular Biology, Faculty of Veterinary Medicine, Universitas Gadjah Mada, Yogyakarta, Indonesia
| | - Wayan T. Artama
- Department of Biochemistry and Molecular Biology, Faculty of Veterinary Medicine, Universitas Gadjah Mada, Yogyakarta, Indonesia
| |
Collapse
|
15
|
Graber JK, Signer-Hasler H, Burren A, Drögemüller C. Evaluation of truncating variants in the LCORL gene in relation to body size of goats from Switzerland. Anim Genet 2022; 53:237-239. [PMID: 35118693 PMCID: PMC9305788 DOI: 10.1111/age.13177] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2022] [Accepted: 01/18/2022] [Indexed: 02/05/2023]
Affiliation(s)
| | - Heidi Signer-Hasler
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
| | - Alexander Burren
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
| | - Cord Drögemüller
- Institute of Genetics, Vetsuisse Faculty, University of Bern, Bern, Switzerland
| |
Collapse
|
16
|
Signer-Hasler H, Henkel J, Bangerter E, Bulut Z, Drögemüller C, Leeb T, Flury C. Runs of homozygosity in Swiss goats reveal genetic changes associated with domestication and modern selection. Genet Sel Evol 2022; 54:6. [PMID: 35073837 PMCID: PMC8785455 DOI: 10.1186/s12711-022-00695-w] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2021] [Accepted: 01/06/2022] [Indexed: 11/16/2022] Open
Abstract
Background The domestication of goat (Capra hircus) started 11,000 years ago in the fertile crescent. Breed formation in the nineteenth century, establishment of herd books, and selection for specific traits resulted in 10 modern goat breeds in Switzerland. We analyzed whole-genome sequencing (WGS) data from 217 modern goats and nine wild Bezoar goats (Capra aegagrus). After quality control, 27,728,288 biallelic single nucleotide variants (SNVs) were used for the identification of runs of homozygosity (ROH) and the detection of ROH islands. Results Across the 226 caprine genomes from 11 populations, we detected 344 ROH islands that harbor 1220 annotated genes. We compared the ROH islands between the modern breeds and the Bezoar goats. As a proof of principle, we confirmed a signature of selection, which contains the ASIP gene that controls several breed-specific coat color patterns. In two other ROH islands, we identified two missense variants, STC1:p.Lys139Arg and TSHR:p.Ala239Thr, which might represent causative functional variants for domestication signatures. Conclusions We have shown that the information from ROH islands using WGS data is suitable for the analysis of signatures of selection and allowed the detection of protein coding variants that may have conferred beneficial phenotypes during goat domestication. We hypothesize that the TSHR:p.Ala239Thr variant may have played a role in changing the seasonality of reproduction in modern domesticated goats. The exact functional significance of the STC1:p.Lys139Arg variant remains unclear and requires further investigation. Nonetheless, STC1 might represent a new domestication gene affecting relevant traits such as body size and/or milk yield in goats. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-022-00695-w.
Collapse
Affiliation(s)
- Heidi Signer-Hasler
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, 3052, Zollikofen, Switzerland.
| | - Jan Henkel
- Institute of Genetics, Vetsuisse Faculty, University of Bern, 3001, Bern, Switzerland
| | - Erika Bangerter
- Swiss Goat Breeding Association SZZV, Schützenstrasse 10, 3052, Zollikofen, Switzerland
| | - Zafer Bulut
- Department of Biochemistry, Faculty of Veterinary Medicine, Selcuk University, Konya, Turkey
| | | | - Cord Drögemüller
- Institute of Genetics, Vetsuisse Faculty, University of Bern, 3001, Bern, Switzerland
| | - Tosso Leeb
- Institute of Genetics, Vetsuisse Faculty, University of Bern, 3001, Bern, Switzerland
| | - Christine Flury
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, 3052, Zollikofen, Switzerland
| |
Collapse
|
17
|
The SNP-Based Profiling of Montecristo Feral Goat Populations Reveals a History of Isolation, Bottlenecks, and the Effects of Management. Genes (Basel) 2022; 13:genes13020213. [PMID: 35205258 PMCID: PMC8872249 DOI: 10.3390/genes13020213] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2021] [Revised: 01/17/2022] [Accepted: 01/17/2022] [Indexed: 02/04/2023] Open
Abstract
The Montecristo wild goat is an endangered feral population that has been on the homonymous island in the Tuscan Archipelago since ancient times. The origins of Montecristo goats are still debated, with authors dating their introduction either back to Neolithic times or between the 6th and 13th century of the Common Era. To investigate the evolutionary history and relationships of this population we assembled a 50K SNP dataset including 55 Mediterranean breeds and two nuclei of Montecristo goats sampled on the island and from an ex situ conservation project. Diversity levels, gene flow, population structure, and genetic relationships were assessed through multiple approaches. The insular population scored the lowest values of both observed and expected heterozygosity, highlighting reduced genetic variation, while the ex situ nucleus highlighted a less severe reduction. Multivariate statistics, network, and population structure analyses clearly separated the insular nucleus from all other breeds, including the population of Montecristo goats from the mainland. Moreover, admixture and gene flow analyses pinpointed possible genetic inputs received by the two Montecristo goat nuclei from different sources, while Runs of Homozygosity (ROHs) indicated an ancient bottleneck/founder effect in the insular population and recent extensive inbreeding in the ex situ one. Overall, our results suggest that Montecristo goats experienced several demographic fluctuations combined with admixture events over time and highlighted a noticeable differentiation between the two nuclei.
Collapse
|
18
|
Ablondi M, Sabbioni A, Stocco G, Cipolat-Gotet C, Dadousis C, van Kaam JT, Finocchiaro R, Summer A. Genetic Diversity in the Italian Holstein Dairy Cattle Based on Pedigree and SNP Data Prior and After Genomic Selection. Front Vet Sci 2022; 8:773985. [PMID: 35097040 PMCID: PMC8792952 DOI: 10.3389/fvets.2021.773985] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2021] [Accepted: 11/30/2021] [Indexed: 01/09/2023] Open
Abstract
Genetic diversity has become an urgent matter not only in small local breeds but also in more specialized ones. While the use of genomic data in livestock breeding programs increased genetic gain, there is increasing evidence that this benefit may be counterbalanced by the potential loss of genetic variability. Thus, in this study, we aimed to investigate the genetic diversity in the Italian Holstein dairy cattle using pedigree and genomic data from cows born between 2002 and 2020. We estimated variation in inbreeding, effective population size, and generation interval and compared those aspects prior to and after the introduction of genomic selection in the breed. The dataset contained 84,443 single-nucleotide polymorphisms (SNPs), and 74,485 cows were analyzed. Pedigree depth based on complete generation equivalent was equal to 10.67. A run of homozygosity (ROH) analysis was adopted to estimate SNP-based inbreeding (FROH). The average pedigree inbreeding was 0.07, while the average FROH was more than double, being equal to 0.17. The pattern of the effective population size based on pedigree and SNP data was similar although different in scale, with a constant decrease within the last five generations. The overall inbreeding rate (ΔF) per year was equal to +0.27% and +0.44% for Fped and FROH throughout the studied period, which corresponded to about +1.35% and +2.2% per generation, respectively. A significant increase in the ΔF was found since the introduction of genomic selection in the breed. This study in the Italian Holstein dairy cattle showed the importance of controlling the loss of genetic diversity to ensure the long-term sustainability of this breed, as well as to guarantee future market demands.
Collapse
Affiliation(s)
- Michela Ablondi
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| | - Alberto Sabbioni
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| | - Giorgia Stocco
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| | - Claudio Cipolat-Gotet
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
- *Correspondence: Claudio Cipolat-Gotet
| | - Christos Dadousis
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| | - Jan-Thijs van Kaam
- Associazione Nazionale Allevatori della Razza Frisona Bruna e Jersey Italiana, Cremona, Italy
| | - Raffaella Finocchiaro
- Associazione Nazionale Allevatori della Razza Frisona Bruna e Jersey Italiana, Cremona, Italy
| | - Andrea Summer
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| |
Collapse
|
19
|
Schiavo G, Bovo S, Ribani A, Moscatelli G, Bonacini M, Prandi M, Mancin E, Mantovani R, Dall'Olio S, Fontanesi L. Comparative analysis of inbreeding parameters and runs of homozygosity islands in 2 Italian autochthonous cattle breeds mainly raised in the Parmigiano-Reggiano cheese production region. J Dairy Sci 2021; 105:2408-2425. [PMID: 34955250 DOI: 10.3168/jds.2021-20915] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2021] [Accepted: 10/25/2021] [Indexed: 01/19/2023]
Abstract
Reggiana and Modenese are autochthonous cattle breeds, reared in the North of Italy, that can be mainly distinguished for their standard coat color (Reggiana is red, whereas Modenese is white with some pale gray shades). Almost all milk produced by these breeds is transformed into 2 mono-breed branded Parmigiano-Reggiano cheeses, from which farmers receive the economic incomes needed for the sustainable conservation of these animal genetic resources. After the setting up of their herd books in 1960s, these breeds experienced a strong reduction in the population size that was subsequently reverted starting in the 1990s (Reggiana) or more recently (Modenese) reaching at present a total of about 2,800 and 500 registered cows, respectively. Due to the small population size of these breeds, inbreeding is a very important cause of concern for their conservation programs. Inbreeding is traditionally estimated using pedigree data, which are summarized in an inbreeding coefficient calculated at the individual level (FPED). However, incompleteness of pedigree information and registration errors can affect the effectiveness of conservation strategies. High-throughput SNP genotyping platforms allow investigation of inbreeding using genome information that can overcome the limits of pedigree data. Several approaches have been proposed to estimate genomic inbreeding, with the use of runs of homozygosity (ROH) considered to be the more appropriate. In this study, several pedigree and genomic inbreeding parameters, calculated using the whole herd book populations or considering genotyping information (GeneSeek GGP Bovine 150K) from 1,684 Reggiana cattle and 323 Modenese cattle, were compared. Average inbreeding values per year were used to calculate effective population size. Reggiana breed had generally lower genomic inbreeding values than Modenese breed. The low correlation between pedigree-based and genomic-based parameters (ranging from 0.187 to 0.195 and 0.319 to 0.323 in the Reggiana and Modenese breeds, respectively) reflected the common problems of local populations in which pedigree records are not complete. The high proportion of short ROH over the total number of ROH indicates no major recent inbreeding events in both breeds. ROH islands spread over the genome of the 2 breeds (15 in Reggiana and 14 in Modenese) identified several signatures of selection. Some of these included genes affecting milk production traits, stature, body conformation traits (with a main ROH island in both breeds on BTA6 containing the ABCG2, NCAPG, and LCORL genes) and coat color (on BTA13 in Modenese containing the ASIP gene). In conclusion, this work provides an extensive comparative analysis of pedigree and genomic inbreeding parameters and relevant genomic information that will be useful in the conservation strategies of these 2 iconic local cattle breeds.
Collapse
Affiliation(s)
- Giuseppina Schiavo
- Department of Agricultural and Food Sciences, Division of Animal Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| | - Samuele Bovo
- Department of Agricultural and Food Sciences, Division of Animal Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| | - Anisa Ribani
- Department of Agricultural and Food Sciences, Division of Animal Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| | - Giulia Moscatelli
- Department of Agricultural and Food Sciences, Division of Animal Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| | - Massimo Bonacini
- Associazione Nazionale Allevatori Bovini di Razza Reggiana (ANABORARE), Via Masaccio 11, 42124 Reggio Emilia, Italy
| | - Marco Prandi
- Associazione Nazionale Allevatori Bovini di Razza Reggiana (ANABORARE), Via Masaccio 11, 42124 Reggio Emilia, Italy
| | - Enrico Mancin
- Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padova, Viale dell'Università 16, 35020 Legnaro (PD), Italy
| | - Roberto Mantovani
- Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padova, Viale dell'Università 16, 35020 Legnaro (PD), Italy
| | - Stefania Dall'Olio
- Department of Agricultural and Food Sciences, Division of Animal Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| | - Luca Fontanesi
- Department of Agricultural and Food Sciences, Division of Animal Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy.
| |
Collapse
|
20
|
Henkel J, Dubacher A, Bangerter E, Herren U, Ammann P, Drögemüller C, Flury C, Leeb T. Introgression of ASIP and TYRP1 Alleles Explains Coat Color Variation in Valais Goats. J Hered 2021; 112:452-457. [PMID: 34050662 DOI: 10.1093/jhered/esab024] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2020] [Accepted: 04/21/2021] [Indexed: 11/13/2022] Open
Abstract
The Valais Blackneck goat is a Swiss goat breed with a characteristic coat color phenotype. Before the revision of the breed standard in 1938, 4 different color varieties of Valais goats were known. Besides Blackneck animals resembling the modern breed standard, the brown and white Copperneck goat, the white Capra Sempione, and the greyish Grüenochte comprised the historic Valais goats. The brown pigmentation of Copperneck goats had previously been traced back to an introgression of a mutant TYRP1 allele from Toggenburg goats. In the present study, we identified additional introgression events of distinct ASIP alleles causing the remaining 2 rare coat color patterns within the Valais Blackneck goat breed. We identified the introgression of the AWt allele from Appenzell or Saanen goats in white Capra Sempione goats. Similarly, introgression of the Apc allele from Peacock goats resulted in the greyish Grüenochte phenotype. These results demonstrate past hybridization events between breeds that are separated today. A perfect genotype-phenotype association in 393 Valais goats supported the causality of the genotyped variants for the different coat color phenotypes. Our study gives insights into the introgression of functionally relevant copy number variant (CNV) alleles controlling pigmentation between goat breeds with strikingly different coat color patterns.
Collapse
Affiliation(s)
- Jan Henkel
- Institute of Genetics, Vetsuisse Faculty, University of Bern, Bremgartenstrasse 109a, 3001 Bern, Switzerland.,DermFocus, University of Bern, 3001 Bern, Switzerland
| | - Alexandra Dubacher
- Institute of Genetics, Vetsuisse Faculty, University of Bern, Bremgartenstrasse 109a, 3001 Bern, Switzerland
| | - Erika Bangerter
- Swiss Goat Breeding Association, 3052 Zollikofen, Switzerland
| | - Ursula Herren
- Swiss Goat Breeding Association, 3052 Zollikofen, Switzerland
| | | | - Cord Drögemüller
- Institute of Genetics, Vetsuisse Faculty, University of Bern, Bremgartenstrasse 109a, 3001 Bern, Switzerland.,DermFocus, University of Bern, 3001 Bern, Switzerland
| | - Christine Flury
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, 3052 Zollikofen, Switzerland
| | - Tosso Leeb
- Institute of Genetics, Vetsuisse Faculty, University of Bern, Bremgartenstrasse 109a, 3001 Bern, Switzerland.,DermFocus, University of Bern, 3001 Bern, Switzerland
| |
Collapse
|
21
|
Abousoliman I, Reyer H, Oster M, Murani E, Mohamed I, Wimmers K. Genome-Wide Analysis for Early Growth-Related Traits of the Locally Adapted Egyptian Barki Sheep. Genes (Basel) 2021; 12:1243. [PMID: 34440417 PMCID: PMC8394750 DOI: 10.3390/genes12081243] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2021] [Revised: 08/09/2021] [Accepted: 08/11/2021] [Indexed: 12/28/2022] Open
Abstract
Sheep play a critical role in the agricultural and livestock sector in Egypt. For sheep meat production, growth traits such as birth and weaning weights are very important and determine the supply and income of local farmers. The Barki sheep originates from the northeastern coastal zone of Africa, and due to its good adaptation to the harsh environmental conditions, it contributes significantly to the meat production in these semi-arid regions. This study aimed to use a genome-wide SNP panel to identify genomic regions that are diversified between groups of individuals of Egyptian Barki sheep with high and low growth performance traits. In this context, from a phenotyped population of 140 lambs of Barki sheep, 69 lambs were considered for a genome-wide scan with the Illumina OvineSNP50 V2 BeadChip. The selected lambs were grouped into divergent subsets with significantly different performance for birth weight and weaning weight. After quality control, 63 animals and 40,383 SNPs were used for analysis. The fixation index (FST) for each SNP was calculated between the groups. The results verified genomic regions harboring some previously proposed candidate genes for traits related to body growth, i.e., EYA2, GDF2, GDF10, MEF2B, SLC16A7, TBX15, TFAP2B, and TNNC2. Moreover, novel candidate genes were proposed with known functional implications on growth processes such as CPXM2 and LRIG3. Subsequent association analysis showed significant effects of the considered SNPs on birth and weaning weights. Results highlight the genetic diversity associated with performance traits and thus the potential to improve growth traits in the Barki sheep breed.
Collapse
Affiliation(s)
- Ibrahim Abousoliman
- Leibniz Institute for Farm Animal Biology, Wilhelm-Stahl-Allee 2, 18196 Dummerstorf, Germany; (I.A.); (H.R.); (M.O.); (E.M.)
- Desert Research Center, Department of Animal and Poultry Breeding, 1 Mathaf El-Matareya St., El-Matareya, Cairo 11753, Egypt;
| | - Henry Reyer
- Leibniz Institute for Farm Animal Biology, Wilhelm-Stahl-Allee 2, 18196 Dummerstorf, Germany; (I.A.); (H.R.); (M.O.); (E.M.)
| | - Michael Oster
- Leibniz Institute for Farm Animal Biology, Wilhelm-Stahl-Allee 2, 18196 Dummerstorf, Germany; (I.A.); (H.R.); (M.O.); (E.M.)
| | - Eduard Murani
- Leibniz Institute for Farm Animal Biology, Wilhelm-Stahl-Allee 2, 18196 Dummerstorf, Germany; (I.A.); (H.R.); (M.O.); (E.M.)
| | - Ismail Mohamed
- Desert Research Center, Department of Animal and Poultry Breeding, 1 Mathaf El-Matareya St., El-Matareya, Cairo 11753, Egypt;
| | - Klaus Wimmers
- Leibniz Institute for Farm Animal Biology, Wilhelm-Stahl-Allee 2, 18196 Dummerstorf, Germany; (I.A.); (H.R.); (M.O.); (E.M.)
- Faculty of Agricultural and Environmental Sciences, University of Rostock, Justus-von-Liebig-Weg 7, 18059 Rostock, Germany
| |
Collapse
|
22
|
Pogorevc N, Simčič M, Khayatzadeh N, Sölkner J, Berger B, Bojkovski D, Zorc M, Dovč P, Medugorac I, Horvat S. Post-genotyping optimization of dataset formation could affect genetic diversity parameters: an example of analyses with alpine goat breeds. BMC Genomics 2021; 22:546. [PMID: 34273960 PMCID: PMC8285797 DOI: 10.1186/s12864-021-07802-z] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2020] [Accepted: 06/13/2021] [Indexed: 12/05/2022] Open
Abstract
Background Local breeds retained unique genetic variability important for adaptive potential especially in light of challenges related to climate change. Our first objective was to perform, for the first time, a genome-wide diversity characterization using Illumina GoatSNP50 BeadChip of autochthonous Drežnica goat breed from Slovenia, and five and one local breeds from neighboring Austria and Italy, respectively. For optimal conservation and breeding programs of endangered local breeds, it is important to detect past admixture events and strive for preservation of purebred representatives of each breed with low or without admixture. In the second objective, we hence investigated the effect of inclusion or exclusion of outliers from datasets on genetic diversity and population structure parameters. Results Distinct genetic origin of the Drežnica goat was demonstrated as having closest nodes to Austrian and Italian breeds. A phylogenetic study of these breeds with other goat breeds having SNP data available in the DRYAD repository positioned them in the alpine, European and global context. Swiss breeds clustered with cosmopolitan alpine breeds and were closer to French and Spanish breeds. On the other hand, the Drežnica goat, Austrian and Italian breeds were closer to Turkish breeds. Datasets where outliers were excluded affected estimates of genetic diversity parameters within the breed and increased the pairwise genetic distances between most of the breeds. Alpine breeds, including Drežnica, Austrian and Italian goats analyzed here, still exhibit relatively high levels of genetic variability, homogeneous genetic structure and strong geographical partitioning. Conclusions Genetic diversity analyses revealed that the Slovenian Drežnica goat has a distinct genetic identity and is closely related to the neighboring Austrian and Italian alpine breeds. These results expand our knowledge on phylogeny of goat breeds from easternmost part of the European Alps. The here employed outlier test and datasets optimization approaches provided an objective and statistically powerful tool for removal of admixed outliers. Importance of this test in selecting the representatives of each breed is warranted to obtain more objective diversity parameters and phylogenetic analysis. Such parameters are often the basis of breeding and management programs and are therefore important for preserving genetic variability and uniqueness of local rare breeds. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07802-z.
Collapse
Affiliation(s)
- Neža Pogorevc
- Department of Animal science, Biotechnical Faculty, University of Ljubljana, Jamnikarjeva 101, SI-1000, Ljubljana, Slovenia
| | - Mojca Simčič
- Department of Animal science, Biotechnical Faculty, University of Ljubljana, Jamnikarjeva 101, SI-1000, Ljubljana, Slovenia
| | - Negar Khayatzadeh
- Division of Livestock Science, Department of Sustainable Agricultural Systems, University of Natural Resources and Life Sciences Vienna, Gregor Mendel Str. 33, A-1180, Vienna, Austria
| | - Johann Sölkner
- Division of Livestock Science, Department of Sustainable Agricultural Systems, University of Natural Resources and Life Sciences Vienna, Gregor Mendel Str. 33, A-1180, Vienna, Austria
| | - Beate Berger
- Department Animal Genetic Resources, AREC Raumberg-Gumpenstein, Institute of Organic Farming and Biodiversity of Farm Animals, 4601 Thalheim b., Wels, Austria
| | - Danijela Bojkovski
- Department of Animal science, Biotechnical Faculty, University of Ljubljana, Jamnikarjeva 101, SI-1000, Ljubljana, Slovenia
| | - Minja Zorc
- Department of Animal science, Biotechnical Faculty, University of Ljubljana, Jamnikarjeva 101, SI-1000, Ljubljana, Slovenia
| | - Peter Dovč
- Department of Animal science, Biotechnical Faculty, University of Ljubljana, Jamnikarjeva 101, SI-1000, Ljubljana, Slovenia
| | - Ivica Medugorac
- Population Genomics Group, Department of Veterinary Sciences, Faculty of Veterinary Medicine, Ludwig-Maximilians-University Munich, Lena-Christ-Straβe 48, 8215, Martinsried/Planegg, Germany
| | - Simon Horvat
- Department of Animal science, Biotechnical Faculty, University of Ljubljana, Jamnikarjeva 101, SI-1000, Ljubljana, Slovenia.
| |
Collapse
|
23
|
An Overview of the Use of Genotyping Techniques for Assessing Genetic Diversity in Local Farm Animal Breeds. Animals (Basel) 2021; 11:ani11072016. [PMID: 34359144 PMCID: PMC8300386 DOI: 10.3390/ani11072016] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2021] [Revised: 06/30/2021] [Accepted: 07/01/2021] [Indexed: 12/18/2022] Open
Abstract
Simple Summary The number of local farm animal breeds is declining worldwide. However, these breeds have different degrees of genetic diversity. Measuring genetic diversity is important for the development of conservation strategies and, therefore, various genomic analysis techniques are available. The aim of the present work was to shed light on the use of these techniques in diversity studies of local breeds. In summary, a total of 133 worldwide studies that examined genetic diversity in local cattle, sheep, goat, chicken and pig breeds were reviewed. The results show that over time, almost all available genomic techniques were used and various diversity parameters were calculated. Therefore, the present results provide a comprehensive overview of the application of these techniques in the field of local breeds. This can provide helpful insights into the advancement of the conservation of breeds with high genetic diversity. Abstract Globally, many local farm animal breeds are threatened with extinction. However, these breeds contribute to the high amount of genetic diversity required to combat unforeseen future challenges of livestock production systems. To assess genetic diversity, various genotyping techniques have been developed. Based on the respective genomic information, different parameters, e.g., heterozygosity, allele frequencies and inbreeding coefficient, can be measured in order to reveal genetic diversity between and within breeds. The aim of the present work was to shed light on the use of genotyping techniques in the field of local farm animal breeds. Therefore, a total of 133 studies across the world that examined genetic diversity in local cattle, sheep, goat, chicken and pig breeds were reviewed. The results show that diversity of cattle was most often investigated with microsatellite use as the main technique. Furthermore, a large variety of diversity parameters that were calculated with different programs were identified. For 15% of the included studies, the used genotypes are publicly available, and, in 6%, phenotypes were recorded. In conclusion, the present results provide a comprehensive overview of the application of genotyping techniques in the field of local breeds. This can provide helpful insights to advance the conservation of breeds.
Collapse
|
24
|
Deniskova TE, Dotsev AV, Selionova MI, Reyer H, Sölkner J, Fornara MS, Aybazov AMM, Wimmers K, Brem G, Zinovieva NA. SNP-Based Genotyping Provides Insight Into the West Asian Origin of Russian Local Goats. Front Genet 2021; 12:708740. [PMID: 34276802 PMCID: PMC8282346 DOI: 10.3389/fgene.2021.708740] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2021] [Accepted: 06/07/2021] [Indexed: 11/13/2022] Open
Abstract
Specific local environmental and sociocultural conditions have led to the creation of various goat populations in Russia. National goat diversity includes breeds that have been selected for down and mohair production traits as well as versatile local breeds for which pastoralism is the main management system. Effective preservation and breeding programs for local goat breeds are missing due to the lack of DNA-based data. In this work, we analyzed the genetic diversity and population structure of Russian local goats, including Altai Mountain, Altai White Downy, Dagestan Downy, Dagestan Local, Karachaev, Orenburg, and Soviet Mohair goats, which were genotyped with the Illumina Goat SNP50 BeadChip. In addition, we addressed genetic relationships between local and global goat populations obtained from the AdaptMap project. Russian goats showed a high level of genetic diversity. Although a decrease in historical effective population sizes was revealed, the recent effective population sizes estimated for three generations ago were larger than 100 in all studied populations. The mean runs of homozygosity (ROH) lengths ranged from 79.42 to 183.94 Mb, and the average ROH number varied from 18 to 41. Short ROH segments (<2 Mb) were predominant in all breeds, while the longest ROH class (>16 Mb) was the least frequent. Principal component analysis, Neighbor-Net graph, and Admixture clustering revealed several patterns in Russian local goats. First, a separation of the Karachaev breed from other populations was observed. Moreover, genetic connections between the Orenburg and Altai Mountain breeds were suggested and the Dagestan breeds were found to be admixed with the Soviet Mohair breed. Neighbor-Net analysis and clustering of local and global breeds demonstrated the close genetic relations between Russian local and Turkish breeds that probably resulted from past admixture events through postdomestication routes. Our findings contribute to the understanding of the genetic relationships of goats originating in West Asia and Eurasia and may be used to design breeding programs for local goats to ensure their effective conservation and proper management.
Collapse
Affiliation(s)
| | - Arsen V Dotsev
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, Russia
| | - Marina I Selionova
- Russian State Agrarian University - Moscow Timiryazev Agricultural Academy, Moscow, Russia
| | - Henry Reyer
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology, Dummerstorf, Germany
| | - Johann Sölkner
- Division of Livestock Sciences, University of Natural Resources and Life Sciences, Vienna, Vienna, Austria
| | | | - Ali-Magomed M Aybazov
- All-Russian Research Institute of Sheep and Goat Breeding - Branch of the Federal State Budgetary Scientific Institution, North Caucasian Agrarian Center, Stavropol, Russia
| | - Klaus Wimmers
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology, Dummerstorf, Germany
| | - Gottfried Brem
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, Russia.,Institute of Animal Breeding and Genetics, University of Veterinary Medicine Vienna, Vienna, Austria
| | | |
Collapse
|
25
|
Addo S, Klingel S, Thaller G, Hinrichs D. Genetic diversity and the application of runs of homozygosity-based methods for inbreeding estimation in German White-headed Mutton sheep. PLoS One 2021; 16:e0250608. [PMID: 33956807 PMCID: PMC8101715 DOI: 10.1371/journal.pone.0250608] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2020] [Accepted: 04/09/2021] [Indexed: 11/19/2022] Open
Abstract
The German White-headed Mutton (GWM) sheep is a monitoring population believed to have been improved through crosses with other breeds, e.g., Texel (TXL) and French Berrichone du Cher (BDC). The primary aim of the study was to analyse genetic diversity and breed composition of GWM sheep. Furthermore, different measures of computing inbreeding from the runs of homozygosity (ROH) were investigated. Data for GWM consisted of pedigree information on 19,000 animals and 40,753 quality filtered SNPs on 46 individuals. Additionally, publicly available genotype data on 209 individuals belonging to nine sheep breeds were included in the analysis. Due to evenness of SNPs spacing and proportionality of the number of SNPs in each autosome to autosome length, a high correlation (rp = 0.99) was found between genomic inbreeding coefficients computed based on the length of ROH (FROH_L) and those computed relative to the number of SNPs in ROH (FROH_N). Total inbreeding was partitioned into values for individual chromosomes revealing the highest levels of inbreeding on chromosomes 1, 2 and 3. Correlations between the ROH-based inbreeding measures and pedigree inbreeding reached 0.82. The observed heterozygosity estimate in GWM was high (0.39), however, the breed suffered low level of effective population size (~50) from a genomic viewpoint. Moreover, effective number of founders (186), and effective number of ancestors (144) implied disequilibrium of founder contribution and a genetic bottleneck in the breed. Multidimensional scaling and network visualisation analyses revealed close connectedness of GWM to BDC and German Texel (GTX). A model-based admixture analysis consistently indicated the flow of genes from other breeds, particularly BDC to GWM. Our analyses highlight the mixed genetic background of GWM sheep and furthermore, suggest a close monitoring of the breed to consolidate its genetic diversity while averting further reduction in the effective population size.
Collapse
Affiliation(s)
- Sowah Addo
- Institute of Animal Breeding and Husbandry, Kiel University, Kiel, Germany
- Department of Animal Breeding, University of Kassel, Witzenhausen, Germany
- * E-mail:
| | - Stefanie Klingel
- Arche Warder, Center for Rare and Endangered Domestic Animals, Warder, Germany
| | - Georg Thaller
- Institute of Animal Breeding and Husbandry, Kiel University, Kiel, Germany
| | - Dirk Hinrichs
- Department of Animal Breeding, University of Kassel, Witzenhausen, Germany
| |
Collapse
|
26
|
Abstract
Runs of homozygosity (ROH) are chromosomal stretches that in a diploid genome appear in a homozygous state and display identical alleles at multiple contiguous loci. This study aimed to systematically compare the genomic distribution of the ROH islands among five populations of wild vs commercial chickens of both layer and broiler type. To this end, we analyzed whole genome sequences of 115 birds including white layer (WL, n = 25), brown layer (BL, n = 25), broiler line A (BRA, n = 20), broiler line B (BRB, n = 20) and Red Junglefowl (RJF, n = 25). The ROH segments varied in size markedly among populations, ranging from 0.3 to 21.83 Mb reflecting their past genealogy. White layers contained the largest portion of the genome in homozygous state with an average ROH length of 432.1 Mb (±18.7) per bird, despite carrying it in short segments (0.3-1 Mb). Population-wise inbreeding measures based on Wright's (Fis) and genomic (FROH) metrics revealed highly inbred genome of layer lines relative to the broilers and Red Junglefowl. We further revealed the ROH islands, among commercial lines overlapped with QTL related to limb development (GREM1, MEOX2), body weight (Meis2a.1, uc_338), eggshell color (GLCCI1, ICA1, UMAD1), antibody response to Newcastle virus (ROBO2), and feather pecking. Comparison of ROH landscape in sequencing resolution demonstrated that a sizable portion of genome of commercial lines segregates in homozygote state, reflecting many generations of assortative mating and intensive selection in their recent history. In contrary, wild birds carry shorter ROH segments, likely suggestive of older evolutionary events.
Collapse
|
27
|
Conservation Assessment of the State Goat Farms by Using SNP Genotyping Data. Genes (Basel) 2020; 11:genes11060652. [PMID: 32545749 PMCID: PMC7349881 DOI: 10.3390/genes11060652] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2020] [Revised: 06/09/2020] [Accepted: 06/11/2020] [Indexed: 01/10/2023] Open
Abstract
Conservation of genetic resources is of great concern globally to maintain genetic diversity for sustainable food security. Comprehensive identification of the breed composition, estimation of inbreeding and effective population size are essential for the effective management of farm animal genetic resources and to prevent the animals from genetic erosion. The Zhongwei male (ZWM), Arbas Cashmere male (ACM) and Jining Grey male (JGM) goats are conserved in three different state goat farms in China but their family information, level of inbreeding and effective population size are unknown. We investigated the genomic relationship, inbreeding coefficient and effective population size in these three breeds from three state goat farms using the Illumina goat SNP50 BeadChip. Genomic relationships and phylogenetic analysis revealed that the breeds are clearly separated and formed separate clusters based on their genetic relationship. We obtained a high proportion of informative SNPs, ranging from 91.8% in the Arbas Cashmere male to 96.2% in the Jining Grey male goat breeds with an average mean of 96.8%. Inbreeding, as measured by FROH, ranged from 1.79% in ZWM to 8.62% in ACM goat populations. High FROH values, elevated genomic coverage of very long ROH (>30 Mb) and severe decline in effective population size were recorded in ACM goat farm. The existence of a high correlation between FHOM and FROH indicates that FROH can be used as an alternative to inbreeding estimates in the absence of pedigree records. The Ne estimates 13 generations ago were 166, 69 and 79 for ZWM, ACM and JGM goat farm, respectively indicating that these goat breeds were strongly affected by selection pressure or genetic drift. This study provides insight into the genomic relationship, levels of inbreeding and effective population size in the studied goat populations conserved in the state goat farms which will be valuable in prioritizing populations for conservation and for developing suitable management practices for further genetic improvement of these Chinese male goats.
Collapse
|
28
|
Henkel J, Saif R, Jagannathan V, Schmocker C, Zeindler F, Bangerter E, Herren U, Posantzis D, Bulut Z, Ammann P, Drögemüller C, Flury C, Leeb T. Selection signatures in goats reveal copy number variants underlying breed-defining coat color phenotypes. PLoS Genet 2019; 15:e1008536. [PMID: 31841508 PMCID: PMC6936872 DOI: 10.1371/journal.pgen.1008536] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2019] [Revised: 12/30/2019] [Accepted: 11/23/2019] [Indexed: 12/16/2022] Open
Abstract
Domestication and human selection have formed diverse goat breeds with characteristic phenotypes. This process correlated with the fixation of causative genetic variants controlling breed-specific traits within regions of reduced genetic diversity, so called selection signatures or selective sweeps. Using whole genome sequencing of DNA pools (pool-seq) from 20 genetically diverse modern goat breeds and bezoars, we identified 2,239 putative selection signatures. In two Pakistani goat breeds, Pak Angora and Barbari, we found selection signatures in a region harboring KIT, a gene involved in melanoblast development, migration, and survival. The search for candidate causative variants responsible for these selective sweeps revealed two different copy number variants (CNVs) downstream of KIT that were exclusively present in white Pak Angora and white-spotted Barbari goats. Several Swiss goat breeds selected for specific coat colors showed selection signatures at the ASIP locus encoding the agouti signaling protein. Analysis of these selective sweeps revealed four different CNVs associated with the white or tan (AWt), Swiss markings (Asm), badgerface (Ab), and the newly proposed peacock (Apc) allele. RNA-seq analyses on skin samples from goats with the different CNV alleles suggest that the identified structural variants lead to an altered expression of ASIP between eumelanistic and pheomelanistic body areas. Our study yields novel insights into the genetic control of pigmentation by identifying six functionally relevant CNVs. It illustrates how structural changes of the genome have contributed to phenotypic evolution in domestic goats. Domestic animals have been selected for hundreds or sometimes even thousands of years for traits that were appreciated by their human owners. This process correlated with the fixation of causative genetic variants controlling breed-specific traits within regions of reduced genetic diversity, so called selection signatures or selective sweeps. We conducted a comprehensive screen for selection signatures in 20 phenotypically and genetically diverse modern goat breeds and identified a total of 2,239 putative selection signatures in our dataset. Follow-up experiments on selection signatures harboring known candidate genes for coat color revealed six different copy number variants (CNVs). Two of these CNVs were located in the 3’-flanking region of KIT and associated with a completely white coat color phenotype in Pak Angora goats and a white-spotted coat color phenotype in Barbari goats, respectively. The other four CNVs were located at the ASIP locus. They were associated with four different types of coat color patterning in seven Swiss goat breeds. Their functional effect is mediated by region-specific quantitative changes in ASIP mRNA expression. Our study illustrates how structural changes of the genome have contributed to phenotypic evolution in domestic goats.
Collapse
Affiliation(s)
- Jan Henkel
- Institute of Genetics, Vetsuisse Faculty, University of Bern, Bern, Switzerland
- DermFocus, University of Bern, Bern, Switzerland
| | - Rashid Saif
- Institute of Genetics, Vetsuisse Faculty, University of Bern, Bern, Switzerland
- Institute of Biotechnology, Gulab Devi Educational Complex, Lahore, Pakistan
| | - Vidhya Jagannathan
- Institute of Genetics, Vetsuisse Faculty, University of Bern, Bern, Switzerland
- DermFocus, University of Bern, Bern, Switzerland
| | - Corinne Schmocker
- Institute of Genetics, Vetsuisse Faculty, University of Bern, Bern, Switzerland
| | - Flurina Zeindler
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
| | | | - Ursula Herren
- Swiss Goat Breeding Association, Zollikofen, Switzerland
| | | | - Zafer Bulut
- Department of Biochemistry, Faculty of Veterinary Medicine, Selcuk University, Konya, Turkey
| | | | - Cord Drögemüller
- Institute of Genetics, Vetsuisse Faculty, University of Bern, Bern, Switzerland
- DermFocus, University of Bern, Bern, Switzerland
| | - Christine Flury
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
| | - Tosso Leeb
- Institute of Genetics, Vetsuisse Faculty, University of Bern, Bern, Switzerland
- DermFocus, University of Bern, Bern, Switzerland
- * E-mail:
| |
Collapse
|
29
|
Michailidou S, Tsangaris GT, Tzora A, Skoufos I, Banos G, Argiriou A, Arsenos G. Analysis of genome-wide DNA arrays reveals the genomic population structure and diversity in autochthonous Greek goat breeds. PLoS One 2019; 14:e0226179. [PMID: 31830089 PMCID: PMC6907847 DOI: 10.1371/journal.pone.0226179] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2018] [Accepted: 11/21/2019] [Indexed: 12/02/2022] Open
Abstract
Goats play an important role in the livestock sector in Greece. The national herd consists mainly of two indigenous breeds, the Eghoria and Skopelos. Here, we report the population structure and genomic profiles of these two native goat breeds using Illumina’s Goat SNP50 BeadChip. Moreover, we present a panel of candidate markers acquired using different genetic models for breed discrimination. Quality control on the initial dataset resulted in 48,841 SNPs kept for downstream analysis. Principal component and admixture analyses were applied to assess population structure. The rate of inbreeding within breed was evaluated based on the distribution of runs of homozygosity in the genome and respective coefficients, the genomic relationship matrix, the patterns of linkage disequilibrium, and the historic effective population size. Results showed that both breeds exhibit high levels of genetic diversity. Level of inbreeding between the two breeds estimated by the Wright’s fixation index FST was low (Fst = 0.04362), indicating the existence of a weak genetic differentiation between them. In addition, grouping of farms according to their geographical locations was observed. This study presents for the first time a genome-based analysis on the genetic structure of the two indigenous Greek goat breeds and identifies markers that can be potentially exploited in future selective breeding programs for traceability purposes, targeted genetic improvement schemes and conservation strategies.
Collapse
Affiliation(s)
- S. Michailidou
- Laboratory of Animal Husbandry, School of Veterinary Medicine, School of Health Sciences, Aristotle University of Thessaloniki, Thessaloniki, Greece
- Institute of Applied Biosciences, Center for Research and Technology Hellas, Thermi, Greece
- * E-mail:
| | - G. Th. Tsangaris
- Proteomics Research Unit, Biomedical Research Foundation of the Academy of Athens, Athens, Greece
| | - A. Tzora
- School of Agriculture, Department of Agriculture, Division of Animal Production, University of Ioannina, Kostakioi Artas, Greece
| | - I. Skoufos
- School of Agriculture, Department of Agriculture, Division of Animal Production, University of Ioannina, Kostakioi Artas, Greece
| | - G. Banos
- Laboratory of Animal Husbandry, School of Veterinary Medicine, School of Health Sciences, Aristotle University of Thessaloniki, Thessaloniki, Greece
- Scotland's Rural College and The Roslin Institute University of Edinburgh, Edinburgh, Scotland, United Kingdom
| | - A. Argiriou
- Institute of Applied Biosciences, Center for Research and Technology Hellas, Thermi, Greece
| | - G. Arsenos
- Laboratory of Animal Husbandry, School of Veterinary Medicine, School of Health Sciences, Aristotle University of Thessaloniki, Thessaloniki, Greece
| |
Collapse
|
30
|
Runs of Homozygosity and NetView analyses provide new insight into the genome-wide diversity and admixture of three German cattle breeds. PLoS One 2019; 14:e0225847. [PMID: 31800604 PMCID: PMC6892555 DOI: 10.1371/journal.pone.0225847] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2019] [Accepted: 11/13/2019] [Indexed: 11/19/2022] Open
Abstract
Angler (RVA) and Red-and-White dual-purpose (RDN) cattle were in the past decades crossed with influential Red Holstein (RH) sires. However, genome-wide diversity studies in these breeds are lacking. The objective of the present study was to elucidate the genome-wide diversity and population structure of the three German cattle breeds. Using 40,851 single nucleotide polymorphism markers scored in 337 individuals, runs of homozygosity (ROH) were analysed in each breed. Clustering and a high-resolution network visualisation analyses were performed on an extended dataset that included 11 additional (outgroup) breeds. Genetic diversity levels were high with observed heterozygosity above 0.35 in all three breeds. Only RVA had a recent past effective population size (Ne) estimate above 100 at 5 generations ago. ROH length distribution followed a similar pattern across breeds and the majority of ROH were found in the length class of >5 to 10 Mb. Estimates of average inbreeding calculated from ROH (FROH) were 0.021 (RVA), 0.045 (RDN) and 0.053 (RH). Moderate to high positive correlations were found between FROH and pedigree inbreeding (FPED) and between FROH and inbreeding derived from the excess of homozygosity (FHOM), while the intercept of the regression of FROH on FPED was above zero. The population structure analysis showed strong evidence of admixture between RVA and RH. Introgression of RDN with RH genes was minimally detected and for the first time, the study uncovered Norwegian Red Cattle ancestry in RVA. Highly heterogeneous genetic background was found for RVA and RH and as expected, the breeds of the extended dataset effectively differentiated mostly based on geographical origin, validating our findings. The results of this study confirm the impact of RH sires on RVA and RDN populations. Furthermore, a close monitoring is suggested to curb further reduction of Ne in the breeds.
Collapse
|
31
|
Muriuki C, Bush SJ, Salavati M, McCulloch ME, Lisowski ZM, Agaba M, Djikeng A, Hume DA, Clark EL. A Mini-Atlas of Gene Expression for the Domestic Goat ( Capra hircus). Front Genet 2019; 10:1080. [PMID: 31749840 PMCID: PMC6844187 DOI: 10.3389/fgene.2019.01080] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2019] [Accepted: 10/09/2019] [Indexed: 12/12/2022] Open
Abstract
Goats (Capra hircus) are an economically important livestock species providing meat and milk across the globe. They are of particular importance in tropical agri-systems contributing to sustainable agriculture, alleviation of poverty, social cohesion, and utilisation of marginal grazing. There are excellent genetic and genomic resources available for goats, including a highly contiguous reference genome (ARS1). However, gene expression information is limited in comparison to other ruminants. To support functional annotation of the genome and comparative transcriptomics, we created a mini-atlas of gene expression for the domestic goat. RNA-Seq analysis of 17 transcriptionally rich tissues and 3 cell-types detected the majority (90%) of predicted protein-coding transcripts and assigned informative gene names to more than 1000 previously unannotated protein-coding genes in the current reference genome for goat (ARS1). Using network-based cluster analysis, we grouped genes according to their expression patterns and assigned those groups of coexpressed genes to specific cell populations or pathways. We describe clusters of genes expressed in the gastro-intestinal tract and provide the expression profiles across tissues of a subset of genes associated with functional traits. Comparative analysis of the goat atlas with the larger sheep gene expression atlas dataset revealed transcriptional similarities between macrophage associated signatures in the sheep and goats sampled in this study. The goat transcriptomic resource complements the large gene expression dataset we have generated for sheep and contributes to the available genomic resources for interpretation of the relationship between genotype and phenotype in small ruminants.
Collapse
Affiliation(s)
- Charity Muriuki
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Edinburgh, United Kingdom
- Centre for Tropical Livestock Genetics and Health (CTLGH), Edinburgh, United Kingdom
| | - Stephen J. Bush
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Edinburgh, United Kingdom
- Nuffield Department of Clinical Medicine, John Radcliffe Hospital, University of Oxford, Oxford, United Kingdom
| | - Mazdak Salavati
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Edinburgh, United Kingdom
- Centre for Tropical Livestock Genetics and Health (CTLGH), Edinburgh, United Kingdom
| | - Mary E.B. McCulloch
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Edinburgh, United Kingdom
| | - Zofia M. Lisowski
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Edinburgh, United Kingdom
| | - Morris Agaba
- Biosciences Eastern and Central Africa - International Livestock Research Institute (BecA - ILRI) Hub, Nairobi, Kenya
| | - Appolinaire Djikeng
- Centre for Tropical Livestock Genetics and Health (CTLGH), Edinburgh, United Kingdom
| | - David A. Hume
- Mater Research Institute-University of Queensland, Woolloongabba, QLD, Australia
| | - Emily L. Clark
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Edinburgh, United Kingdom
- Centre for Tropical Livestock Genetics and Health (CTLGH), Edinburgh, United Kingdom
| |
Collapse
|
32
|
E GX, Duan XH, Zhang JH, Huang YF, Zhao YJ, Na RS, Zhao ZQ, Ma YH, Chu MX, Basang WD, Zhu YB, An TW, Luo XL, Han YG, Zeng Y. Genome-wide selection signatures analysis of litter size in Dazu black goats using single-nucleotide polymorphism. 3 Biotech 2019; 9:336. [PMID: 31475088 DOI: 10.1007/s13205-019-1869-3] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2019] [Accepted: 08/12/2019] [Indexed: 11/24/2022] Open
Abstract
Litter size is considered to be the most important index for estimating domestic animal productivity. The number of indigenous goats in China with higher litter sizes than those of commercial breeds in other countries may be helpful for accelerating genetic improvements in goat breeding. We performed a genome-wide selective sweep analysis of 31 Dazu black goats with extreme standard deviation in litter size within the third fetus to identify significant genomic regions and candidate genes through different analyses. The analysis identified a total of 33,917,703 variants, including 32,262,179 SNPs and 1,655,524 indels. In addition, two novel candidate genes (LRP1B and GLRB), which are related to litter size, were obtained with π, Tajima's D, πA/πB, and F ST at the individual level with a 95% threshold for each parameter. These two genes were annotated in five GO terms (localization, binding, macromolecular complex, membrane part, and membrane) and two pathways (long-term depression and neuroactive ligand-receptor interaction pathway). Regarding the result of linkage disequilibrium (LD) analysis, in LRP1B and GRID2, the high-yield Dazu black goats exhibit significantly different LD patterns from low-yield goats. Litter size variability has low heritability and is related to multiple complex factors found in domestic animals. Obtaining a clear explanation and significant signal by genome-wide selective sweep analysis with a small sample size is difficult. However, we investigated some candidate genes, particularly LRP1B and GLRB, which may provide useful information for further research.
Collapse
Affiliation(s)
- Guang-Xin E
- 1Chongqing Key Laboratory of Forage & Herbivore, Chongqing Engineering Research Centre for Herbivores Resource Protection and Utilization, College of Animal Science and Technology, Southwest University, Chongqing, China
| | - Xing-Hai Duan
- 1Chongqing Key Laboratory of Forage & Herbivore, Chongqing Engineering Research Centre for Herbivores Resource Protection and Utilization, College of Animal Science and Technology, Southwest University, Chongqing, China
| | - Jia-Hua Zhang
- 1Chongqing Key Laboratory of Forage & Herbivore, Chongqing Engineering Research Centre for Herbivores Resource Protection and Utilization, College of Animal Science and Technology, Southwest University, Chongqing, China
| | - Yong-Fu Huang
- 1Chongqing Key Laboratory of Forage & Herbivore, Chongqing Engineering Research Centre for Herbivores Resource Protection and Utilization, College of Animal Science and Technology, Southwest University, Chongqing, China
| | - Yong-Ju Zhao
- 1Chongqing Key Laboratory of Forage & Herbivore, Chongqing Engineering Research Centre for Herbivores Resource Protection and Utilization, College of Animal Science and Technology, Southwest University, Chongqing, China
| | - Ri-Su Na
- 1Chongqing Key Laboratory of Forage & Herbivore, Chongqing Engineering Research Centre for Herbivores Resource Protection and Utilization, College of Animal Science and Technology, Southwest University, Chongqing, China
| | - Zhong-Quan Zhao
- 1Chongqing Key Laboratory of Forage & Herbivore, Chongqing Engineering Research Centre for Herbivores Resource Protection and Utilization, College of Animal Science and Technology, Southwest University, Chongqing, China
| | - Yue-Hui Ma
- 2Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193 China
| | - Ming-Xing Chu
- 2Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193 China
| | - Wang-Dui Basang
- Institute of Animal Husbandry and Veterinary Medicine, Tibet Academy of Agriculture and Animal Husandry Science, Lasa, 850009 China
| | - Yan-Bin Zhu
- Institute of Animal Husbandry and Veterinary Medicine, Tibet Academy of Agriculture and Animal Husandry Science, Lasa, 850009 China
| | - Tian-Wu An
- 4Sichuan Academy of Grassland Sciences, Chengdu, 611731 Sichuan China
| | - Xiao-Lin Luo
- 4Sichuan Academy of Grassland Sciences, Chengdu, 611731 Sichuan China
| | - Yan-Guo Han
- 1Chongqing Key Laboratory of Forage & Herbivore, Chongqing Engineering Research Centre for Herbivores Resource Protection and Utilization, College of Animal Science and Technology, Southwest University, Chongqing, China
| | - Yan Zeng
- 1Chongqing Key Laboratory of Forage & Herbivore, Chongqing Engineering Research Centre for Herbivores Resource Protection and Utilization, College of Animal Science and Technology, Southwest University, Chongqing, China
| |
Collapse
|
33
|
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: 39] [Impact Index Per Article: 7.8] [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
|
34
|
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
|
35
|
Oget C, Servin B, Palhière I. Genetic diversity analysis of French goat populations reveals selective sweeps involved in their differentiation. Anim Genet 2018; 50:54-63. [PMID: 30549070 PMCID: PMC6590323 DOI: 10.1111/age.12752] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/19/2018] [Indexed: 12/31/2022]
Abstract
After domestication 11 000 years ago in Asia Minor, the goat followed human migration to Europe and Asia. It was then introduced in Africa and is now raised all over the world. In this study, we exploited a dataset composed of 54 000 SNPs (Illumina goat DNA chip) to analyze the genetic diversity of 223 individuals belonging to eight French breeds (Alpine, Angora, Corse, Fossés, Poitevine, Provençale, Pyrénées and Saanen). Analyses carried out included individual-based approaches (principal component analysis and population structure) and population-based approaches (phylogenetic tree constructions). The results of the genetic diversity analyses revealed that French breeds are clearly differentiated, in particular, the Angora breed that originates from south west Asia. The Provençale breed shows a very original genetic pattern that could be the result of ancient admixture. Then, selection signatures were detected by identifying regions of outlying genetic differentiation between populations. Five genomic regions were detected under selection on chromosomes 5, 6, 11, 13 and 20, revealing mainly soft selective sweeps and a few hard selective sweeps and highlighting candidate genes that had been selected for during the evolutionary history of these breeds. Among them, two coat coloration genes (ADAMTS20 and ASIP) and one gene related to milk composition (CSN1S1) were involved.
Collapse
Affiliation(s)
- C Oget
- GenPhySE, INRA, Université de Toulouse, INPT, ENVT, 31326, Castanet Tolosan, France
| | - B Servin
- GenPhySE, INRA, Université de Toulouse, INPT, ENVT, 31326, Castanet Tolosan, France
| | - I Palhière
- GenPhySE, INRA, Université de Toulouse, INPT, ENVT, 31326, Castanet Tolosan, France
| |
Collapse
|
36
|
Bertolini F, Servin B, Talenti A, Rochat E, Kim ES, Oget C, Palhière I, Crisà A, Catillo G, Steri R, Amills M, Colli L, Marras G, Milanesi M, Nicolazzi E, Rosen BD, Van Tassell CP, Guldbrandtsen B, Sonstegard TS, Tosser-Klopp G, Stella A, Rothschild MF, Joost S, Crepaldi P. Signatures of selection and environmental adaptation across the goat genome post-domestication. Genet Sel Evol 2018; 50:57. [PMID: 30449276 PMCID: PMC6240954 DOI: 10.1186/s12711-018-0421-y] [Citation(s) in RCA: 95] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2017] [Accepted: 10/15/2018] [Indexed: 01/03/2023] Open
Abstract
BACKGROUND Since goat was domesticated 10,000 years ago, many factors have contributed to the differentiation of goat breeds and these are classified mainly into two types: (i) adaptation to different breeding systems and/or purposes and (ii) adaptation to different environments. As a result, approximately 600 goat breeds have developed worldwide; they differ considerably from one another in terms of phenotypic characteristics and are adapted to a wide range of climatic conditions. In this work, we analyzed the AdaptMap goat dataset, which is composed of data from more than 3000 animals collected worldwide and genotyped with the CaprineSNP50 BeadChip. These animals were partitioned into groups based on geographical area, production uses, available records on solid coat color and environmental variables including the sampling geographical coordinates, to investigate the role of natural and/or artificial selection in shaping the genome of goat breeds. RESULTS Several signatures of selection on different chromosomal regions were detected across the different breeds, sub-geographical clusters, phenotypic and climatic groups. These regions contain genes that are involved in important biological processes, such as milk-, meat- or fiber-related production, coat color, glucose pathway, oxidative stress response, size, and circadian clock differences. Our results confirm previous findings in other species on adaptation to extreme environments and human purposes and provide new genes that could explain some of the differences between goat breeds according to their geographical distribution and adaptation to different environments. CONCLUSIONS These analyses of signatures of selection provide a comprehensive first picture of the global domestication process and adaptation of goat breeds and highlight possible genes that may have contributed to the differentiation of this species worldwide.
Collapse
Affiliation(s)
- Francesca Bertolini
- Department of Animal Science, Iowa State University, Ames, IA 50011 USA
- National Institute of Aquatic Resources, Technical University of Denmark (DTU), 2800 Lyngby, Denmark
| | - Bertrand Servin
- GenPhySE, INRA, Université de Toulouse, INPT, ENVT, 31326 Castanet Tolosan, France
| | - Andrea Talenti
- Dipartimento di Medicina Veterinaria, Università degli Studi di Milano, 20133 Milan, Italy
| | - Estelle Rochat
- Laboratory of Geographic Information Systems (LASIG), School of Architecture, Civil and Environmental Engineering (ENAC), Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | | | - Claire Oget
- GenPhySE, INRA, Université de Toulouse, INPT, ENVT, 31326 Castanet Tolosan, France
| | - Isabelle Palhière
- GenPhySE, INRA, Université de Toulouse, INPT, ENVT, 31326 Castanet Tolosan, France
| | - Alessandra Crisà
- Consiglio per la ricerca in agricoltura e l’analisi dell’economia agraria (CREA) - Research Centre for Animal Production and Acquaculture, 00015 Monterotondo, Roma, Italy
| | - Gennaro Catillo
- Consiglio per la ricerca in agricoltura e l’analisi dell’economia agraria (CREA) - Research Centre for Animal Production and Acquaculture, 00015 Monterotondo, Roma, Italy
| | - Roberto Steri
- Consiglio per la ricerca in agricoltura e l’analisi dell’economia agraria (CREA) - Research Centre for Animal Production and Acquaculture, 00015 Monterotondo, Roma, Italy
| | - Marcel Amills
- Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus Universitat Autonoma de Barcelona, Bellaterra, 08193 Barcelona, Spain
| | - Licia Colli
- DIANA Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti, Università Cattolica del S. Cuore, 29100 Piacenza, Italy
- BioDNA Centro di Ricerca sulla Biodiversità e sul DNA Antico, Università Cattolica del S. Cuore, 29100 Piacenza, Italy
| | - Gabriele Marras
- Fondazione Parco Tecnologico Padano (PTP), 26900 Lodi, Italy
| | - Marco Milanesi
- DIANA Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti, Università Cattolica del S. Cuore, 29100 Piacenza, Italy
- Department of Support, Production and Animal Health, School of Veterinary Medicine, São Paulo State University (UNESP), Araçatuba, Brazil
| | | | - Benjamin D. Rosen
- Animal Genomics and Improvement Laboratory, ARS USDA, Beltsville, MD 20705 USA
| | | | - Bernt Guldbrandtsen
- Center for Quantitative Genetics and Genomics, Aarhus University, 8830 Tjele, Denmark
| | | | - Gwenola Tosser-Klopp
- GenPhySE, INRA, Université de Toulouse, INPT, ENVT, 31326 Castanet Tolosan, France
| | - Alessandra Stella
- BioDNA Centro di Ricerca sulla Biodiversità e sul DNA Antico, Università Cattolica del S. Cuore, 29100 Piacenza, Italy
| | - Max F. Rothschild
- Department of Animal Science, Iowa State University, Ames, IA 50011 USA
| | - Stéphane Joost
- Laboratory of Geographic Information Systems (LASIG), School of Architecture, Civil and Environmental Engineering (ENAC), Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | - Paola Crepaldi
- Dipartimento di Medicina Veterinaria, Università degli Studi di Milano, 20133 Milan, Italy
| | - the AdaptMap consortium
- Department of Animal Science, Iowa State University, Ames, IA 50011 USA
- National Institute of Aquatic Resources, Technical University of Denmark (DTU), 2800 Lyngby, Denmark
- GenPhySE, INRA, Université de Toulouse, INPT, ENVT, 31326 Castanet Tolosan, France
- Dipartimento di Medicina Veterinaria, Università degli Studi di Milano, 20133 Milan, Italy
- Laboratory of Geographic Information Systems (LASIG), School of Architecture, Civil and Environmental Engineering (ENAC), Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
- Recombinetics Inc, St Paul, 55104 MN USA
- Consiglio per la ricerca in agricoltura e l’analisi dell’economia agraria (CREA) - Research Centre for Animal Production and Acquaculture, 00015 Monterotondo, Roma, Italy
- Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus Universitat Autonoma de Barcelona, Bellaterra, 08193 Barcelona, Spain
- DIANA Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti, Università Cattolica del S. Cuore, 29100 Piacenza, Italy
- BioDNA Centro di Ricerca sulla Biodiversità e sul DNA Antico, Università Cattolica del S. Cuore, 29100 Piacenza, Italy
- Fondazione Parco Tecnologico Padano (PTP), 26900 Lodi, Italy
- Department of Support, Production and Animal Health, School of Veterinary Medicine, São Paulo State University (UNESP), Araçatuba, Brazil
- Animal Genomics and Improvement Laboratory, ARS USDA, Beltsville, MD 20705 USA
- Center for Quantitative Genetics and Genomics, Aarhus University, 8830 Tjele, Denmark
| |
Collapse
|
37
|
Bertolini F, Cardoso TF, Marras G, Nicolazzi EL, Rothschild MF, Amills M. Genome-wide patterns of homozygosity provide clues about the population history and adaptation of goats. Genet Sel Evol 2018; 50:59. [PMID: 30449279 PMCID: PMC6241033 DOI: 10.1186/s12711-018-0424-8] [Citation(s) in RCA: 56] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2017] [Accepted: 10/15/2018] [Indexed: 01/15/2023] Open
Abstract
Background Patterns of homozygosity can be influenced by several factors, such as demography, recombination, and selection. Using the goat SNP50 BeadChip, we genotyped 3171 goats belonging to 117 populations with a worldwide distribution. Our objectives were to characterize the number and length of runs of homozygosity (ROH) and to detect ROH hotspots in order to gain new insights into the consequences of neutral and selection processes on the genome-wide homozygosity patterns of goats. Results The proportion of the goat genome covered by ROH is, in general, less than 15% with an inverse relationship between ROH length and frequency i.e. short ROH (< 3 Mb) are the most frequent ones. Our data also indicate that ~ 60% of the breeds display low FROH coefficients (< 0.10), while ~ 30 and ~ 10% of the goat populations show moderate (0.10 < FROH < 0.20) or high (> 0.20) FROH values. For populations from Asia, the average number of ROH is smaller and their coverage is lower in goats from the Near East than in goats from Central Asia, which is consistent with the role of the Fertile Crescent as the primary centre of goat domestication. We also observed that local breeds with small population sizes tend to have a larger fraction of the genome covered by ROH compared to breeds with tens or hundreds of thousands of individuals. Five regions on three goat chromosomes i.e. 11, 12 and 18, contain ROH hotspots that overlap with signatures of selection. Conclusions Patterns of homozygosity (average number of ROH of 77 and genome coverage of 248 Mb; FROH < 0.15) are similar in goats from different geographic areas. The increased homozygosity in local breeds is the consequence of their small population size and geographic isolation as well as of founder effects and recent inbreeding. The existence of three ROH hotspots that co-localize with signatures of selection demonstrates that selection has also played an important role in increasing the homozygosity of specific regions in the goat genome. Finally, most of the goat breeds analysed in this work display low levels of homozygosity, which is favourable for their genetic management and viability. Electronic supplementary material The online version of this article (10.1186/s12711-018-0424-8) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Francesca Bertolini
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA. .,National Institute of Aquatic Resources, Technical University of Denmark (DTU), Lyngby, 2800, Denmark.
| | - Tainã Figueiredo Cardoso
- Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus, Universitat Autonoma de Barcelona, 08193, Bellaterra, Barcelona, Spain
| | - Gabriele Marras
- Fondazione Parco Tecnologico Padano (PTP), 26900, Lodi, Italy
| | | | - Max F Rothschild
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA
| | - Marcel Amills
- Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus, Universitat Autonoma de Barcelona, 08193, Bellaterra, Barcelona, Spain
| | | |
Collapse
|
38
|
Kumar C, Song S, Dewani P, Kumar M, Parkash O, Ma Y, Malhi KK, Yang N, Mwacharo JM, He X, Jiang L. Population structure, genetic diversity and selection signatures within seven indigenous Pakistani goat populations. Anim Genet 2018; 49:592-604. [PMID: 30229969 DOI: 10.1111/age.12722] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/25/2018] [Indexed: 12/20/2022]
Abstract
Goat farming in Pakistan depends on indigenous breeds that have adapted to specific agro-ecological conditions. Pakistan has a rich resource of goat breeds, and the genetic diversity of these goat breeds is largely unknown. In this study, genetic diversity and population structure were characterized from seven indigenous goat breeds using the goat 50K SNP chip. The genetic diversity analysis showed that Bugi toori goats have the highest inbreeding level, consistent with the highest linkage disequilibrium, lowest diversity and long run of heterozygosity segments. This indicates that this breed should be prioritized in future conservation activities. The population structure analysis revealed four fairly distinct clusters (including Bugi toori, Bari, Black Tapri and some Kamori) and three other breeds that are seemingly the results of admixture between these or related groups (some Kamori, Pateri, Tapri and White Tapri). The selection signatures were evaluated in each breed. A total of 2508 putative selection signals were reported. The 26 significant windows were identified in more than four breeds, and selection signatures spanned several genes that directly or indirectly influence traits included coat colour variation (KIT), reproduction (BMPR1B, GNRHR, INSL6, JAK2 and EGR4), body size (SOCS2), ear size (MSRB3) and milk composition (ABCG2, SPP1, CSN1S2, CSN2, CSN3 and PROLACTIN).
Collapse
Affiliation(s)
- C Kumar
- Institute of Animal Science (IAS), Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China.,Directorate of Veterinary Research and Diagnostic Central Veterinary Diagnostic Laboratory, Tando Jam, 70050, Sindh, Pakistan.,Department of Animal Breeding and Genetics, Faculty of Animal Husbandry and Veterinary Sciences, Sindh Agricultural University, Tando Jam, 70060, Sindh, Pakistan
| | - S Song
- Institute of Animal Science (IAS), Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China.,Department of Animal Genetics and Breeding, China Agricultural University, Beijing, 100094, China
| | - P Dewani
- Directorate of Veterinary Research and Diagnostic Central Veterinary Diagnostic Laboratory, Tando Jam, 70050, Sindh, Pakistan
| | - M Kumar
- Department of Animal Breeding and Genetics, Faculty of Animal Husbandry and Veterinary Sciences, Sindh Agricultural University, Tando Jam, 70060, Sindh, Pakistan
| | - O Parkash
- Directorate of Veterinary Research and Diagnostic Central Veterinary Diagnostic Laboratory, Tando Jam, 70050, Sindh, Pakistan
| | - Y Ma
- Institute of Animal Science (IAS), Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China
| | - K K Malhi
- Department of Animal Breeding and Genetics, Faculty of Animal Husbandry and Veterinary Sciences, Sindh Agricultural University, Tando Jam, 70060, Sindh, Pakistan
| | - N Yang
- Department of Animal Genetics and Breeding, China Agricultural University, Beijing, 100094, China
| | - J M Mwacharo
- Small Ruminant Genomics Group, International Center for Agricultural Research in the Dry Areas (ICARDA), P.O. Box 5689, Addis Ababa, Ethiopia
| | - X He
- Institute of Animal Science (IAS), Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China
| | - L Jiang
- Institute of Animal Science (IAS), Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China
| |
Collapse
|
39
|
Onzima RB, Upadhyay MR, Doekes HP, Brito LF, Bosse M, Kanis E, Groenen MAM, Crooijmans RPMA. Genome-Wide Characterization of Selection Signatures and Runs of Homozygosity in Ugandan Goat Breeds. Front Genet 2018; 9:318. [PMID: 30154830 PMCID: PMC6102322 DOI: 10.3389/fgene.2018.00318] [Citation(s) in RCA: 70] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2018] [Accepted: 07/25/2018] [Indexed: 01/06/2023] Open
Abstract
Both natural and artificial selection are among the main driving forces shaping genetic variation across the genome of livestock species. Selection typically leaves signatures in the genome, which are often characterized by high genetic differentiation across breeds and/or a strong reduction in genetic diversity in regions associated with traits under intense selection pressure. In this study, we evaluated selection signatures and genomic inbreeding coefficients, FROH, based on runs of homozygosity (ROH), in six Ugandan goat breeds: Boer (n = 13), and the indigenous breeds Karamojong (n = 15), Kigezi (n = 29), Mubende (n = 29), Small East African (n = 29), and Sebei (n = 29). After genotyping quality control, 45,294 autosomal single nucleotide polymorphisms (SNPs) remained for further analyses. A total of 394 and 6 breed-specific putative selection signatures were identified across all breeds, based on marker-specific fixation index (FST-values) and haplotype differentiation (hapFLK), respectively. These regions were enriched with genes involved in signaling pathways associated directly or indirectly with environmental adaptation, such as immune response (e.g., IL10RB and IL23A), growth and fatty acid composition (e.g., FGF9 and IGF1), and thermo-tolerance (e.g., MTOR and MAPK3). The study revealed little overlap between breeds in genomic regions under selection and generally did not display the typical classic selection signatures as expected due to the complex nature of the traits. In the Boer breed, candidate genes associated with production traits, such as body size and growth (e.g., GJB2 and GJA3) were also identified. Furthermore, analysis of ROH in indigenous goat breeds showed very low levels of genomic inbreeding (with the mean FROH per breed ranging from 0.8% to 2.4%), as compared to higher inbreeding in Boer (mean FROH = 13.8%). Short ROH were more frequent than long ROH, except in Karamojong, providing insight in the developmental history of these goat breeds. This study provides insights into the effects of long-term selection in Boer and indigenous Ugandan goat breeds, which are relevant for implementation of breeding programs and conservation of genetic resources, as well as their sustainable use and management.
Collapse
Affiliation(s)
- Robert B. Onzima
- Animal Breeding and Genomics, Wageningen University and Research, Wageningen, Netherlands
- National Agricultural Research Organization (NARO), Entebbe, Uganda
| | - Maulik R. Upadhyay
- Animal Breeding and Genomics, Wageningen University and Research, Wageningen, Netherlands
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Harmen P. Doekes
- Animal Breeding and Genomics, Wageningen University and Research, Wageningen, Netherlands
| | - Luiz. F. Brito
- Department of Animal Biosciences, Centre for Genetic Improvement of Livestock (CGIL), University of Guelph, Guelph, ON, Canada
| | - Mirte Bosse
- Animal Breeding and Genomics, Wageningen University and Research, Wageningen, Netherlands
| | - Egbert Kanis
- Animal Breeding and Genomics, Wageningen University and Research, Wageningen, Netherlands
| | - Martien A. M. Groenen
- Animal Breeding and Genomics, Wageningen University and Research, Wageningen, Netherlands
| | | |
Collapse
|
40
|
Browett S, McHugo G, Richardson IW, Magee DA, Park SDE, Fahey AG, Kearney JF, Correia CN, Randhawa IAS, MacHugh DE. Genomic Characterisation of the Indigenous Irish Kerry Cattle Breed. Front Genet 2018. [PMID: 29520297 PMCID: PMC5827531 DOI: 10.3389/fgene.2018.00051] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Kerry cattle are an endangered landrace heritage breed of cultural importance to Ireland. In the present study we have used genome-wide SNP array data to evaluate genomic diversity within the Kerry population and between Kerry cattle and other European breeds. Patterns of genetic differentiation and gene flow among breeds using phylogenetic trees with ancestry graphs highlighted historical gene flow from the British Shorthorn breed into the ancestral population of modern Kerry cattle. Principal component analysis (PCA) and genetic clustering emphasised the genetic distinctiveness of Kerry cattle relative to comparator British and European cattle breeds. Modelling of genetic effective population size (Ne) revealed a demographic trend of diminishing Ne over time and that recent estimated Ne values for the Kerry breed may be less than the threshold for sustainable genetic conservation. In addition, analysis of genome-wide autozygosity (FROH) showed that genomic inbreeding has increased significantly during the 20 years between 1992 and 2012. Finally, signatures of selection revealed genomic regions subject to natural and artificial selection as Kerry cattle adapted to the climate, physical geography and agro-ecology of southwest Ireland.
Collapse
Affiliation(s)
- Sam Browett
- Ecosystems and Environment Research Centre, School of Environment and Life Sciences, University of Salford, Salford, United Kingdom
| | - Gillian McHugo
- Animal Genomics Laboratory, UCD School of Agriculture and Food Science, University College Dublin, Dublin, Ireland
| | | | - David A Magee
- Animal Genomics Laboratory, UCD School of Agriculture and Food Science, University College Dublin, Dublin, Ireland
| | | | - Alan G Fahey
- 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
| | - Imtiaz A S Randhawa
- Sydney School of Veterinary Science, The University of Sydney, Camden, NSW, Australia
| | - 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
|
41
|
Druml T, Neuditschko M, Grilz-Seger G, Horna M, Ricard A, Mesarič M, Cotman M, Pausch H, Brem G. Population Networks Associated with Runs of Homozygosity Reveal New Insights into the Breeding History of the Haflinger Horse. J Hered 2017; 109:384-392. [DOI: 10.1093/jhered/esx114] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2017] [Accepted: 12/23/2017] [Indexed: 02/04/2023] Open
Affiliation(s)
- Thomas Druml
- Institute of Animal Breeding and Genetics, University of Veterinary Sciences Vienna, Vienna, Austria
| | | | | | - Michaela Horna
- Department of Animal Husbandry, Slovak University of Agriculture in Nitra, Nitra-Chrenová, Slovak Republic
| | - Anne Ricard
- Institut National de la Recherche Agronomique, UMR 1313 Génétique Animale et Biologie Intégrative, Jouy-en-Josas, France
- Institut Français du Cheval et de l’Equitation, Recherche et Innovation, Exmes, France
| | - Matjaz Mesarič
- Clinic for Reproduction and Large Animals, Veterinary Faculty, University of Lubljana, Cesta v Mestni log, Ljubljana, Slovenia
| | - Marco Cotman
- Institute of Preclinical Sciences, Veterinary Faculty, University of Ljubljana, Cesta v Mestni log, Ljubljana, Slovenia
| | | | - Gottfried Brem
- Institute of Animal Breeding and Genetics, University of Veterinary Sciences Vienna, Vienna, Austria
| |
Collapse
|
42
|
Blant A, Kwong M, Szpiech ZA, Pemberton TJ. Weighted likelihood inference of genomic autozygosity patterns in dense genotype data. BMC Genomics 2017; 18:928. [PMID: 29191164 PMCID: PMC5709839 DOI: 10.1186/s12864-017-4312-3] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2017] [Accepted: 11/16/2017] [Indexed: 12/14/2022] Open
Abstract
Background Genomic regions of autozygosity (ROA) arise when an individual is homozygous for haplotypes inherited identical-by-descent from ancestors shared by both parents. Over the past decade, they have gained importance for understanding evolutionary history and the genetic basis of complex diseases and traits. However, methods to infer ROA in dense genotype data have not evolved in step with advances in genome technology that now enable us to rapidly create large high-resolution genotype datasets, limiting our ability to investigate their constituent ROA patterns. Methods We report a weighted likelihood approach for inferring ROA in dense genotype data that accounts for autocorrelation among genotyped positions and the possibilities of unobserved mutation and recombination events, and variability in the confidence of individual genotype calls in whole genome sequence (WGS) data. Results Forward-time genetic simulations under two demographic scenarios that reflect situations where inbreeding and its effect on fitness are of interest suggest this approach is better powered than existing state-of-the-art methods to infer ROA at marker densities consistent with WGS and popular microarray genotyping platforms used in human and non-human studies. Moreover, we present evidence that suggests this approach is able to distinguish ROA arising via consanguinity from ROA arising via endogamy. Using subsets of The 1000 Genomes Project Phase 3 data we show that, relative to WGS, intermediate and long ROA are captured robustly with popular microarray platforms, while detection of short ROA is more variable and improves with marker density. Worldwide ROA patterns inferred from WGS data are found to accord well with those previously reported on the basis of microarray genotype data. Finally, we highlight the potential of this approach to detect genomic regions enriched for autozygosity signals in one group relative to another based upon comparisons of per-individual autozygosity likelihoods instead of inferred ROA frequencies. Conclusions This weighted likelihood ROA inference approach can assist population- and disease-geneticists working with a wide variety of data types and species to explore ROA patterns and to identify genomic regions with differential ROA signals among groups, thereby advancing our understanding of evolutionary history and the role of recessive variation in phenotypic variation and disease. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-4312-3) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Alexandra Blant
- Department of Biochemistry and Medical Genetics, University of Manitoba, Winnipeg, MB, Canada
| | - Michelle Kwong
- Department of Biochemistry and Medical Genetics, University of Manitoba, Winnipeg, MB, Canada
| | - Zachary A Szpiech
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, CA, USA
| | - Trevor J Pemberton
- Department of Biochemistry and Medical Genetics, University of Manitoba, Winnipeg, MB, Canada.
| |
Collapse
|
43
|
Vidal O, Drögemüller C, Obexer-Ruff G, Reber I, Jordana J, Martínez A, Bâlteanu VA, Delgado JV, Eghbalsaied S, Landi V, Goyache F, Traoré A, Pazzola M, Vacca GM, Badaoui B, Pilla F, D'Andrea M, Álvarez I, Capote J, Sharaf A, Pons À, Amills M. Differential distribution of Y-chromosome haplotypes in Swiss and Southern European goat breeds. Sci Rep 2017; 7:16161. [PMID: 29170508 PMCID: PMC5701018 DOI: 10.1038/s41598-017-15593-1] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2017] [Accepted: 10/27/2017] [Indexed: 11/09/2022] Open
Abstract
The analysis of Y-chromosome variation has provided valuable clues about the paternal history of domestic animal populations. The main goal of the current work was to characterize Y-chromosome diversity in 31 goat populations from Central Eastern (Switzerland and Romania) and Southern Europe (Spain and Italy) as well as in reference populations from Africa and the Near East. Towards this end, we have genotyped seven single nucleotide polymorphisms (SNPs), mapping to the SRY, ZFY, AMELY and DDX3Y Y-linked loci, in 275 bucks from 31 populations. We have observed a low level of variability in the goat Y-chromosome, with just five haplotypes segregating in the whole set of populations. We have also found that Swiss bucks carry exclusively Y1 haplotypes (Y1A: 24%, Y1B1: 15%, Y1B2: 43% and Y1C: 18%), while in Italian and Spanish bucks Y2A is the most abundant haplotype (77%). Interestingly, in Carpathian goats from Romania the Y2A haplotype is also frequent (42%). The high Y-chromosome differentiation between Swiss and Italian/Spanish breeds might be due to the post-domestication spread of two different Near Eastern genetic stocks through the Danubian and Mediterranean corridors. Historical gene flow between Southern European and Northern African goats might have also contributed to generate such pattern of genetic differentiation.
Collapse
Affiliation(s)
- Oriol Vidal
- Departament de Biologia, Universitat de Girona, 17003, Girona, Spain.
| | - Cord Drögemüller
- Institute of Genetics, University of Bern, Bern, 3001, Switzerland
| | | | - Irene Reber
- Institute of Genetics, University of Bern, Bern, 3001, Switzerland
| | - Jordi Jordana
- Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Amparo Martínez
- Departamento de Genética, Universidad de Córdoba, 14071, Córdoba, Spain
| | - Valentin Adrian Bâlteanu
- Institute of Life Sciences, Faculty of Animal Science and Biotechnologies, University of Agricultural Sciences and Veterinary Medicine, 400372, Cluj-Napoca, Romania
| | | | - Shahin Eghbalsaied
- Transgenesis Center of Excellence, Isfahan (Khorasgan) Branch, Islamic Azad University, Isfahan, Iran
| | - Vincenzo Landi
- Departamento de Genética, Universidad de Córdoba, 14071, Córdoba, Spain
| | - Felix Goyache
- Área de Genética y Reproducción Animal, SERIDA-Deva, Camino de Rioseco 1225, Gijón, 33394, Spain
| | - Amadou Traoré
- Institut de l'Environnement et Recherches Agricoles, 04 BP 8645, Ouagadougou, 04, Burkina Faso
| | - Michele Pazzola
- Department of Veterinary Medicine, University of Sassari, 07100, Sassari, Italy
| | | | - Bouabid Badaoui
- University Mohammed V, Agdal, Faculty of Sciences, 4 Av. Ibn Battota, Rabat, Morocco
| | - Fabio Pilla
- Dipartimento Agricoltura, Ambiente e Alimenti, Università Degli Studi Del Molise, Campobasso, Italy
| | - Mariasilvia D'Andrea
- Dipartimento Agricoltura, Ambiente e Alimenti, Università Degli Studi Del Molise, Campobasso, Italy
| | - Isabel Álvarez
- Área de Genética y Reproducción Animal, SERIDA-Deva, Camino de Rioseco 1225, Gijón, 33394, Spain
| | - Juan Capote
- Instituto Canario de Investigaciones Agrarias, Canary Islands, Tenerife, La Laguna 38108, Spain
| | - Abdoallah Sharaf
- Genetic Department, Faculty of Agriculture, Ain Shams University, Cairo, 11241, Egypt.,Institute of Parasitology, Biology Centre, Czech Academy of Sciences, 37005, České Budějovice, Czechia.,Department of Animal Genetics, Center for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus Universitat Autònoma de Barcelona, Bellaterra, 08193, Spain
| | - Àgueda Pons
- Unitat de Races Autòctones, Servei de Millora Agrària, (SEMILLA-SAU), Son Ferriol, 07198, Spain
| | - Marcel Amills
- Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain.,Department of Animal Genetics, Center for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus Universitat Autònoma de Barcelona, Bellaterra, 08193, Spain
| |
Collapse
|
44
|
Signer-Hasler H, Burren A, Neuditschko M, Frischknecht M, Garrick D, Stricker C, Gredler B, Bapst B, Flury C. Population structure and genomic inbreeding in nine Swiss dairy cattle populations. Genet Sel Evol 2017; 49:83. [PMID: 29115934 PMCID: PMC5674839 DOI: 10.1186/s12711-017-0358-6] [Citation(s) in RCA: 44] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2016] [Accepted: 10/26/2017] [Indexed: 12/25/2022] Open
Abstract
BACKGROUND Domestication, breed formation and intensive selection have resulted in divergent cattle breeds that likely exhibit their own genomic signatures. In this study, we used genotypes from 27,612 autosomal single nucleotide polymorphisms to characterize population structure based on 9214 sires representing nine Swiss dairy cattle populations: Brown Swiss (BS), Braunvieh (BV), Original Braunvieh (OB), Holstein (HO), Red Holstein (RH), Swiss Fleckvieh (SF), Simmental (SI), Eringer (ER) and Evolèner (EV). Genomic inbreeding (F ROH) and signatures of selection were determined by calculating runs of homozygosity (ROH). The results build the basis for a better understanding of the genetic development of Swiss dairy cattle populations and highlight differences between the original populations (i.e. OB, SI, ER and EV) and those that have become more popular in Switzerland as currently reflected by their larger populations (i.e. BS, BV, HO, RH and SF). RESULTS The levels of genetic diversity were highest and lowest in the SF and BS breeds, respectively. Based on F ST values, we conclude that, among all pairwise comparisons, BS and HO (0.156) differ more than the other pairs of populations. The original Swiss cattle populations OB, SI, ER, and EV are clearly genetically separated from the Swiss cattle populations that are now more common and represented by larger numbers of cows. Mean levels of F ROH ranged from 0.027 (ER) to 0.091 (BS). Three of the original Swiss cattle populations, ER (F ROH: 0.027), OB (F ROH: 0.029), and SI (F ROH: 0.039), showed low levels of genomic inbreeding, whereas it was much higher in EV (F ROH: 0.074). Private signatures of selection for the original Swiss cattle populations are reported for BTA4, 5, 11 and 26. CONCLUSIONS The low levels of genomic inbreeding observed in the original Swiss cattle populations ER, OB and SI compared to the other breeds are explained by a lesser use of artificial insemination and greater use of natural service. Natural service results in more sires having progeny at each generation and thus this breeding practice is likely the major reason for the remarkable levels of genetic diversity retained within these populations. The fact that the EV population is regionally restricted and its small census size of herd-book cows explain its high level of genomic inbreeding.
Collapse
Affiliation(s)
- Heidi Signer-Hasler
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
| | - Alexander Burren
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
| | | | - Mirjam Frischknecht
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
- Qualitas AG, Zug, Switzerland
| | | | | | | | | | - Christine Flury
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
| |
Collapse
|
45
|
The Effect of Botulinum Toxin Type A on Expression Profiling of Long Noncoding RNAs in Human Dermal Fibroblasts. BIOMED RESEARCH INTERNATIONAL 2017; 2017:2957941. [PMID: 28265570 PMCID: PMC5318640 DOI: 10.1155/2017/2957941] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/04/2016] [Revised: 12/04/2016] [Accepted: 12/20/2016] [Indexed: 01/13/2023]
Abstract
Objective. This study was aimed at analyzing the expressions of long noncoding RNAs (lncRNAs) in Botulinum Toxin Type A (BoNTA) treated human dermal fibroblasts (HDFs) in vitro. Methods. We used RNA sequencing to characterize the lncRNAs and mRNAs transcriptome in the control and BoNTA treated group, in conjunction with application of GO (gene ontology) analysis and KEGG (kyoto encyclopedia of genes and genomes) analysis to delineate the alterations in gene expression. We also obtained quantitative real time polymerase chain reaction (qRT-PCR) to confirm some differentially expressed genes. Results. Numerous differentially expressed genes were observed by microarrays between the two groups. qRT-PCR confirmed the changes of six lncRNAs (RP11-517C16.2-001, FR271872, LOC283352, RP11-401E9.3, FGFR3P, and XXbac-BPG16N22.5) and nine mRNAs (NOS2, C13orf15, FOS, FCN2, SPINT1, PLAC8, BIRC5, NOS2, and COL19A1). Farther studies indicated that the downregulating effect of BoNTA on the expression of FGFR3P was time-related and the dosage of BoNTA at a range from 2.5 U/106 cells to 7.5 U/106 cells increased the expression of FGFR3P and COL19A1 in HDFs as well. Conclusion. The expression profiling of lncRNAs was visibly changed in BoNTA treated HDFs. Further studies should focus on several lncRNAs to investigate their functions in BoNTA treated HDFs and the underlying mechanisms.
Collapse
|
46
|
Ginja C, Gama LT, Martínez A, Sevane N, Martin-Burriel I, Lanari MR, Revidatti MA, Aranguren-Méndez JA, Bedotti DO, Ribeiro MN, Sponenberg P, Aguirre EL, Alvarez-Franco LA, Menezes MPC, Chacón E, Galarza A, Gómez-Urviola N, Martínez-López OR, Pimenta-Filho EC, da Rocha LL, Stemmer A, Landi V, Delgado-Bermejo JV. Genetic diversity and patterns of population structure in Creole goats from the Americas. Anim Genet 2017; 48:315-329. [PMID: 28094449 DOI: 10.1111/age.12529] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/02/2016] [Indexed: 01/03/2023]
Abstract
Biodiversity studies are more efficient when large numbers of breeds belonging to several countries are involved, as they allow for an in-depth analysis of the within- and between-breed components of genetic diversity. A set of 21 microsatellites was used to investigate the genetic composition of 24 Creole goat breeds (910 animals) from 10 countries to estimate levels of genetic variability, infer population structure and understand genetic relationships among populations across the American continent. Three commercial transboundary breeds were included in the analyses to investigate admixture with Creole goats. Overall, the genetic diversity of Creole populations (mean number of alleles = 5.82 ± 1.14, observed heterozygosity = 0.585 ± 0.074) was moderate and slightly lower than what was detected in other studies with breeds from other regions. The Bayesian clustering analysis without prior information on source populations identified 22 breed clusters. Three groups comprised more than one population, namely from Brazil (Azul and Graúna; Moxotó and Repartida) and Argentina (Long and shorthair Chilluda, Pampeana Colorada and Angora-type goat). Substructure was found in Criolla Paraguaya. When prior information on sample origin was considered, 92% of the individuals were assigned to the source population (threshold q ≥ 0.700). Creole breeds are well-differentiated entities (mean coefficient of genetic differentiation = 0.111 ± 0.048, with the exception of isolated island populations). Dilution from admixture with commercial transboundary breeds appears to be negligible. Significant levels of inbreeding were detected (inbreeding coefficient > 0 in most Creole goat populations, P < 0.05). Our results provide a broad perspective on the extant genetic diversity of Creole goats, however further studies are needed to understand whether the observed geographical patterns of population structure may reflect the mode of goat colonization in the Americas.
Collapse
Affiliation(s)
- C Ginja
- CIBIO-InBIO - Centro de Investigação em Biodiversidade e Recursos Genéticos, Universidade do Porto, Campus Agrário de Vairão, Rua Padre Armando Quintas n. 7, 4485-661, Vairão, Portugal
| | - L T Gama
- CIISA, Faculdade de Medicina Veterinária, Universidade de Lisboa, Lisboa, Portugal
| | - A Martínez
- Departamento de Genética, Universidad de Córdoba, Campus de Excelencia Internacional Agroalimentario ceiA3, Córdoba, Spain
| | - N Sevane
- Departamento de Producción Animal, Universidad Complutense de Madrid, Madrid, Spain
| | - I Martin-Burriel
- Martin-Burriel, Laboratorio de Genética Bioquímica, Facultad de Veterinaria, Universidad de Zaragoza, Zaragoza, Spain
| | - M R Lanari
- Area de Producción Animal, Instituto Nacional de Tecnología Agropecuaria EEA, Bariloche, Argentina
| | - M A Revidatti
- Facultad de Ciencias Veterinarias, Universidad Nacional del Nordeste, Corrientes, Argentina
| | - J A Aranguren-Méndez
- Facultad de Ciencias Veterinarias, Universidad de Zulia, Maracaibo-Zulia, Venezuela
| | - D O Bedotti
- Instituto Nacional de Tecnología Agropecuaria EEA Anguil "Ing. Agr. Guillermo Covas", Bariloche, Argentina
| | - M N Ribeiro
- Departamento de Zootecnia, Universidade Federal Rural de Pernambuco, Recife, PE, Brazil
| | - P Sponenberg
- Virginia-Maryland Regional College of Veterinary Medicine, Virginia Tech, Blacksburg, VA, USA
| | - E L Aguirre
- Universidad Nacional de Loja, Loja, Ecuador.,Grupo de Melhoramento Animal e Biotecnologias GMAB-FZEA-USP, Brazil
| | | | | | - E Chacón
- Universidad Técnica de Cotopaxi, La Maná, Ecuador
| | - A Galarza
- Universidad Mayor de San Simón, Cochabamba, Bolivia
| | - N Gómez-Urviola
- Universidad Nacional Micaela Bastidas de Apurímac, Abancay, Perú
| | - O R Martínez-López
- Centro Multidisciplinario de Investigaciones Tecnológicas, Dirección General de Investigación Científica y Tecnológica, Universidad Nacional de Asunción, San Lorenzo, Paraguay
| | | | - L L da Rocha
- Departamento de Zootecnia, Universidade Federal Rural de Pernambuco, Recife, PE, Brazil
| | - A Stemmer
- Universidad Mayor de San Simón, Cochabamba, Bolivia
| | - V Landi
- Departamento de Genética, Universidad de Córdoba, Campus de Excelencia Internacional Agroalimentario ceiA3, Córdoba, Spain
| | - J V Delgado-Bermejo
- Departamento de Genética, Universidad de Córdoba, Campus de Excelencia Internacional Agroalimentario ceiA3, Córdoba, Spain
| |
Collapse
|