1
|
Zhang X, Yang L, Xie Z, Gan J, Zhu P, Song J, Kang H, Zhang Z, Liu L, Xiang H, Li H. Genetic assessment and candidate genes identification for breed-specific characteristics of Qingyuan partridge chicken based on runs of homozygosity. BMC Genomics 2024; 25:577. [PMID: 38858651 PMCID: PMC11163754 DOI: 10.1186/s12864-024-10492-y] [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: 10/17/2023] [Accepted: 06/03/2024] [Indexed: 06/12/2024] Open
Abstract
BACKGROUND Several core breeding and supporting lines of the Qingyuan partridge chicken, a representative local chicken breed in China, have been developed over 20 years. Consequently, its economic traits related to growth and reproduction have been significantly improved by breeding selection and commercial utilization, but some characteristic traits, such as partridge feathers, high meat quality and sufficient flavor, have always been retained. However, effective methods for genetic assessment and functional gene exploration of similar trait groups are lacking. The presence of identical haplotype fragments transmitted from parent to offspring results in runs of homozygosity (ROH), which offer an efficient solution. In this study, genomes of 134 Qingyuan partridge chickens representing two breeding populations and one preserved population were re-sequenced to evaluate the genetic diversity and explore functional genes by analyzing the diversity, distribution, and frequency of ROH. RESULTS The results showed a low level of genomic linkage and degree of inbreeding within both the bred and preserved populations, suggesting abundant genetic diversity and an adequate genetic potential of the Qingyuan partridge chicken. Throughout the long-term selection process, 21 genes, including GLI3, ANO5, BLVRA, EFNB2, SLC5A12, and SVIP, associated with breed-specific characteristics were accumulated within three ROH islands, whereas another 21 genes associated with growth traits including IRX1, IRX2, EGFR, TPK1, NOVA1, BDNF and so on were accumulated within five ROH islands. CONCLUSIONS These findings provide new insights into the genetic assessment and identification of genes with breed-specific and selective characteristics, offering a solid genetic basis for breeding and protection of Qingyuan partridge chickens.
Collapse
Affiliation(s)
- Xing Zhang
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China
| | - Liu Yang
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China
- Key Laboratory of National Forestry and Grassland Administration on Conservation Biology of Rare Animals in the Giant, China Conservation and Research Centre for the Giant Panda, Panda National Park, Chengdu, 611830, China
| | - Zhuojun Xie
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China
| | - Jiankang Gan
- Guangdong Tinoo's Foods Group Co., Ltd, Qingyuan, 511827, China
| | - Piao Zhu
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China
| | - Jiani Song
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China
| | - Huimin Kang
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China
| | - Zhengfen Zhang
- Guangdong Tinoo's Foods Group Co., Ltd, Qingyuan, 511827, China
| | - Lingbin Liu
- College of Animal Science and Technology, Southwest University, Chongqing, 400715, China
| | - Hai Xiang
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China.
| | - Hua Li
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, School of Life Science and Engineering, Foshan University, Foshan, 528225, China.
- Guangdong Tinoo's Foods Group Co., Ltd, Qingyuan, 511827, China.
| |
Collapse
|
2
|
Barmentlo NWG, Meirmans PG, Stiver WH, Yarkovich JG, McCann BE, Piaggio AJ, Wright D, Smyser TJ, Bosse M. Natural selection on feralization genes contributed to the invasive spread of wild pigs throughout the United States. Mol Ecol 2024; 33:e17383. [PMID: 38747342 DOI: 10.1111/mec.17383] [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: 01/30/2024] [Revised: 04/23/2024] [Accepted: 04/30/2024] [Indexed: 06/07/2024]
Abstract
Despite a long presence in the contiguous United States (US), the distribution of invasive wild pigs (Sus scrofa × domesticus) has expanded rapidly since the 1980s, suggesting a more recent evolutionary shift towards greater invasiveness. Contemporary populations of wild pigs represent exoferal hybrid descendants of domestic pigs and European wild boar, with such hybridization expected to enrich genetic diversity and increase the adaptive potential of populations. Our objective was to characterize how genetic enrichment through hybridization increases the invasiveness of populations by identifying signals of selection and the ancestral origins of selected loci. Our study focused on invasive wild pigs within Great Smoky Mountains National Park, which represents a hybrid population descendent from the admixture of established populations of feral pigs and an introduction of European wild boar to North America. Accordingly, we genotyped 881 wild pigs with multiple high-density single-nucleotide polymorphism (SNP) arrays. We found 233 markers under putative selection spread over 79 regions across 16 out of 18 autosomes, which contained genes involved in traits affecting feralization. Among these, genes were found to be related to skull formation and neurogenesis, with two genes, TYRP1 and TYR, also encoding for crucial melanogenesis enzymes. The most common haplotypes associated with regions under selection for the Great Smoky Mountains population were also common among other populations throughout the region, indicating a key role of putatively selective variants in the fitness of invasive populations. Interestingly, many of these haplotypes were absent among European wild boar reference genotypes, indicating feralization through genetic adaptation.
Collapse
Affiliation(s)
- Niek W G Barmentlo
- Section Ecology & Evolution, Amsterdam Institute for Life and Environment (A-LIFE), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
- Institute for Biodiversity and Ecosystem Dynamics, University of Amsterdam, Amsterdam, The Netherlands
| | - Patrick G Meirmans
- Institute for Biodiversity and Ecosystem Dynamics, University of Amsterdam, Amsterdam, The Netherlands
| | | | | | - Blake E McCann
- Theodore Roosevelt National Park, Medora, North Dakota, USA
| | | | - Dominic Wright
- Department of Physics, Chemistry and Biology, Linköping University, Linköping, Sweden
| | - Timothy J Smyser
- USDA APHIS WS National Wildlife Research Center, Fort Collins, Colorado, USA
| | - Mirte Bosse
- Section Ecology & Evolution, Amsterdam Institute for Life and Environment (A-LIFE), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
- Wageningen University & Research - Animal Breeding and Genomics, Wageningen, The Netherlands
| |
Collapse
|
3
|
Kar D, Ganguly I, Singh S, Bhatia AK, Dixit SP. Genome-wide runs of homozygosity signatures in diverse Indian goat breeds. 3 Biotech 2024; 14:81. [PMID: 38375512 PMCID: PMC10874352 DOI: 10.1007/s13205-024-03921-y] [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/22/2023] [Accepted: 01/05/2024] [Indexed: 02/21/2024] Open
Abstract
The present study analyzed ROH and consensus ROH regions in 102 animals of eleven diverse Indian goat (Capra hircus) breeds using whole genome sequencing. A total of 51,705 ROH and 21,271 consensus regions were identified. The mean number of ROH per animal was highest in the meat breed, Jharkhand Black (2693) and lowest in the pashmina breed, Changthangi (60). The average length of ROH (ALROH) was maximum in Kanniadu (974.11 Kb) and minimum in Tellicherry (146.98 Kb). Long ROH is typically associated with more recent inbreeding, whereas short ROH is connected to more ancient inbreeding. The overall ROH-based genomic inbreeding (FROH) was highest for Jharkhand Black (0.602) followed by Kanniadu (0.120) and Sangamneri (0.108) among all breeds. FROH of Jharkhand Black was higher than Kanniadu up to 5 Mb ROH length category. However, in > 20 Mb ROH length category, Kanniadu (0.98) exhibited significantly higher FROH than Jharkhand Black (0.46). This implies that Kanniadu had higher levels of recent inbreeding than Jharkhand Black. Despite this, due to the presence of both recent and ancient inbreeding, Jharkhand Black demonstrated higher overall FROH compared to Kanniadu. ROH patterns revealed dual purpose (meat and dairy) and pashmina breeds as less consanguineous while recent inbreeding was apparent in meat breeds. Analysis of ROH consensus regions identified selection sweeps in key genes governing intramuscular fat deposition, meat tenderisation, lean meat production and carcass weight (CDK4, ALOX15, CASP9, PRDM16, DVL1) in meat breeds; milk fat percentage and mammary gland development (POLD1, NOTCH2, ARHGAP35) in dual purpose (meat and dairy) breeds; while cold adaptation and hair follicle development (APOBEC1, DNAJC3, F2RL1, FGF9) in pashmina breed. MAPK, RAS, BMP and Wnt signaling pathways associated with hair follicle morphogenesis in Changthangi were also identified. PCA analysis based on ROH consensus regions revealed that meat breeds are more diverse than other goat breeds/populations. Supplementary Information The online version contains supplementary material available at 10.1007/s13205-024-03921-y.
Collapse
Affiliation(s)
- Dibyasha Kar
- Division of Animal Genetics, ICAR-National Bureau of Animal Genetic Resources, Karnal, Haryana 132001 India
- Division of Animal Genetics and Breeding, ICAR-National Dairy Research Institute, Karnal, Haryana 132001 India
| | - Indrajit Ganguly
- Division of Animal Genetics, ICAR-National Bureau of Animal Genetic Resources, Karnal, Haryana 132001 India
| | - Sanjeev Singh
- Division of Animal Genetics, ICAR-National Bureau of Animal Genetic Resources, Karnal, Haryana 132001 India
| | - Avnish Kumar Bhatia
- Division of Animal Genetics, ICAR-National Bureau of Animal Genetic Resources, Karnal, Haryana 132001 India
| | - S. P. Dixit
- Division of Animal Genetics, ICAR-National Bureau of Animal Genetic Resources, Karnal, Haryana 132001 India
| |
Collapse
|
4
|
Yuan J, Kitchener AC, Lackey LB, Sun T, Jiangzuo Q, Tuohetahong Y, Zhao L, Yang P, Wang G, Huang C, Wang J, Hou W, Liu Y, Chen W, Mi D, Murphy WJ, Li G. The genome of the black-footed cat: Revealing a rich natural history and urgent conservation priorities for small felids. Proc Natl Acad Sci U S A 2024; 121:e2310763120. [PMID: 38165928 PMCID: PMC10786289 DOI: 10.1073/pnas.2310763120] [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: 06/27/2023] [Accepted: 10/31/2023] [Indexed: 01/04/2024] Open
Abstract
Habitat degradation and loss of genetic diversity are common threats faced by almost all of today's wild cats. Big cats, such as tigers and lions, are of great concern and have received considerable conservation attention through policies and international actions. However, knowledge of and conservation actions for small wild cats are lagging considerably behind. The black-footed cat, Felis nigripes, one of the smallest felid species, is experiencing increasing threats with a rapid reduction in population size. However, there is a lack of genetic information to assist in developing effective conservation actions. A de novo assembly of a high-quality chromosome-level reference genome of the black-footed cat was made, and comparative genomics and population genomics analyses were carried out. These analyses revealed that the most significant genetic changes in the evolution of the black-footed cat are the rapid evolution of sensory and metabolic-related genes, reflecting genetic adaptations to its characteristic nocturnal hunting and a high metabolic rate. Genomes of the black-footed cat exhibit a high level of inbreeding, especially for signals of recent inbreeding events, which suggest that they may have experienced severe genetic isolation caused by habitat fragmentation. More importantly, inbreeding associated with two deleterious mutated genes may exacerbate the risk of amyloidosis, the dominant disease that causes mortality of about 70% of captive individuals. Our research provides comprehensive documentation of the evolutionary history of the black-footed cat and suggests that there is an urgent need to investigate genomic variations of small felids worldwide to support effective conservation actions.
Collapse
Affiliation(s)
- Jiaqing Yuan
- College of Life Sciences, Shaanxi Normal University, Xi’an710119, China
| | - Andrew C. Kitchener
- Department of Natural Sciences, National Museums Scotland, EdinburghEH1 1JF, United Kingdom
- School of Geosciences, University of Edinburgh, EdinburghEH9 3PX, United Kingdom
| | | | - Ting Sun
- College of Life Sciences, Shaanxi Normal University, Xi’an710119, China
| | - Qigao Jiangzuo
- Key Laboratory of Vertebrate Evolution and Human Origins of Chinese Academy of Sciences, Institute of Vertebrate Paleontology and Paleoanthropology, Chinese Academy of Sciences, Beijing100044, China
| | | | - Le Zhao
- College of Life Sciences, Shaanxi Normal University, Xi’an710119, China
- QinLing-Bashan Mountains Bioresources Comprehensive Development Collaborative Innovation Center, School of Bioscience and Engineering, Shaanxi University of Technology, Hanzhong723099, China
| | - Peng Yang
- College of Life Sciences, Shaanxi Normal University, Xi’an710119, China
| | - Guiqiang Wang
- College of Life Sciences, Shaanxi Normal University, Xi’an710119, China
| | - Chen Huang
- College of Life Sciences, Shaanxi Normal University, Xi’an710119, China
| | - Jinhong Wang
- College of Life Sciences, Shaanxi Normal University, Xi’an710119, China
| | - Wenhui Hou
- College of Life Sciences, Shaanxi Normal University, Xi’an710119, China
| | - Yang Liu
- College of Life Sciences, Shaanxi Normal University, Xi’an710119, China
| | - Wu Chen
- Guangzhou Zoo, Guangzhou Wildlife Research Center, Guangzhou510070, China
| | - Da Mi
- Xi’an Haorui Genomics Technology Co., Ltd., Xi’an710116, China
- School of Life Science and Technology, Xi’an Jiaotong University, Xi’an710049, China
| | - William J. Murphy
- Veterinary Integrative Biosciences, Texas A&M University, College Station, TX77843
| | - Gang Li
- College of Life Sciences, Shaanxi Normal University, Xi’an710119, China
- Guangzhou Zoo, Guangzhou Wildlife Research Center, Guangzhou510070, China
| |
Collapse
|
5
|
Xu Z, Wu J, Zhang Y, Qiao M, Zhou J, Feng Y, Li Z, Sun H, Lin R, Song Z, Zhao H, Li L, Chen N, Li Y, Oyelami FO, Peng X, Mei S. Genome-wide detection of selection signatures in Jianli pigs reveals novel cis-regulatory haplotype in EDNRB associated with two-end black coat color. BMC Genomics 2024; 25:23. [PMID: 38166718 PMCID: PMC10763394 DOI: 10.1186/s12864-023-09943-9] [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: 10/07/2023] [Accepted: 12/26/2023] [Indexed: 01/05/2024] Open
Abstract
BACKGROUND Jianli pig, a renowned indigenous breed in China, has the characteristics of a two-end black (TEB) coat color, excellent meat quality, strong adaptability and increased prolificacy. However, there is limited information available regarding the genetic diversity, population structure and genomic regions under selection of Jianli pig. On the other hand, the genetic mechanism of TEB coat color has remained largely unknown. RESULTS In this study, the whole genome resequencing of 30 Jianli pigs within a context of 153 individuals representing 13 diverse breeds was performed. The population structure analysis revealed that Jianli pigs have close genetic relationships with the Tongcheng pig breed, their geographical neighbors. Three methods (observed heterozygosity, expected heterozygosity, and runs of homozygosity) implied a relatively high level of genetic diversity and, a low inbreeding coefficient in Jianli compared with other pigs. We used Fst and XP-EHH to detect the selection signatures in Jianli pigs compared with Asian wild boar. A total of 451 candidate genes influencing meat quality (CREBBP, ADCY9, EEPD1 and HDAC9), reproduction (ESR1 and FANCA), and coat color (EDNRB, MITF and MC1R), were detected by gene annotation analysis. Finally, to fine-map the genomic region for the two-end black (TEB) coat color phenotype in Jianli pigs, we performed three signature selection methods between the TEB coat color and no-TEB coat color pig breeds. The current study, further confirmed that the EDNRB gene is a candidate gene for TEB color phenotype found in Chinese pigs, including Jinhua pigs, and the haplotype harboring 25 SNPs in the EDNRB gene may promote the formation of TEB coat color. Further ATAC-seq and luciferase reporter assays of these regions suggest that the 25-SNPs region was a strong candidate causative mutation that regulates the TEB coat color phenotype by altering enhancer function. CONCLUSION Our results advanced the understanding of the genetic mechanism behind artificial selection, and provided further resources for the protection and breeding improvement of Jianli pigs.
Collapse
Affiliation(s)
- Zhong Xu
- Hubei Key Laboratory of Animal Embryo and Molecular Breeding, Institute of Animal Husbandry and Veterinary, Hubei Provincial Academy of Agricultural Sciences, Wuhan, 430064, China
| | - Junjing Wu
- Hubei Key Laboratory of Animal Embryo and Molecular Breeding, Institute of Animal Husbandry and Veterinary, Hubei Provincial Academy of Agricultural Sciences, Wuhan, 430064, China
| | - Yu Zhang
- Hubei Key Laboratory of Animal Embryo and Molecular Breeding, Institute of Animal Husbandry and Veterinary, Hubei Provincial Academy of Agricultural Sciences, Wuhan, 430064, China
| | - Mu Qiao
- Hubei Key Laboratory of Animal Embryo and Molecular Breeding, Institute of Animal Husbandry and Veterinary, Hubei Provincial Academy of Agricultural Sciences, Wuhan, 430064, China
| | - Jiawei Zhou
- Hubei Key Laboratory of Animal Embryo and Molecular Breeding, Institute of Animal Husbandry and Veterinary, Hubei Provincial Academy of Agricultural Sciences, Wuhan, 430064, China
| | - Yue Feng
- Hubei Key Laboratory of Animal Embryo and Molecular Breeding, Institute of Animal Husbandry and Veterinary, Hubei Provincial Academy of Agricultural Sciences, Wuhan, 430064, China
| | - Zipeng Li
- Hubei Key Laboratory of Animal Embryo and Molecular Breeding, Institute of Animal Husbandry and Veterinary, Hubei Provincial Academy of Agricultural Sciences, Wuhan, 430064, China
| | - Hua Sun
- Hubei Key Laboratory of Animal Embryo and Molecular Breeding, Institute of Animal Husbandry and Veterinary, Hubei Provincial Academy of Agricultural Sciences, Wuhan, 430064, China
| | - Ruiyi Lin
- (College of Animal Sciences, College of Bee Science), Fujian Agriculture and Forestry University, Fuzhou, 350002, China
| | - Zhongxu Song
- Hubei Key Laboratory of Animal Embryo and Molecular Breeding, Institute of Animal Husbandry and Veterinary, Hubei Provincial Academy of Agricultural Sciences, Wuhan, 430064, China
| | - Haizhong Zhao
- Hubei Key Laboratory of Animal Embryo and Molecular Breeding, Institute of Animal Husbandry and Veterinary, Hubei Provincial Academy of Agricultural Sciences, Wuhan, 430064, China
| | - Lianghua Li
- Hubei Key Laboratory of Animal Embryo and Molecular Breeding, Institute of Animal Husbandry and Veterinary, Hubei Provincial Academy of Agricultural Sciences, Wuhan, 430064, China
| | - Nanqi Chen
- Hubei Key Laboratory of Animal Embryo and Molecular Breeding, Institute of Animal Husbandry and Veterinary, Hubei Provincial Academy of Agricultural Sciences, Wuhan, 430064, China
| | - Yujie Li
- Hubei Key Laboratory of Animal Embryo and Molecular Breeding, Institute of Animal Husbandry and Veterinary, Hubei Provincial Academy of Agricultural Sciences, Wuhan, 430064, China
| | | | - Xianwen Peng
- Hubei Key Laboratory of Animal Embryo and Molecular Breeding, Institute of Animal Husbandry and Veterinary, Hubei Provincial Academy of Agricultural Sciences, Wuhan, 430064, China.
- Hubei Hongshan Laboratory, Wuhan, 430064, China.
| | - Shuqi Mei
- Hubei Key Laboratory of Animal Embryo and Molecular Breeding, Institute of Animal Husbandry and Veterinary, Hubei Provincial Academy of Agricultural Sciences, Wuhan, 430064, China.
- Hubei Hongshan Laboratory, Wuhan, 430064, China.
| |
Collapse
|
6
|
Saif-Ur-Rehman M, Hassan FU, Reecy J, Deng T. Whole-genome SNP markers reveal runs of homozygosity in indigenous cattle breeds of Pakistan. Anim Biotechnol 2023; 34:1384-1396. [PMID: 35044288 DOI: 10.1080/10495398.2022.2026369] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/01/2022]
Abstract
The runs of homozygosity (ROH) were identified in 14 Pakistani cattle breeds (n = 105) by genotyping with the Illumina 50 K SNP BeadChip. These breeds were categorized into Dairy, Dual, and Draft breeds based on their utility and production performance. We identified a total of 10,936 ROHs which mainly consisted of a high number of shorter segments (1-4 Mb). Dairy group exhibited the highest level of inbreeding (FROH: 0.078 ± 0.028) while the lowest (FROH: 0.002 ± 0.008) was observed in Dual group. In 48 genomic regions identified with a high frequency of ROH, 207 genes were detected in the three breed groups. A substantially higher number of ROH islands detected in dairy breeds indicated the impact of the positive selection pressure over the years. Important candidate genes and QTL were detected in the ROH islands associated with economic traits like milk production, reproduction, meat, carcass, and health traits in dairy cattle.
Collapse
Affiliation(s)
| | - Faiz-Ul Hassan
- Institute of Animal and Dairy Sciences, University of Agriculture, Faisalabad, Pakistan
| | - James Reecy
- Department of Animal Science, Iowa State University, Ames, IA, USA
| | - Tingxian Deng
- Guangxi Provincial Key Laboratory of Buffalo Genetics, Breeding and Reproduction Technology, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, China
| |
Collapse
|
7
|
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
|
8
|
Pacheco HA, Rossoni A, Cecchinato A, Peñagaricano F. Identification of runs of homozygosity associated with male fertility in Italian Brown Swiss cattle. Front Genet 2023; 14:1227310. [PMID: 37485336 PMCID: PMC10356982 DOI: 10.3389/fgene.2023.1227310] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Accepted: 06/23/2023] [Indexed: 07/25/2023] Open
Abstract
Intensive selection for improved productivity has been accompanied by an increase in inbreeding rates and a reduction in genetic diversity. The increase in inbreeding tends to impact performance, especially fitness-related traits such as male fertility. Inbreeding can be monitored using runs of homozygosity (ROH), defined as contiguous lengths of homozygous genotypes observed in an individual's chromosome. The goal of this study was to evaluate the presence of ROH in Italian Brown Swiss cattle and assess its association with bull fertility. First, we evaluated the association between ROH and male fertility using 1,102 Italian Brown Swiss bulls with sire conception rate records and 572 K SNPs spanning the entire genome. Second, we split the entire population into 100 high-fertility and 100 low-fertility bulls to investigate the potential enrichment of ROH segments in the low-fertility group. Finally, we mapped the significant ROH regions to the bovine genome to identify candidate genes associated with sperm biology and male fertility. Notably, there was a negative association between bull fertility and the amount of homozygosity. Four different ROH regions located in chromosomes 6, 10, 11, and 24 were significantly overrepresented in low-fertility bulls (Fisher's exact test, p-value <0.01). Remarkably, these four genomic regions harbor many genes such as WDR19, RPL9, LIAS, UBE2K, DPF3, 5S-rRNA, 7SK, U6, and WDR7 that are related to sperm biology and male fertility. Overall, our findings suggest that inbreeding and increased homozygosity have a negative impact on male fertility in Italian Brown Swiss cattle. The quantification of ROH can contribute to minimizing the inbreeding rate and avoid its negative effect on fitness-related traits, such as male fertility.
Collapse
Affiliation(s)
- Hendyel A. Pacheco
- Department of Animal and Dairy Sciences, University of Wisconsin-Madison, Madison, WI, United States
| | | | - Alessio Cecchinato
- Department of Agronomy, Food, Natural Resources, Animals and Environment, University of Padova, Padua, Italy
| | - Francisco Peñagaricano
- Department of Animal and Dairy Sciences, University of Wisconsin-Madison, Madison, WI, United States
| |
Collapse
|
9
|
Chen C, Zhu B, Tang X, Chen B, Liu M, Gao N, Li S, Gu J. Genome-Wide Assessment of Runs of Homozygosity by Whole-Genome Sequencing in Diverse Horse Breeds Worldwide. Genes (Basel) 2023; 14:1211. [PMID: 37372391 DOI: 10.3390/genes14061211] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2023] [Revised: 05/29/2023] [Accepted: 05/30/2023] [Indexed: 06/29/2023] Open
Abstract
In the genomes of diploid organisms, runs of homozygosity (ROH), consecutive segments of homozygosity, are extended. ROH can be applied to evaluate the inbreeding situation of individuals without pedigree data and to detect selective signatures via ROH islands. We sequenced and analyzed data derived from the whole-genome sequencing of 97 horses, investigated the distribution of genome-wide ROH patterns, and calculated ROH-based inbreeding coefficients for 16 representative horse varieties from around the world. Our findings indicated that both ancient and recent inbreeding occurrences had varying degrees of impact on various horse breeds. However, recent inbreeding events were uncommon, particularly among indigenous horse breeds. Consequently, the ROH-based genomic inbreeding coefficient could aid in monitoring the level of inbreeding. Using the Thoroughbred population as a case study, we discovered 24 ROH islands containing 72 candidate genes associated with artificial selection traits. We found that the candidate genes in Thoroughbreds were involved in neurotransmission (CHRNA6, PRKN, and GRM1), muscle development (ADAMTS15 and QKI), positive regulation of heart rate and heart contraction (HEY2 and TRDN), regulation of insulin secretion (CACNA1S, KCNMB2, and KCNMB3), and spermatogenesis (JAM3, PACRG, and SPATA6L). Our findings provide insight into horse breed characteristics and future breeding strategies.
Collapse
Affiliation(s)
- Chujie Chen
- Hunan Provincial Key Laboratory for Genetic Improvement of Domestic Animal, College of Animal Science and Technology, Hunan Agricultural University, Changsha 410128, China
| | - Bo Zhu
- Novogene Bioinformatics Institute, Beijing 100015, China
| | - Xiangwei Tang
- Hunan Provincial Key Laboratory for Genetic Improvement of Domestic Animal, College of Animal Science and Technology, Hunan Agricultural University, Changsha 410128, China
| | - Bin Chen
- Hunan Provincial Key Laboratory for Genetic Improvement of Domestic Animal, College of Animal Science and Technology, Hunan Agricultural University, Changsha 410128, China
| | - Mei Liu
- Hunan Provincial Key Laboratory for Genetic Improvement of Domestic Animal, College of Animal Science and Technology, Hunan Agricultural University, Changsha 410128, China
| | - Ning Gao
- Hunan Provincial Key Laboratory for Genetic Improvement of Domestic Animal, College of Animal Science and Technology, Hunan Agricultural University, Changsha 410128, China
| | - Sheng Li
- Maxun Biotechnology Institute, Changsha 410024, China
| | - Jingjing Gu
- Hunan Provincial Key Laboratory for Genetic Improvement of Domestic Animal, College of Animal Science and Technology, Hunan Agricultural University, Changsha 410128, China
| |
Collapse
|
10
|
Martinez R, Bejarano D, Ramírez J, Ocampo R, Polanco N, Perez JE, Onofre HG, Rocha JF. Genomic variability and population structure of six Colombian cattle breeds. Trop Anim Health Prod 2023; 55:185. [PMID: 37130925 PMCID: PMC10154261 DOI: 10.1007/s11250-023-03574-8] [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: 09/03/2022] [Accepted: 03/29/2023] [Indexed: 05/04/2023]
Abstract
Analyses of the genetic diversity of indigenous cattle are essential for implementing conservation programs, promoting their sustainable use and maintaining productive advantages offered by these breeds in local conditions. The aim of this study was to investigate the genetic diversity and population structure of six Colombian cattle breeds: Blanco Orejinegro (BON), Costeño con Cuernos (CCC), Romosinuano (ROM), Sanmartinero (SAM), Casanareño (CAS), and Hartón del Valle (HDV). Two additional breed groups were included for comparison: Zebu (CEB) and a crossbreed of Colombian cattle breeds × Zebu. Genetic diversity within breeds was analyzed using expected heterozygosity (He), inbreeding coefficient (f), and runs of homozygosity (ROH). Population structure was assessed using model-based clustering (ADMIXTURE) and principal components analysis (PCA). Zebu cattle showed the lowest genetic diversity (He = 0.240). Breeds with the highest genetic diversity level were HDV and BON (He = 0.350 and 0.340, respectively). Inbreeding was lower for Colombian cattle breeds ranging between 0.005 and 0.045. Overall, the largest average genetic distance was found among the group of Colombian cattle breeds and Zebu, while the smallest was found between ROM and CCC. Model-based clustering revealed some level of admixture among HDV and CAS cattle which is consistent with their recent history. The results of the present study provide a useful insight on the genetic structure of Colombian cattle breeds.
Collapse
Affiliation(s)
- Rodrigo Martinez
- Corporación Colombiana de Investigación Agropecuaria - Agrosavia, Mosquera, Cundinamarca, Colombia
| | - Diego Bejarano
- Corporación Colombiana de Investigación Agropecuaria - Agrosavia, Mosquera, Cundinamarca, Colombia
| | - Julián Ramírez
- Corporación Colombiana de Investigación Agropecuaria - Agrosavia, Mosquera, Cundinamarca, Colombia
| | - Ricardo Ocampo
- Corporación Colombiana de Investigación Agropecuaria - Agrosavia, Mosquera, Cundinamarca, Colombia.
| | - Nelson Polanco
- Corporación Colombiana de Investigación Agropecuaria - Agrosavia, Mosquera, Cundinamarca, Colombia
| | - Juan Esteban Perez
- Corporación Colombiana de Investigación Agropecuaria - Agrosavia, Mosquera, Cundinamarca, Colombia
| | - Héctor Guillermo Onofre
- Corporación Colombiana de Investigación Agropecuaria - Agrosavia, Mosquera, Cundinamarca, Colombia
| | - Juan Felipe Rocha
- Corporación Colombiana de Investigación Agropecuaria - Agrosavia, Mosquera, Cundinamarca, Colombia
| |
Collapse
|
11
|
Hewett AM, Stoffel MA, Peters L, Johnston SE, Pemberton JM. Selection, recombination and population history effects on runs of homozygosity (ROH) in wild red deer (Cervus elaphus). Heredity (Edinb) 2023; 130:242-250. [PMID: 36801920 PMCID: PMC10076382 DOI: 10.1038/s41437-023-00602-z] [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: 05/11/2022] [Revised: 02/07/2023] [Accepted: 02/07/2023] [Indexed: 02/19/2023] Open
Abstract
The distribution of runs of homozygosity (ROH) may be shaped by a number of interacting processes such as selection, recombination and population history, but little is known about the importance of these mechanisms in shaping ROH in wild populations. We combined an empirical dataset of >3000 red deer genotyped at >35,000 genome-wide autosomal SNPs and evolutionary simulations to investigate the influence of each of these factors on ROH. We assessed ROH in a focal and comparison population to investigate the effect of population history. We investigated the role of recombination using both a physical map and a genetic linkage map to search for ROH. We found differences in ROH distribution between both populations and map types indicating that population history and local recombination rate have an effect on ROH. Finally, we ran forward genetic simulations with varying population histories, recombination rates and levels of selection, allowing us to further interpret our empirical data. These simulations showed that population history has a greater effect on ROH distribution than either recombination or selection. We further show that selection can cause genomic regions where ROH is common, only when the effective population size (Ne) is large or selection is particularly strong. In populations having undergone a population bottleneck, genetic drift can outweigh the effect of selection. Overall, we conclude that in this population, genetic drift resulting from a historical population bottleneck is most likely to have resulted in the observed ROH distribution, with selection possibly playing a minor role.
Collapse
Affiliation(s)
- Anna M Hewett
- Institute of Ecology and Evolution, School of Biological Sciences, University of Edinburgh, Edinburgh, EH9 3FL, UK.
| | - Martin A Stoffel
- Institute of Ecology and Evolution, School of Biological Sciences, University of Edinburgh, Edinburgh, EH9 3FL, UK
| | - Lucy Peters
- Institute of Ecology and Evolution, School of Biological Sciences, University of Edinburgh, Edinburgh, EH9 3FL, UK
- GenPhySE, Université de Toulouse, INRAE, ENVT, 31326, Castanet-Tolosan, France
| | - Susan E Johnston
- Institute of Ecology and Evolution, School of Biological Sciences, University of Edinburgh, Edinburgh, EH9 3FL, UK
| | - Josephine M Pemberton
- Institute of Ecology and Evolution, School of Biological Sciences, University of Edinburgh, Edinburgh, EH9 3FL, UK
| |
Collapse
|
12
|
Rocha RDFB, Garcia AO, Otto PI, da Silva MVB, Martins MF, Machado MA, Panetto JCDC, Guimarães SEF. Runs of homozygosity and signatures of selection for number of oocytes and embryos in the Gir Indicine cattle. Mamm Genome 2023:10.1007/s00335-023-09989-w. [PMID: 37000236 DOI: 10.1007/s00335-023-09989-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2022] [Accepted: 03/11/2023] [Indexed: 04/01/2023]
Abstract
Runs of homozygosity (ROH) and signatures of selection are the results of selection processes in livestock species that have been shown to affect several traits in cattle. The aim of the current work was to verify the profile of ROH and inbreeding depression in the number of total (TO) and viable oocytes (VO) and the number of embryos (EMBR) in Gir Indicine cattle. In addition, we aim to identify signatures of selection, genes, and enriched regions between Gir subpopulations sorted by breeding value for these traits. The genotype file contained 2093 animals and 420,718 SNP markers. Breeding values used to sort Gir animals were previously obtained. ROH and signature of selection analyses were performed using PLINK software, followed by ROH-based (FROH) and pedigree-based inbreeding (Fped) and a search for genes and their functions. An average of 50 ± 8.59 ROHs were found per animal. ROHs were separated into classes according to size, ranging from 1 to 2 Mb (ROH1-2Mb: 58.17%), representing ancient inbreeding, ROH2-4Mb (22.74%), ROH4-8Mb (11.34%), ROH8-16Mb (5.51%), and ROH>16Mb (2.24%). Combining our results, we conclude that the increase in general FROH and Fped significantly decreases TO and VO; however, in different chromosomes traits can increase or decrease with FROH. In the analysis for signatures of selection, we identified 15 genes from 47 significant genomic regions, indicating differences in populations with high and low breeding value for the three traits.
Collapse
Affiliation(s)
| | | | - Pamela Itajara Otto
- Department of Animal Science, Universidade Federal de Santa Maria, Santa Maria, Rio Grande Do Sul, Brazil
| | | | | | | | | | | |
Collapse
|
13
|
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
|
14
|
Lavanchy E, Goudet J. Effect of reduced genomic representation on using runs of homozygosity for inbreeding characterization. Mol Ecol Resour 2023; 23:787-802. [PMID: 36626297 DOI: 10.1111/1755-0998.13755] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2022] [Revised: 12/22/2022] [Accepted: 01/05/2023] [Indexed: 01/11/2023]
Abstract
Genomic measures of inbreeding based on identical-by-descent (IBD) segments are increasingly used to measure inbreeding and mostly estimated on SNP arrays and whole-genome sequencing (WGS) data. However, some softwares recurrently used for their estimation assume that genomic positions which have not been genotyped are nonvariant. This might be true for WGS data, but not for reduced genomic representations and can lead to spurious IBD segments estimation. In this project, we simulated the outputs of WGS, two SNP arrays of different sizes and RAD-sequencing for three populations with different sizes and histories. We compare the results of IBD segments estimation with two softwares: runs of homozygosity (ROHs) estimated with PLINK and homozygous-by-descent (HBD) segments estimated with RZooRoH. We demonstrate that to obtain meaningful estimates of inbreeding, RZooRoH requires a SNPs density 11 times smaller compared to PLINK: ranks of inbreeding coefficients were conserved among individuals above 22 SNPs/Mb for PLINK and 2 SNPs/Mb for RZooRoH. We also show that in populations with simple demographic histories, distribution of ROHs and HBD segments are correctly estimated with both SNP arrays and WGS. PLINK correctly estimated distribution of ROHs with SNP densities above 22 SNPs/Mb, while RZooRoH correctly estimated distribution of HBD segments with SNPs densities above 11 SNPs/Mb. However, in a population with a more complex demographic history, RZooRoH resulted in better distribution of IBD segments estimation compared to PLINK even with WGS data. Consequently, we advise researchers to use either methods relying on excess homozygosity averaged across SNPs or model-based HBD segments calling methods for inbreeding estimations.
Collapse
Affiliation(s)
- Eléonore Lavanchy
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland.,Swiss Institute of Bioinformatics, University of Lausanne, Lausanne, Switzerland
| | - Jérôme Goudet
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland.,Swiss Institute of Bioinformatics, University of Lausanne, Lausanne, Switzerland
| |
Collapse
|
15
|
Ablondi M, Summer A, Stocco G, Finocchiaro R, van Kaam JT, Cassandro M, Dadousis C, Sabbioni A, Cipolat-Gotet C. The role of inbreeding depression on productive performance in the Italian Holstein breed. J Anim Sci 2023; 101:skad382. [PMID: 37983004 PMCID: PMC10693289 DOI: 10.1093/jas/skad382] [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: 05/31/2023] [Accepted: 11/18/2023] [Indexed: 11/21/2023] Open
Abstract
Inbreeding depression has become an urgent issue in cosmopolitan breeds where the massive genetic progress achieved in the latest generations is counterbalanced by a dramatic loss of genetic diversity causing increased health issues. Thus, the aim of this study was to estimate inbreeding depression on productive traits in Holstein dairy cattle. More precisely, we aimed to i) determine the level of inbreeding in 27,735 Italian Holstein dairy cows using pedigree and genotype data, ii) quantify the effect of inbreeding on 305-d in milk yield (MY; kg), fat yield (FY; kg), and protein yield (PY; kg) based on different statistical approaches, iii) determine if recent inbreeding has a more harmful impact than ancestral ones, and iv) quantify chromosomal homozygosity effect on productive traits. Quality control was performed on the autosomal chromosomes resulting in a final dataset of 84,443 single nucleotide polymorphisms. Four statistical models were used to evaluate the presence of inbreeding depression, which included linear regression analysis and division of FPED and FROH into percentile classes. Moreover, FROH was partitioned into i) length classes to assess the role of recent and ancestral inbreeding and ii) chromosome-specific contributions (FROH-CHR). Results evidenced that inbreeding negatively impacted the productive performance of Italian Holstein Friesian cows. However, differences between the estimated FPED and FROH coefficients resulted in different estimates of inbreeding depression. For instance, a 1% increase in FPED and FROH was associated with a decrease in MY of about 44 and 61 kg (P < 0.01). Further, when considering the extreme inbreeding percentile classes moving from the 5th lowest to the 95th highest, there was a reduction of -263 kg and -561 kg per lactation for FPED and FROH. Increased inbreeding, estimated by FPED and FROH, had also a negative effect on PY and FY, either fit as a regressor or percentile classes. When evaluating the impact of inbreeding based on runs of homozygosity (ROH) length classes, longer ROH (over 8 Mb) had a negative effect in all traits, indicating that recent inbreeding might be more harmful than the ancestral one. Finally, results within chromosome homozygosity highlighted specific chromosomes with a more deleterious effect on productive traits.
Collapse
Affiliation(s)
- Michela Ablondi
- Department of Veterinary Science, University of Parma, 43126 Parma, Italy
| | - Andrea Summer
- Department of Veterinary Science, University of Parma, 43126 Parma, Italy
| | - Giorgia Stocco
- Department of Veterinary Science, University of Parma, 43126 Parma, Italy
| | - Raffaella Finocchiaro
- Associazione Nazionale Allevatori della Razza Frisona Bruna e Jersey Italiana (ANAFIBJ), 26100 Cremona, Italy
| | - Jan-Thijs van Kaam
- Associazione Nazionale Allevatori della Razza Frisona Bruna e Jersey Italiana (ANAFIBJ), 26100 Cremona, Italy
| | - Martino Cassandro
- Associazione Nazionale Allevatori della Razza Frisona Bruna e Jersey Italiana (ANAFIBJ), 26100 Cremona, Italy
- Department of Agronomy, Food, Natural resources, Animals and Environment, University of Padova, 35020 Legnaro, Italy
| | - Christos Dadousis
- Department of Veterinary Science, University of Parma, 43126 Parma, Italy
| | - Alberto Sabbioni
- Department of Veterinary Science, University of Parma, 43126 Parma, Italy
| | | |
Collapse
|
16
|
Chessari G, Criscione A, Tolone M, Bordonaro S, Rizzuto I, Riggio S, Macaluso V, Moscarelli A, Portolano B, Sardina MT, Mastrangelo S. High-density SNP markers elucidate the genetic divergence and population structure of Noticiana sheep breed in the Mediterranean context. Front Vet Sci 2023; 10:1127354. [PMID: 37205231 PMCID: PMC10185747 DOI: 10.3389/fvets.2023.1127354] [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: 12/19/2022] [Accepted: 04/13/2023] [Indexed: 05/21/2023] Open
Abstract
Among livestock species, sheep have played an early major role in the Mediterranean area. Italy has a long history of sheep breeding and, despite a dramatic contraction in numbers, still raise several local populations that may represent a unique source of genetic diversity. The Noticiana is a breed of the south-eastern part of Sicily appreciated both for its dairy products and for its resistance to harsh environment. In this study, the high-density Illumina Ovine SNP600K BeadChip array was used for the first genome-wide characterization of 48 individuals of Noticiana sheep to investigate its diversity, the genome structure and the relationship within the context of worldwide and Italian breeds. Moreover, the runs of homozygosity (ROH) pattern and the pairwise FST-outliers were examined. Noticiana reported moderate levels of genetic diversity. The high percentage of short and medium length ROH segments (93% under 4 Mb) is indicative of a within breed relatedness dating back to ancient times, despite the absence of management for the mating plans and the reduced population size. In the worldwide context, the Southern Italian, Spanish and Albanian breeds overlapped in a macro cluster which also included the Noticiana sheep. The results highlighted ancestral genetic components of Noticiana shared with Comisana breed, and showed the clear separation from the other Italian sheep. This is likely the consequence of the combined effects of genetic drift, small population size and reproductive isolation. ROH islands and FST-outliers approaches in Noticiana identified genes and QTLs involved in milk and meat production, as well as related to the local adaptation, and therefore are consistent with the phenotypic traits of the studied breed. Although a wider sampling could be useful to deepen the genomic survey on Noticiana, these results represent a crucial starting point for the characterization of an important local genetic resource, with a view of supporting the local economy and preserving the biodiversity of the sheep species.
Collapse
Affiliation(s)
- Giorgio Chessari
- Dipartimento Agricoltura, Alimentazione e Ambiente, University of Catania, Catania, Italy
| | - Andrea Criscione
- Dipartimento Agricoltura, Alimentazione e Ambiente, University of Catania, Catania, Italy
| | - Marco Tolone
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Salvatore Bordonaro
- Dipartimento Agricoltura, Alimentazione e Ambiente, University of Catania, Catania, Italy
| | - Ilaria Rizzuto
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Silvia Riggio
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Vito Macaluso
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Angelo Moscarelli
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Baldassare Portolano
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Maria Teresa Sardina
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Salvatore Mastrangelo
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
- *Correspondence: Salvatore Mastrangelo,
| |
Collapse
|
17
|
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
|
18
|
Bosse M, van Loon S. Challenges in quantifying genome erosion for conservation. Front Genet 2022; 13:960958. [PMID: 36226192 PMCID: PMC9549127 DOI: 10.3389/fgene.2022.960958] [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: 06/03/2022] [Accepted: 08/09/2022] [Indexed: 11/18/2022] Open
Abstract
Massive defaunation and high extinction rates have become characteristic of the Anthropocene. Genetic effects of population decline can lead populations into an extinction vortex, where declining populations show lower genetic fitness, in turn leading to lower populations still. The lower genetic fitness in a declining population due to a shrinking gene pool is known as genetic erosion. Three different types of genetic erosion are highlighted in this review: overall homozygosity, genetic load and runs of homozygosity (ROH), which are indicative of inbreeding. The ability to quantify genetic erosion could be a very helpful tool for conservationists, as it can provide them with an objective, quantifiable measure to use in the assessment of species at risk of extinction. The link between conservation status and genetic erosion should become more apparent. Currently, no clear correlation can be observed between the current conservation status and genetic erosion. However, the high quantities of genetic erosion in wild populations, especially in those species dealing with habitat fragmentation and habitat decline, may be early signs of deteriorating populations. Whole genome sequencing data is the way forward to quantify genetic erosion. Extra screening steps for genetic load and hybridization can be included, since they could potentially have great impact on population fitness. This way, the information yielded from genetic sequence data can provide conservationists with an objective genetic method in the assessment of species at risk of extinction. However, the great complexity of genome erosion quantification asks for consensus and bridging science and its applications, which remains challenging.
Collapse
Affiliation(s)
- Mirte Bosse
- Amsterdam Institute for Life and Environment (A-LIFE), Section Ecology and Evolution, Vrije Universiteit Amsterdam, Amsterdam, Netherlands
- Animal Breeding and Genomics, Wageningen University and Research, Wageningen, Netherlands
- *Correspondence: Mirte Bosse,
| | - Sam van Loon
- Amsterdam Institute for Life and Environment (A-LIFE), Section Ecology and Evolution, Vrije Universiteit Amsterdam, Amsterdam, Netherlands
| |
Collapse
|
19
|
Lozada-Soto EA, Tiezzi F, Jiang J, Cole JB, VanRaden PM, Maltecca C. Genomic characterization of autozygosity and recent inbreeding trends in all major breeds of US dairy cattle. J Dairy Sci 2022; 105:8956-8971. [PMID: 36153159 DOI: 10.3168/jds.2022-22116] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Accepted: 06/20/2022] [Indexed: 11/19/2022]
Abstract
Maintaining a genetically diverse dairy cattle population is critical to preserving adaptability to future breeding goals and avoiding declines in fitness. This study characterized the genomic landscape of autozygosity and assessed trends in genetic diversity in 5 breeds of US dairy cattle. We analyzed a sizable genomic data set containing 4,173,679 pedigreed and genotyped animals of the Ayrshire, Brown Swiss, Guernsey, Holstein, and Jersey breeds. Runs of homozygosity (ROH) of 2 Mb or longer in length were identified in each animal. The within-breed means for number and the combined length of ROH were highest in Jerseys (62.66 ± 8.29 ROH and 426.24 ± 83.40 Mb, respectively; mean ± SD) and lowest in Ayrshires (37.24 ± 8.27 ROH and 265.05 ± 85.00 Mb, respectively). Short ROH were the most abundant, but moderate to large ROH made up the largest proportion of genome autozygosity in all breeds. In addition, we identified ROH islands in each breed. This revealed selection patterns for milk production, productive life, health, and reproduction in most breeds and evidence for parallel selective pressure for loci on chromosome 6 between Ayrshire and Brown Swiss and for loci on chromosome 20 between Holstein and Jersey. We calculated inbreeding coefficients using 3 different approaches, pedigree-based (FPED), marker-based using a genomic relationship matrix (FGRM), and segment-based using ROH (FROH). The average inbreeding coefficient ranged from 0.06 in Ayrshires and Brown Swiss to 0.08 in Jerseys and Holsteins using FPED, from 0.22 in Holsteins to 0.29 in Guernsey and Jerseys using FGRM, and from 0.11 in Ayrshires to 0.17 in Jerseys using FROH. In addition, the effective population size at past generations (5-100 generations ago), the yearly rate of inbreeding, and the effective population size in 3 recent periods (2000-2009, 2010-2014, and 2015-2018) were determined in each breed to ascertain current and historical trends of genetic diversity. We found a historical trend of decreasing effective population size in the last 100 generations in all breeds and breed differences in the effect of the recent implementation of genomic selection on inbreeding accumulation.
Collapse
Affiliation(s)
| | - Francesco Tiezzi
- Department of Agriculture, Food, Environment and Forestry (DAGRI), University of Florence, 50144 Florence, Italy
| | - Jicai Jiang
- Department of Animal Science, North Carolina State University, Raleigh 27607
| | | | - Paul M VanRaden
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Service, USDA, Beltsville, MD 20705
| | - Christian Maltecca
- Department of Animal Science, North Carolina State University, Raleigh 27607
| |
Collapse
|
20
|
Cortes-Hernández JG, Ruiz-López FJ, Vásquez-Peláez CG, García-Ruiz A. Runs of homocigosity and its association with productive traits in Mexican Holstein cattle. PLoS One 2022; 17:e0274743. [PMID: 36121861 PMCID: PMC9484644 DOI: 10.1371/journal.pone.0274743] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Accepted: 09/03/2022] [Indexed: 11/21/2022] Open
Abstract
The objective of this study was to describe the runs of homozygosity (ROH) detected in the Mexican Holstein population and to associate them with milk, fat and protein yields, and conformation final score. After imputation and genomic quality control, 4,227 genotyped animals with 100,806 SNPs markers each were used. ROH with a minimum length of 1 Mb and a minimum of 10 SNPs were included in the analysis. One heterozygous SNP marker and five missing genotypes per ROH were allowed. A total of 425,098 ROH were found in the studied population (71.83 ± 10.73 ROH per animal), with an average length and coverage of 4.80 ± 0.77 Mb, and 276.89 Mb, respectively. The average chromosome length covered by ROH was 10.40 ± 3.70 Mb. ROH between 1 and 2 Mb were the most frequent in the population (51.33%) while those between 14 and 16 Mb were the least frequent (1.20%). Long chromosomes showed a larger number of ROH. Chromosomes 10 and 20, had a greater percentage of their length covered by ROH because they presented a largest number of long ROH (>8 Mb). From the total ROH, 17 were detected in 1,847 animals and distributed among different chromosomes, and were associated with milk, fat and protein yield and percentage, and conformation final score. Of the ROH with effects on production traits, the majority were found with a length between 1 and 4 Mb. These results show evidence of genomic regions preserved by genetic selection and associated with the improvement of the productivity and functionality of dairy cattle.
Collapse
Affiliation(s)
- José G. Cortes-Hernández
- Programa de Maestría y Doctorado en Ciencias de la Producción y de la Salud Animal, Universidad Nacional Autónoma de México, Ciudad de México, México
| | - Felipe J. Ruiz-López
- Centro Nacional de Investigación Disciplinaria en Fisiología y Mejoramiento Animal, Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias, Ajuchitlán Colón Querétaro, México
| | - Carlos G. Vásquez-Peláez
- Departamento de Genética y Bioestadística de la Facultad de Medicina Veterinaria y Zootecnia, Universidad Nacional Autónoma de México, Ciudad de México, México
| | - Adriana García-Ruiz
- Centro Nacional de Investigación Disciplinaria en Fisiología y Mejoramiento Animal, Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias, Ajuchitlán Colón Querétaro, México
- * E-mail:
| |
Collapse
|
21
|
|
22
|
Yuan J, Zhou X, Xu G, Xu S, Liu B. Genetic diversity and population structure of Tongcheng pigs in China using whole-genome SNP chip. Front Genet 2022; 13:910521. [PMID: 36092902 PMCID: PMC9455598 DOI: 10.3389/fgene.2022.910521] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2022] [Accepted: 07/29/2022] [Indexed: 11/13/2022] Open
Abstract
Tongcheng (TC) pigs, distinguished by their superior meat quality, are a Chinese indigenous pig breed. Recently, the genetic resources of TC pigs are under tremendous threat due to the introduction of cosmopolitan pig breeds and African swine fever disease. To promote their management and conservation, the present study assessed genetic diversity and population structure of TC pigs using single nucleotide polymorphism (SNP) markers. A total of 26, 999 SNPs were screened from 51, 315 SNPs in 68 TC pigs. The multi-dimensional scaling (MDS) analysis and neighbor-joining tree revealed that all 68 pigs were from a purebred population. The effective population size decreased over time, and it was 96 prior to generation 20. Both linkage disequilibrium (LD) and neutrality test indicated a low selection of TC pigs with average LD value of 0.15 ± 0.23. Genetic diversity results exhibited a minor allele frequency (MAF) of 0.23, observed heterozygosity (HO) of 0.32, expected heterozygosity (He) of 0.31, and nucleotide diversity (Pi) of 0.31. All these parameters indicated a remarkably high genetic diversity of TC pigs. Additionally, 184 runs of homozygosity (ROH) segments were detected from the whole genome of TC pigs with an average ROH length of 23.71Mb, ranging from 11.26Mb to 69.02 Mb. The highest ROH coverage was found on chromosome 1 (10.12%), while the lowest was on chromosome 18 (1.49%). The average inbreeding coefficients based on ROH (FROH) was 0.04%. Fourteen ROH islands containing 240 genes were detected on 9 different autosomes. Some of these 240 genes were overlapped with the genes related to biological processes such as immune function, reproduction, muscular development, and fat deposition, including FFAR2, FFAR4, MAPK8, NPY5R, KISS1, and these genes might be associated with such traits as meat quality and disease resistance in TC pigs. Taken together, population structure and genetic diversity results suggested that the TC pig represented a valuable genetic resource. However, TC pig breed conservation program remains to be further optimized to ensure adequate genetic diversity and avoid inbreeding depression. Our findings provide theoretical basis for formulating management and conservation strategies for TC pigs.
Collapse
Affiliation(s)
- Jiao Yuan
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, China
| | - Xiang Zhou
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, China
- The Cooperative Innovation Center for Sustainable Pig Production, Wuhan, China
- The Engineering Technology Research Center of Local Pig Breed Improvement of Hubei Province, Wuhan, China
| | - Guoqiang Xu
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, China
| | - Sanping Xu
- Department of Agricultural and Rural Bureau, Xianning, China
| | - Bang Liu
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, China
- The Cooperative Innovation Center for Sustainable Pig Production, Wuhan, China
- The Engineering Technology Research Center of Local Pig Breed Improvement of Hubei Province, Wuhan, China
- *Correspondence: Bang Liu,
| |
Collapse
|
23
|
Genomic inbreeding and runs of homozygosity analysis of indigenous cattle populations in southern China. PLoS One 2022; 17:e0271718. [PMID: 36006904 PMCID: PMC9409551 DOI: 10.1371/journal.pone.0271718] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Accepted: 07/06/2022] [Indexed: 11/19/2022] Open
Abstract
Runs of homozygosity (ROH) are continuous homozygous segments from the common ancestor of parents. Evaluating ROH pattern can help to understand inbreeding level and genetic basis of important traits. In this study, three representative cattle populations including Leiqiong cattle (LQC), Lufeng cattle (LFC) and Hainan cattle (HNC) were genotyped using the Illumina BovineHD SNPs array (770K) to assess ROH pattern at genome wide level. Totally, we identified 26,537 ROH with an average of 153 ROH per individual. The sizes of ROH ranged from 0.5 to 53.26Mb, and the average length was 1.03Mb. The average of FROH ranged from 0.10 (LQC) to 0.15 (HNC). Moreover, we identified 34 ROH islands (with frequency > 0.5) across genome. Based on these regions, we observed several breed-specific candidate genes related to adaptive traits. Several common genes related to immunity (TMEM173, MZB1 and SIL1), and heat stress (DNAJC18) were identified in all three populations. Three genes related to immunity (UGP2), development (PURA) and reproduction (VPS54) were detected in both HNC and LQC. Notably, we identified several breed-specific genes related to sperm development (BRDT and SPAG6) and heat stress (TAF7) in HNC, and immunity (CDC23 and NME5) and development (WNT87) in LFC. Our findings provided valuable insights into understanding the genomic homozygosity pattern and promoting the conservation of genetic resources of Chinese indigenous cattle.
Collapse
|
24
|
Chokoe TC, Hadebe K, Muchadeyi FC, Nephawe KA, Dzomba EF, Mphahlele TD, Matelele TC, Mtileni BJ. Conservation status and historical relatedness of South African communal indigenous goat populations using a genome-wide single-nucleotide polymorphism marker. Front Genet 2022; 13:909472. [PMID: 36017496 PMCID: PMC9395594 DOI: 10.3389/fgene.2022.909472] [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/31/2022] [Accepted: 06/27/2022] [Indexed: 11/18/2022] Open
Abstract
Indigenous goats form the majority of populations in smallholder, low input, low output production systems and are considered an important genetic resource due to their adaptability to different production environments and support of communal farming. Effective population size (Ne), inbreeding levels, and the runs of homozygosity (ROHs) are effective tools for exploring the genetic diversity and understanding the demographic history in efforts to support breeding strategies to use and conserve genetic resources. Across populations, the current Ne of Gauteng was the lowest at 371 animals, while the historical Ne across populations suggests that the ancestor Ne has decreased by 53.86%, 44.58%, 42.16%, and 41.16% in Free State (FS), North West (NW), Limpopo (LP), and Gauteng (GP), respectively, over the last 971 generations. Genomic inbreeding levels related to ancient kinship (FROH > 5 Mb) were highest in FS (0.08 ± 0.09) and lowest in the Eastern Cape (EC) (0.02 ± 0.02). A total of 871 ROH island regions which include important environmental adaptation and hermo-tolerance genes such as IL10RB, IL23A, FGF9, IGF1, EGR1, MTOR, and MAPK3 were identified (occurring in over 20% of the samples) in FS (n = 37), GP (n = 42), and NW (n = 2) populations only. The mean length of ROH across populations was 7.76 Mb and ranged from 1.61 Mb in KwaZulu-Natal (KZN) to 98.05 Mb (GP and NW). The distribution of ROH according to their size showed that the majority (n = 1949) of the detected ROH were > 5 Mb in length compared to the other categories. Assuming two hypothetical ancestral populations, the populations from KZN and LP are revealed, supporting PC 1. The genomes of KZN and LP share a common origin but have substantial admixture from the EC and NW populations. The findings revealed that the occurrence of high Ne and autozygosity varied largely across breeds in communal indigenous goat populations at recent and ancient events when a genome-wide single-nucleotide polymorphism (SNP) marker was used. The use of Illumina goat SNP50K BeadChip shows that there was a migration route of communal indigenous goat populations from the northern part (LP) of South Africa to the eastern areas of the KZN that confirmed their historical relatedness and coincides with the migration periods of the Bantu nation.
Collapse
Affiliation(s)
- T. C. Chokoe
- Farm Animal Genetic Resources, Department of Agriculture, Land Reform and Rural Development, Pretoria, South Africa
- School of Agriculture & Environmental Sciences, University of Limpopo, Polokwane, South Africa
- *Correspondence: T. C. Chokoe,
| | - K. Hadebe
- Biotechnology Platform, Agricultural Research Council, Pretoria, South Africa
| | - F. C. Muchadeyi
- Biotechnology Platform, Agricultural Research Council, Pretoria, South Africa
| | - K. A. Nephawe
- Department of Animal Sciences, Tshwane University of Technology, Pretoria, South Africa
| | - E. F. Dzomba
- Discipline of Genetics, School of Life Sciences, University of Kwazulu-Natal, Scottsville, South African
| | - T. D. Mphahlele
- Farm Animal Genetic Resources, Department of Agriculture, Land Reform and Rural Development, Pretoria, South Africa
| | - T. C. Matelele
- Farm Animal Genetic Resources, Department of Agriculture, Land Reform and Rural Development, Pretoria, South Africa
| | - B. J. Mtileni
- Department of Animal Sciences, Tshwane University of Technology, Pretoria, South Africa
| |
Collapse
|
25
|
Berghöfer J, Khaveh N, Mundlos S, Metzger J. Simultaneous testing of rule- and model-based approaches for runs of homozygosity detection opens up a window into genomic footprints of selection in pigs. BMC Genomics 2022; 23:564. [PMID: 35933356 PMCID: PMC9357325 DOI: 10.1186/s12864-022-08801-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2022] [Accepted: 07/29/2022] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Past selection events left footprints in the genome of domestic animals, which can be traced back by stretches of homozygous genotypes, designated as runs of homozygosity (ROHs). The analysis of common ROH regions within groups or populations displaying potential signatures of selection requires high-quality SNP data as well as carefully adjusted ROH-defining parameters. In this study, we used a simultaneous testing of rule- and model-based approaches to perform strategic ROH calling in genomic data from different pig populations to detect genomic regions under selection for specific phenotypes. RESULTS Our ROH analysis using a rule-based approach offered by PLINK, as well as a model-based approach run by RZooRoH demonstrated a high efficiency of both methods. It underlined the importance of providing a high-quality SNP set as input as well as adjusting parameters based on dataset and population for ROH calling. Particularly, ROHs ≤ 20 kb were called in a high frequency by both tools, but to some extent covered different gene sets in subsequent analysis of ROH regions common for investigated pig groups. Phenotype associated ROH analysis resulted in regions under potential selection characterizing heritage pig breeds, known to harbour a long-established breeding history. In particular, the selection focus on fitness-related traits was underlined by various ROHs harbouring disease resistance or tolerance-associated genes. Moreover, we identified potential selection signatures associated with ear morphology, which confirmed known candidate genes as well as uncovered a missense mutation in the ABCA6 gene potentially supporting ear cartilage formation. CONCLUSIONS The results of this study highlight the strengths and unique features of rule- and model-based approaches as well as demonstrate their potential for ROH analysis in animal populations. We provide a workflow for ROH detection, evaluating the major steps from filtering for high-quality SNP sets to intersecting ROH regions. Formula-based estimations defining ROHs for rule-based method show its limits, particularly for efficient detection of smaller ROHs. Moreover, we emphasize the role of ROH detection for the identification of potential footprints of selection in pigs, displaying their breed-specific characteristics or favourable phenotypes.
Collapse
Affiliation(s)
- Jan Berghöfer
- Research Group Veterinary Functional Genomics, Max Planck Institute for Molecular Genetics, Berlin, Germany.,Department of Biology, Chemistry and Pharmacy, Institute of Chemistry and Biochemistry, Freie Universität Berlin, Berlin, Germany
| | - Nadia Khaveh
- Research Group Veterinary Functional Genomics, Max Planck Institute for Molecular Genetics, Berlin, Germany.,Institute of Animal Breeding and Genetics, University of Veterinary Medicine Hannover, Hannover, Germany
| | - Stefan Mundlos
- Research Group Development & Disease, Max Planck Institute for Molecular Genetics, Berlin, Germany.,Institute of Medical Genetics and Human Genetics, Charité-Universitätsmedizin Berlin, Berlin, Germany.,Charité-Universitätsmedizin Berlin, BCRT, Berlin Institute of Health Center for Regenerative Therapies, Berlin, Germany
| | - Julia Metzger
- Research Group Veterinary Functional Genomics, Max Planck Institute for Molecular Genetics, Berlin, Germany. .,Institute of Animal Breeding and Genetics, University of Veterinary Medicine Hannover, Hannover, Germany.
| |
Collapse
|
26
|
Runs of Homozygosity and Quantitative Trait Locus/Association for Semen Parameters in Selected Chinese and South African Beef Cattle. Animals (Basel) 2022; 12:ani12121546. [PMID: 35739882 PMCID: PMC9219517 DOI: 10.3390/ani12121546] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2022] [Revised: 05/31/2022] [Accepted: 06/06/2022] [Indexed: 12/01/2022] Open
Abstract
In this study, runs of homozygosity (ROH) and quantitative trait locus/association (QTL) for semen parameters in selected Chinese and South African beef cattle breed were estimated. The computed results showed 7516 ROH were observed between classes 0−5 Mb with no ROH observed in classes >40 Mb. Distribution of ROH showed high level of genomic coverage for ANG, NGU, CSI, and BEL breeds. Approximately 13 genomic regions with QTL were controlling sperm motility, sperm concentration, semen volume, sperm count, sperm head abnormalities, sperm tail abnormalities, sperm integrity, and percentage of abnormal sperm traits. Nine candidate genes, CDF9, MARCH1, WDR19, SLOICI, ST7, DOP1B, CFAF9, INHBA, and ADAMTS1, were suggested to be associated with above mentioned QTL traits. The results for inbreeding coefficient showed moderate correlation between FROH vs FHOM at 0.603 and high correlation between FROH 0−5 Mb 0.929, and lowest correlation for 0−>40 Mb 0.400. This study suggested recent inbreeding in CSI, BEL, ANG, BON, SIM, and NGU breeds. Furthermore, it highlighted varied inbreeding levels and identified QTL for semen traits and genes of association. These results can assist in implementation of genetic improvement strategies for bulls and provide awareness and proper guidelines in developing breeding programs.
Collapse
|
27
|
Zhang Q, Zhang Q, Jensen J. Association Studies and Genomic Prediction for Genetic Improvements in Agriculture. FRONTIERS IN PLANT SCIENCE 2022; 13:904230. [PMID: 35720549 PMCID: PMC9201771 DOI: 10.3389/fpls.2022.904230] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/25/2022] [Accepted: 05/16/2022] [Indexed: 06/15/2023]
Abstract
To feed the fast growing global population with sufficient food using limited global resources, it is urgent to develop and utilize cutting-edge technologies and improve efficiency of agricultural production. In this review, we specifically introduce the concepts, theories, methods, applications and future implications of association studies and predicting unknown genetic value or future phenotypic events using genomics in the area of breeding in agriculture. Genome wide association studies can identify the quantitative genetic loci associated with phenotypes of importance in agriculture, while genomic prediction utilizes individual genetic value to rank selection candidates to improve the next generation of plants or animals. These technologies and methods have improved the efficiency of genetic improvement programs for agricultural production via elite animal breeds and plant varieties. With the development of new data acquisition technologies, there will be more and more data collected from high-through-put technologies to assist agricultural breeding. It will be crucial to extract useful information among these large amounts of data and to face this challenge, more efficient algorithms need to be developed and utilized for analyzing these data. Such development will require knowledge from multiple disciplines of research.
Collapse
Affiliation(s)
- Qianqian Zhang
- Institute of Biotechnology, Beijing Academy of Agricultural and Forestry Sciences, Beijing, China
| | - Qin Zhang
- College of Animal Science and Technology, Shandong Agricultural University, Taian, China
- College of Animal Science and Technology, China Agricultural University, BeijingChina
| | - Just Jensen
- Centre for Quantitative Genetics and Genomics, Aarhus University, Aarhus, Denmark
| |
Collapse
|
28
|
Jiang Y, Li X, Liu J, Zhang W, Zhou M, Wang J, Liu L, Su S, Zhao F, Chen H, Wang C. Genome-wide detection of genetic structure and runs of homozygosity analysis in Anhui indigenous and Western commercial pig breeds using PorcineSNP80k data. BMC Genomics 2022; 23:373. [PMID: 35581549 PMCID: PMC9115978 DOI: 10.1186/s12864-022-08583-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2021] [Accepted: 04/22/2022] [Indexed: 11/25/2022] Open
Abstract
Background Runs of homozygosity (ROH) are continuous homozygous regions typically located in the DNA sequence of diploid organisms. Identifications of ROH that lead to reduced performance can provide valuable insight into the genetic architecture of complex traits. Here, we systematically investigated the population genetic structure of five Anhui indigenous pig breeds (AHIPs), and compared them to those of five Western commercial pig breeds (WECPs). Furthermore, we examined the occurrence and distribution of ROHs in the five AHIPs and estimated the inbreeding coefficients based on the ROHs (FROH) and homozygosity (FHOM). Finally, we identified genomic regions with high frequencies of ROHs and annotated candidate genes contained therein. Results The WECPs and AHIPs were clearly differentiated into two separate clades consistent with their geographical origins, as revealed by the population structure and principal component analysis. We identified 13,530 ROHs across all individuals, of which 4,555 and 8,975 ROHs were unique to AHIPs and WECPs, respectively. Most ROHs identified in our study were short (< 10 Mb) or medium (10–20 Mb) in length. WECPs had significantly higher numbers of short ROHs, and AHIPs generally had longer ROHs. FROH values were significantly lower in AHIPs than in WECPs, indicating that breed improvement and conservation programmes were successful in AHIPs. On average, FROH and FHOM values were highly correlated (0.952–0.991) in AHIPs and WECPs. A total of 27 regions had a high frequency of ROHs and contained 17 key candidate genes associated with economically important traits in pigs. Among these, nine candidate genes (CCNT2, EGR2, MYL3, CDH13, PROX1, FLVCR1, SETD2, FGF18, and FGF20) found in WECPs were related to muscular and skeletal development, whereas eight candidate genes (CSN1S1, SULT1E1, TJP1, ZNF366, LIPC, MCEE, STAP1, and DUSP) found in AHIPs were associated with health, reproduction, and fatness traits. Conclusion Our findings provide a useful reference for the selection and assortative mating of pig breeds, laying the groundwork for future research on the population genetic structures of AHIPs, ultimately helping protect these local varieties. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08583-9.
Collapse
Affiliation(s)
- Yao Jiang
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China
| | - Xiaojin Li
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China
| | - Jiali Liu
- Key Laboratory of Animal Genetics, Breeding and Reproduction (Poultry) of Ministry of Agriculture, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Wei Zhang
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China
| | - Mei Zhou
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China
| | - Jieru Wang
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China
| | - Linqing Liu
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China
| | - Shiguang Su
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China
| | - Fuping Zhao
- Key Laboratory of Animal Genetics, Breeding and Reproduction (Poultry) of Ministry of Agriculture, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Hongquan Chen
- College of Animal Science and Technology, Anhui Agricultural University, Hefei, 230036, China
| | - Chonglong Wang
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China.
| |
Collapse
|
29
|
Nguyen TN, Chen N, Cosgrove EJ, Bowman R, Fitzpatrick JW, Clark AG. Dynamics of reduced genetic diversity in increasingly fragmented populations of Florida
Scrub‐Jays
,
Aphelocoma coerulescens. Evol Appl 2022; 15:1018-1027. [PMID: 35782006 PMCID: PMC9234620 DOI: 10.1111/eva.13421] [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/16/2021] [Revised: 04/04/2022] [Accepted: 05/04/2022] [Indexed: 11/29/2022] Open
Abstract
Understanding the genomic consequences of population decline is important for predicting species' vulnerability to intensifying global change. Empirical information about genomic changes in populations in the early stages of decline, especially for those still experiencing immigration, remains scarce. We used 7834 autosomal SNPs and demographic data for 288 Florida scrub jays (Aphelocoma coerulescens; FSJ) sampled in 2000 and 2008 to compare levels of genetic diversity, inbreeding, relatedness, and lengths of runs of homozygosity (ROH) between two subpopulations within dispersal distance of one another but have experienced contrasting demographic trajectories. At Archbold Biological Station (ABS), the FSJ population has been stable because of consistent habitat protection and management, while at nearby Placid Lakes Estates (PLE), the population declined precipitously due to suburban development. By the onset of our sampling in 2000, birds in PLE were already less heterozygous, more inbred, and on average more related than birds in ABS. No significant changes occurred in heterozygosity or inbreeding across the 8‐year sampling interval, but average relatedness among individuals decreased in PLE, thus by 2008 average relatedness did not differ between sites. PLE harbored a similar proportion of short ROH but a greater proportion of long ROH than ABS, suggesting one continuous population of shared demographic history in the past, which is now experiencing more recent inbreeding. These results broadly uphold the predictions of simple population genetic models based on inferred effective population sizes and rates of immigration. Our study highlights how, in just a few generations, formerly continuous populations can diverge in heterozygosity and levels of inbreeding with severe local population decline despite ongoing gene flow.
Collapse
Affiliation(s)
- Tram N. Nguyen
- Department of Ecology and Evolutionary Biology Cornell University Ithaca NY USA
- Cornell Lab of Ornithology Ithaca NY USA
| | - Nancy Chen
- Department of Biology University of Rochester Rochester NY USA
| | - Elissa J. Cosgrove
- Department of Molecular Biology and Genetics Cornell University Ithaca NY USA
| | - Reed Bowman
- Avian Ecology Lab, Archbold Biological Station, Venus FL USA
| | - John W. Fitzpatrick
- Department of Ecology and Evolutionary Biology Cornell University Ithaca NY USA
- Cornell Lab of Ornithology Ithaca NY USA
| | - Andrew. G. Clark
- Department of Ecology and Evolutionary Biology Cornell University Ithaca NY USA
- Avian Ecology Lab, Archbold Biological Station, Venus FL USA
| |
Collapse
|
30
|
Li G, Tang J, Huang J, Jiang Y, Fan Y, Wang X, Ren J. Genome-Wide Estimates of Runs of Homozygosity, Heterozygosity, and Genetic Load in Two Chinese Indigenous Goat Breeds. Front Genet 2022; 13:774196. [PMID: 35559012 PMCID: PMC9086400 DOI: 10.3389/fgene.2022.774196] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2021] [Accepted: 03/14/2022] [Indexed: 11/13/2022] Open
Abstract
Runs of homozygosity (ROH) and heterozygosity (ROHet) are windows into population demographic history and adaptive evolution. Numerous studies have shown that deleterious mutations are enriched in the ROH of humans, pigs, cattle, and chickens. However, the relationship of deleterious variants to ROH and the pattern of ROHet in goats have been largely understudied. Here, 240 Guangfeng and Ganxi goats from Jiangxi Province, China, were genotyped using the Illumina GoatSNP50 BeadChip and genome-wide ROH, ROHet, and genetic load analyses were performed in the context of 32 global goat breeds. The classes with the highest percentage of ROH and ROHet were 0.5–2 Mb and 0.5–1 Mb, respectively. The results of inbreeding coefficients (based on SNP and ROH) and ROHet measurements showed that Guangfeng goats had higher genetic variability than most Chinese goats, while Ganxi goats had a high degree of inbreeding, even exceeding that of commercial goat breeds. Next, the predicted damaging homozygotes were more enriched in long ROHs, especially in Guangfeng goats. Therefore, we suggest that information on damaging alleles should also be incorporated into the design of breeding and conservation programs. A list of genes related to fecundity, growth, and environmental adaptation were identified in the ROH hotspots of two Jiangxi goats. A sense-related ROH hotspot (chromosome 12: 50.55–50.81 Mb) was shared across global goat breeds and may have undergone selection prior to goat domestication. Furthermore, an identical ROHet hotspot (chromosome 1: 132.21–132.54 Mb) containing two genes associated with embryonic development (STAG1 and PCCB) was detected in domestic goat breeds worldwide. Tajima’s D and BetaScan2 statistics indicated that this region may be caused by long-term balancing selection. These findings not only provide guidance for the design of conservation strategies for Jiangxi goat breeds but also enrich our understanding of the adaptive evolution of goats.
Collapse
Affiliation(s)
- Guixin Li
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou, China
| | - Jianhong Tang
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou, China.,Laboratory Animal Engineering Research Center of Ganzhou, Gannan Medical University, Ganzhou, China
| | - Jinyan Huang
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou, China
| | - Yongchuang Jiang
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou, China
| | - Yin Fan
- Department of Animal Science, Jiangxi Biotech Vocational College, Nanchang, China
| | - Xiaopeng Wang
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou, China
| | - Jun Ren
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou, China
| |
Collapse
|
31
|
Yıldız B, Megens H, Hvilsom C, Bosse M. Genomic consequences of a century of inbreeding and isolation in the Danish wild boar population. Evol Appl 2022; 15:954-966. [PMID: 35782012 PMCID: PMC9234630 DOI: 10.1111/eva.13385] [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: 10/25/2021] [Revised: 02/26/2022] [Accepted: 03/07/2022] [Indexed: 11/30/2022] Open
Abstract
Demographic events such as series of bottlenecks impact the genetic variation and adaptive potential of populations. European megafauna, such as wild boars (Sus scrofa), have experienced severe climatic and size fluctuations that have shaped their genetic variation. Habitat fragmentation and human‐mediated translocations have further contributed to the complex demographic history of European wild boar. Danish wild boars represent an extreme case of a small and isolated population founded by four wild boars from Germany. Here, we explore the genetic composition of the Danish wild boar population in Klelund. We genotyped all 21 Danish wild boars that were recently transferred from the source population in Lille Vildmose into the Klelund Plantation to establish a novel wild boar population. We compared the Danish wild boars with high‐density single‐nucleotide polymorphism genotypes from a comprehensive reference set of 1263 wild and domesticated pigs, including 11 individuals from Ulm, one of two presumed founder locations in Germany. Our findings support the European wild background of the Danish population, and no traces of gene flow with wild or domesticated pigs were found. The narrow genetic origin of the Danish wild boars is illustrated by extremely long and frequent runs of homozygous stretches in their genomes, indicative of recent inbreeding. This study provides the first insights into one of the most inbred wild boar populations globally established a century ago from a narrow base of only four founders.
Collapse
Affiliation(s)
- Beril Yıldız
- Wageningen University & Research ‐ Animal Breeding and Genomics Wageningen The Netherlands
- Netherlands Institute of Ecology (NIOO‐KNAW) – Department of Animal Ecology Wageningen the Netherlands
| | - Hendrik‐Jan Megens
- Wageningen University & Research ‐ Animal Breeding and Genomics Wageningen The Netherlands
| | | | - Mirte Bosse
- Wageningen University & Research ‐ Animal Breeding and Genomics Wageningen The Netherlands
- Department of Ecological Science Animal Ecology Group Vrije Universiteit Amsterdam Amsterdam the Netherlands
| |
Collapse
|
32
|
Wang Y, Dong R, Li X, Cui C, Yu G. Analysis of the Genetic Diversity and Family Structure of the Licha Black Pig Population on Jiaodong Peninsula, Shandong Province, China. Animals (Basel) 2022; 12:ani12081045. [PMID: 35454291 PMCID: PMC9026534 DOI: 10.3390/ani12081045] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2022] [Revised: 04/14/2022] [Accepted: 04/16/2022] [Indexed: 12/31/2022] Open
Abstract
Simple Summary This study investigated the current conservation status, including the genetic diversity, the family structure, and inbreeding, of the Licha black pig population on Jiaodong Peninsula (Shandong Province, China). The results provide insights into the management and conservation of a local pig breed. Breeders are encouraged to utilize genomic data to improve mating schemes based on the family information obtained in this study, such as keeping an equivalent number of boars and sows in each family and selecting individuals with a kinship coefficient of less than 0.1 for mating. Abstract The Licha black pig, a popular indigenous Chinese pig breed, is known for its multi-vertebral trait and higher lean meat rate. Understanding the current conservation status, family structure, and degree of inbreeding of the Licha black pig population will be useful to maintain a sufficient level of genetic diversity in these animal resources. In the present study, the genetic diversity, population structure, and inbreeding coefficient of this conserved population were analyzed using SNP genotyping data from 209 Licha black pigs. Based on the genomic information, this population was divided into eight different families with boars. The effective population size (Ne), polymorphic marker ratio (PN), expected heterozygosity (He), and observed heterozygosity (Ho) of this population were 8.7, 0.827, 0.3576, and 0.3512, respectively. In addition, a total of 5976 runs of homozygosity (ROHs) were identified, and most of the ROHs (54.9%) were greater than 5 Mb. The genomic inbreeding coefficient of each individual was estimated based on ROHs (FROH) with an average inbreeding coefficient of 0.11 for the population. Five statistics (Ne, PN, Ho, He, and FROH) showed a decrease in the level of genetic diversity and a high degree of inbreeding in this population. Thus, special preservation programs need to be implemented in the future, such as introducing new individuals or improving the mating plan. Altogether, our study provides the first genomic overview of the genetic diversity and population structure of Licha black pigs, which will be useful for the management and long-term preservation of this breed.
Collapse
Affiliation(s)
- Yuan Wang
- College of Animal Science and Technology, Qingdao Agricultural University, Qingdao 266109, China; (Y.W.); (R.D.); (X.L.)
| | - Ruilan Dong
- College of Animal Science and Technology, Qingdao Agricultural University, Qingdao 266109, China; (Y.W.); (R.D.); (X.L.)
| | - Xiao Li
- College of Animal Science and Technology, Qingdao Agricultural University, Qingdao 266109, China; (Y.W.); (R.D.); (X.L.)
| | - Chao Cui
- Bureau of Agriculture and Rural Affairs of Jiaozhou, Jiaozhou 266300, China;
| | - Guanghui Yu
- College of Animal Science and Technology, Qingdao Agricultural University, Qingdao 266109, China; (Y.W.); (R.D.); (X.L.)
- Correspondence:
| |
Collapse
|
33
|
Kim YM, Seong HS, Kim YS, Hong JK, Sa SJ, Lee J, Lee JH, Cho KH, Chung WH, Choi JW, Cho ES. Genome-Wide Assessment of a Korean Composite Pig Breed, Woori-Heukdon. Front Genet 2022; 13:779152. [PMID: 35186025 PMCID: PMC8847790 DOI: 10.3389/fgene.2022.779152] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2021] [Accepted: 01/07/2022] [Indexed: 02/03/2023] Open
Abstract
A Korean synthetic pig breed, Woori-Heukdon (WRH; F3), was developed by crossing parental breeds (Korean native pig [KNP] and Korean Duroc [DUC]) with their crossbred populations (F1 and F2). This study in genome-wide assessed a total of 2,074 pigs which include the crossbred and the parental populations using the Illumina PorcineSNP60 BeadChip. After quality control of the initial datasets, we performed population structure, genetic diversity, and runs of homozygosity (ROH) analyses. Population structure analyses showed that crossbred populations were genetically influenced by the parental breeds according to their generation stage in the crossbreeding scheme. Moreover, principal component analysis showed the dispersed cluster of WRH, which might reflect introducing a new breeding group into the previous one. Expected heterozygosity values, which were used to assess genetic diversity, were .365, .349, .336, .330, and .211 for WRH, F2, F1, DUC, and KNP, respectively. The inbreeding coefficient based on ROH was the highest in KNP (.409), followed by WRH (.186), DUC (.178), F2 (.107), and F1 (.035). Moreover, the frequency of short ROH decreased according to the crossing stage (from F1 to WRH). Alternatively, the frequency of medium and long ROH increased, which indicated recent inbreeding in F2 and WRH. Furthermore, gene annotation of the ROH islands in WRH that might be inherited from their parental breeds revealed several interesting candidate genes that may be associated with adaptation, meat quality, production, and reproduction traits in pigs.
Collapse
Affiliation(s)
- Yong-Min Kim
- Swine Science Division, National Institute of Animal Science, Rural Development Administration, Cheonan, South Korea.,Department of Animal Science, College of Animal Life Sciences, Kangwon National University, Chuncheon, South Korea
| | - Ha-Seung Seong
- Swine Science Division, National Institute of Animal Science, Rural Development Administration, Cheonan, South Korea.,Department of Animal Science, College of Animal Life Sciences, Kangwon National University, Chuncheon, South Korea
| | - Young-Sin Kim
- Swine Science Division, National Institute of Animal Science, Rural Development Administration, Cheonan, South Korea
| | - Joon-Ki Hong
- Swine Science Division, National Institute of Animal Science, Rural Development Administration, Cheonan, South Korea
| | - Soo-Jin Sa
- Swine Science Division, National Institute of Animal Science, Rural Development Administration, Cheonan, South Korea
| | - Jungjae Lee
- Department of Animal Science and Technology, College of Biotechnology and Natural Resources, Chung-Ang University, Anseong, South Korea
| | - Jun-Hee Lee
- Institute of Agriculture and Life Science, College of Agriculture and Life Sciences, Gyeongsang National University, Jinju, South Korea
| | - Kyu-Ho Cho
- Swine Science Division, National Institute of Animal Science, Rural Development Administration, Cheonan, South Korea
| | - Won-Hyong Chung
- Research Group of Healthcare, Korea Food Research Institute, Wanju, South Korea
| | - Jung-Woo Choi
- Department of Animal Science, College of Animal Life Sciences, Kangwon National University, Chuncheon, South Korea
| | - Eun-Seok Cho
- Swine Science Division, National Institute of Animal Science, Rural Development Administration, Cheonan, South Korea
| |
Collapse
|
34
|
Genetic Diversity, Admixture and Analysis of Homozygous-by-Descent (HBD) Segments of Russian Wild Boar. BIOLOGY 2022; 11:biology11020203. [PMID: 35205070 PMCID: PMC8869248 DOI: 10.3390/biology11020203] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/19/2021] [Revised: 01/25/2022] [Accepted: 01/25/2022] [Indexed: 02/04/2023]
Abstract
The wild boar is the wild ancestor of the domestic pig and one of the most common species of ungulates. At the beginning of the 20th century, the wild boar was practically exterminated in the European part of Russia. In the period 1935-1988, 7705 boars were caught in various regions of the European part of Russia, the Far East, Ukraine, Belarus, Kyrgyzstan, Kazakhstan, Latvia, Lithuania, Estonia, Tajikistan and resettled in the territory of Russia. Asian and European wild boars dwell the territory of Russia. The aim of our research was to study the genetic diversity and structure of wild boar populations in different regions of Russia using genome-wide genotyping. We have determined the genetic distances, population structure, parameters of genetic diversity and significantly expanded our understanding of the genetic state of the Russian wild boar. For the first time, we calculated autozygosity of the wild boar of the European and Asian subspecies using Homozygous-by-Descent (HBD) Segments analysis, which is important in terms of population recovery. We also found evidence of hybridization between Russian wild boar and domestic pigs. A group of European wild boars showed introgression of the Asian boar into population. The mean level of the inbreeding coefficient in European wild boar was higher than in Asian wild boar, and combined groups of the European boar had higher inbreeding coefficient than Russian wild boars. These results obtained can be used in population management.
Collapse
|
35
|
De los Ríos-Pérez L, Druet T, Goldammer T, Wittenburg D. Analysis of Autozygosity Using Whole-Genome Sequence Data of Full-Sib Families in Pikeperch (Sander lucioperca). Front Genet 2022; 12:786934. [PMID: 35111201 PMCID: PMC8801746 DOI: 10.3389/fgene.2021.786934] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2021] [Accepted: 12/23/2021] [Indexed: 11/13/2022] Open
Abstract
Pikeperch (Sander lucioperca) has emerged as a high value species to the aquaculture industry. However, its farming techniques are at an early stage and its production is often performed without a selective breeding program, potentially leading to high levels of inbreeding. In this study, we identified and characterized autozygosity based on genome-wide runs of homozygosity (ROH) on a sample of parental and offspring individuals, determined effective population size (Ne), and assessed relatedness among parental individuals. A mean of 2,235 ± 526 and 1,841 ± 363 ROH segments per individual, resulting in a mean inbreeding coefficient of 0.33 ± 0.06 and 0.25 ± 0.06 were estimated for the progeny and parents, respectively. Ne was about 12 until four generations ago and at most 106 for 63 generations in the past, with varying genetic relatedness amongst the parents. This study shows the importance of genomic information when family relationships are unknown and the need of selective breeding programs for reproductive management decisions in the aquaculture industry.
Collapse
Affiliation(s)
- Lidia De los Ríos-Pérez
- Institute of Genetics and Biometry, Research Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
| | - Tom Druet
- Unit of Animal Genomics, GIGA-R & Faculty of Veterinary Medicine, University of Liège, Liège, Belgium
| | - Tom Goldammer
- Institute of Genome Biology, Research Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
- Molecular Biology and Fish Genetics, Faculty of Agriculture and Environmental Sciences, University of Rostock, Rostock, Germany
| | - Dörte Wittenburg
- Institute of Genetics and Biometry, Research Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
- *Correspondence: Dörte Wittenburg,
| |
Collapse
|
36
|
Wolfsberger WW, Ayala NM, Castro-Marquez SO, Irizarry-Negron VM, Potapchuk A, Shchubelka K, Potish L, Majeske AJ, Oliver LF, Lameiro AD, Martínez-Cruzado JC, Lindgren G, Oleksyk TK. Genetic diversity and selection in Puerto Rican horses. Sci Rep 2022; 12:515. [PMID: 35017609 PMCID: PMC8752667 DOI: 10.1038/s41598-021-04537-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2021] [Accepted: 12/23/2021] [Indexed: 11/21/2022] Open
Abstract
Since the first Spanish settlers brought horses to America centuries ago, several local varieties and breeds have been established in the New World. These were generally a consequence of the admixture of the different breeds arriving from Europe. In some instances, local horses have been selectively bred for specific traits, such as appearance, endurance, strength, and gait. We looked at the genetics of two breeds, the Puerto Rican Non-Purebred (PRNPB) (also known as the "Criollo") horses and the Puerto Rican Paso Fino (PRPF), from the Caribbean Island of Puerto Rico. While it is reasonable to assume that there was a historic connection between the two, the genetic link between them has never been established. In our study, we started by looking at the genetic ancestry and diversity of current Puerto Rican horse populations using a 668 bp fragment of the mitochondrial DNA D-loop (HVR1) in 200 horses from 27 locations on the island. We then genotyped all 200 horses in our sample for the "gait-keeper" DMRT3 mutant allele previously associated with the paso gait especially cherished in this island breed. We also genotyped a subset of 24 samples with the Illumina Neogen Equine Community genome-wide array (65,000 SNPs). This data was further combined with the publicly available PRPF genomes from other studies. Our analysis show an undeniable genetic connection between the two varieties in Puerto Rico, consistent with the hypothesis that PRNPB horses represent the descendants of the original genetic pool, a mix of horses imported from the Iberian Peninsula and elsewhere in Europe. Some of the original founders of PRNRB population must have carried the "gait-keeper" DMRT3 allele upon arrival to the island. From this admixture, the desired traits were selected by the local people over the span of centuries. We propose that the frequency of the mutant "gait-keeper" allele originally increased in the local horses due to the selection for the smooth ride and other characters, long before the PRPF breed was established. To support this hypothesis, we demonstrate that PRNPB horses, and not the purebred PRPF, carry a signature of selection in the genomic region containing the DMRT3 locus to this day. The lack of the detectable signature of selection associated with the DMRT3 in the PRPF would be expected if this native breed was originally derived from the genetic pool of PRNPB horses established earlier and most of the founders already had the mutant allele. Consequently, selection specific to PRPF later focused on allels in other genes (including CHRM5, CYP2E1, MYH7, SRSF1, PAM, PRN and others) that have not been previously associated with the prized paso gait phenotype in Puerto Rico or anywhere else.
Collapse
Affiliation(s)
- Walter W Wolfsberger
- Department of Biological Sciences, Oakland University, Rochester, MI, USA
- Biology Department, University of Puerto Rico at Mayaguez, Mayaguez, Puerto Rico
- Biology Department, Uzhhorod National University, Uzhhorod, Ukraine
| | - Nikole M Ayala
- Biology Department, University of Puerto Rico at Mayaguez, Mayaguez, Puerto Rico
| | - Stephanie O Castro-Marquez
- Department of Biological Sciences, Oakland University, Rochester, MI, USA
- Biology Department, University of Puerto Rico at Mayaguez, Mayaguez, Puerto Rico
| | | | - Antoliy Potapchuk
- Biology Department, University of Puerto Rico at Mayaguez, Mayaguez, Puerto Rico
| | - Khrystyna Shchubelka
- Department of Biological Sciences, Oakland University, Rochester, MI, USA
- Biology Department, University of Puerto Rico at Mayaguez, Mayaguez, Puerto Rico
- Biology Department, Uzhhorod National University, Uzhhorod, Ukraine
| | - Ludvig Potish
- Department of Forestry, Uzhhorod National University, Uzhhorod, Ukraine
| | - Audrey J Majeske
- Department of Biological Sciences, Oakland University, Rochester, MI, USA
- Biology Department, University of Puerto Rico at Mayaguez, Mayaguez, Puerto Rico
| | - Luis Figueroa Oliver
- Biology Department, University of Puerto Rico at Mayaguez, Mayaguez, Puerto Rico
| | - Alondra Diaz Lameiro
- Biology Department, University of Puerto Rico at Mayaguez, Mayaguez, Puerto Rico
| | | | - Gabriella Lindgren
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden
- Livestock Genetics, Department of Biosystems, KU Leuven, Leuven, Belgium
| | - Taras K Oleksyk
- Department of Biological Sciences, Oakland University, Rochester, MI, USA.
- Biology Department, University of Puerto Rico at Mayaguez, Mayaguez, Puerto Rico.
- Biology Department, Uzhhorod National University, Uzhhorod, Ukraine.
| |
Collapse
|
37
|
Machado AP, Topaloudis A, Cumer T, Lavanchy E, Bontzorlos V, Ceccherelli R, Charter M, Kassinis N, Lymberakis P, Manzia F, Ducrest AL, Dupasquier M, Guex N, Roulin A, Goudet J. Genomic consequences of colonisation, migration and genetic drift in barn owl insular populations of the eastern Mediterranean. Mol Ecol 2021; 31:1375-1388. [PMID: 34894026 PMCID: PMC9305133 DOI: 10.1111/mec.16324] [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/17/2021] [Revised: 11/01/2021] [Accepted: 11/17/2021] [Indexed: 01/25/2023]
Abstract
The study of insular populations was key in the development of evolutionary theory. The successful colonisation of an island depends on the geographic context, and specific characteristics of the organism and the island, but also on stochastic processes. As a result, apparently identical islands may harbour populations with contrasting histories. Here, we use whole genome sequences of 65 barn owls to investigate the patterns of inbreeding and genetic diversity of insular populations in the eastern Mediterranean Sea. We focus on Crete and Cyprus, islands with similar size, climate and distance to mainland, that provide natural replicates for a comparative analysis of the impacts of microevolutionary processes on isolated populations. We show that barn owl populations from each island have a separate origin, Crete being genetically more similar to other Greek islands and mainland Greece, and Cyprus more similar to the Levant. Further, our data show that their respective demographic histories following colonisation were also distinct. On the one hand, Crete harbours a small population and maintains very low levels of gene flow with neighbouring populations. This has resulted in low genetic diversity, strong genetic drift, increased relatedness in the population and remote inbreeding. Cyprus, on the other hand, appears to maintain enough gene flow with the mainland to avoid such an outcome. Our study provides a comparative population genomic analysis of the effects of neutral processes on a classical island‐mainland model system. It provides empirical evidence for the role of stochastic processes in determining the fate of diverging isolated populations.
Collapse
Affiliation(s)
- Ana Paula Machado
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
| | | | - Tristan Cumer
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
| | - Eléonore Lavanchy
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
| | - Vasileios Bontzorlos
- Green Fund, Kifisia, Athens, Greece.,"TYTO" - Organization for the Management and Conservation of Biodiversity in Agricultural Ecosystems, Larisa, Greece
| | | | - Motti Charter
- Shamir Research Institute, University of Haifa, Katzrin, Israel.,Department of Geography and Environmental Sciences, University of Haifa, Haifa, Israel
| | | | - Petros Lymberakis
- Natural History Museum of Crete, University of Crete, Herakleio, Greece
| | | | - Anne-Lyse Ducrest
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
| | | | - Nicolas Guex
- Bioinformatics Competence Centre, University of Lausanne, Lausanne, Switzerland
| | - Alexandre Roulin
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
| | - Jérôme Goudet
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland.,Swiss Institute of Bioinformatics, Lausanne, Switzerland
| |
Collapse
|
38
|
Sanglard LP, Huang Y, Gray KA, Linhares DCL, Dekkers JCM, Niederwerder MC, Fernando RL, Serão NVL. Further host-genomic characterization of total antibody response to PRRSV vaccination and its relationship with reproductive performance in commercial sows: genome-wide haplotype and zygosity analyses. Genet Sel Evol 2021; 53:91. [PMID: 34875996 PMCID: PMC8650375 DOI: 10.1186/s12711-021-00676-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2021] [Accepted: 10/19/2021] [Indexed: 11/16/2022] Open
Abstract
Background The possibility of using antibody response (S/P ratio) to PRRSV vaccination measured in crossbred commercial gilts as a genetic indicator for reproductive performance in vaccinated crossbred sows has motivated further studies of the genomic basis of this trait. In this study, we investigated the association of haplotypes and runs of homozygosity (ROH) and heterozygosity (ROHet) with S/P ratio and their impact on reproductive performance. Results There was no association (P-value ≥ 0.18) of S/P ratio with the percentage of ROH or ROHet, or with the percentage of heterozygosity across the whole genome or in the major histocompatibility complex (MHC) region. However, specific ROH and ROHet regions were significantly associated (P-value ≤ 0.01) with S/P ratio on chromosomes 1, 4, 5, 7, 10, 11, 13, and 17 but not (P-value ≥ 0.10) with reproductive performance. With the haplotype-based genome-wide association study (GWAS), additional genomic regions associated with S/P ratio were identified on chromosomes 4, 7, and 9. These regions harbor immune-related genes, such as SLA-DOB, TAP2, TAPBP, TMIGD3, and ADORA. Four haplotypes at the identified region on chromosome 7 were also associated with multiple reproductive traits. A haplotype significantly associated with S/P ratio that is located in the MHC region may be in stronger linkage disequilibrium (LD) with the quantitative trait loci (QTL) than the previously identified single nucleotide polymorphism (SNP) (H3GA0020505) given the larger estimate of genetic variance explained by the haplotype than by the SNP. Conclusions Specific ROH and ROHet regions were significantly associated with S/P ratio. The haplotype-based GWAS identified novel QTL for S/P ratio on chromosomes 4, 7, and 9 and confirmed the presence of at least one QTL in the MHC region. The chromosome 7 region was also associated with reproductive performance. These results narrow the search for causal genes in this region and suggest SLA-DOB and TAP2 as potential candidate genes associated with S/P ratio on chromosome 7. These results provide additional opportunities for marker-assisted selection and genomic selection for S/P ratio as genetic indicator for litter size in commercial pig populations. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-021-00676-5.
Collapse
Affiliation(s)
- Leticia P Sanglard
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA
| | - Yijian Huang
- Smithfield Premium Genetic, Rose Hill, NC, 28458, USA
| | - Kent A Gray
- Smithfield Premium Genetic, Rose Hill, NC, 28458, USA
| | - Daniel C L Linhares
- Department of Veterinary Diagnostic & Production Animal Medicine, Iowa State University, Ames, IA, 50011, USA
| | - Jack C M Dekkers
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA
| | - Megan C Niederwerder
- Department of Diagnostic Medicine/Pathobiology, Kansas State University, Manhattan, KS, 66506, USA
| | - Rohan L Fernando
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA
| | - Nick V L Serão
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA.
| |
Collapse
|
39
|
Tsartsianidou V, Sánchez-Molano E, Kapsona VV, Basdagianni Z, Chatziplis D, Arsenos G, Triantafyllidis A, Banos G. A comprehensive genome-wide scan detects genomic regions related to local adaptation and climate resilience in Mediterranean domestic sheep. Genet Sel Evol 2021; 53:90. [PMID: 34856922 PMCID: PMC8641236 DOI: 10.1186/s12711-021-00682-7] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2021] [Accepted: 11/09/2021] [Indexed: 02/07/2023] Open
Abstract
Background The management of farm animal genetic resources and the adaptation of animals to climate change will probably have major effects on the long-term sustainability of the livestock sector. Genomic data harbour useful relevant information that needs to be harnessed for effectively managing genetic resources. In this paper, we report the genome characterization of the highly productive Mediterranean Chios dairy sheep and focus on genetic diversity measures related with local adaptation and selection and the genetic architecture of animal resilience to weather fluctuations as a novel adaptative trait linked to climate change. Results We detected runs of homozygosity (ROH) and heterozygosity (ROHet) that revealed multiple highly homozygous and heterozygous hotspots across the Chios sheep genome. A particularly highly homozygous region was identified on chromosome 13 as a candidate of directional genetic selection associated with milk traits, which includes annotated genes that were previously shown to be linked to local adaptation to harsh environmental conditions. Favourable heterozygosity related with a potentially protective role against livestock diseases and enhanced overall fitness was revealed in heterozygous-rich regions on sheep chromosomes 3, 10, 13 and 19. Furthermore, genomic analyses were conducted on sheep resilience phenotypes that display changes in milk production in response to weather variation. Sheep resilience to heat stress was a significantly heritable trait (h2 = 0.26) and genetically antagonistic to milk production. Genome-wide association and regional heritability mapping analyses revealed novel genomic markers and regions on chromosome 5 that were significantly associated with sheep resilience to climate change. Subsequently, an annotation analysis detected a set of genes on chromosome 5 that were associated with olfactory receptor complexes that could participate in heat stress mitigation through changes in respiration rate and respiratory evaporation. Other genes were grouped in previously reported biological processes relevant to livestock heat dissipation, including stress and immune response. Conclusions Our results may contribute to the optimal management of sheep genetic resources and inform modern selective breeding programmes that aim at mitigating future environmental challenges towards sustainable farming, while better balancing animal adaptation and productivity. Our results are directly relevant to the studied breed and the respective environmental conditions; however, the methodology may be extended to other livestock species of interest. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-021-00682-7.
Collapse
Affiliation(s)
- Valentina Tsartsianidou
- Department of Genetics, Development & Molecular Biology, School of Biology, Aristotle University of Thessaloniki, 54124, Thessaloniki, Greece.
| | - Enrique Sánchez-Molano
- Division of Genetics and Genomics, School of Veterinary Studies, The Roslin Institute and Royal (Dick), University of Edinburgh, Easter Bush, Midlothian, EH25 9RG, UK
| | - Vanessa Varvara Kapsona
- Department of Animal and Veterinary Sciences, Scotland's Rural College, Roslin Institute Building, Easter Bush, Midlothian, EH25 9RG, UK
| | - Zoitsa Basdagianni
- Department of Animal Production, School of Agriculture, Aristotle University of Thessaloniki, 54124, Thessaloniki, Greece
| | - Dimitrios Chatziplis
- Laboratory of Agrobiotechnology and Inspection of Agricultural Products, Department of Agriculture, International Hellenic University, Alexander Campus, 57400, Sindos, Greece
| | - Georgios Arsenos
- Laboratory of Animal Husbandry, School of Veterinary Medicine, Aristotle University of Thessaloniki, 54124, Thessaloniki, Greece
| | - Alexandros Triantafyllidis
- Department of Genetics, Development & Molecular Biology, School of Biology, Aristotle University of Thessaloniki, 54124, Thessaloniki, Greece
| | - Georgios Banos
- Department of Animal and Veterinary Sciences, Scotland's Rural College, Roslin Institute Building, Easter Bush, Midlothian, EH25 9RG, UK.,Laboratory of Animal Husbandry, School of Veterinary Medicine, Aristotle University of Thessaloniki, 54124, Thessaloniki, Greece
| |
Collapse
|
40
|
Weng Z, Xu Y, Zhong M, Li W, Chen J, Zhong F, Du B, Zhang B, Huang X. Runs of homozygosity analysis reveals population characteristics of yellow-feathered chickens using re-sequencing data. Br Poult Sci 2021; 63:307-315. [PMID: 34747677 DOI: 10.1080/00071668.2021.2003752] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Abstract
1. To evaluate the inbreeding of yellow-feathered chickens (YFCs) and identify genes related to their unique characteristics, whole-genome re-sequencing data were applied to detect runs of homozygosity (ROH) in the genomes of ten YFC breeds. The number, length, distribution of ROH, and inbreeding coefficient in different YFC populations were calculated. Genomic regions with high frequency in ROH were annotated.2. In total, 25,547 ROH with an average length of 335 kb were detected, with most being <1 Mb. The domination of short ROH reflected the long breeding history of this chicken. The number, length, frequency, and distribution of ROH varied among chicken populations, and high genetic diversity was maintained.3. Numerous genes related to YFC characteristics were identified in the high-frequency ROH regions. Among these, IFNA, IFNB, IL11RA, IL22RA1, IFNLR1, and TRIF genes were involved in disease resistance. The AMY, G6PC, SDHB, GCNT4, and ACO genes were associated with energy material metabolism; and FABPL, AQP7, ACAA2, and RYR2 were related to meat quality and flavour. The KITLG, CREB3, RYR2, and LGR4 genes, related to pigmentation, were detected.4. This ROH-based inbreeding evaluation laid the foundation for breeding and conservation of YFC populations, and the candidate genes identified can be used for marker-assisted selection.
Collapse
Affiliation(s)
- Zhuoxian Weng
- Jiaying University/Guangdong Provincial Key Laboratory of Conservation and Precision Utilization of Characteristic Agricultural Resources in Mountainous Areas, Meizhou 514015, China.,College of Animal Science and Technology, Hunan Agricultural University, Changsha 410128, China.,Guangdong Innovation Centre for Science and Technology of Wuhua Yellow Chicken, Meizhou, 514015, China
| | - Yongjie Xu
- Jiaying University/Guangdong Provincial Key Laboratory of Conservation and Precision Utilization of Characteristic Agricultural Resources in Mountainous Areas, Meizhou 514015, China.,Guangdong Innovation Centre for Science and Technology of Wuhua Yellow Chicken, Meizhou, 514015, China
| | - Ming Zhong
- Jiaying University/Guangdong Provincial Key Laboratory of Conservation and Precision Utilization of Characteristic Agricultural Resources in Mountainous Areas, Meizhou 514015, China.,Guangdong Innovation Centre for Science and Technology of Wuhua Yellow Chicken, Meizhou, 514015, China
| | - Weina Li
- Jiaying University/Guangdong Provincial Key Laboratory of Conservation and Precision Utilization of Characteristic Agricultural Resources in Mountainous Areas, Meizhou 514015, China.,Guangdong Innovation Centre for Science and Technology of Wuhua Yellow Chicken, Meizhou, 514015, China
| | - Jiebo Chen
- Guangdong Innovation Centre for Science and Technology of Wuhua Yellow Chicken, Meizhou, 514015, China
| | - Fusheng Zhong
- Jiaying University/Guangdong Provincial Key Laboratory of Conservation and Precision Utilization of Characteristic Agricultural Resources in Mountainous Areas, Meizhou 514015, China.,Guangdong Innovation Centre for Science and Technology of Wuhua Yellow Chicken, Meizhou, 514015, China
| | - Bingwang Du
- Jiaying University/Guangdong Provincial Key Laboratory of Conservation and Precision Utilization of Characteristic Agricultural Resources in Mountainous Areas, Meizhou 514015, China.,Guangdong Innovation Centre for Science and Technology of Wuhua Yellow Chicken, Meizhou, 514015, China
| | - Bin Zhang
- College of Animal Science and Technology, Hunan Agricultural University, Changsha 410128, China
| | - Xunhe Huang
- Jiaying University/Guangdong Provincial Key Laboratory of Conservation and Precision Utilization of Characteristic Agricultural Resources in Mountainous Areas, Meizhou 514015, China.,Guangdong Innovation Centre for Science and Technology of Wuhua Yellow Chicken, Meizhou, 514015, China
| |
Collapse
|
41
|
Kava R, Peripolli E, Berton MP, Lemos M, Lobo RB, Stafuzza NB, Pereira AS, Baldi F. Genome-wide structural variations in Brazilian Senepol cattle, a tropically adapted taurine breed. Livest Sci 2021. [DOI: 10.1016/j.livsci.2021.104708] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
42
|
Runs of homozygosity analysis reveals consensus homozygous regions affecting production traits in Chinese Simmental beef cattle. BMC Genomics 2021; 22:678. [PMID: 34548021 PMCID: PMC8454143 DOI: 10.1186/s12864-021-07992-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2021] [Accepted: 09/07/2021] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Genomic regions with a high frequency of runs of homozygosity (ROH) are related to important traits in farm animals. We carried out a comprehensive analysis of ROH and evaluated their association with production traits using the BovineHD (770 K) SNP array in Chinese Simmental beef cattle. RESULTS We detected a total of 116,953 homozygous segments with 2.47Gb across the genome in the studied population. The average number of ROH per individual was 99.03 and the average length was 117.29 Mb. Notably, we detected 42 regions with a frequency of more than 0.2. We obtained 17 candidate genes related to body size, meat quality, and reproductive traits. Furthermore, using Fisher's exact test, we found 101 regions were associated with production traits by comparing high groups with low groups in terms of production traits. Of those, we identified several significant regions for production traits (P < 0.05) by association analysis, within which candidate genes including ECT2, GABRA4, and GABRB1 have been previously reported for those traits in beef cattle. CONCLUSIONS Our study explored ROH patterns and their potential associations with production traits in beef cattle. These results may help to better understand the association between production traits and genome homozygosity and offer valuable insights into managing inbreeding by designing reasonable breeding programs in farm animals.
Collapse
|
43
|
Investigating inbreeding in the turkey (Meleagris gallopavo) genome. Poult Sci 2021; 100:101366. [PMID: 34525446 PMCID: PMC8445901 DOI: 10.1016/j.psj.2021.101366] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2021] [Revised: 06/02/2021] [Accepted: 06/24/2021] [Indexed: 02/06/2023] Open
Abstract
The detrimental effects of increased homozygosity due to inbreeding have prompted the development of methods to reduce inbreeding. The detection of runs of homozygosity (ROH), or contiguous stretches of homozygous marker genotypes, can be used to describe and quantify the level of inbreeding in an individual. The estimation of inbreeding coefficients can be calculated based on pedigree information, ROH, or the genomic relationship matrix. The aim of this study was to detect and describe ROH in the turkey genome and compare estimates of pedigree-based inbreeding coefficients (FPED) with genomic-based inbreeding coefficients estimated from ROH (FROH) and the genomic relationship matrix (FGRM). A total of 2,616,890 pedigree records were available. Of these records, 6,371 genotyped animals from three purebred turkey (Meleagris gallopavo) lines between 2013 and 2019 were available, and these were obtained using a dense single nucleotide polymorphism array (56,452 SNPs). The overall mean length of detected ROH was 2.87 ± 0.29 Mb with a mean number of 84.87 ± 8.79 ROH per animal. Short ROH with lengths of 1 to 2 Mb long were the most abundant throughout the genome. Mean ROH coverage differed greatly between chromosomes and lines. Considering inbreeding coefficient means across all lines, genomic derived inbreeding coefficients (FROH = 0.27; FGRM = 0.32) were higher than coefficients estimated from pedigree records (FPED = 0.14). Correlations between FROH and FPED, FROH and FGRM, and FPED and FGRM ranged between 0.19 to 0.31, 0.68 to 0.73, and 0.17 to 0.30, respectively. Additionally, correlations between FROH from different lengths and FPED substantially increased with ROH length from -0.06 to 0.33. Results of the current research, including the distribution of ROH throughout the genome and ROH-derived inbreeding estimates, can provide a more comprehensive description of inbreeding in the turkey genome. This knowledge can be used to evaluate genetic diversity, a requirement for genetic improvement, and develop methods to minimize inbreeding in turkey breeding programs.
Collapse
|
44
|
Selli A, Ventura RV, Fonseca PAS, Buzanskas ME, Andrietta LT, Balieiro JCC, Brito LF. Detection and Visualization of Heterozygosity-Rich Regions and Runs of Homozygosity in Worldwide Sheep Populations. Animals (Basel) 2021; 11:2696. [PMID: 34573664 PMCID: PMC8472390 DOI: 10.3390/ani11092696] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2021] [Revised: 09/11/2021] [Accepted: 09/13/2021] [Indexed: 12/25/2022] Open
Abstract
In this study, we chose 17 worldwide sheep populations of eight breeds, which were intensively selected for different purposes (meat, milk, or wool), or locally-adapted breeds, in order to identify and characterize factors impacting the detection of runs of homozygosity (ROH) and heterozygosity-rich regions (HRRs) in sheep. We also applied a business intelligence (BI) tool to integrate and visualize outputs from complementary analyses. We observed a prevalence of short ROH, and a clear distinction between the ROH profiles across populations. The visualizations showed a fragmentation of medium and long ROH segments. Furthermore, we tested different scenarios for the detection of HRR and evaluated the impact of the detection parameters used. Our findings suggest that HRRs are small and frequent in the sheep genome; however, further studies with higher density SNP chips and different detection methods are suggested for future research. We also defined ROH and HRR islands and identified common regions across the populations, where genes related to a variety of traits were reported, such as body size, muscle development, and brain functions. These results indicate that such regions are associated with many traits, and thus were under selective pressure in sheep breeds raised for different purposes. Interestingly, many candidate genes detected within the HRR islands were associated with brain integrity. We also observed a strong association of high linkage disequilibrium pattern with ROH compared with HRR, despite the fact that many regions in linkage disequilibrium were not located in ROH regions.
Collapse
Affiliation(s)
- Alana Selli
- Department of Nutrition and Animal Production, School of Veterinary Medicine and Animal Science (FMVZ), University of São Paulo (USP), Pirassununga 13635-900, São Paulo, Brazil; (L.T.A.); (J.C.C.B.)
| | - Ricardo V. Ventura
- Department of Nutrition and Animal Production, School of Veterinary Medicine and Animal Science (FMVZ), University of São Paulo (USP), Pirassununga 13635-900, São Paulo, Brazil; (L.T.A.); (J.C.C.B.)
| | - Pablo A. S. Fonseca
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON N1G 2W1, Canada;
| | - Marcos E. Buzanskas
- Department of Animal Science, Federal University of Paraíba, João Pessoa 58051-900, Paraiba, Brazil;
| | - Lucas T. Andrietta
- Department of Nutrition and Animal Production, School of Veterinary Medicine and Animal Science (FMVZ), University of São Paulo (USP), Pirassununga 13635-900, São Paulo, Brazil; (L.T.A.); (J.C.C.B.)
| | - Júlio C. C. Balieiro
- Department of Nutrition and Animal Production, School of Veterinary Medicine and Animal Science (FMVZ), University of São Paulo (USP), Pirassununga 13635-900, São Paulo, Brazil; (L.T.A.); (J.C.C.B.)
| | - Luiz F. Brito
- Department of Animal Sciences, Purdue University, West Lafayette, IN 47907, USA;
| |
Collapse
|
45
|
Ochoa A, Gibbs HL. Genomic signatures of inbreeding and mutation load in a threatened rattlesnake. Mol Ecol 2021; 30:5454-5469. [PMID: 34448259 DOI: 10.1111/mec.16147] [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: 04/27/2021] [Revised: 08/04/2021] [Accepted: 08/16/2021] [Indexed: 11/28/2022]
Abstract
Theory predicts that threatened species living in small populations will experience high levels of inbreeding that will increase their genetic load, but recent work suggests that the impact of load may be minimized by purging resulting from long-term population bottlenecks. Empirical studies that examine this idea using genome-wide estimates of inbreeding and genetic load in threatened species are limited. Here we use individual genome resequencing data to compare levels of inbreeding, levels of genetic load (estimated as mutation load) and population history in threatened Eastern massasauga rattlesnakes (Sistrurus catenatus), which exist in small isolated populations, and closely related yet outbred Western massasauga rattlesnakes (Sistrurus tergeminus). In terms of inbreeding, S. catenatus genomes had a greater number of runs of homozygosity of varying sizes, indicating sustained inbreeding through repeated bottlenecks when compared to S. tergeminus. At the species level, outbred S. tergeminus had higher genome-wide levels of mutation load in the form of greater numbers of derived deleterious mutations compared to S. catenatus, presumably due to long-term purging of deleterious mutations in S. catenatus. In contrast, mutations that escaped species-level drift effects within S. catenatus populations were in general more frequent and more often found in homozygous genotypes than in S. tergeminus, suggesting a reduced efficiency of purifying selection in smaller S. catenatus populations for most mutations. Our results support an emerging idea that the historical demography of a threatened species has a significant impact on the type of genetic load present, which impacts implementation of conservation actions such as genetic rescue.
Collapse
Affiliation(s)
- Alexander Ochoa
- Department of Evolution, Ecology, and Organismal Biology, Ohio Biodiversity Conservation Partnership, Ohio State University, Columbus, Ohio, USA
| | - H Lisle Gibbs
- Department of Evolution, Ecology, and Organismal Biology, Ohio Biodiversity Conservation Partnership, Ohio State University, Columbus, Ohio, USA
| |
Collapse
|
46
|
Sumreddee P, Hay EH, Toghiani S, Roberts A, Aggrey SE, Rekaya R. Grid search approach to discriminate between old and recent inbreeding using phenotypic, pedigree and genomic information. BMC Genomics 2021; 22:538. [PMID: 34256689 PMCID: PMC8278650 DOI: 10.1186/s12864-021-07872-z] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2020] [Accepted: 07/05/2021] [Indexed: 12/02/2022] Open
Abstract
Background Although inbreeding caused by the mating of animals related through a recent common ancestor is expected to have more harmful effects on phenotypes than ancient inbreeding (old inbreeding), estimating these effects requires a clear definition of recent (new) and ancient (old) inbreeding. Several methods have been proposed to classify inbreeding using pedigree and genomic data. Unfortunately, these methods are largely based on heuristic criteria such as the number of generations from a common ancestor or length of runs of homozygosity (ROH) segments. To mitigate these deficiencies, this study aimed to develop a method to classify pedigree and genomic inbreeding into recent and ancient classes based on a grid search algorithm driven by the assumption that new inbreeding tends to have a more pronounced detrimental effect on traits. The proposed method was tested using a cattle population characterized by a deep pedigree. Results Effects of recent and ancient inbreeding were assessed on four growth traits (birth, weaning and yearling weights and average daily gain). Thresholds to classify inbreeding into recent and ancient classes were trait-specific and varied across traits and sources of information. Using pedigree information, inbreeding generated in the last 10 to 11 generations was considered as recent. When genomic information (ROH) was used, thresholds ranged between four to seven generations, indicating, in part, the ability of ROH segments to characterize the harmful effects of inbreeding in shorter periods of time. Nevertheless, using the proposed classification method, the discrimination between new and old inbreeding was less robust when ROH segments were used compared to pedigree. Using several model comparison criteria, the proposed approach was generally better than existing methods. Recent inbreeding appeared to be more harmful across the growth traits analyzed. However, both new and old inbreeding were found to be associated with decreased yearling weight and average daily gain. Conclusions The proposed method provided a more objective quantitative approach for the classification of inbreeding. The proposed method detected a clear divergence in the effects of old and recent inbreeding using pedigree data and it was superior to existing methods for all analyzed traits. Using ROH data, the discrimination between old and recent inbreeding was less clear and the proposed method was superior to existing approaches for two out of the four analyzed traits. Deleterious effects of recent inbreeding were detected sooner (fewer generations) using genomic information than pedigree. Difference in the results using genomic and pedigree information could be due to the dissimilarity in the number of generations to a common ancestor. Additionally, the uncertainty associated with the identification of ROH segments and associated inbreeding could have an effect on the results. Potential biases in the estimation of inbreeding effects may occur when new and old inbreeding are discriminated based on arbitrary thresholds. To minimize the impact of inbreeding, mating designs should take the different inbreeding origins into consideration. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07872-z.
Collapse
Affiliation(s)
- Pattarapol Sumreddee
- Department of Animal and Dairy Science, The University of Georgia, Athens, GA, 30602, USA
| | - El Hamidi Hay
- USDA Agricultural Research Service, Fort Keogh Livestock and Range Research Laboratory, Miles City, MT, 59301, USA.
| | - Sajjad Toghiani
- USDA Agricultural Research Service, Beltsville Agricultural Research Center, Beltsville, MD, 20705, USA
| | - Andrew Roberts
- USDA Agricultural Research Service, Fort Keogh Livestock and Range Research Laboratory, Miles City, MT, 59301, USA
| | - Samuel E Aggrey
- Department of Poultry Science, The University of Georgia, Athens, GA, 30602, USA.,Institute of Bioinformatics, The University of Georgia, Athens, GA, 30602, USA
| | - Romdhane Rekaya
- Department of Animal and Dairy Science, The University of Georgia, Athens, GA, 30602, USA.,Institute of Bioinformatics, The University of Georgia, Athens, GA, 30602, USA.,Department of Statistics, The University of Georgia, Athens, GA, 30602, USA
| |
Collapse
|
47
|
Duntsch L, Whibley A, Brekke P, Ewen JG, Santure AW. Genomic data of different resolutions reveal consistent inbreeding estimates but contrasting homozygosity landscapes for the threatened Aotearoa New Zealand hihi. Mol Ecol 2021; 30:6006-6020. [PMID: 34242449 DOI: 10.1111/mec.16068] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2020] [Accepted: 07/01/2021] [Indexed: 12/19/2022]
Abstract
Inbreeding can lead to a loss of heterozygosity in a population and when combined with genetic drift may reduce the adaptive potential of a species. However, there is uncertainty about whether resequencing data can provide accurate and consistent inbreeding estimates. Here, we performed an in-depth inbreeding analysis for hihi (Notiomystis cincta), an endemic and nationally vulnerable passerine bird of Aotearoa New Zealand. We first focused on subsampling variants from a reference genome male, and found that low-density data sets tend to miss runs of homozygosity (ROH) in some places and overestimate ROH length in others, resulting in contrasting homozygosity landscapes. Low-coverage resequencing and 50 K SNP array densities can yield comparable inbreeding results to high-coverage resequencing approaches, but the results for all data sets are highly dependent on the software settings employed. Second, we extended our analysis to 10 hihi where low-coverage whole genome resequencing, RAD-seq and SNP array genotypes are available. We inferred ROH and individual inbreeding to evaluate the relative effects of sequencing depth versus SNP density on estimating inbreeding coefficients and found that high rates of missingness downwardly bias both the number and length of ROH. In summary, when using genomic data to evaluate inbreeding, studies must consider that ROH estimates are heavily dependent on analysis parameters, data set density and individual sequencing depth.
Collapse
Affiliation(s)
- Laura Duntsch
- School of Biological Sciences, University of Auckland, Auckland, New Zealand
| | - Annabel Whibley
- School of Biological Sciences, University of Auckland, Auckland, New Zealand
| | - Patricia Brekke
- Institute of Zoology, Zoological Society of London, London, UK
| | - John G Ewen
- Institute of Zoology, Zoological Society of London, London, UK
| | - Anna W Santure
- School of Biological Sciences, University of Auckland, Auckland, New Zealand
| |
Collapse
|
48
|
Tarekegn GM, Khayatzadeh N, Liu B, Osama S, Haile A, Rischkowsky B, Zhang W, Tesfaye K, Dessie T, Mwai OA, Djikeng A, Mwacharo JM. Ethiopian indigenous goats offer insights into past and recent demographic dynamics and local adaptation in sub-Saharan African goats. Evol Appl 2021; 14:1716-1731. [PMID: 34295359 PMCID: PMC8287980 DOI: 10.1111/eva.13118] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2020] [Revised: 07/15/2020] [Accepted: 08/13/2020] [Indexed: 12/13/2022] Open
Abstract
Knowledge on how adaptive evolution and human socio-cultural and economic interests shaped livestock genomes particularly in sub-Saharan Africa remains limited. Ethiopia is in a geographic region that has been critical in the history of African agriculture with ancient and diverse human ethnicity and bio-climatic conditions. Using 52K genome-wide data analysed in 646 individuals from 13 Ethiopian indigenous goat populations, we observed high levels of genetic variation. Although runs of homozygosity (ROH) were ubiquitous genome-wide, there were clear differences in patterns of ROH length and abundance and in effective population sizes illustrating differences in genome homozygosity, evolutionary history, and management. Phylogenetic analysis incorporating patterns of genetic differentiation and gene flow with ancestry modelling highlighted past and recent intermixing and possible two deep ancient genetic ancestries that could have been brought by humans with the first introduction of goats in Africa. We observed four strong selection signatures that were specific to Arsi-Bale and Nubian goats. These signatures overlapped genomic regions with genes associated with morphological, adaptation, reproduction and production traits due possibly to selection under environmental constraints and/or human preferences. The regions also overlapped uncharacterized genes, calling for a comprehensive annotation of the goat genome. Our results provide insights into mechanisms leading to genome variation and differentiation in sub-Saharan Africa indigenous goats.
Collapse
Affiliation(s)
- Getinet M. Tarekegn
- Department of Animal Production and TechnologySchool of Animal Sciences and Veterinary MedicineBahir Dar UniversityBahir DarEthiopia
- Department of Animal Breeding and GeneticsSwedish University of Agricultural Sciences (SLU)UppsalaSweden
| | - Negar Khayatzadeh
- Department of Sustainable Agricultural SystemsDivision of Livestock SciencesUniversity of Natural Resources and Life SciencesViennaAustria
| | - Bin Liu
- Inner Mongolia Agricultural UniversityHohhotChina
| | - Sarah Osama
- The University of QueenslandSaint LuciaQLDAustralia
| | - Aynalem Haile
- Small Ruminant GenomicsInternational Centre for Agricultural Research in the Dry Areas (ICARDA)Addis AbabaEthiopia
| | - Barbara Rischkowsky
- Small Ruminant GenomicsInternational Centre for Agricultural Research in the Dry Areas (ICARDA)Addis AbabaEthiopia
| | | | - Kassahun Tesfaye
- Department of Microbial, Cellular and Molecular BiologyAddis Ababa UniversityAddis AbabaEthiopia
| | - Tadelle Dessie
- International Livestock Research Institute (ILRI)Addis AbabaEthiopia
| | - Okeyo A. Mwai
- International Livestock Research Institute (ILRI)NairobiKenya
| | - Appolinaire Djikeng
- Animal and Veterinary Sciences Group, SRUC and Centre for Tropical Livestock Genetics and Health (CTLGH)The Roslin InstituteEaster BushMidlothianUK
| | - Joram M. Mwacharo
- Small Ruminant GenomicsInternational Centre for Agricultural Research in the Dry Areas (ICARDA)Addis AbabaEthiopia
- Animal and Veterinary Sciences Group, SRUC and Centre for Tropical Livestock Genetics and Health (CTLGH)The Roslin InstituteEaster BushMidlothianUK
| |
Collapse
|
49
|
Druet T, Oleński K, Flori L, Bertrand AR, Olech W, Tokarska M, Kaminski S, Gautier M. Genomic Footprints of Recovery in the European Bison. J Hered 2021; 111:194-203. [PMID: 32027367 DOI: 10.1093/jhered/esaa002] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/05/2020] [Indexed: 02/04/2023] Open
Abstract
After extinction in the wild in the beginning of the 20th century, the European bison has been successfully recovered in 2 distinct genetic lines from only 12 and 7 captive founders. We here aimed at characterizing the levels of realized inbreeding in these 2 restored lines to provide empirical insights into the genomic footprints left by population recovery from a small number of founders. To that end, we genotyped 183 European bison born over the last 40 years with the Illumina BovineHD beadchip that contained 22 602 informative autosomal single-nucleotide polymorphisms after data filtering. We then identified homozygous-by-descent (HBD) segments and classified them into different age-related classes relying on a model-based approach. As expected, we observed that the strong and recent founder effect experienced by the 2 lines resulted in very high levels of recent inbreeding and in the presence of long HBD tracks (up to 120 Mb). These long HBD tracks were associated with ancestors living approximately from 4 to 32 generations in the past, suggesting that inbreeding accumulated over multiple generations after the bottleneck. The contribution to inbreeding of the most recent groups of ancestors was however found to be decreasing in both lines. In addition, comparison of Lowland individuals born at different time periods showed that the levels of inbreeding tended to stabilize, HBD segments being shorter in animals born more recently which indicates efficient control of inbreeding. Monitoring HBD segment lengths over generations may thus be viewed as a valuable genomic diagnostic tool for populations in conservation or recovery programs.
Collapse
Affiliation(s)
- Tom Druet
- Unit of Animal Genomics, GIGA-R and Faculty of Veterinary Medicine, University of Liège, Liège, Belgium
| | - Kamil Oleński
- Department of Animal Genetics, University of Warmia and Mazury in Olsztyn, Olsztyn, Poland
| | - Laurence Flori
- SELMET, INRAE, CIRAD, Montpellier Supagro, University of Montpellier, Montpellier, France
| | - Amandine R Bertrand
- Unit of Animal Genomics, GIGA-R and Faculty of Veterinary Medicine, University of Liège, Liège, Belgium
| | - Wanda Olech
- Department of Animal Genetics and Breeding, Warsaw University of Life Sciences-SGGW, Warsaw, Poland
| | | | - Stanislaw Kaminski
- Department of Animal Genetics, University of Warmia and Mazury in Olsztyn, Olsztyn, Poland
| | - Mathieu Gautier
- INRAE, UMR CBGP (INRAE-IRD-Cirad-Montpellier SupAgro), Montferrier-sur-Lez, France
| |
Collapse
|
50
|
The Nero Lucano Pig Breed: Recovery and Variability. Animals (Basel) 2021; 11:ani11051331. [PMID: 34067067 PMCID: PMC8150585 DOI: 10.3390/ani11051331] [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: 02/16/2021] [Revised: 05/01/2021] [Accepted: 05/04/2021] [Indexed: 01/14/2023] Open
Abstract
Simple Summary The reduction of biodiversity determines the loss of species and breeds, with the consequent disappearance of production systems, knowledge, cultures and local traditions. The Nero Lucano pig is a native breed of Southern Italy (Basilicata region) recovered, starting from 2001, because of the high quality of its cured meat products. This study gives a picture of the low genetic variability of this breed. Knowledge of individual inbreeding levels allows for planning of interventions to reduce the negative effects of the low effective population size and, then, improve the efficiency of the actual recovery project. Abstract The Nero Lucano (NL) pig is a black coat colored breed characterized by a remarkable ability to adapt to the difficult territory and climatic conditions of Basilicata region in Southern Italy. In the second half of the twentieth century, technological innovation, agricultural evolution, new breeding methods and the demand for increasingly lean meat brought the breed almost to extinction. Only in 2001, thanks to local institutions such as: the Basilicata Region, the University of Basilicata, the Regional Breeders Association and the Medio Basento mountain community, the NL pig returned to populate the area with the consequent possibility to appreciate again its specific cured meat products. We analyzed the pedigrees recorded by the breeders and the Illumina Porcine SNP60 BeadChip genotypes in order to obtain the genetic structure of the NL pig. Results evidenced that this population is characterized by long mean generation intervals (up to 3.5 yr), low effective population size (down to 7.2) and high mean inbreeding coefficients (FMOL = 0.53, FROH = 0.39). This picture highlights the low level of genetic variability and the critical issues to be faced for the complete recovery of this population.
Collapse
|