1
|
Bayraktar M. Analysing the genetic diversity of three sheep breeds in Turkey and nearby countries using 50 K SNPs data. Anim Biotechnol 2024; 35:2329106. [PMID: 38497403 DOI: 10.1080/10495398.2024.2329106] [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] [Indexed: 03/19/2024]
Abstract
This study analysed the genetic diversity and population structure of eight sheep breeds in Turkey and nearby countries. Moderate genetic diversity was observed, with the Sakiz (SKZ) exhibiting the highest diversity based on heterozygosity and allelic richness (AR) values. Genetic distances revealed differentiation between the populations, with the most significant divergence between the Cyprus Fat Tail (CFT) and SKZ breeds. PCA demonstrated SKZ and Chios (CHI) clustering together, indicating genetic similarity. Karakas (KRS), Norduz (NDZ), Afshari (AFS), Moghani (MOG) and others showed overlap, reflecting genetic relationships. Ancestry analysis found that KRS was predominantly inherited from the second ancestral population, while SKZ and NDZ were primarily derived from the first and second ancestral lineages. This illustrated the populations' diverse origins. Most genetic variation (96.84%) was within, not between, populations. The phi-statistic (PhiPT) indicated moderate differentiation overall. Phylogenetic analysis further demonstrated the genetic distinctiveness of the SKZ breed. ROH and FROH analyses showed that SKZ exhibited the highest homozygosity and inbreeding, while KRS displayed the lowest. This study elucidates these breeds' genetic diversity, structure and relationships. Key findings include moderate diversity, evidence of differentiation between breeds, diverse ancestral origins and distinct ROH patterns. This provides insights into the population's genetic characteristics and conservation requirements.
Collapse
Affiliation(s)
- Mervan Bayraktar
- Department of Animal Science, Faculty of Agriculture, Çukurova University, Adana, Turkey
| |
Collapse
|
2
|
Klawatsch J, Papachristou D, Koutsouli P, Upadhyay M, Seichter D, Russ I, Mioč B, Simčič M, Bizelis I, Medugorac I. Genetic basis of ear length in sheep breeds sampled across the region from the Middle East to the Alps. Anim Genet 2024; 55:123-133. [PMID: 38069488 DOI: 10.1111/age.13387] [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/28/2023] [Revised: 11/28/2023] [Accepted: 11/29/2023] [Indexed: 01/04/2024]
Abstract
Ear length in sheep (Ovis aries) shows a wide range of natural variation, from the absence of an outer ear structure (anotia), to small outer ears (microtia), to regular ear length. Up until now, the underlying genetics of this phenotype has been studied in four sheep breeds from China, Jordan and Italy. These studies revealed a broad range of genes significantly associated with ear length, potentially indicating genetic heterogeneity across breeds or geographic regions. In the current study, we performed genome-wide SNP genotyping and haplotype-based mapping, in a population of 340 individuals, to identify loci influencing ear length variation in additional sheep breeds from Slovenia, Croatia, Cyprus and Greece. Additionally, two previously described candidate variants were also genotyped in our mapping population. The mapping model without candidate variant genotypes revealed only one genome-wide significant signal, which was located next to HMX1 on OAR6. This region was previously described as being associated with ear length variation in the Altay and Awassi sheep breeds. The mapping model including the candidate duplication genotype near HMX1 as a fixed effect explained the phenotypic variance on OAR6 and revealed an additional genome-wide significant locus on OAR13 associated with ear length. Our results, combined with published evidence, suggest that a duplication in the evolutionarily conserved region near HMX1 is the major regulator of ear length in sheep breeds descended from a larger region from Central Asia, to the Middle East, Cyprus, Greece and to the Alps. This distribution suggests an ancient origin of the derived allele.
Collapse
Affiliation(s)
- Jürgen Klawatsch
- Population Genomics Group, Department of Veterinary Sciences, Ludwig Maximillians University Munich, Munich, Germany
- Tierzuchtforschung e.V. Munich, Grub, Germany
| | | | - Panagiota Koutsouli
- Department of Animal Science, Agricultural University of Athens, Athens, Greece
| | - Maulik Upadhyay
- Population Genomics Group, Department of Veterinary Sciences, Ludwig Maximillians University Munich, Munich, Germany
| | | | - Ingolf Russ
- Tierzuchtforschung e.V. Munich, Grub, Germany
| | - Boro Mioč
- Department of Animal Science and Technology, Faculty of Agriculture, University of Zagreb, Zagreb, Croatia
| | - Mojca Simčič
- Department of Animal Science, Biotechnical Faculty, University of Ljubljana, Ljubljana, Slovenia
| | - Iosif Bizelis
- Department of Animal Science, Agricultural University of Athens, Athens, Greece
| | - Ivica Medugorac
- Population Genomics Group, Department of Veterinary Sciences, Ludwig Maximillians University Munich, Munich, Germany
| |
Collapse
|
3
|
Jung H, Choi CH, Lee M, Kim SY, Aknazarov B, Nyrgaziev R, Atabekova N, Jetigenov E, Chung YS, Lee HI. Molecular Detection and Phylogenetic Analysis of Tick-Borne Encephalitis Virus from Ticks Collected from Cattle in Kyrgyzstan, 2023. Viruses 2024; 16:107. [PMID: 38257807 PMCID: PMC10821214 DOI: 10.3390/v16010107] [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: 12/15/2023] [Revised: 01/05/2024] [Accepted: 01/10/2024] [Indexed: 01/24/2024] Open
Abstract
Ticks are important vectors of the tick-borne encephalitis virus (TBEV). In Kyrgyzstan, the livestock farming trade and nomadic lifestyle enable tick-borne diseases to be imported from neighboring countries, but there are few relevant studies. In this study, we collected 40 ticks from cattle in Kyrgyzstan. Molecular marker analysis identified the ticks as Ixodes persulcatus (97.5%; n = 39) and Haemaphysalis punctata (2.5%; n = 1). Real-time PCR screening revealed two ticks to be positive for TBEV, but only one tick was amplified using nested PCR targeting the TBEV envelope (E) and non-structure 5 (NS5) gene. The obtained sequences belonged to the TBEV Siberian subtype and phylogenetic tree analysis results confirmed that the virus was related to the Bosnia strain. We also performed next-generation sequencing, which confirmed the TBEV Siberian subtype. Continuous research and surveillance of TBEV in Kyrgyzstan are required to provide further information on tick-borne diseases.
Collapse
Affiliation(s)
- Haneul Jung
- Division of Vectors and Parasitic Diseases, Korea Disease Control and Prevention Agency (KDCA), Cheongju 28159, Republic of Korea; (H.J.); (S.-Y.K.)
| | - Chi-Hwan Choi
- Division of High-Risk Pathogens, Korea Disease Control and Prevention Agency (KDCA), Cheongju 28159, Republic of Korea; (C.-H.C.); (M.L.); (Y.-S.C.)
| | - Minji Lee
- Division of High-Risk Pathogens, Korea Disease Control and Prevention Agency (KDCA), Cheongju 28159, Republic of Korea; (C.-H.C.); (M.L.); (Y.-S.C.)
| | - Seong-Yoon Kim
- Division of Vectors and Parasitic Diseases, Korea Disease Control and Prevention Agency (KDCA), Cheongju 28159, Republic of Korea; (H.J.); (S.-Y.K.)
| | - Bekbolsun Aknazarov
- Faculty of Veterinary Medicine, Kyrgyz National Agrarian University Named after K. I. Skryabin, Bishkek 720005, Kyrgyzstan; (B.A.); (R.N.); (N.A.); (E.J.)
| | - Rysbek Nyrgaziev
- Faculty of Veterinary Medicine, Kyrgyz National Agrarian University Named after K. I. Skryabin, Bishkek 720005, Kyrgyzstan; (B.A.); (R.N.); (N.A.); (E.J.)
| | - Nurzina Atabekova
- Faculty of Veterinary Medicine, Kyrgyz National Agrarian University Named after K. I. Skryabin, Bishkek 720005, Kyrgyzstan; (B.A.); (R.N.); (N.A.); (E.J.)
| | - Elmurat Jetigenov
- Faculty of Veterinary Medicine, Kyrgyz National Agrarian University Named after K. I. Skryabin, Bishkek 720005, Kyrgyzstan; (B.A.); (R.N.); (N.A.); (E.J.)
| | - Yoon-Seok Chung
- Division of High-Risk Pathogens, Korea Disease Control and Prevention Agency (KDCA), Cheongju 28159, Republic of Korea; (C.-H.C.); (M.L.); (Y.-S.C.)
| | - Hee-Il Lee
- Division of Vectors and Parasitic Diseases, Korea Disease Control and Prevention Agency (KDCA), Cheongju 28159, Republic of Korea; (H.J.); (S.-Y.K.)
| |
Collapse
|
4
|
Zhang CL, Zhang J, Tuersuntuoheti M, Zhou W, Han Z, Li X, Yang R, Zhang L, Zheng L, Liu S. Landscape genomics reveals adaptive divergence of indigenous sheep in different ecological environments of Xinjiang, China. THE SCIENCE OF THE TOTAL ENVIRONMENT 2023; 904:166698. [PMID: 37683864 DOI: 10.1016/j.scitotenv.2023.166698] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/22/2023] [Revised: 08/22/2023] [Accepted: 08/28/2023] [Indexed: 09/10/2023]
Abstract
Sheep are important livestock animals that have evolved under various ecological pressures. Xinjiang is a region with diverse and harsh environments that have shaped many local sheep breeds with unique characteristics and environmental adaptability. However, these breeds are losing ecological flexibility due to the promotion of intensive farming practices. Here we sequenced 14 local sheep breeds from Xinjiang and analyzed their genetic structure and gene flow with other sheep breeds from neighboring regions. The Tibetan Plateau was the geographic origin of Xinjiang native sheep evolution. We performed genome-environment association analysis and identified Bio9: Mean Temperature of Driest Quarter and Bio15: Precipitation Seasonality as the key environmental factors affecting Xinjiang local sheep and the key genes involved in their survival and adaptation. We classified Xinjiang native sheep breeds into six groups based on their differential genes by pairwise selective sweep analysis and Community Network Analysis. We analyzed transcriptome expression data of 832 sheep tissues and detected tissue-specific enrichment of six group-specific genes in different biological systems. Our results revealed the genetic basis of year-round estrus, drought tolerance, hypoxia resistance, and cold tolerance traits of Xinjiang sheep breeds. Moreover, we proposed conservation strategies for Xinjiang local sheep breeds and provided theoretical guidance for breeding new sheep breeds under global extreme environments.
Collapse
Affiliation(s)
- Cheng-Long Zhang
- College of Animal Science and Technology, Tarim University, Xingfu Road, Alar 843300, Xinjiang, China; Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Xingfu Road, Alar 843300, Xinjiang, China
| | - Jihu Zhang
- College of Animal Science and Technology, Tarim University, Xingfu Road, Alar 843300, Xinjiang, China; Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Xingfu Road, Alar 843300, Xinjiang, China
| | - Mirenisa Tuersuntuoheti
- College of Animal Science and Technology, Tarim University, Xingfu Road, Alar 843300, Xinjiang, China; Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Xingfu Road, Alar 843300, Xinjiang, China
| | - Wen Zhou
- College of Animal Science and Technology, Tarim University, Xingfu Road, Alar 843300, Xinjiang, China; Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Xingfu Road, Alar 843300, Xinjiang, China
| | - Zhipeng Han
- College of Animal Science and Technology, Tarim University, Xingfu Road, Alar 843300, Xinjiang, China; Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Xingfu Road, Alar 843300, Xinjiang, China
| | - Xiaopeng Li
- College of Animal Science and Technology, Tarim University, Xingfu Road, Alar 843300, Xinjiang, China; Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Xingfu Road, Alar 843300, Xinjiang, China
| | - Ruizhi Yang
- College of Animal Science and Technology, Tarim University, Xingfu Road, Alar 843300, Xinjiang, China; Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Xingfu Road, Alar 843300, Xinjiang, China
| | - Lulu Zhang
- College of Animal Science and Technology, Tarim University, Xingfu Road, Alar 843300, Xinjiang, China; Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Xingfu Road, Alar 843300, Xinjiang, China
| | - Langman Zheng
- College of Animal Science and Technology, Tarim University, Xingfu Road, Alar 843300, Xinjiang, China; Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Xingfu Road, Alar 843300, Xinjiang, China
| | - Shudong Liu
- College of Animal Science and Technology, Tarim University, Xingfu Road, Alar 843300, Xinjiang, China; Key Laboratory of Tarim Animal Husbandry Science and Technology, Xinjiang Production and Construction Corps, Xingfu Road, Alar 843300, Xinjiang, China.
| |
Collapse
|
5
|
Amandykova M, Akhatayeva Z, Kozhakhmet A, Kapassuly T, Orazymbetova Z, Yergali K, Khamzin K, Iskakov K, Dossybayev K. Distribution of Runs of Homozygosity and Their Relationship with Candidate Genes for Productivity in Kazakh Meat-Wool Sheep Breed. Genes (Basel) 2023; 14:1988. [PMID: 38002931 PMCID: PMC10671688 DOI: 10.3390/genes14111988] [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: 09/09/2023] [Revised: 10/17/2023] [Accepted: 10/23/2023] [Indexed: 11/26/2023] Open
Abstract
Increasing the fertility of sheep remains one of the crucial issues of modern sheep breeding. The Kazakh meat-wool sheep is an excellent breed with high meat and wool productivity and well adapted to harsh conditions. Nowadays, runs of homozygosity (ROHs) are considered a suitable approach for studying the genetic characteristics of farm animals. The aims of the study were to analyze the distribution of ROHs, describe autozygosity, and detect genomic regions with high ROH islands. In this study, we genotyped a total of 281 Kazakh meat-wool sheep using the Illumina iScan® system (EquipNet, Canton, MA, USA) via Ovine SNP50 BeadChip array. As a results, a total of 15,069 ROHs were found in the three Kazakh meat-wool sheep populations. The mean number of ROH per animal across populations varied from 40.3 (POP1) to 42.2 (POP2) in the category 1+ Mb. Furthermore, the number of ROH per animal in ROH1-2 Mb were much higher than ROH2-4 Mb and ROH8-16 Mb in the three sheep populations. Most of individuals had small number of ROH>16 Mb. The highest and lowest genomic inbreeding coefficient values were observed in POP2 and POP3, respectively. The estimated FROH presented the impact that recent inbreeding has had in all sheep populations. Furthermore, a set of interesting candidate genes (BMP2, BMPR2, BMPRIB, CLOCK, KDM2B, TIAM1, TASP1, MYBPC1, MYOM1, and CACNA2D1), which are related to the productive traits, were found. Collectively, these findings will contribute to the breeding and conservation strategies of the Kazakh meat-wool sheep breed.
Collapse
Affiliation(s)
- Makpal Amandykova
- Laboratory of Animal Genetics and Cytogenetics, Institute of Genetics and Physiology SC MSHE RK, Al-Farabi Ave. 93, Almaty 050060, Kazakhstan; (M.A.); (Z.A.); (A.K.); (T.K.); (Z.O.); (K.Y.); (K.I.)
- Department of Molecular Biology and Genetics, Faculty of Biology and Biotechnology, Al-Farabi Kazakh National University, Al-Farabi Ave. 71, Almaty 050042, Kazakhstan
| | - Zhanerke Akhatayeva
- Laboratory of Animal Genetics and Cytogenetics, Institute of Genetics and Physiology SC MSHE RK, Al-Farabi Ave. 93, Almaty 050060, Kazakhstan; (M.A.); (Z.A.); (A.K.); (T.K.); (Z.O.); (K.Y.); (K.I.)
- Laboratory of Molecular Genetics, Kazakh Research Institute of Livestock and Fodder Production, Zhandosov Str. 51, Almaty 050035, Kazakhstan;
| | - Altynay Kozhakhmet
- Laboratory of Animal Genetics and Cytogenetics, Institute of Genetics and Physiology SC MSHE RK, Al-Farabi Ave. 93, Almaty 050060, Kazakhstan; (M.A.); (Z.A.); (A.K.); (T.K.); (Z.O.); (K.Y.); (K.I.)
- Department of Molecular Biology and Genetics, Faculty of Biology and Biotechnology, Al-Farabi Kazakh National University, Al-Farabi Ave. 71, Almaty 050042, Kazakhstan
- Laboratory of Molecular Genetics, Kazakh Research Institute of Livestock and Fodder Production, Zhandosov Str. 51, Almaty 050035, Kazakhstan;
| | - Tilek Kapassuly
- Laboratory of Animal Genetics and Cytogenetics, Institute of Genetics and Physiology SC MSHE RK, Al-Farabi Ave. 93, Almaty 050060, Kazakhstan; (M.A.); (Z.A.); (A.K.); (T.K.); (Z.O.); (K.Y.); (K.I.)
- Department of Molecular Biology and Genetics, Faculty of Biology and Biotechnology, Al-Farabi Kazakh National University, Al-Farabi Ave. 71, Almaty 050042, Kazakhstan
| | - Zarina Orazymbetova
- Laboratory of Animal Genetics and Cytogenetics, Institute of Genetics and Physiology SC MSHE RK, Al-Farabi Ave. 93, Almaty 050060, Kazakhstan; (M.A.); (Z.A.); (A.K.); (T.K.); (Z.O.); (K.Y.); (K.I.)
| | - Kanagat Yergali
- Laboratory of Animal Genetics and Cytogenetics, Institute of Genetics and Physiology SC MSHE RK, Al-Farabi Ave. 93, Almaty 050060, Kazakhstan; (M.A.); (Z.A.); (A.K.); (T.K.); (Z.O.); (K.Y.); (K.I.)
- Laboratory of Molecular Genetics, Kazakh Research Institute of Livestock and Fodder Production, Zhandosov Str. 51, Almaty 050035, Kazakhstan;
| | - Kadyrzhan Khamzin
- Laboratory of Molecular Genetics, Kazakh Research Institute of Livestock and Fodder Production, Zhandosov Str. 51, Almaty 050035, Kazakhstan;
| | - Kairat Iskakov
- Laboratory of Animal Genetics and Cytogenetics, Institute of Genetics and Physiology SC MSHE RK, Al-Farabi Ave. 93, Almaty 050060, Kazakhstan; (M.A.); (Z.A.); (A.K.); (T.K.); (Z.O.); (K.Y.); (K.I.)
- Laboratory of Molecular Genetics, Kazakh Research Institute of Livestock and Fodder Production, Zhandosov Str. 51, Almaty 050035, Kazakhstan;
| | - Kairat Dossybayev
- Laboratory of Animal Genetics and Cytogenetics, Institute of Genetics and Physiology SC MSHE RK, Al-Farabi Ave. 93, Almaty 050060, Kazakhstan; (M.A.); (Z.A.); (A.K.); (T.K.); (Z.O.); (K.Y.); (K.I.)
- Department of Molecular Biology and Genetics, Faculty of Biology and Biotechnology, Al-Farabi Kazakh National University, Al-Farabi Ave. 71, Almaty 050042, Kazakhstan
- Laboratory of Molecular Genetics, Kazakh Research Institute of Livestock and Fodder Production, Zhandosov Str. 51, Almaty 050035, Kazakhstan;
| |
Collapse
|
6
|
Barani S, Nejati-Javaremi A, Moradi MH, Moradi-Sharbabak M, Gholizadeh M, Esfandyari H. Genome-wide study of linkage disequilibrium, population structure, and inbreeding in Iranian indigenous sheep breeds. PLoS One 2023; 18:e0286463. [PMID: 37267244 DOI: 10.1371/journal.pone.0286463] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2022] [Accepted: 05/16/2023] [Indexed: 06/04/2023] Open
Abstract
Knowledge of linkage disequilibrium (LD), genetic structure and genetic diversity are some key parameters to study the breeding history of indigenous small ruminants. In this study, the OvineSNP50 Bead Chip array was used to estimate and compare LD, genetic diversity, effective population size (Ne) and genomic inbreeding in 186 individuals, from three Iranian indigenous sheep breeds consisting of Baluchi (n = 96), Lori-Bakhtiari (n = 47) and Zel (n = 47). The results of principal component analysis (PCA) revealed that all animals were allocated to the groups that they sampled and the admixture analysis revealed that the structure within the populations is best explained when separated into three groups (K = 3). The average r2 values estimated between adjacent single nucleotide polymorphisms (SNPs) at distances up to 10Kb, were 0.388±0.324, 0.353±0.311, and 0.333±0.309 for Baluchi, Lori-Bakhtiari and Zel, respectively. Estimation of genetic diversity and effective population size (Ne) showed that the Zel breed had the highest heterozygosity and Ne, whereas the lowest value was found in Baluchi breed. Estimation of genomic inbreeding using FROH (based on the long stretches of consecutive homozygous genotypes) showed the highest inbreeding coefficient in Baluchi and the lowest in Zel breed that could be due to higher pressure of artificial selection on Baluchi breed. The results of genomic inbreeding and Ne showed an increase in sharing haplotypes in Baluchi, leading to the enlargement of LD and the consequences of linkage disequilibrium and haplotype blocks confirmed this point. Also, the persistence of the LD phase between Zel and Lori-Bakhtiari was highest indicating that these two breeds would be combined in a multi-breed training population in genomic selection studies.
Collapse
Affiliation(s)
- S Barani
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, Iran
| | - A Nejati-Javaremi
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, Iran
| | - M H Moradi
- Department of Animal Science, Faculty of Agriculture and Natural Resources, Arak University, Arak, Iran
| | - M Moradi-Sharbabak
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, Iran
| | - M Gholizadeh
- Department of Animal Science, Sari Agricultural Sciences and Natural Resources University, Sari, Mazandaran, Iran
| | | |
Collapse
|
7
|
Ceccobelli S, Landi V, Senczuk G, Mastrangelo S, Sardina MT, Ben-Jemaa S, Persichilli C, Karsli T, Bâlteanu VA, Raschia MA, Poli MA, Ciappesoni G, Muchadeyi FC, Dzomba EF, Kunene NW, Lühken G, Deniskova TE, Dotsev AV, Zinovieva NA, Zsolnai A, Anton I, Kusza S, Carolino N, Santos-Silva F, Kawęcka A, Świątek M, Niżnikowski R, Špehar M, Anaya G, Granero A, Perloiro T, Cardoso P, Grande S, de Los Santos BL, Danchin-Burge C, Pasquini M, Martínez Martínez A, Delgado Bermejo JV, Lasagna E, Ciani E, Sarti FM, Pilla F. A comprehensive analysis of the genetic diversity and environmental adaptability in worldwide Merino and Merino-derived sheep breeds. Genet Sel Evol 2023; 55:24. [PMID: 37013467 PMCID: PMC10069132 DOI: 10.1186/s12711-023-00797-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2022] [Accepted: 03/24/2023] [Indexed: 04/05/2023] Open
Abstract
BACKGROUND To enhance and extend the knowledge about the global historical and phylogenetic relationships between Merino and Merino-derived breeds, 19 populations were genotyped with the OvineSNP50 BeadChip specifically for this study, while an additional 23 populations from the publicly available genotypes were retrieved. Three complementary statistical tests, Rsb (extended haplotype homozygosity between-populations), XP-EHH (cross-population extended haplotype homozygosity), and runs of homozygosity (ROH) islands were applied to identify genomic variants with potential impact on the adaptability of Merino genetic type in two contrasting climate zones. RESULTS The results indicate that a large part of the Merino's genetic relatedness and admixture patterns are explained by their genetic background and/or geographic origin, followed by local admixture. Multi-dimensional scaling, Neighbor-Net, Admixture, and TREEMIX analyses consistently provided evidence of the role of Australian, Rambouillet and German strains in the extensive gene introgression into the other Merino and Merino-derived breeds. The close relationship between Iberian Merinos and other South-western European breeds is consistent with the Iberian origin of the Merino genetic type, with traces from previous contributions of other Mediterranean stocks. Using Rsb and XP-EHH approaches, signatures of selection were detected spanning four genomic regions located on Ovis aries chromosomes (OAR) 1, 6 and 16, whereas two genomic regions on OAR6, that partially overlapped with the previous ones, were highlighted by ROH islands. Overall, the three approaches identified 106 candidate genes putatively under selection. Among them, genes related to immune response were identified via the gene interaction network. In addition, several candidate genes were found, such as LEKR1, LCORL, GHR, RBPJ, BMPR1B, PPARGC1A, and PRKAA1, related to morphological, growth and reproductive traits, adaptive thermogenesis, and hypoxia responses. CONCLUSIONS To the best of our knowledge, this is the first comprehensive dataset that includes most of the Merino and Merino-derived sheep breeds raised in different regions of the world. The results provide an in-depth picture of the genetic makeup of the current Merino and Merino-derived breeds, highlighting the possible selection pressures associated with the combined effect of anthropic and environmental factors. The study underlines the importance of Merino genetic types as invaluable resources of possible adaptive diversity in the context of the occurring climate changes.
Collapse
Affiliation(s)
- Simone Ceccobelli
- Department of Agricultural, Food and Environmental Sciences, Università Politecnica delle Marche, 60131, Ancona, Italy.
| | - Vincenzo Landi
- Department of Veterinary Medicine, University of Bari ''Aldo Moro", 70010, Valenzano, Italy
| | - Gabriele Senczuk
- Department of Agricultural, Environmental and Food Sciences, University of Molise, 86100, Campobasso, Italy
| | - Salvatore Mastrangelo
- Department of Agricultural, Food and Forest Sciences, University of Palermo, 90128, Palermo, Italy
| | - Maria Teresa Sardina
- Department of Agricultural, Food and Forest Sciences, University of Palermo, 90128, Palermo, Italy
| | - Slim Ben-Jemaa
- Laboratoire des Productions Animales et Fourragères, Institut National de la Recherche Agronomique de Tunisie, Université de Carthage, 2049, Ariana, Tunisia
| | - Christian Persichilli
- Department of Agricultural, Environmental and Food Sciences, University of Molise, 86100, Campobasso, Italy
| | - Taki Karsli
- Department of Animal Science, Faculty of Agriculture, Eskisehir Osmangazi University, 26040, Eskisehir, Turkey
| | - Valentin-Adrian Bâlteanu
- Laboratory of Genomics, Biodiversity, Animal Breeding and Molecular Pathology, Institute of Life Sciences, University of Agricultural Sciences and Veterinary Medicine of Cluj-Napoca, 400372, Cluj-Napoca, Romania
| | - María Agustina Raschia
- Instituto de Genética "Ewald A. Favret", Instituto Nacional de Tecnología Agropecuaria, CICVyA-CNIA, B1686, Hurlingham, Buenos Aires, Argentina
| | - Mario Andrés Poli
- Instituto de Genética "Ewald A. Favret", Instituto Nacional de Tecnología Agropecuaria, CICVyA-CNIA, B1686, Hurlingham, Buenos Aires, Argentina
| | - Gabriel Ciappesoni
- Instituto Nacional de Investigación Agropecuaria, 90200, Canelones, Uruguay
| | | | - Edgar Farai Dzomba
- Discipline of Genetics, School of Life Sciences, University of KwaZulu-Natal, 3209, Scottsville, Pietermaritzburg, South Africa
| | | | - Gesine Lühken
- Institute of Animal Breeding and Genetics, Justus Liebig University, 35390, Giessen, Germany
| | | | | | | | - Attila Zsolnai
- Department of Animal Breeding, Institute of Animal Science, Hungarian University of Agriculture and Life Sciences, Kaposvár Campus, 2053, Herceghalom, Hungary
| | - István Anton
- Department of Animal Breeding, Institute of Animal Science, Hungarian University of Agriculture and Life Sciences, Kaposvár Campus, 2053, Herceghalom, Hungary
| | - Szilvia Kusza
- Centre for Agricultural Genomics and Biotechnology, Faculty of Agricultural and Food Sciences and Environmental Management, University of Debrecen, 4032, Debrecen, Hungary
| | - Nuno Carolino
- Instituto Nacional de Investigação Agrária e Veterinária, 2005-048, Vale de Santarém, Portugal
| | - Fátima Santos-Silva
- Instituto Nacional de Investigação Agrária e Veterinária, 2005-048, Vale de Santarém, Portugal
| | - Aldona Kawęcka
- Department of Sheep and Goat Breeding, National Research Institute of Animal Production, 32-083, Kraków, Poland
| | - Marcin Świątek
- Department of Animal Breeding, Institute of Animal Sciences, Warsaw University of Life Sciences-SGGW, 02-786, Warsaw, Poland
| | - Roman Niżnikowski
- Department of Animal Breeding, Institute of Animal Sciences, Warsaw University of Life Sciences-SGGW, 02-786, Warsaw, Poland
| | - Marija Špehar
- Croatian Agency for Agriculture and Food, 10000, Zagreb, Croatia
| | - Gabriel Anaya
- MERAGEM Group, Department of Genetics, University of Córdoba, 14071, Córdoba, Spain
| | - Antonio Granero
- Asociación Nacional de Criadores de Ganado Merino (ACME), 28028, Madrid, Spain
| | - Tiago Perloiro
- Associação Nacional de Criadores de Ovinos da Raça Merina (ANCORME), 7005-665, Évora, Portugal
| | - Pedro Cardoso
- Associação de Produtores Agropecuários (OVIBEIRA), 6000-244, Castelo Branco, Portugal
| | - Silverio Grande
- Associazione Nazionale della Pastorizia (ASSONAPA), 00187, Rome, Italy
| | | | | | - Marina Pasquini
- Department of Agricultural, Food and Environmental Sciences, Università Politecnica delle Marche, 60131, Ancona, Italy
| | | | | | - Emiliano Lasagna
- Department of Agricultural, Food and Environmental Sciences, University of Perugia, 06121, Perugia, Italy
| | - Elena Ciani
- Department of Bioscience, Biotechnology and Biopharmaceutics, University of Bari "Aldo Moro", 70124, Bari, Italy
| | - Francesca Maria Sarti
- Department of Agricultural, Food and Environmental Sciences, University of Perugia, 06121, Perugia, Italy
| | - Fabio Pilla
- Department of Agricultural, Environmental and Food Sciences, University of Molise, 86100, Campobasso, Italy
| |
Collapse
|
8
|
Machová K, Marina H, Arranz JJ, Pelayo R, Rychtářová J, Milerski M, Vostrý L, Suárez-Vega A. Genetic diversity of two native sheep breeds by genome-wide analysis of single nucleotide polymorphisms. Animal 2023; 17:100690. [PMID: 36566708 DOI: 10.1016/j.animal.2022.100690] [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/31/2022] [Revised: 11/22/2022] [Accepted: 11/24/2022] [Indexed: 12/05/2022] Open
Abstract
Wallachian and Sumava sheep are autochthonous breeds that have undergone a significant bottleneck effect and subsequent restoration efforts. The first objective of this study was to evaluate the degree of genetic variability of both breeds and, therefore, the current management of the breeding. The second was to determine whether these two breeds still retain their genetic uniqueness in relation to each other and other breeds, despite regenerative interventions. Our data consisted of 48 individuals of Sumava and 37 individuals of Wallachian sheep. The comparison data contained 25 other breeds (primarily European) from the HapMap dataset generated by the International Sheep Genomics Consortium. When comparing all 27 breeds, the Czech breeds clustered with 15 other breeds and formed a single branch with them according to Nei's distances. At the same time, however, the clusters of both breeds were integral and easily distinguishable from the others when displayed with principal component analysis (PCA). Population substructure analysis did not show any common genetic ancestry of the Czech national breeds and breeds used for regeneration or, eventually, breeds whose ancestral population was used for regeneration. The average values of FST were higher in Wallachian sheep (FST = 0.14) than in Sumava sheep (FST = 0.08). The linkage disequilibrium (LD) extension per autosome was higher in Wallachian than in Sumava sheep. Consequently, the Ne estimates five generations ago were 68 for Sumava versus 34 for Wallachian sheep. Both native Czech breeds exhibit a wide range of inbreeding based on the excess of homozygosity (FHOM) among individuals, from -0.04 to 0.16 in Sumava and from -0.13 to 0.12 in Wallachian. Average inbreeding based on runs of homozygosity was 0.21 in Sumava and 0.27 in Wallachian. Most detected runs of homozygosity (ROH) were less than 5 Mb long for both breeds. ROH segments longer than 15 Mb were absent in Wallachian sheep. Concerning putative selection signatures, a total of 471 candidate genes in Wallachian sheep within 11 hotspots and 653 genes within 13 hotspots in Sumava sheep were identified. Czech breeds appear to be well differentiated from each other and other European breeds. Their genetic diversity is low, especially in the case of the Wallachian breed. Sumava is not so threatened by low diversity but has a larger share of the non-native gene pool.
Collapse
Affiliation(s)
- Karolína Machová
- Department of Genetics and Breeding, Czech University of Life Sciences (CZU), Prague, Czech Republic, Faculty of Agrobiology, Food and Natural Resources, Kamýcká 129, 165 00 Praha, Suchdol, Czech Republic.
| | - Héctor Marina
- Department of Animal Production, University of Leon, Leon, Spain, Veterinary Faculty, Campus de Vegazana, s/n, 24071 Leon, Spain
| | - Juan Jose Arranz
- Department of Animal Production, University of Leon, Leon, Spain, Veterinary Faculty, Campus de Vegazana, s/n, 24071 Leon, Spain
| | - Rocío Pelayo
- Department of Animal Production, University of Leon, Leon, Spain, Veterinary Faculty, Campus de Vegazana, s/n, 24071 Leon, Spain
| | - Jana Rychtářová
- Institute of Animal Science, Prague, Czech republic, Přátelství 815, 104 00 Praha, Uhříněves, Czech Republic
| | - Michal Milerski
- Institute of Animal Science, Prague, Czech republic, Přátelství 815, 104 00 Praha, Uhříněves, Czech Republic
| | - Luboš Vostrý
- Department of Genetics and Breeding, Czech University of Life Sciences (CZU), Prague, Czech Republic, Faculty of Agrobiology, Food and Natural Resources, Kamýcká 129, 165 00 Praha, Suchdol, Czech Republic
| | - Aroa Suárez-Vega
- Department of Animal Production, University of Leon, Leon, Spain, Veterinary Faculty, Campus de Vegazana, s/n, 24071 Leon, Spain
| |
Collapse
|
9
|
Wang H, Wang Q, Tan X, Wang J, Zhang J, Zheng M, Zhao G, Wen J. Estimation of genetic variability and identification of regions under selection based on runs of homozygosity in Beijing-You Chickens. Poult Sci 2022; 102:102342. [PMID: 36470032 PMCID: PMC9719870 DOI: 10.1016/j.psj.2022.102342] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2021] [Revised: 11/08/2022] [Accepted: 11/09/2022] [Indexed: 11/16/2022] Open
Abstract
The genetic composition of populations is the result of a long-term process of selection and adaptation to specific environments and ecosystems. Runs of homozygosity (ROHs) are homozygous segments of the genome where the 2 haplotypes inherited from the parents are identical. The detection of ROH can be used to describe the genetic variability and quantify the level of inbreeding in an individual. Here, we investigated the occurrence and distribution of ROHs in 40 Beijing-You Chickens from the random breeding population (BJY_C) and 40 Beijing-You Chickens from the intramuscular fat (IMF) selection population (BJY_S). Principal component analysis (PCA) and maximum likelihood (ML) analyses showed that BJY_C was completely separated from the BJY_S. The nucleotide diversity of BJY_C was higher than that of BJY_S, and the decay rate of LD of BJY_C was faster. The ROHs were identified for a total of 7,101 in BJY_C and 9,273 in BJY_S, respectively. The ROH-based inbreeding estimate (FROH) of BJY_C was 0.079, which was significantly lower than that of BJY_S (FROH = 0.114). The results were the same as the estimates of the inbreeding coefficients calculated based on homozygosity (FHOM), the correlation between uniting gametes (FUNI), and the genomic relationship matrix (FGRM). Additionally, the distribution and number of ROH islands in chromosomes of BJY_C and BJY_S were significantly different. The ROH islands of BJY_S that included genes associated with lipid metabolism and fat deposition, such as CIDEA and S1PR1, were absent in BJY_C. However, GPR161 was detected in both populations, which is a candidate gene for the formation of the unique five-finger trait in Beijing-You chickens. Our findings contributed to the understanding of the genetic diversity of random or artificially selected populations, and allowed the accurate monitoring of population inbreeding using genomic information, as well as the detection of genomic regions that affect traits under selection.
Collapse
Affiliation(s)
- Hailong Wang
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China
| | - Qiao Wang
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China
| | - Xiaodong Tan
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China
| | - Jie Wang
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China
| | - Jin Zhang
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China
| | - Maiqing Zheng
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China
| | - Guiping Zhao
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China
| | - Jie Wen
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China.
| |
Collapse
|
10
|
Qiao R, Zhang M, Zhang B, Li X, Han X, Wang K, Li X, Yang F, Hu P. Population genetic structure analysis and identification of backfat thickness loci of Chinese synthetic Yunan pigs. Front Genet 2022; 13:1039838. [PMID: 36437945 PMCID: PMC9681789 DOI: 10.3389/fgene.2022.1039838] [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: 09/08/2022] [Accepted: 10/26/2022] [Indexed: 11/10/2022] Open
Abstract
Yunan is a crossed lean meat pig breed in China. Backfat thickness is the gold standard for carcass quality grading. However, over 14 years after breed registration, the backfat of Yunan thickened and the consistency of backfat thickness decreased. Meanwhile, no genetic study has been ever performed on Yunan population. So, in this study we collected all the 120 nucleus individuals of Yunan and recorded six backfat traits of them, carried out population genetic structure analysis, selection signals analysis and genome-wide association study of Yunan pigs with the help of their founder population Duroc and Chinese native Huainan pigs, to determine the genomic loci on backfat of Yunan. Genetic diversity indexes suggested Yunan pigs had no inbreeding risk while population genetic structure showed they had few molecular pedigrees and were stratified. A total of 71 common selection signals affecting growth and fat deposition were detected by FST and XP-CLR methods. 34 significant loci associated with six backfat traits were detected, among which a 1.40 Mb region on SSC4 (20.03–21.43 Mb) were outstanding as the strong region underlying backfat. This region was common with the results of selection signature analysis, former reported QTLs for backfat and was common for different kinds of backfat traits at different development stage. ENPP2, EXT1 and SLC30A8 genes around were fat deposition related genes and were of Huainan pig’s origin, among which Type 2 diabetes related gene SLC30A8 was the most reasonable for being in a 193.21 Kb haplotype block of the 1.40 Mb region. Our results had application value for conservation, mating and breeding improvement of backfat thickness of Yunan pigs and provided evidence for a human function gene might be reproduced in pigs.
Collapse
Affiliation(s)
- Ruimin Qiao
- *Correspondence: Ruimin Qiao, ; Panyang Hu, hpy9809.@163.com
| | | | | | | | | | | | | | | | - Panyang Hu
- *Correspondence: Ruimin Qiao, ; Panyang Hu, hpy9809.@163.com
| |
Collapse
|
11
|
Marković M, Radonjić D, Zorc M, Đokić M, Marković B. Genetic Diversity of Montenegrin Local Sheep Breeds Based on Microsatellite Markers. Animals (Basel) 2022; 12:3029. [PMID: 36359153 PMCID: PMC9653887 DOI: 10.3390/ani12213029] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Revised: 10/31/2022] [Accepted: 11/01/2022] [Indexed: 10/29/2023] Open
Abstract
The Montenegrin sheep population mostly consists of local breeds and their crossbreeds that are very valuable from their genome preservation point of view. The aim of this study was the investigation of the genetic diversity of seven Montenegrin sheep breeds (Jezeropivska-JP, Sora-SOR, Zetska zuja-ZZ, Bardoka-BAR, Sjenička-SJ, Ljaba-Lj, and Piperska zuja-PIP) using 18 microsatellite sets of markers. The genotyping was done for 291 samples from seven populations using the multiplex amplification of sequences with polymerase chain reaction (PCR). The parameters of genetic diversity were estimated using several software tools. In total, 243 alleles were found, with a range of 6 to 25 by locus. The mean observed heterozygosity (Ho), polymorphism information content (PIC), and Fis values (fixation index) per marker were 0.728, 0.781, and -0.007, respectively. The mean number of alleles per breed varied from 4.889 in ZZ to 10.056 in JP. The highest Ho was estimated for JP (0.763) and the lowest for ZZ (0.640). The genetic structure showed close relations between SOR and JP, and both of them with SJ, while ZZ, LJ, and PIP were more distanced. This study provides useful indicators for the development of further in-depth studies and the creation of appropriate conservation programs.
Collapse
Affiliation(s)
- Milan Marković
- Biotechnical Faculty, University of Montenegro, Mihaila Lalića 15, 81000 Podgorica, Montenegro
| | - Dušica Radonjić
- Biotechnical Faculty, University of Montenegro, Mihaila Lalića 15, 81000 Podgorica, Montenegro
| | - Minja Zorc
- Biotechnical Faculty, University of Ljubljana, Groblje 3, 1230 Domžale, Slovenia
| | - Milena Đokić
- Biotechnical Faculty, University of Montenegro, Mihaila Lalića 15, 81000 Podgorica, Montenegro
| | - Božidarka Marković
- Biotechnical Faculty, University of Montenegro, Mihaila Lalića 15, 81000 Podgorica, Montenegro
| |
Collapse
|
12
|
Lin R, Li J, Yang Y, Yang Y, Chen J, Zhao F, Xiao T. Genome-Wide Population Structure Analysis and Genetic Diversity Detection of Four Chinese Indigenous Duck Breeds from Fujian Province. Animals (Basel) 2022; 12:ani12172302. [PMID: 36078022 PMCID: PMC9454422 DOI: 10.3390/ani12172302] [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: 07/18/2022] [Revised: 08/29/2022] [Accepted: 09/02/2022] [Indexed: 12/03/2022] Open
Abstract
Simple Summary The aim of this study was to conduct a genome-wide comparative analysis of four indigenous Chinese duck breeds (Jinding, Liancheng white, Putian black, and Shanma ducks) from Fujian Province, to understand their genetic diversity and population structure. Population parameters showed that the four indigenous breeds were separated groups. Five genomic regions are presented as hotspots of autozygosity among these indigenous duck breeds, with candidate genes involved in muscle growth, pigmentation, and neuroregulation. Genomic information may play a vital role in the improvement of conservation strategies. Abstract The assessment of population genetic structure is the basis for understanding the genetic information of indigenous breeds and is important for the protection and management of indigenous breeds. However, the population genetic differentiation of many local breeds still remains unclear. Here, we performed a genome-wide comparative analysis of Jinding, Liancheng white, Putian black, and Shanma ducks based on the genomic sequences using RAD sequencing to understand their population structure and genetic diversity. The population parameters showed that there were obvious genetic differences among the four indigenous breeds, which were separated groups. Among them, Liancheng white and Shanma ducks may come from the same ancestor because the phylogenetic tree forms three tree trunks. In addition, during the runs of homozygosity (ROH), we found that the average inbreeding coefficient of Liancheng white and Putian black ducks was the lowest and the highest, respectively. Five genomic regions were considered to be the hotspots of autozygosity among these indigenous duck breeds, and the candidate genes involved a variety of potential variations, such as muscle growth, pigmentation, and neuroregulation. These findings provide insights into the further improvement and conservation of Fujian duck breeds.
Collapse
|
13
|
Salek Ardestani S, Zandi MB, Vahedi SM, Janssens S. Population structure and genomic footprints of selection in five major Iranian horse breeds. Anim Genet 2022; 53:627-639. [PMID: 35919961 DOI: 10.1111/age.13243] [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: 03/03/2022] [Revised: 06/08/2022] [Accepted: 07/04/2022] [Indexed: 11/28/2022]
Abstract
The genetic structure and characteristics of Iranian native breeds are yet to be comprehensibly investigated and studied. Therefore, we employed genomic information of 364 Iranian native horses representing the Asil (n = 109), Caspian (n = 40), Dareshuri (n = 44), Kurdish (n = 95), and Turkoman (n = 76) breeds to reveal the genetic structure and characteristics. For these and 19 other horse breeds, principal component analysis, Bayesian model-based, Neighbor-Net, and bootstrap-based TreeMix approaches were applied to investigate and compare their genetic structure. Additionally, three haplotype-based methods including haplotype homozygosity pooled, integrated haplotype score, and number of segregating sites by length were applied to trace genomic footprints of selection of Asil, Caspian, Dareshuri, Kurdish, and Turkoman groups. Then, the Mahalanobis distance based on the negative-log10 rank-based P-values was estimated based on the haplotype homozygosity pooled, integrated haplotype score, and number of segregating sites by length values. Asil, Caspian, Dareshuri, Kurdish, and Turkoman can be categorized into five different genetic clusters. Based on the top 1% of Mahalanobis distance based on the negative-log10 rank-based P-values of SNPs, we identified 24 SNPs formerly reported to be associated with different traits and >100 genes undergoing selection pressures in Asil, Caspian, Dareshuri, Kurdish, and Turkoman. The detected QTL undergoing selection pressures were associated with withers height, equine metabolic syndrome, overall body size, insect bite hypersensitivity, guttural pouch tympany, white markings, Rhodococcus equi infection, jumping test score, alternate gaits, and body weight traits. Our findings will aid to have a better perspective of the genetic characteristics and population structure of Asil, Caspian, Dareshuri, Kurdish, and Turkoman horses as Iranian native horse breeds.
Collapse
Affiliation(s)
| | | | - Seyed Milad Vahedi
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, Nova Scotia, Canada
| | - Steven Janssens
- Department Biosystems, Center Animal Breeding and Genetics, KU Leuven, Leuven, Belgium
| |
Collapse
|
14
|
Drzaic I, Curik I, Lukic B, Shihabi M, Li MH, Kantanen J, Mastrangelo S, Ciani E, Lenstra JA, Cubric-Curik V. High-Density Genomic Characterization of Native Croatian Sheep Breeds. Front Genet 2022; 13:940736. [PMID: 35910220 PMCID: PMC9337876 DOI: 10.3389/fgene.2022.940736] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2022] [Accepted: 06/16/2022] [Indexed: 11/13/2022] Open
Abstract
A recent comprehensive genomic analysis based on 50K SNP profiles has shown that the regional Balkan sheep populations have considerable genetic overlap but are distinctly different from surrounding breeds. All eight Croatian sheep breeds were represented by a small number of individuals per breed. Here, we genotyped 220 individuals representing the native Croatian sheep breeds (Istrian Sheep, Krk Island Sheep, Cres Island Sheep, Rab Island Sheep, Lika Pramenka, Pag Island Sheep, Dalmatian Pramenka, Dubrovnik Sheep) and mouflon using the Ovine Infinium® HD SNP BeadChip (606,006 SNPs). In addition, we included publicly available Balkan Pramenka and other Mediterranean sheep breeds. Our analyses revealed the complex population structure of Croatian sheep breeds and their origin and geographic barriers (island versus mainland). Migration patterns confirmed the historical establishment of breeds and the pathways of gene flow. Inbreeding coefficients (FROH>2 Mb) between sheep populations ranged from 0.025 to 0.070, with lower inbreeding coefficients observed in Dalmatian Pramenka and Pag Island Sheep and higher inbreeding in Dubrovnik sheep. The estimated effective population size ranged from 61 to 1039 for Krk Island Sheep and Dalmatian Pramenka, respectively. Higher inbreeding levels and lower effective population size indicate the need for improved conservation management to maintain genetic diversity in some breeds. Our results will contribute to breeding and conservation strategies of native Croatian sheep breeds.
Collapse
Affiliation(s)
- Ivana Drzaic
- Department of Animal Science, University of Zagreb Faculty of Agriculture, Zagreb, Croatia
- *Correspondence: Ivana Drzaic, ; Vlatka Cubric-Curik,
| | - Ino Curik
- Department of Animal Science, University of Zagreb Faculty of Agriculture, Zagreb, Croatia
| | - Boris Lukic
- Department of Animal Production and Biotechnology, Faculty of Agrobiotechnical Sciences Osijek, Chair for Domestic Animal Breeding and Genetics, J. J. Strossmayer University of Osijek, Osijek, Croatia
| | - Mario Shihabi
- Department of Animal Science, University of Zagreb Faculty of Agriculture, Zagreb, Croatia
| | - Meng-Hua Li
- College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Juha Kantanen
- Production Systems, Natural Resources Institute Finland (Luke), Jokioinen, Finland
| | - Salvatore Mastrangelo
- Dipartimento di Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Elena Ciani
- Dipartimento di Bioscienze, Biotecnologie e Biofarmaceutica, Universita Degli Studi di Bari “Aldo Moro”, Bari, Italy
| | | | - Vlatka Cubric-Curik
- Department of Animal Science, University of Zagreb Faculty of Agriculture, Zagreb, Croatia
- *Correspondence: Ivana Drzaic, ; Vlatka Cubric-Curik,
| |
Collapse
|
15
|
Zhumadillayev N, Dossybayev K, Khamzina A, Kapasuly T, Khamzina Z, Tlevlesov N. SNP Genotyping Characterizes the Genome Composition of the New Baisary Fat-Tailed Sheep Breed. Animals (Basel) 2022; 12:ani12111468. [PMID: 35681932 PMCID: PMC9179407 DOI: 10.3390/ani12111468] [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: 02/11/2022] [Revised: 05/27/2022] [Accepted: 05/30/2022] [Indexed: 12/04/2022] Open
Abstract
Lamb meat has become increasingly popular in several nations during the last few decades, especially in Kazakhstan. Due to the rising demand for lamb meat, our sheep breeders developed a new fat-tailed sheep and named the breed Baisary. Animals of the Baisary breed are characterized by a large physique, strong constitution, stretched body, deep and wide chest, medium or large-sized fat tail, long legs (height at the withers of adult rams 85–100 cm, sheep 75–90 cm), long lanceolate ears and strong hooves. Lambs of the Baisary breed surpass their peers of the original parent breeds by 15–20% in live weight at the weaning period. To characterize the genetic structure of Baisary sheep and compare it with the ancestral breeds, we genotyped 247 individuals from five sheep breeds with Ovine SNP50K. The estimated private allelic richness ranged from 0.0030 to 0.0047, with the minimum and maximum provided by the Gissar (Giss1) and Kazakh meat-wool breeds, respectively. The highest and lowest FIS values, meanwhile, were observed in the Afghan fat-tailed population and Baisary sheep, respectively. The calculated inbreeding coefficient showed that Edilbay and Baisary sheep have excess heterozygosity. According to principal components analysis, Baisary are close to Gissar populations, the Afghan fat-tailed breed and Edilbay sheep. These results were consistent with the Admixture and phylogenetic analysis. Overall, our results indicated that Baisary sheep differ genetically from their progenitors.
Collapse
Affiliation(s)
- Narzhan Zhumadillayev
- Test Center, Kazakh Scientific Research Institute of Animal Husbandry and Forage Production, Zhandosov, 51, Almaty 050035, Kazakhstan; (N.Z.); (Z.K.); (N.T.)
| | - Kairat Dossybayev
- Laboratory of Genetics and Cytogenetics, RSE “Institute of Genetics and Physiology” CS MES RK, Al-Farabi Avenue, 93, Almaty 050060, Kazakhstan;
- Department of Molecular Biology and Genetics, Faculty of Biology and Biotechnology, Al-Farabi Kazakh National University, Almaty 050040, Kazakhstan
- Correspondence:
| | - Aigerim Khamzina
- Green Biotechnology and Cell Engineering Laboratory, Kazakh National Agrarian Research University, Almaty 050010, Kazakhstan;
| | - Tilek Kapasuly
- Laboratory of Genetics and Cytogenetics, RSE “Institute of Genetics and Physiology” CS MES RK, Al-Farabi Avenue, 93, Almaty 050060, Kazakhstan;
| | - Zhangylsyn Khamzina
- Test Center, Kazakh Scientific Research Institute of Animal Husbandry and Forage Production, Zhandosov, 51, Almaty 050035, Kazakhstan; (N.Z.); (Z.K.); (N.T.)
- Department of Molecular Biology and Genetics, Faculty of Biology and Biotechnology, Al-Farabi Kazakh National University, Almaty 050040, Kazakhstan
| | - Nurlan Tlevlesov
- Test Center, Kazakh Scientific Research Institute of Animal Husbandry and Forage Production, Zhandosov, 51, Almaty 050035, Kazakhstan; (N.Z.); (Z.K.); (N.T.)
| |
Collapse
|
16
|
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
|
17
|
Araujo AC, Carneiro PLS, Oliveira HR, Schenkel FS, Veroneze R, Lourenco DAL, Brito LF. A Comprehensive Comparison of Haplotype-Based Single-Step Genomic Predictions in Livestock Populations With Different Genetic Diversity Levels: A Simulation Study. Front Genet 2021; 12:729867. [PMID: 34721524 PMCID: PMC8551834 DOI: 10.3389/fgene.2021.729867] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2021] [Accepted: 09/07/2021] [Indexed: 11/13/2022] Open
Abstract
The level of genetic diversity in a population is inversely proportional to the linkage disequilibrium (LD) between individual single nucleotide polymorphisms (SNPs) and quantitative trait loci (QTLs), leading to lower predictive ability of genomic breeding values (GEBVs) in high genetically diverse populations. Haplotype-based predictions could outperform individual SNP predictions by better capturing the LD between SNP and QTL. Therefore, we aimed to evaluate the accuracy and bias of individual-SNP- and haplotype-based genomic predictions under the single-step-genomic best linear unbiased prediction (ssGBLUP) approach in genetically diverse populations. We simulated purebred and composite sheep populations using literature parameters for moderate and low heritability traits. The haplotypes were created based on LD thresholds of 0.1, 0.3, and 0.6. Pseudo-SNPs from unique haplotype alleles were used to create the genomic relationship matrix ( G ) in the ssGBLUP analyses. Alternative scenarios were compared in which the pseudo-SNPs were combined with non-LD clustered SNPs, only pseudo-SNPs, or haplotypes fitted in a second G (two relationship matrices). The GEBV accuracies for the moderate heritability-trait scenarios fitting individual SNPs ranged from 0.41 to 0.55 and with haplotypes from 0.17 to 0.54 in the most (Ne ≅ 450) and less (Ne < 200) genetically diverse populations, respectively, and the bias fitting individual SNPs or haplotypes ranged between -0.14 and -0.08 and from -0.62 to -0.08, respectively. For the low heritability-trait scenarios, the GEBV accuracies fitting individual SNPs ranged from 0.24 to 0.32, and for fitting haplotypes, it ranged from 0.11 to 0.32 in the more (Ne ≅ 250) and less (Ne ≅ 100) genetically diverse populations, respectively, and the bias ranged between -0.36 and -0.32 and from -0.78 to -0.33 fitting individual SNPs or haplotypes, respectively. The lowest accuracies and largest biases were observed fitting only pseudo-SNPs from blocks constructed with an LD threshold of 0.3 (p < 0.05), whereas the best results were obtained using only SNPs or the combination of independent SNPs and pseudo-SNPs in one or two G matrices, in both heritability levels and all populations regardless of the level of genetic diversity. In summary, haplotype-based models did not improve the performance of genomic predictions in genetically diverse populations.
Collapse
Affiliation(s)
- Andre C Araujo
- Postgraduate Program in Animal Sciences, State University of Southwestern Bahia, Itapetinga, Brazil.,Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
| | - Paulo L S Carneiro
- Department of Biology, State University of Southwestern Bahia, Jequié, Brazil
| | - Hinayah R Oliveira
- Department of Animal Sciences, Purdue University, West Lafayette, IN, United States.,Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada
| | - Flavio S Schenkel
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada
| | - Renata Veroneze
- Department of Animal Sciences, Federal University of Viçosa, Viçosa, Brazil
| | - Daniela A L Lourenco
- Department of Animal and Dairy Science, University of Georgia, Athens, GA, United States
| | - Luiz F Brito
- Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
| |
Collapse
|
18
|
Biodiversity of Russian Local Sheep Breeds Based on Pattern of Runs of Homozygosity. DIVERSITY 2021. [DOI: 10.3390/d13080360] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Russian sheep breeds traditionally raised in specific environments are valuable parts of sociocultural heritage and economic component of the regions. However, the import of commercial breeds negatively influences the population sizes of local sheep populations and might lead to biodiversity loss. Estimation of the runs of homozygosity (ROH) in local sheep genomes is an informative tool to address their current genetic state. In this work, we aimed to address the ROH distribution and to estimate genome inbreeding based on SNP data to evaluate genetic diversity in Russian local sheep breeds. Materials for this study included SNP-genotypes from twenty-seven Russian local sheep breeds which were generated using the Illumina OvineSNP50 BeadChip (n = 391) or the Illumina Ovine Infinium HD BeadChip (n = 315). A consecutive runs method was used to calculate ROH which were estimated for each animal and then categorized in the ROH length classes. The ROH were found in all breeds. The mean ROH length varied from 86 to 280 Mb, while the ROH number ranged from 37 to 123. The genomic inbreeding coefficient varied from 0.033 to 0.106. Our findings provide evidence of low to moderate genomic inbreeding in major local sheep populations.
Collapse
|
19
|
Moosanezhad Khabisi M, Asadi Foozi M, Lv FH, Esmailizadeh A. Genome-wide DNA arrays profiling unravels the genetic structure of Iranian sheep and pattern of admixture with worldwide coarse-wool sheep breeds. Genomics 2021; 113:3501-3511. [PMID: 34293474 DOI: 10.1016/j.ygeno.2021.07.019] [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: 02/19/2021] [Revised: 05/18/2021] [Accepted: 07/17/2021] [Indexed: 10/20/2022]
Abstract
Archaeological and genetic evidence show that sheep were originally domesticated in area around the North of Zagros mountains, North-west of Iran. The Persian plateau exhibits a variety of native sheep breeds with a common characteristic of coarse-wool production. Therefore, knowledge about the genetic structure and diversity of Iranian sheep and genetic connections with other sheep breeds is of great interest. To this end, we genotyped 154 samples from 11 sheep breeds distributed across Iran with the Ovine Infinium HD SNP 600 K BeadChip array, and analyzed this dataset combined with the retrieved data of 558 samples from 19 worldwide coarse-wool sheep breeds. The average genetic diversity ranged from 0.315 to 0.354, while the FST values ranged from 0.016 to 0.177 indicating a low differentiation of Iranian sheep. Analysis of molecular variance showed that 90.21 and 9.79% of the source of variation were related to differences within and between populations, respectively. Our results indicated that the coarse-wool sheep from Europe were clearly different from those of the Asia. Accordingly, the Asiatic mouflon was positioned between Asian and European countries. In addition, we found that the genetic background of Iranian sheep is present in sheep from China and Kyrgyzstan, as well as India. The revealed admixture patterns of the Iranian sheep and other coarse-wool sheep breeds probably resulted from the expansion of nomads and through the Silk Road trade network.
Collapse
Affiliation(s)
- Mozhdeh Moosanezhad Khabisi
- Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, PB 76169-133 Kerman, Iran
| | - Masood Asadi Foozi
- Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, PB 76169-133 Kerman, Iran
| | - Feng-Hua Lv
- College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Ali Esmailizadeh
- Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, PB 76169-133 Kerman, Iran.
| |
Collapse
|
20
|
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
|
21
|
British Sheep Breeds as a Part of World Sheep Gene Pool Landscape: Looking into Genomic Applications. Animals (Basel) 2021; 11:ani11040994. [PMID: 33916207 PMCID: PMC8103502 DOI: 10.3390/ani11040994] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2021] [Revised: 03/27/2021] [Accepted: 03/29/2021] [Indexed: 01/18/2023] Open
Abstract
Sheep farming has been an important sector of the UK's economy and rural life for many centuries. It is the favored source of wool, meat and milk products. In the era of exponential progress in genomic technologies, we can now address the questions of what is special about UK sheep breed genotypes and how they differ genetically form one another and from other countries. We can reflect how their natural history has been determined at the level of their genetic code and what traces have been left in their genomes because of selection for phenotypic traits. These include adaptability to certain environmental conditions and management, as well as resistance to disease. Application of these advancements in genetics and genomics to study sheep breeds of British domestic selection has begun and will continue in order to facilitate conservation solutions and production improvement.
Collapse
|
22
|
Calderon M, More MJ, Gutierrez GA, Ponce de León FA. Development of a 76k Alpaca ( Vicugna pacos) Single Nucleotide Polymorphisms (SNPs) Microarray. Genes (Basel) 2021; 12:genes12020291. [PMID: 33669871 PMCID: PMC7923280 DOI: 10.3390/genes12020291] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2020] [Revised: 02/13/2021] [Accepted: 02/16/2021] [Indexed: 11/16/2022] Open
Abstract
Small farm producers’ sustenance depends on their alpaca herds and the production of fiber. Genetic improvement of fiber characteristics would increase their economic benefits and quality of life. The incorporation of molecular marker technology could overcome current limitations for the implementation of genetic improvement programs. Hence, the aim of this project was the generation of an alpaca single nucleotide polymorphism (SNP) microarray. A sample of 150 Huacaya alpacas from four farms, two each in Puno and Cerro de Pasco were used for SNP discovery by genotyping by sequencing (GBS). Reduced representation libraries, two per animal, were produced after DNA digestion with ApeK1 and double digestion with Pst1-Msp1. Ten alpaca genomes, sequenced at depths between 12× to 30×, and the VicPac3.1 reference genome were used for read alignments. Bioinformatics analysis discovered 76,508 SNPs included in the microarray. Candidate genes SNPs (302) for fiber quality and color are also included. The microarray SNPs cover 90.5% of the genome length with a density of about 39 ± 2.51 SNPs/Mb of DNA at an average interval of 26.45 ± 18.57 kbp. The performance was evaluated by genotyping 30 family trios and comparing them to their pedigrees, as well as comparing microarray to GBS genotypes. Concordance values of 0.93 and 0.94 for ApeK1 and Pst1-Msp1 generated SNPs were observed. Similarly, 290 fiber quality and color candidate gene SNPs were validated. Availability of this microarray will facilitate genome-wide association studies, marker-assisted selection and, in time, genomic selection.
Collapse
Affiliation(s)
- Marcos Calderon
- Facultad de Zootecnia, Universidad Nacional Agraria La Molina, Lima 15024, Peru; (M.C.); (M.J.M.); (G.A.G.)
- Escuela de Formación Profesional de Zootecnia, Facultad de Ciencias Agropecuarias, Universidad Nacional Daniel Alcídes Carrión, Cerro de Pasco 19001, Peru
| | - Manuel J. More
- Facultad de Zootecnia, Universidad Nacional Agraria La Molina, Lima 15024, Peru; (M.C.); (M.J.M.); (G.A.G.)
| | - Gustavo A. Gutierrez
- Facultad de Zootecnia, Universidad Nacional Agraria La Molina, Lima 15024, Peru; (M.C.); (M.J.M.); (G.A.G.)
| | - Federico Abel Ponce de León
- Department of Animal Science, University of Minnesota, Minneapolis, MN 55108, USA
- Correspondence: ; Tel.: +1-612-419-7870
| |
Collapse
|
23
|
Eydivandi S, Roudbar MA, Karimi MO, Sahana G. Genomic scans for selective sweeps through haplotype homozygosity and allelic fixation in 14 indigenous sheep breeds from Middle East and South Asia. Sci Rep 2021; 11:2834. [PMID: 33531649 PMCID: PMC7854752 DOI: 10.1038/s41598-021-82625-2] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2020] [Accepted: 01/22/2021] [Indexed: 01/30/2023] Open
Abstract
The performance and productivity of livestock have consistently improved by natural and artificial selection over the centuries. Both these selections are expected to leave patterns on the genome and lead to changes in allele frequencies, but natural selection has played the major role among indigenous populations. Detecting selective sweeps in livestock may assist in understanding the processes involved in domestication, genome evolution and discovery of genomic regions associated with economically important traits. We investigated population genetic diversity and selection signals in this study using SNP genotype data of 14 indigenous sheep breeds from Middle East and South Asia, including six breeds from Iran, namely Iranian Balochi, Afshari, Moghani, Qezel, Zel, and Lori-Bakhtiari, three breeds from Afghanistan, namely Afghan Balochi, Arabi, and Gadik, three breeds from India, namely Indian Garole, Changthangi, and Deccani, and two breeds from Bangladesh, namely Bangladeshi Garole and Bangladesh East. The SNP genotype data were generated by the Illumina OvineSNP50 Genotyping BeadChip array. To detect genetic diversity and population structure, we used principal component analysis (PCA), admixture, phylogenetic analyses, and Runs of homozygosity. We applied four complementary statistical tests, FST (fixation index), xp-EHH (cross-population extended haplotype homozygosity), Rsb (extended haplotype homozygosity between-populations), and FLK (the extension of the Lewontin and Krakauer) to detect selective sweeps. Our results not only confirm the previous studies but also provide a suite of novel candidate genes involved in different traits in sheep. On average, FST, xp-EHH, Rsb, and FLK detected 128, 207, 222, and 252 genomic regions as candidates for selective sweeps, respectively. Furthermore, nine overlapping candidate genes were detected by these four tests, especially TNIK, DOCK1, USH2A, and TYW1B which associate with resistance to diseases and climate adaptation. Knowledge of candidate genomic regions in sheep populations may facilitate the identification and potential exploitation of the underlying genes in sheep breeding.
Collapse
Affiliation(s)
- Sirous Eydivandi
- Department of Animal Science, Behbahan Branch, Islamic Azad University, Behbahan, Iran.
- Center for Quantitative Genetics and Genomics, Faculty of Technical Sciences, Aarhus University, 8830, Tjele, Denmark.
| | - Mahmoud Amiri Roudbar
- Department of Animal Science, Safiabad-Dezful Agricultural and Natural Resources Research and Education Center, Agricultural Research, Education and Extension Organization (AREEO), Dezful, Iran
| | - Mohammad Osman Karimi
- Department of Animal Science, Faculty of Agriculture, Herat University, Herat, Afghanistan
| | - Goutam Sahana
- Center for Quantitative Genetics and Genomics, Faculty of Technical Sciences, Aarhus University, 8830, Tjele, Denmark
| |
Collapse
|
24
|
Saravanan K, Panigrahi M, Kumar H, Bhushan B, Dutt T, Mishra B. Genome-wide analysis of genetic diversity and selection signatures in three Indian sheep breeds. Livest Sci 2021. [DOI: 10.1016/j.livsci.2020.104367] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|
25
|
Wu F, Sun H, Lu S, Gou X, Yan D, Xu Z, Zhang Z, Qadri QR, Zhang Z, Wang Z, Chen Q, Li M, Wang X, Dong X, Wang Q, Pan Y. Genetic Diversity and Selection Signatures Within Diannan Small-Ear Pigs Revealed by Next-Generation Sequencing. Front Genet 2020; 11:733. [PMID: 32849777 PMCID: PMC7406676 DOI: 10.3389/fgene.2020.00733] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2020] [Accepted: 06/16/2020] [Indexed: 12/26/2022] Open
Abstract
Genetic characterization of Chinese indigenous pig breeds is essential to promote scientific conservation and sustainable development of pigs. Here, we systematically surveyed the genomes of 75 unrelated Diannan small-ear (DSE) pigs from three diverse regions (Yingjiang County, Jinping County, and Sipsongpanna in Yunnan Province) to describe their population structures, genetic diversity, inbreeding coefficients, and selection signatures. First, these individuals were sequenced and genotyped using the genome reducing and sequencing (GGRS) protocol. A total of 438,038 autosomal single-nucleotide polymorphisms (SNPs) were obtained and used for subsequent statistical analysis. The results showed that these DSE pigs were clearly differentiated into three separate clades revealed by the population structure and principal component analysis, which is consistent with their geographical origins. Diannan small-ear pigs owned lower genetic diversity when compared with some other pig breeds, which demonstrated the need to strengthen the conservation strategies for DSE pigs. In addition, the inbreeding coefficients based on runs of homozygosity (ROH) length (F ROH) were calculated in each ROH length categories, respectively. And the results indicated that the ancient (up to 50 generations ago) inbreeding had greater impacts than recent (within the last five generations) inbreeding within DSE pigs. Some candidate selection signatures within the DSE pig population were detected through the ROH islands and integrated haplotype homozygosity score (iHS) methods. And genes associated with meat quality (COL15A1, RPL3L, and SLC9A3R2), body size (PALM2-AKAP2, NANS, TRAF7, and PACSIN1), adaptability (CLDN9 and E4F1), and appetite (GRM4) were identified. These findings can help to understand the genetic characteristics and provide insights into the molecular background of special phenotypes of DSE pigs to promote conservation and sustainability of the breed.
Collapse
Affiliation(s)
- Fen Wu
- Department of Animal Breeding and Reproduction, College of Animal Sciences, Zhejiang University, Hangzhou, China
| | - Hao Sun
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China
| | - Shaoxiong Lu
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China
| | - Xiao Gou
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China
| | - Dawei Yan
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China
| | - Zhong Xu
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China
| | - Zhenyang Zhang
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China
| | - Qamar Raza Qadri
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China
| | - Zhe Zhang
- Department of Animal Breeding and Reproduction, College of Animal Sciences, Zhejiang University, Hangzhou, China
| | - Zhen Wang
- Department of Animal Breeding and Reproduction, College of Animal Sciences, Zhejiang University, Hangzhou, China
| | - Qiang Chen
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China
| | - Mingli Li
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China
| | - Xiaoyi Wang
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China
| | - Xinxing Dong
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China
| | - Qishan Wang
- Department of Animal Breeding and Reproduction, College of Animal Sciences, Zhejiang University, Hangzhou, China
| | - Yuchun Pan
- Department of Animal Breeding and Reproduction, College of Animal Sciences, Zhejiang University, Hangzhou, China
| |
Collapse
|