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
|
Giovannini S, Chessari G, Riggio S, Marletta D, Sardina MT, Mastrangelo S, Sarti FM. Insight into the current genomic diversity, conservation status and population structure of Tunisian Barbarine sheep breed. Front Genet 2024; 15:1379086. [PMID: 38881792 PMCID: PMC11176520 DOI: 10.3389/fgene.2024.1379086] [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: 01/30/2024] [Accepted: 05/07/2024] [Indexed: 06/18/2024] Open
Abstract
Local livestock breeds play a crucial role in global biodiversity, connecting natural and human-influenced environments and contributing significantly to ecosystem services. While commercial breeds dominate industrial systems, local livestock breeds in developing countries, like Barbarine sheep in Tunisia, are vital for food security and community maintenance. The Tunisian Barbarine sheep, known for its adaptability and distinctive fat-tailed morphology, faces challenges due to historical crossbreeding. In this study, the Illumina Ovine SNP50K BeadChip array was used to perform a genome-wide characterization of Tunisian Barbarine sheep to investigate its genetic diversity, the genome structure, and the relationship within the context of Mediterranean breeds. The results show moderate genetic diversity and low inbreeding. Runs of Homozygosity analysis find genomic regions linked to important traits, including fat tail characteristics. Genomic relationship analysis shows proximity to Algerian thin-tailed breeds, suggesting crossbreeding impacts. Admixture analysis reveals unique genetic patterns, emphasizing the Tunisian Barbarine's identity within the Mediterranean context and its closeness to African breeds. Current results represent a starting point for the creation of monitoring and conservation plans. In summary, despite genetic dilution due to crossbreeding, the identification of genomic regions offers crucial insights for conservation. The study confirms the importance of preserving unique genetic characteristics of local breeds, particularly in the face of ongoing crossbreeding practices and environmental challenges. These findings contribute valuable insights for the sustainable management of this unique genetic reservoir, supporting local economies and preserving sheep species biodiversity.
Collapse
Affiliation(s)
- Samira Giovannini
- Dipartimento di Scienze Agrarie, Alimentari e Ambientali, University of Perugia, Perugia, Italy
| | - Giorgio Chessari
- Dipartimento Agricoltura, Alimentazione e Ambiente, University of Catania, Catania, Italy
| | - Silvia Riggio
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Donata Marletta
- Dipartimento Agricoltura, Alimentazione e Ambiente, University of Catania, Catania, 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
| | - Francesca Maria Sarti
- Dipartimento di Scienze Agrarie, Alimentari e Ambientali, University of Perugia, Perugia, Italy
| |
Collapse
|
3
|
Tan X, Liu L, Dong J, Huang M, Zhang J, Li Q, Wang H, Bai L, Cui M, Zhou Z, Wu D, Xiang Y, Li W, Wang D. Genome-wide detections for runs of homozygosity and selective signatures reveal novel candidate genes under domestication in chickens. BMC Genomics 2024; 25:485. [PMID: 38755540 PMCID: PMC11097469 DOI: 10.1186/s12864-024-10349-4] [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: 03/01/2024] [Accepted: 04/25/2024] [Indexed: 05/18/2024] Open
Abstract
BACKGROUND Indigenous chickens were developed through a combination of natural and artificial selection; essentially, changes in genomes led to the formation of these modern breeds via admixture events. However, their confusing genetic backgrounds include a genomic footprint regulating complex traits, which is not conducive to modern animal breeding. RESULTS To better evaluate the candidate regions under domestication in indigenous chickens, we considered both runs of homozygosity (ROHs) and selective signatures in 13 indigenous chickens. The genomes of Silkie feather chickens presented the highest heterozygosity, whereas the highest inbreeding status and ROH number were found in Luhua chickens. Short ROH (< 1 Mb), were the principal type in all chickens. A total of 291 ROH islands were detected, and QTLdb mapping results indicated that body weight and carcass traits were the most important traits. An ROH on chromosome 2 covering VSTM2A gene was detected in 12 populations. Combined analysis with the Tajima's D index revealed that 18 genes (e.g., VSTM2A, BBOX1, and RYR2) were under selection and covered by ROH islands. Transcriptional analysis results showed that RYR2 and BBOX1 were specifically expressed in the heart and muscle tissue, respectively. CONCLUSION Based on genome-wide scanning for ROH and selective signatures, we evaluated the genomic characteristics and detected significant candidate genes covered by ROH islands and selective signatures. The findings in this study facilitated the understanding of genetic diversity and provided valuable insights for chicken breeding and conservation strategies.
Collapse
Affiliation(s)
- Xiaodong Tan
- Institute of Animal Husbandry and Veterinary Science, Zhejiang Academy of Agricultural Sciences, Hangzhou, 310021, China
| | - Lu Liu
- College of Animal Sciences, Zhejiang University, Hangzhou, Zhejiang, 310058, China
- Jinhua Jinfan Feed Co., Ltd, Jinhua, Zhejiang, 321000, China
| | - Jie Dong
- Institute of Animal Husbandry and Veterinary Science, Zhejiang Academy of Agricultural Sciences, Hangzhou, 310021, China
| | - Minjie Huang
- Institute of Animal Husbandry and Veterinary Science, Zhejiang Academy of Agricultural Sciences, Hangzhou, 310021, China
| | - Jiawen Zhang
- Institute of Animal Husbandry and Veterinary Science, Zhejiang Academy of Agricultural Sciences, Hangzhou, 310021, China
| | - Qinghai Li
- Animal Husbandry Institute, Hangzhou Academy of Agricultural Sciences, Hangzhou, 310024, China
| | - Huanhuan Wang
- Animal Husbandry Institute, Hangzhou Academy of Agricultural Sciences, Hangzhou, 310024, China
| | - Lijuan Bai
- Zhejiang Animal Husbandry Technology Extension and Breeding Livestock and Poultry Monitoring Station, Hangzhou, 310020, China
| | - Ming Cui
- Zhejiang Animal Husbandry Technology Extension and Breeding Livestock and Poultry Monitoring Station, Hangzhou, 310020, China
| | - Zhenzhen Zhou
- Institute of Animal Husbandry and Veterinary Science, Zhejiang Academy of Agricultural Sciences, Hangzhou, 310021, China
| | - De Wu
- Postdoctoral Research Station, Jinhua Development Zone, Jinhua, Zhejiang, 321000, China
| | - Yun Xiang
- Jinhua Jinfan Feed Co., Ltd, Jinhua, Zhejiang, 321000, China.
| | - Weifen Li
- College of Animal Sciences, Zhejiang University, Hangzhou, Zhejiang, 310058, China.
| | - Deqian Wang
- Institute of Animal Husbandry and Veterinary Science, Zhejiang Academy of Agricultural Sciences, Hangzhou, 310021, China.
| |
Collapse
|
4
|
Illa SK, Mumtaz S, Nath S, Mukherjee S, Mukherjee A. Characterization of runs of Homozygosity revealed genomic inbreeding and patterns of selection in indigenous sahiwal cattle. J Appl Genet 2024; 65:167-180. [PMID: 38110827 DOI: 10.1007/s13353-023-00816-1] [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: 02/10/2023] [Revised: 11/29/2023] [Accepted: 12/05/2023] [Indexed: 12/20/2023]
Abstract
Runs of homozygosity (ROH) are contiguous genomic regions, homozygous across all sites which arise in an individual due to the parents transmitting identical haplotypes to their offspring. The genetic improvement program of Sahiwal cattle after decades of selection needs re-assessment of breeding strategy and population phenomena. Hence, the present study was carried out to optimize input parameters in PLINK for ROH estimates, to explore ROH islands and assessment of pedigree and genome-based inbreeding in Sahiwal cattle. The sliding window approach with parameters standardized to define ROH for the specific population under study was used for the identification of runs. The optimum maximum gap, density, window-snp and window-threshold were 250 Kb, 120 Kb/SNP, 10, 0.05 respectively and ROH patterns were also characterized. ROH islands were defined as the short homozygous genomic regions shared by a large proportion of individuals in a population, containing significantly higher occurrences of ROH than the population specific threshold level. These were identified using the -homozyg-group function of the PLINK v1.9 program. Our results indicated that the Islands of ROH harbor a few candidate genes, ACAD11, RFX4, BANP, UBA5 that are associated with major economic traits. The average FPED (Pedigree based inbreeding coefficient), FROH (Genomic inbreeding coefficient), FHOM (Inbreeding estimated as the ratio of observed and expected homozygous genotypes), FGRM (Inbreeding estimated on genomic relationship method) and FGRM0.5 (Inbreeding estimated from the diagonal of a GRM with allele frequencies near to 0.5) were 0.009, 0.091, 0.035, -0.104 and -0.009, respectively. Our study revealed the optimum parameter setting in PLINK viz. maximal gaps between two SNPs, minimal density of SNPs in a segment (in kb/SNP) and scanning window size to identify ROH segments, which will enable ROH estimation more efficient and comparable across various SNP genotyping-based studies. The result further emphasized the significant role of genomics in unraveling population diversity, selection signatures and inbreeding in the ongoing Sahiwal breed improvement programs.
Collapse
Affiliation(s)
- Satish Kumar Illa
- Livestock Research Station, Garividi, Sri Venkateswara Veterinary University, Tirupati, Andhra Pradesh State, India
| | - Shabahat Mumtaz
- Animal Husbandry Department, Kolkata, West Bengal State, India
| | - Sapna Nath
- College of Veterinary Science, Garividi, Sri Venkateswara Veterinary University, Tirupati, Andhra Pradesh State, India
| | - Sabyasachi Mukherjee
- Animal Genetics & Breeding Division, Indian Council of Agricultural Research (ICAR)-National Dairy Research Institute (NDRI), Karnal, Haryana State, India.
| | - Anupama Mukherjee
- Animal Genetics & Breeding Division, Indian Council of Agricultural Research (ICAR)-National Dairy Research Institute (NDRI), Karnal, Haryana State, India.
| |
Collapse
|
5
|
Sun X, Guo J, Li R, Zhang H, Zhang Y, Liu GE, Emu Q, Zhang H. Whole-Genome Resequencing Reveals Genetic Diversity and Wool Trait-Related Genes in Liangshan Semi-Fine-Wool Sheep. Animals (Basel) 2024; 14:444. [PMID: 38338087 PMCID: PMC10854784 DOI: 10.3390/ani14030444] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2023] [Revised: 01/12/2024] [Accepted: 01/25/2024] [Indexed: 02/12/2024] Open
Abstract
Understanding the genetic makeup of local sheep breeds is essential for their scientific conservation and sustainable utilization. The Liangshan semi-fine-wool sheep (LSS), a Chinese semi-fine-wool breed renowned for its soft wool, was analyzed using whole-genome sequencing data including 35 LSS, 84 sheep from other domestic breeds, and 20 Asiatic mouflons. We investigated the genetic composition of LSS by conducting analyses of the population structure, runs of homozygosity, genomic inbreeding coefficients, and selection signature. Our findings indicated that LSS shares greater genetic similarity with Border Leicester and Romney sheep than with Tibetan (TIB), Yunnan (YNS), and Chinese Merino sheep. Genomic analysis indicated low to moderate inbreeding coefficients, ranging from 0.014 to 0.154. In identifying selection signals across the LSS genome, we pinpointed 195 candidate regions housing 74 annotated genes (e.g., IRF2BP2, BVES, and ALOX5). We also found the overlaps between the candidate regions and several known quantitative trait loci related to wool traits, such as the wool staple length and wool fiber diameter. A selective sweep region, marked by the highest value of cross-population extended haplotype homozygosity, encompassed IRF2BP2-an influential candidate gene affecting fleece fiber traits. Furthermore, notable differences in genotype frequency at a mutation site (c.1051 + 46T > C, Chr25: 6,784,190 bp) within IRF2BP2 were observed between LSS and TIB and YNS sheep (Fisher's exact test, p < 2.2 × 10-16). Taken together, these findings offer insights crucial for the conservation and breeding enhancement of LSS.
Collapse
Affiliation(s)
- Xueliang Sun
- 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; (X.S.); (J.G.)
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, 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; (X.S.); (J.G.)
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 611130, China
| | - Ran Li
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
| | - Huanhuan Zhang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
| | - Yifei 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; (X.S.); (J.G.)
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 611130, China
| | - George E. Liu
- Animal Genomics and Improvement Laboratory, Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, MD 20705, USA
| | - Quzhe Emu
- Animal Genetics and Breeding Key Laboratory of Sichuan Province, Sichuan Animal Science Academy, No. 7, Niusha Road, Chengdu 610066, 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; (X.S.); (J.G.)
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 611130, China
| |
Collapse
|
6
|
Saravanan KA, Rajawat D, Kumar H, Nayak SS, Bhushan B, Dutt T, Panigrahi M. Signatures of selection in riverine buffalo populations revealed by genome-wide SNP data. Anim Biotechnol 2023; 34:3343-3354. [PMID: 36384399 DOI: 10.1080/10495398.2022.2145292] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The detection of selection signatures assists in understanding domestication, evolution, and the identification of genomic regions related to adaptation and production traits in buffaloes. The emergence of high-throughput technologies like Next Generation Sequencing and SNP genotyping had expanded our ability to detect these signatures of selection. In this study, we sought to identify signatures of selection in five buffalo populations (Brazilian Murrah, Bulgarian Murrah, Indian Murrah, Nili-Ravi, and Kundi) using Axiom Buffalo 90 K Genotyping Array data. Using seven different methodologies (Tajima's D, CLR, ROH, iHS, FST, FLK and hapFLK), we identified selection signatures in 374 genomic regions, spanning a total of 381 genes and 350 quantitative trait loci (QTLs). Among these, several candidate genes were associated with QTLs for milk production, reproduction, growth and carcass traits. The genes and QTLs reported in this study provide insight into selection signals shaping the genome of buffalo breeds. Our findings can aid in further genomic association studies, genomic prediction, and the implementation of breeding programmes in Indian buffaloes.
Collapse
Affiliation(s)
- K A Saravanan
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, India
| | - Divya Rajawat
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, India
| | - Harshit Kumar
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, India
| | - Sonali Sonejita Nayak
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, India
| | - Bharat Bhushan
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, India
| | - Triveni Dutt
- Livestock Production and Management Section, Indian Veterinary Research Institute, Izatnagar, Bareilly, India
| | - Manjit Panigrahi
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, India
| |
Collapse
|
7
|
Ghoreishifar M, Vahedi SM, Salek Ardestani S, Khansefid M, Pryce JE. Genome-wide assessment and mapping of inbreeding depression identifies candidate genes associated with semen traits in Holstein bulls. BMC Genomics 2023; 24:230. [PMID: 37138201 PMCID: PMC10157977 DOI: 10.1186/s12864-023-09298-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2023] [Accepted: 04/05/2023] [Indexed: 05/05/2023] Open
Abstract
BACKGROUND The reduction in phenotypic performance of a population due to mating between close relatives is called inbreeding depression. The genetic background of inbreeding depression for semen traits is poorly understood. Thus, the objectives were to estimate the effect of inbreeding and to identify genomic regions underlying inbreeding depression of semen traits including ejaculate volume (EV), sperm concentration (SC), and sperm motility (SM). The dataset comprised ~ 330 K semen records from ~ 1.5 K Holstein bulls genotyped with 50 K single nucleotide polymorphism (SNP) BeadChip. Genomic inbreeding coefficients were estimated using runs of homozygosity (i.e., FROH > 1 Mb) and excess of SNP homozygosity (FSNP). The effect of inbreeding was estimated by regressing phenotypes of semen traits on inbreeding coefficients. Associated variants with inbreeding depression were also detected by regressing phenotypes on ROH state of the variants. RESULTS Significant inbreeding depression was observed for SC and SM (p < 0.01). A 1% increase in FROH reduced SM and SC by 0.28% and 0.42% of the population mean, respectively. By splitting FROH into different lengths, we found significant reduction in SC and SM due to longer ROH, which is indicative of more recent inbreeding. A genome-wide association study revealed two signals positioned on BTA 8 associated with inbreeding depression of SC (p < 0.00001; FDR < 0.02). Three candidate genes of GALNTL6, HMGB2, and ADAM29, located in these regions, have established and conserved connections with reproduction and/or male fertility. Moreover, six genomic regions on BTA 3, 9, 21 and 28 were associated with SM (p < 0.0001; FDR < 0.08). These genomic regions contained genes including PRMT6, SCAPER, EDC3, and LIN28B with established connections to spermatogenesis or fertility. CONCLUSIONS Inbreeding depression adversely affects SC and SM, with evidence that longer ROH, or more recent inbreeding, being especially detrimental. There are genomic regions associated with semen traits that seems to be especially sensitive to homozygosity, and evidence to support some from other studies. Breeding companies may wish to consider avoiding homozygosity in these regions for potential artificial insemination sires.
Collapse
Affiliation(s)
- Mohammad Ghoreishifar
- Agriculture Victoria Research, AgriBio, Centre for AgriBioscience, 5 Ring Road, Bundoora, Victoria, 3083, Australia.
- School of Applied Systems Biology, La Trobe University, Bundoora, Victoria, 3083, Australia.
| | - Seyed Milad Vahedi
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, B2N5E3, Canada
| | | | - Majid Khansefid
- Agriculture Victoria Research, AgriBio, Centre for AgriBioscience, 5 Ring Road, Bundoora, Victoria, 3083, Australia
- School of Applied Systems Biology, La Trobe University, Bundoora, Victoria, 3083, Australia
| | - Jennie E Pryce
- Agriculture Victoria Research, AgriBio, Centre for AgriBioscience, 5 Ring Road, Bundoora, Victoria, 3083, Australia
- School of Applied Systems Biology, La Trobe University, Bundoora, Victoria, 3083, Australia
| |
Collapse
|
8
|
Ceccobelli S, Landi V, Senczuk G, Mastrangelo S, Sardina MT, Ben-Jemaa S, Persichilli C, Karsli T, Bâlteanu VA, Raschia MA, Poli MA, Ciappesoni G, Muchadeyi FC, Dzomba EF, Kunene NW, Lühken G, Deniskova TE, Dotsev AV, Zinovieva NA, Zsolnai A, Anton I, Kusza S, Carolino N, Santos-Silva F, Kawęcka A, Świątek M, Niżnikowski R, Špehar M, Anaya G, Granero A, Perloiro T, Cardoso P, Grande S, de Los Santos BL, Danchin-Burge C, Pasquini M, Martínez Martínez A, Delgado Bermejo JV, Lasagna E, Ciani E, Sarti FM, Pilla F. A comprehensive analysis of the genetic diversity and environmental adaptability in worldwide Merino and Merino-derived sheep breeds. Genet Sel Evol 2023; 55:24. [PMID: 37013467 PMCID: PMC10069132 DOI: 10.1186/s12711-023-00797-z] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2022] [Accepted: 03/24/2023] [Indexed: 04/05/2023] Open
Abstract
BACKGROUND To enhance and extend the knowledge about the global historical and phylogenetic relationships between Merino and Merino-derived breeds, 19 populations were genotyped with the OvineSNP50 BeadChip specifically for this study, while an additional 23 populations from the publicly available genotypes were retrieved. Three complementary statistical tests, Rsb (extended haplotype homozygosity between-populations), XP-EHH (cross-population extended haplotype homozygosity), and runs of homozygosity (ROH) islands were applied to identify genomic variants with potential impact on the adaptability of Merino genetic type in two contrasting climate zones. RESULTS The results indicate that a large part of the Merino's genetic relatedness and admixture patterns are explained by their genetic background and/or geographic origin, followed by local admixture. Multi-dimensional scaling, Neighbor-Net, Admixture, and TREEMIX analyses consistently provided evidence of the role of Australian, Rambouillet and German strains in the extensive gene introgression into the other Merino and Merino-derived breeds. The close relationship between Iberian Merinos and other South-western European breeds is consistent with the Iberian origin of the Merino genetic type, with traces from previous contributions of other Mediterranean stocks. Using Rsb and XP-EHH approaches, signatures of selection were detected spanning four genomic regions located on Ovis aries chromosomes (OAR) 1, 6 and 16, whereas two genomic regions on OAR6, that partially overlapped with the previous ones, were highlighted by ROH islands. Overall, the three approaches identified 106 candidate genes putatively under selection. Among them, genes related to immune response were identified via the gene interaction network. In addition, several candidate genes were found, such as LEKR1, LCORL, GHR, RBPJ, BMPR1B, PPARGC1A, and PRKAA1, related to morphological, growth and reproductive traits, adaptive thermogenesis, and hypoxia responses. CONCLUSIONS To the best of our knowledge, this is the first comprehensive dataset that includes most of the Merino and Merino-derived sheep breeds raised in different regions of the world. The results provide an in-depth picture of the genetic makeup of the current Merino and Merino-derived breeds, highlighting the possible selection pressures associated with the combined effect of anthropic and environmental factors. The study underlines the importance of Merino genetic types as invaluable resources of possible adaptive diversity in the context of the occurring climate changes.
Collapse
Affiliation(s)
- Simone Ceccobelli
- Department of Agricultural, Food and Environmental Sciences, Università Politecnica delle Marche, 60131, Ancona, Italy.
| | - Vincenzo Landi
- Department of Veterinary Medicine, University of Bari ''Aldo Moro", 70010, Valenzano, Italy
| | - Gabriele Senczuk
- Department of Agricultural, Environmental and Food Sciences, University of Molise, 86100, Campobasso, Italy
| | - Salvatore Mastrangelo
- Department of Agricultural, Food and Forest Sciences, University of Palermo, 90128, Palermo, Italy
| | - Maria Teresa Sardina
- Department of Agricultural, Food and Forest Sciences, University of Palermo, 90128, Palermo, Italy
| | - Slim Ben-Jemaa
- Laboratoire des Productions Animales et Fourragères, Institut National de la Recherche Agronomique de Tunisie, Université de Carthage, 2049, Ariana, Tunisia
| | - Christian Persichilli
- Department of Agricultural, Environmental and Food Sciences, University of Molise, 86100, Campobasso, Italy
| | - Taki Karsli
- Department of Animal Science, Faculty of Agriculture, Eskisehir Osmangazi University, 26040, Eskisehir, Turkey
| | - Valentin-Adrian Bâlteanu
- Laboratory of Genomics, Biodiversity, Animal Breeding and Molecular Pathology, Institute of Life Sciences, University of Agricultural Sciences and Veterinary Medicine of Cluj-Napoca, 400372, Cluj-Napoca, Romania
| | - María Agustina Raschia
- Instituto de Genética "Ewald A. Favret", Instituto Nacional de Tecnología Agropecuaria, CICVyA-CNIA, B1686, Hurlingham, Buenos Aires, Argentina
| | - Mario Andrés Poli
- Instituto de Genética "Ewald A. Favret", Instituto Nacional de Tecnología Agropecuaria, CICVyA-CNIA, B1686, Hurlingham, Buenos Aires, Argentina
| | - Gabriel Ciappesoni
- Instituto Nacional de Investigación Agropecuaria, 90200, Canelones, Uruguay
| | | | - Edgar Farai Dzomba
- Discipline of Genetics, School of Life Sciences, University of KwaZulu-Natal, 3209, Scottsville, Pietermaritzburg, South Africa
| | | | - Gesine Lühken
- Institute of Animal Breeding and Genetics, Justus Liebig University, 35390, Giessen, Germany
| | | | | | | | - Attila Zsolnai
- Department of Animal Breeding, Institute of Animal Science, Hungarian University of Agriculture and Life Sciences, Kaposvár Campus, 2053, Herceghalom, Hungary
| | - István Anton
- Department of Animal Breeding, Institute of Animal Science, Hungarian University of Agriculture and Life Sciences, Kaposvár Campus, 2053, Herceghalom, Hungary
| | - Szilvia Kusza
- Centre for Agricultural Genomics and Biotechnology, Faculty of Agricultural and Food Sciences and Environmental Management, University of Debrecen, 4032, Debrecen, Hungary
| | - Nuno Carolino
- Instituto Nacional de Investigação Agrária e Veterinária, 2005-048, Vale de Santarém, Portugal
| | - Fátima Santos-Silva
- Instituto Nacional de Investigação Agrária e Veterinária, 2005-048, Vale de Santarém, Portugal
| | - Aldona Kawęcka
- Department of Sheep and Goat Breeding, National Research Institute of Animal Production, 32-083, Kraków, Poland
| | - Marcin Świątek
- Department of Animal Breeding, Institute of Animal Sciences, Warsaw University of Life Sciences-SGGW, 02-786, Warsaw, Poland
| | - Roman Niżnikowski
- Department of Animal Breeding, Institute of Animal Sciences, Warsaw University of Life Sciences-SGGW, 02-786, Warsaw, Poland
| | - Marija Špehar
- Croatian Agency for Agriculture and Food, 10000, Zagreb, Croatia
| | - Gabriel Anaya
- MERAGEM Group, Department of Genetics, University of Córdoba, 14071, Córdoba, Spain
| | - Antonio Granero
- Asociación Nacional de Criadores de Ganado Merino (ACME), 28028, Madrid, Spain
| | - Tiago Perloiro
- Associação Nacional de Criadores de Ovinos da Raça Merina (ANCORME), 7005-665, Évora, Portugal
| | - Pedro Cardoso
- Associação de Produtores Agropecuários (OVIBEIRA), 6000-244, Castelo Branco, Portugal
| | - Silverio Grande
- Associazione Nazionale della Pastorizia (ASSONAPA), 00187, Rome, Italy
| | | | | | - Marina Pasquini
- Department of Agricultural, Food and Environmental Sciences, Università Politecnica delle Marche, 60131, Ancona, Italy
| | | | | | - Emiliano Lasagna
- Department of Agricultural, Food and Environmental Sciences, University of Perugia, 06121, Perugia, Italy
| | - Elena Ciani
- Department of Bioscience, Biotechnology and Biopharmaceutics, University of Bari "Aldo Moro", 70124, Bari, Italy
| | - Francesca Maria Sarti
- Department of Agricultural, Food and Environmental Sciences, University of Perugia, 06121, Perugia, Italy
| | - Fabio Pilla
- Department of Agricultural, Environmental and Food Sciences, University of Molise, 86100, Campobasso, Italy
| |
Collapse
|
9
|
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
|
10
|
Machová K, Marina H, Arranz JJ, Pelayo R, Rychtářová J, Milerski M, Vostrý L, Suárez-Vega A. Genetic diversity of two native sheep breeds by genome-wide analysis of single nucleotide polymorphisms. Animal 2023; 17:100690. [PMID: 36566708 DOI: 10.1016/j.animal.2022.100690] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2022] [Revised: 11/22/2022] [Accepted: 11/24/2022] [Indexed: 12/05/2022] Open
Abstract
Wallachian and Sumava sheep are autochthonous breeds that have undergone a significant bottleneck effect and subsequent restoration efforts. The first objective of this study was to evaluate the degree of genetic variability of both breeds and, therefore, the current management of the breeding. The second was to determine whether these two breeds still retain their genetic uniqueness in relation to each other and other breeds, despite regenerative interventions. Our data consisted of 48 individuals of Sumava and 37 individuals of Wallachian sheep. The comparison data contained 25 other breeds (primarily European) from the HapMap dataset generated by the International Sheep Genomics Consortium. When comparing all 27 breeds, the Czech breeds clustered with 15 other breeds and formed a single branch with them according to Nei's distances. At the same time, however, the clusters of both breeds were integral and easily distinguishable from the others when displayed with principal component analysis (PCA). Population substructure analysis did not show any common genetic ancestry of the Czech national breeds and breeds used for regeneration or, eventually, breeds whose ancestral population was used for regeneration. The average values of FST were higher in Wallachian sheep (FST = 0.14) than in Sumava sheep (FST = 0.08). The linkage disequilibrium (LD) extension per autosome was higher in Wallachian than in Sumava sheep. Consequently, the Ne estimates five generations ago were 68 for Sumava versus 34 for Wallachian sheep. Both native Czech breeds exhibit a wide range of inbreeding based on the excess of homozygosity (FHOM) among individuals, from -0.04 to 0.16 in Sumava and from -0.13 to 0.12 in Wallachian. Average inbreeding based on runs of homozygosity was 0.21 in Sumava and 0.27 in Wallachian. Most detected runs of homozygosity (ROH) were less than 5 Mb long for both breeds. ROH segments longer than 15 Mb were absent in Wallachian sheep. Concerning putative selection signatures, a total of 471 candidate genes in Wallachian sheep within 11 hotspots and 653 genes within 13 hotspots in Sumava sheep were identified. Czech breeds appear to be well differentiated from each other and other European breeds. Their genetic diversity is low, especially in the case of the Wallachian breed. Sumava is not so threatened by low diversity but has a larger share of the non-native gene pool.
Collapse
Affiliation(s)
- Karolína Machová
- Department of Genetics and Breeding, Czech University of Life Sciences (CZU), Prague, Czech Republic, Faculty of Agrobiology, Food and Natural Resources, Kamýcká 129, 165 00 Praha, Suchdol, Czech Republic.
| | - Héctor Marina
- Department of Animal Production, University of Leon, Leon, Spain, Veterinary Faculty, Campus de Vegazana, s/n, 24071 Leon, Spain
| | - Juan Jose Arranz
- Department of Animal Production, University of Leon, Leon, Spain, Veterinary Faculty, Campus de Vegazana, s/n, 24071 Leon, Spain
| | - Rocío Pelayo
- Department of Animal Production, University of Leon, Leon, Spain, Veterinary Faculty, Campus de Vegazana, s/n, 24071 Leon, Spain
| | - Jana Rychtářová
- Institute of Animal Science, Prague, Czech republic, Přátelství 815, 104 00 Praha, Uhříněves, Czech Republic
| | - Michal Milerski
- Institute of Animal Science, Prague, Czech republic, Přátelství 815, 104 00 Praha, Uhříněves, Czech Republic
| | - Luboš Vostrý
- Department of Genetics and Breeding, Czech University of Life Sciences (CZU), Prague, Czech Republic, Faculty of Agrobiology, Food and Natural Resources, Kamýcká 129, 165 00 Praha, Suchdol, Czech Republic
| | - Aroa Suárez-Vega
- Department of Animal Production, University of Leon, Leon, Spain, Veterinary Faculty, Campus de Vegazana, s/n, 24071 Leon, Spain
| |
Collapse
|
11
|
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
|
12
|
Population genomics reveals moderate genetic differentiation between populations of endangered Forest Musk Deer located in Shaanxi and Sichuan. BMC Genomics 2022; 23:668. [PMID: 36138352 PMCID: PMC9503231 DOI: 10.1186/s12864-022-08896-9] [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: 12/15/2021] [Accepted: 09/12/2022] [Indexed: 11/30/2022] Open
Abstract
Background Many endangered species exist in small, genetically depauperate, or inbred populations, hence promoting genetic differentiation and reducing long-term population viability. Forest Musk Deer (Moschus berezovskii) has been subject to illegal hunting for hundreds of years due to the medical and commercial values of musk, resulting in a significant decline in population size. However, it is still unclear to what extent the genetic exchange and inbreeding levels are between geographically isolated populations. By using whole-genome data, we reconstructed the demographic history, evaluated genetic diversity, and characterized the population genetic structure of Forest Musk Deer from one wild population in Sichuan Province and two captive populations from two ex-situ centers in Shaanxi Province. Results SNP calling by GATK resulted in a total of 44,008,662 SNPs. Principal component analysis (PCA), phylogenetic tree (NJ tree), ancestral component analysis (ADMIXTURE) and the ABBA-BABA test separated Sichuan and Shaanxi Forest Musk Deer as two genetic clusters, but no obvious genetic differentiation was observed between the two captive populations. The average pairwise FST value between the populations in Sichuan and Shaanxi ranged from 0.05–0.07, suggesting a low to moderate genetic differentiation. The mean heterozygous SNPs rate was 0.14% (0.11%—0.15%) for Forest Musk Deer at the genomic scale, and varied significantly among three populations (Chi-square = 1.22, p < 0.05, Kruskal–Wallis Test), with the Sichuan population having the lowest (0.11%). The nucleotide diversity of three populations varied significantly (p < 0.05, Kruskal–Wallis Test), with the Sichuan population having the lowest genetic θπ (1.69 × 10–3). Conclusions Genetic diversity of Forest Musk Deer was moderate at the genomic scale compared with other endangered species. Genetic differentiation between populations in Sichuan and Shaanxi may not only result from historical biogeographical factors but also be associated with contemporary human disturbances. Our findings provide scientific aid for the conservation and management of Forest Musk Deer. They can extend the proposed measures at the genomic level to apply to other musk deer species worldwide. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08896-9.
Collapse
|
13
|
Mulim HA, Brito LF, Batista Pinto LF, Moletta JL, Da Silva LR, Pedrosa VB. Genetic and Genomic Characterization of a New Beef Cattle Composite Breed (Purunã) Developed for Production in Pasture-Based Systems. Front Genet 2022; 13:858970. [PMID: 35923708 PMCID: PMC9341487 DOI: 10.3389/fgene.2022.858970] [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: 01/20/2022] [Accepted: 06/10/2022] [Indexed: 11/13/2022] Open
Abstract
Purunã is a composite beef cattle breed, developed in Southern Brazil by crossing the Angus, Charolais, Canchim, and Caracu breeds. The goal of this study was to perform the first genetic characterization of the Purunã breed, based on both pedigree and genomic information. For this, 100 randomly selected animals were genotyped, and 11,205 animals born from 1997 to 2019 had pedigree information. The genetic analyses performed were principal component analysis, admixture, phylogenic tree, pedigree and genomic inbreeding, linkage disequilibrium (LD), effective population size (Ne), consistency of the gametic phase, runs of homozygosity (ROH), heterozygosity-enriched regions (HERs), and functional analyses of the ROH and HER regions identified. Our findings indicate that Purunã is more genetically related to the Charolais, Canchim, and Angus breeds than Caracu or Nellore. The levels of inbreeding were shown to be small based on all the metrics evaluated and ranged from −0.009 to 0.029. A low (−0.12–0.31) correlation of the pedigree-based inbreeding compared to all the genomic inbreeding coefficients evaluated was observed. The LD average was 0.031 (±0.0517), and the consistency of the gametic phase was shown to be low for all the breed pairs, ranging from 0.42 to 0.27 to the distance of 20 Mb. The Ne values based on pedigree and genomic information were 158 and 115, respectively. A total of 1,839 ROHs were found, and the majority of them are of small length (<4 Mb). An important homozygous region was identified on BTA5 with pathways related to behavioral traits (sensory perception, detection of stimulus, and others), as well as candidate genes related to heat tolerance (MY O 1A), feed conversion rate (RDH5), and reproduction (AMDHD1). A total of 1,799 HERs were identified in the Purunã breed with 92.3% of them classified within the 0.5–1 Mb length group, and 19 HER islands were identified in the autosomal genome. These HER islands harbor genes involved in growth pathways, carcass weight (SDCBP), meat and carcass quality (MT2A), and marbling deposition (CISH). Despite the genetic relationship between Purunã and the founder breeds, a multi-breed genomic evaluation is likely not feasible due to their population structure and low consistency of the gametic phase among them.
Collapse
Affiliation(s)
| | - Luiz F. Brito
- Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
| | | | | | | | - Victor Breno Pedrosa
- Department of Animal Science, Federal University of Bahia, Salvador, Brazil
- Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
- Department of Animal Science, State University of Ponta Grossa, Ponta Grossa, Brazil
- *Correspondence: Victor Breno Pedrosa,
| |
Collapse
|
14
|
Chakraborty D, Sharma N, Kour S, Sodhi SS, Gupta MK, Lee SJ, Son YO. Applications of Omics Technology for Livestock Selection and Improvement. Front Genet 2022; 13:774113. [PMID: 35719396 PMCID: PMC9204716 DOI: 10.3389/fgene.2022.774113] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2021] [Accepted: 05/16/2022] [Indexed: 12/16/2022] Open
Abstract
Conventional animal selection and breeding methods were based on the phenotypic performance of the animals. These methods have limitations, particularly for sex-limited traits and traits expressed later in the life cycle (e.g., carcass traits). Consequently, the genetic gain has been slow with high generation intervals. With the advent of high-throughput omics techniques and the availability of multi-omics technologies and sophisticated analytic packages, several promising tools and methods have been developed to estimate the actual genetic potential of the animals. It has now become possible to collect and access large and complex datasets comprising different genomics, transcriptomics, proteomics, metabolomics, and phonemics data as well as animal-level data (such as longevity, behavior, adaptation, etc.,), which provides new opportunities to better understand the mechanisms regulating animals’ actual performance. The cost of omics technology and expertise of several fields like biology, bioinformatics, statistics, and computational biology make these technology impediments to its use in some cases. The population size and accurate phenotypic data recordings are other significant constraints for appropriate selection and breeding strategies. Nevertheless, omics technologies can estimate more accurate breeding values (BVs) and increase the genetic gain by assisting the section of genetically superior, disease-free animals at an early stage of life for enhancing animal productivity and profitability. This manuscript provides an overview of various omics technologies and their limitations for animal genetic selection and breeding decisions.
Collapse
Affiliation(s)
- Dibyendu Chakraborty
- Division of Animal Genetics and Breeding, Faculty of Veterinary Sciences and Animal Husbandry, Sher-e-Kashmir University of Agricultural Sciences and Technology of Jammu, Ranbir Singh Pura, India
| | - Neelesh Sharma
- Division of Veterinary Medicine, Faculty of Veterinary Sciences and Animal Husbandry, Sher-e-Kashmir University of Agricultural Sciences and Technology of Jammu, Ranbir Singh Pura, India
- *Correspondence: Neelesh Sharma, ; Young Ok Son,
| | - Savleen Kour
- Division of Veterinary Medicine, Faculty of Veterinary Sciences and Animal Husbandry, Sher-e-Kashmir University of Agricultural Sciences and Technology of Jammu, Ranbir Singh Pura, India
| | - Simrinder Singh Sodhi
- Department of Animal Biotechnology, College of Animal Biotechnology, Guru Angad Dev Veterinary and Animal Sciences University, Ludhiana, India
| | - Mukesh Kumar Gupta
- Department of Biotechnology and Medical Engineering, National Institute of Technology, Rourkela, India
| | - Sung Jin Lee
- Department of Animal Biotechnology, College of Animal Life Sciences, Kangwon National University, Chuncheon-si, South Korea
| | - Young Ok Son
- Department of Animal Biotechnology, Faculty of Biotechnology, College of Applied Life Sciences and Interdisciplinary Graduate Program in Advanced Convergence Technology and Science, Jeju National University, Jeju, South Korea
- *Correspondence: Neelesh Sharma, ; Young Ok Son,
| |
Collapse
|
15
|
Liu SH, Ma XY, Hassan FU, Gao TY, Deng TX. Genome-wide analysis of runs of homozygosity in Italian Mediterranean buffalo. J Dairy Sci 2022; 105:4324-4334. [PMID: 35307184 DOI: 10.3168/jds.2021-21543] [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: 11/07/2021] [Accepted: 02/07/2022] [Indexed: 11/19/2022]
Abstract
Runs of homozygosity (ROH) are a powerful tool to explore patterns of genomic inbreeding in animal populations and detect signatures of selection. The present study used ROH analysis to evaluate the genome-wide patterns of homozygosity, inbreeding levels, and distribution of ROH islands using the SNP data sets from 899 Mediterranean buffaloes. A total of 42,433 ROH segments were identified, with an average of 47.20 segments per individual. The ROH comprising mostly shorter segments (1-4 Mb) accounted for approximately 72.29% of all ROH. In contrast, the larger ROH (>8 Mb) class accounted for only 7.97% of all ROH segments. Estimated inbreeding coefficients from ROH (FROH) ranged from 0.0201 to 0.0371. Pearson correlations between FROH and genomic relationship matrix increased with the increase of ROH length. We identified ROH hotspots in 12 genomic regions, located on chromosomes 1, 2, 3, 5, 17, and 19, harboring a total of 122 genes. Protein-protein interaction (PPI) analysis revealed the clustering of these genes into 7 PPI networks. Many genes located in these regions were associated with different production traits. In addition, 5 ROH islands overlapped with cattle quantitative trait loci that were mainly associated with milk traits. These findings revealed the genome-wide autozygosity patterns and inbreeding levels in Mediterranean buffalo. Our study identified many candidate genes related to production traits that could be used to assist in selective breeding for genetic improvement of buffalo.
Collapse
Affiliation(s)
- Shen-He Liu
- College of Animal Science and Technology, Henan Agricultural University, Zhengzhou 450046, China.
| | - Xiao-Ya Ma
- Guangxi Provincial Key Laboratory of Buffalo Genetics, Breeding and Reproduction Technology, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning 530001, China
| | - Faiz-Ul Hassan
- Institute of Animal and Dairy Sciences, University of Agriculture, Faisalabad 38040, Pakistan
| | - Teng-Yun Gao
- College of Animal Science and Technology, Henan Agricultural University, Zhengzhou 450046, China
| | - Ting-Xian Deng
- Guangxi Provincial Key Laboratory of Buffalo Genetics, Breeding and Reproduction Technology, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning 530001, China.
| |
Collapse
|
16
|
Mulim HA, Brito LF, Pinto LFB, Ferraz JBS, Grigoletto L, Silva MR, Pedrosa VB. Characterization of runs of homozygosity, heterozygosity-enriched regions, and population structure in cattle populations selected for different breeding goals. BMC Genomics 2022; 23:209. [PMID: 35291953 PMCID: PMC8925140 DOI: 10.1186/s12864-022-08384-0] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Accepted: 02/10/2022] [Indexed: 01/12/2023] Open
Abstract
Background A decline in the level of genetic diversity in livestock can result in reduced response to selection, greater incidence of genetic defects, and inbreeding depression. In this context, various metrics have been proposed to assess the level of genetic diversity in selected populations. Therefore, the main goals of this study were to: 1) investigate the population structure of 16 cattle populations from 15 different pure breeds or composite populations, which have been selected for different breeds goals; and, 2) identify and compare runs of homozygosity (ROH) and heterozygosity-enriched regions (HER) based on different single nucleotide polymorphism (SNP) panels and whole-genome sequence data (WGS), followed by functional genomic analyses. Results A total of 24,187 ROH were found across all cattle populations, with 55% classified in the 2-4 Mb size group. Fourteen homozygosity islands were found in five populations, where four ROH islands located on BTA1, BTA5, BTA16, and BTA19 overlapped between the Brahman (BRM) and Gyr (GIR) breeds. A functional analysis of the genes found in these islands revealed candidate genes known to play a role in the melanogenesis, prolactin signaling, and calcium signaling pathways. The correlations between inbreeding metrics ranged from 0.02 to 0.95, where the methods based on homozygous genotypes (FHOM), uniting of gametes (FUNI), and genotype additive variance (FGRM) showed strong correlations among them. All methods yielded low to moderate correlations with the inbreeding coefficients based on runs of homozygosity (FROH). For the HER, 3576 runs and 26 islands, distributed across all autosomal chromosomes, were found in regions containing genes mainly related to the immune system, indicating potential balancing selection. Although the analyses with WGS did not enable detection of the same island patterns, it unraveled novel regions not captured when using SNP panel data. Conclusions The cattle populations that showed the largest amount of ROH and HER were Senepol (SEN) and Montana (MON), respectively. Overlapping ROH islands were identified between GIR and BRM breeds, indicating a possible historical connection between the populations. The distribution and pattern of ROH and HER are population specific, indicating that different breeds have experienced divergent selection processes or different genetic processes. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08384-0.
Collapse
Affiliation(s)
| | - Luiz F Brito
- Department of Animal Science, Purdue University, West Lafayette, Indiana, USA
| | | | - José Bento Sterman Ferraz
- Department of Animal Sciences, College of Animal Sciences and Food Engineering, University of São Paulo, Pirassununga, São Paulo, Brazil
| | - Lais Grigoletto
- Department of Animal Science, Purdue University, West Lafayette, Indiana, USA.,Department of Animal Sciences, College of Animal Sciences and Food Engineering, University of São Paulo, Pirassununga, São Paulo, Brazil
| | | | - Victor Breno Pedrosa
- Department of Animal Science, Federal University of Bahia, Salvador, Bahia, Brazil. .,Department of Animal Science, State University of Ponta Grossa, Av. General Carlos Cavalcanti, 4748 - Uvaranas, Ponta Grossa, PR, 84030-900, Brazil.
| |
Collapse
|
17
|
Nazari F, Seyedabadi HR, Noshary A, Emamjomeh-Kashan N, Banabazi MH. A genome-wide scan for signatures of selection in Kurdish horse breed. J Equine Vet Sci 2022; 113:103916. [PMID: 35218903 DOI: 10.1016/j.jevs.2022.103916] [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: 11/10/2021] [Revised: 02/20/2022] [Accepted: 02/21/2022] [Indexed: 11/17/2022]
Abstract
The genetic diversity and genomic regions being under putative natural selection in Kurdish horse population were studied. The samples from 72 horses were genotyped by using GGP Equine 70K SNP arrays. The Ne Slope (NeS) analyses revealed that a sharp decline in Ne has probably occurred around four generations ago, and high frequency of ROH with 2-4 Mbp in length suggested that the inbreeding has probably occurred around 20 generations ago. The effective population size (Ne) was 104 horses up to three generations ago and the average inbreeding (FROH) was 0.047(±0.045). Using de-correlated composite of multiple selection signals (DCMS) and runs of homozygosity (ROH) analyses the genomic regions being under putative selection were detected. By using DCMS, a total of 148 significant SNP (FDR < 0.05) were identified, 40% of which were located on ECA9, where the greatest peak was observed. This genomic region harbors several known QTL which are associated with withers height (body size). Also significant genomic regions (FDR<0.05), harboring QTL associated with insect bite hypersensitivity (IBH), hair density and coat texture, alternate gaits, guttural pouch tympany and temperament were identified. By using outputs of ROH analyses, two hotspot regions (i.e. 30% of individuals was considered as threshold), were identified on ECA7 (50.11-54.36 Mbp) and ECA11 (26.10-29.07 Mbp) harboring QTL associated with withers height, alternate gait and IBH. In summary, the genomic regions being under putative natural selection which harbors known QTL associated with body size and IBH, among others, were introduced. Nevertheless, additional functional and comparative studies are necessary to corroborate their effect on the observed genetic and phenotypic diversity of the Kurdish horses.
Collapse
Affiliation(s)
- Faezeh Nazari
- Department of Animal Science, Science and Research Branch, Islamic Azad University, Tehran, Iran
| | - Hamid-Reza Seyedabadi
- Animal Science Research Institute of Iran (ASRI), Agricultural Research, Education and Extension Organization (AREEO), Karaj, Iran.
| | - Alireza Noshary
- Department of Animal Science, Faculty of Agriculture, Karaj Branch, Islamic Azad University, Karaj, Iran
| | - Nasser Emamjomeh-Kashan
- Department of Animal Science, Science and Research Branch, Islamic Azad University, Tehran, Iran
| | - Mohamad-Hosein Banabazi
- Animal Science Research Institute of Iran (ASRI), Agricultural Research, Education and Extension Organization (AREEO), Karaj, Iran; Department of animal breeding and genetics (HGEN), Centre for Veterinary Medicine and Animal Science (VHC), Swedish University of Agricultural Sciences (SLU), Uppsala 75007, Sweden
| |
Collapse
|
18
|
Rahimmadar S, Ghaffari M, Mokhber M, Williams JL. Linkage Disequilibrium and Effective Population Size of Buffalo Populations of Iran, Turkey, Pakistan, and Egypt Using a Medium Density SNP Array. Front Genet 2021; 12:608186. [PMID: 34950186 PMCID: PMC8689148 DOI: 10.3389/fgene.2021.608186] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2020] [Accepted: 11/03/2021] [Indexed: 11/21/2022] Open
Abstract
Linkage disequilibrium (LD) across the genome provides information to identify the genes and variations related to quantitative traits in genome-wide association studies (GWAS) and for the implementation of genomic selection (GS). LD can also be used to evaluate genetic diversity and population structure and reveal genomic regions affected by selection. LD structure and Ne were assessed in a set of 83 water buffaloes, comprising Azeri (AZI), Khuzestani (KHU), and Mazandarani (MAZ) breeds from Iran, Kundi (KUN) and Nili-Ravi (NIL) from Pakistan, Anatolian (ANA) buffalo from Turkey, and buffalo from Egypt (EGY). The values of corrected r2 (defined as the correlation between two loci) of adjacent SNPs for three pooled Iranian breeds (IRI), ANA, EGY, and two pooled Pakistani breeds (PAK) populations were 0.24, 0.28, 0.27, and 0.22, respectively. The corrected r2 between SNPs decreased with increasing physical distance from 100 Kb to 1 Mb. The LD values for IRI, ANA, EGY, and PAK populations were 0.16, 0.23, 0.24, and 0.21 for less than 100Kb, respectively, which reduced rapidly to 0.018, 0.042, 0.059, and 0.024, for a distance of 1 Mb. In all the populations, the decay rate was low for distances greater than 2Mb, up to the longest studied distance (15 Mb). The r2 values for adjacent SNPs in unrelated samples indicated that the Affymetrix Axiom 90 K SNP genomic array was suitable for GWAS and GS in these populations. The persistency of LD phase (PLDP) between populations was assessed, and results showed that PLPD values between the populations were more than 0.9 for distances of less than 100 Kb. The Ne in the recent generations has declined to the extent that breeding plans are urgently required to ensure that these buffalo populations are not at risk of being lost. We found that results are affected by sample size, which could be partially corrected for; however, additional data should be obtained to be confident of the results.
Collapse
Affiliation(s)
- Shirin Rahimmadar
- Department of Animal Science, Faculty of Agricultural Science, Urmia University, Urmia, Iran
| | - Mokhtar Ghaffari
- Department of Animal Science, Faculty of Agricultural Science, Urmia University, Urmia, Iran
| | - Mahdi Mokhber
- Department of Animal Science, Faculty of Agricultural Science, Urmia University, Urmia, Iran
| | - John L Williams
- Davies Research Centre, School of Animal and Veterinary Sciences, University of Adelaide, Roseworthy, SA, Australia.,Department of Animal Science, Food and Nutrition, Università Cattolica Del Sacro Cuore, Piacenza, Italy
| |
Collapse
|
19
|
Rostamzadeh Mahdabi E, Esmailizadeh A, Ayatollahi Mehrgardi A, Asadi Fozi M. A genome-wide scan to identify signatures of selection in two Iranian indigenous chicken ecotypes. Genet Sel Evol 2021; 53:72. [PMID: 34503452 PMCID: PMC8428137 DOI: 10.1186/s12711-021-00664-9] [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: 01/11/2021] [Accepted: 08/25/2021] [Indexed: 11/10/2022] Open
Abstract
Background Various regions of the chicken genome have been under natural and artificial selection for thousands of years. The substantial diversity that exits among chickens from different geographic regions provides an excellent opportunity to investigate the genomic regions under selection which, in turn, will increase our knowledge about the mechanisms that underlie chicken diversity and adaptation. Several statistics have been developed to detect genomic regions that are under selection. In this study, we applied approaches based on differences in allele or haplotype frequencies (FST and hapFLK, respectively) between populations, differences in long stretches of consecutive homozygous sequences (ROH), and differences in allele frequencies within populations (composite likelihood ratio (CLR)) to identify inter- and intra-populations traces of selection in two Iranian indigenous chicken ecotypes, the Lari fighting chicken and the Khazak or creeper (short-leg) chicken. Results Using whole-genome resequencing data of 32 individuals from the two chicken ecotypes, approximately 11.9 million single nucleotide polymorphisms (SNPs) were detected and used in genomic analyses after quality processing. Examination of the distribution of ROH in the two populations indicated short to long ROH, ranging from 0.3 to 5.4 Mb. We found 90 genes that were detected by at least two of the four applied methods. Gene annotation of the detected putative regions under selection revealed candidate genes associated with growth (DCN, MEOX2 and CACNB1), reproduction (ESR1 and CALCR), disease resistance (S1PR1, ALPK1 and MHC-B), behavior pattern (AGMO, GNAO1 and PSEN1), and morphological traits (IHH and NHEJ1). Conclusions Our findings show that these two phenotypically different indigenous chicken populations have been under selection for reproduction, immune, behavioral, and morphology traits. The results illustrate that selection can play an important role in shaping signatures of differentiation across the genomic landscape of two chicken populations. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-021-00664-9.
Collapse
Affiliation(s)
- Elaheh Rostamzadeh Mahdabi
- Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, 22 Bahman Blvd, Kerman, Iran
| | - Ali Esmailizadeh
- Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, 22 Bahman Blvd, Kerman, Iran
| | - Ahmad Ayatollahi Mehrgardi
- Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, 22 Bahman Blvd, Kerman, Iran
| | - Masood Asadi Fozi
- Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, 22 Bahman Blvd, Kerman, Iran.
| |
Collapse
|
20
|
Strillacci MG, Moradi-Shahrbabak H, Davoudi P, Ghoreishifar SM, Mokhber M, Masroure AJ, Bagnato A. A genome-wide scan of copy number variants in three Iranian indigenous river buffaloes. BMC Genomics 2021; 22:305. [PMID: 33902439 PMCID: PMC8077898 DOI: 10.1186/s12864-021-07604-3] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2020] [Accepted: 04/11/2021] [Indexed: 01/29/2023] Open
Abstract
BACKGROUND In Iran, river buffalo is of great importance. It plays an important role in the economy of the Country, because its adaptation to harsh climate conditions and long productive lifespan permitting its farming across the Country and to convert low-quality feed into valuable milk. The genetic variability in Iranian buffalo breeds have been recently studied using SNPs genotyping data, but a whole genome Copy Number Variants (CNVs) mapping was not available. The aim of this study was to perform a genome wide CNV scan in 361 buffaloes of the three Iranian river breeds (Azeri, Khuzestani and Mazandarani) through the analysis of data obtained using the Axiom® Buffalo Genotyping Array 90 K. RESULTS CNVs detection resulted in a total of 9550 CNVs and 302 CNVRs identified in at least 5% of samples within breed, covering around 1.97% of the buffalo genome. and A total of 22 CNVRs were identified in all breeds and a different proportion of regions were in common among the three populations. Within the more represented CNVRs (n = 302) mapped a total of 409 buffalo genes, some of which resulted associated with morphological, healthy, milk, meat and reproductive traits, according to Animal Genome Cattle database. CONCLUSIONS This work provides a step forward in the interpretation of genomic variation within and among the buffalo populations, releasing a first map of CNVs and providing insights about their recent selection and adaptation to environment. The presence of the set of genes and QTL traits harbored in the CNVRs could be possibly linked with the buffalo's natural adaptive history together to a recent selection for milk used as primary food source from this species.
Collapse
Affiliation(s)
- Maria G. Strillacci
- Department of Veterinary Medicine, Università degli Studi di Milano, Via dell’Università 6, 26900 Lodi, Italy
| | - Hossein Moradi-Shahrbabak
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167 Iran
| | - Pourya Davoudi
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS B2N5E3 Canada
| | - Seyed Mohammad Ghoreishifar
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167 Iran
| | - Mahdi Mokhber
- Department of Animal Science, Faculty of Agriculture and Natural resources, Urmia University, 11Km Sero Road, P. O. Box: 165, Urmia, 57561-51818 Iran
| | - Anoar Jamai Masroure
- Department of Veterinary Medicine, Università degli Studi di Milano, Via dell’Università 6, 26900 Lodi, Italy
| | - Alessandro Bagnato
- Department of Veterinary Medicine, Università degli Studi di Milano, Via dell’Università 6, 26900 Lodi, Italy
| |
Collapse
|
21
|
Cesarani A, Biffani S, Garcia A, Lourenco D, Bertolini G, Neglia G, Misztal I, Macciotta NPP. Genomic investigation of milk production in Italian buffalo. ITALIAN JOURNAL OF ANIMAL SCIENCE 2021. [DOI: 10.1080/1828051x.2021.1902404] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Alberto Cesarani
- Department of Animal and Dairy Science, University of Georgia, Athens, GA, USA
| | - Stefano Biffani
- Consiglio Nazionale delle Ricerche (CNR), Istituto di biologia e biotecnologia agraria (IBBA), Milano, Italy
| | - Andre Garcia
- Department of Animal and Dairy Science, University of Georgia, Athens, GA, USA
| | - Daniela Lourenco
- Department of Animal and Dairy Science, University of Georgia, Athens, GA, USA
| | - Giacomo Bertolini
- Associazione Nazionale Allevatori Specie Bufalina (ANASB), Caserta, Italy
| | - Gianluca Neglia
- Dipartimento di Medicina Veterinaria e Produzioni Animali, II University of Naples, Napoli, Italy
| | - Ignacy Misztal
- Department of Animal and Dairy Science, University of Georgia, Athens, GA, USA
| | | |
Collapse
|
22
|
Ghoreishifar SM, Rochus CM, Moghaddaszadeh-Ahrabi S, Davoudi P, Salek Ardestani S, Zinovieva NA, Deniskova TE, Johansson AM. Shared Ancestry and Signatures of Recent Selection in Gotland Sheep. Genes (Basel) 2021; 12:genes12030433. [PMID: 33802939 PMCID: PMC8002741 DOI: 10.3390/genes12030433] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2021] [Revised: 03/01/2021] [Accepted: 03/10/2021] [Indexed: 12/13/2022] Open
Abstract
Gotland sheep, a breed native to Gotland, Sweden (an island in the Baltic Sea), split from the Gute sheep breed approximately 100 years ago, and since, has probably been crossed with other breeds. This breed has recently gained popularity, due to its pelt quality. This study estimates the shared ancestors and identifies recent selection signatures in Gotland sheep using 600 K single nucleotide polymorphism (SNP) genotype data. Admixture analysis shows that the Gotland sheep is a distinct breed, but also has shared ancestral genomic components with Gute (~50%), Karakul (~30%), Romanov (~20%), and Fjällnäs (~10%) sheep breeds. Two complementary methods were applied to detect selection signatures: A Bayesian population differentiation FST and an integrated haplotype homozygosity score (iHS). Our results find that seven significant SNPs (q-value < 0.05) using the FST analysis and 55 significant SNPs (p-value < 0.0001) using the iHS analysis. Of the candidate genes that contain significant markers, or are in proximity to them, we identify several belongings to the keratin genes, RXFP2, ADCY1, ENOX1, USF2, COX7A1, ARHGAP28, CRYBB2, CAPNS1, FMO3, and GREB1. These genes are involved in wool quality, polled and horned phenotypes, fertility, twining rate, meat quality, and growth traits. In summary, our results provide shared founders of Gotland sheep and insight into genomic regions maintained under selection after the breed was formed. These results contribute to the detection of candidate genes and QTLs underlying economic traits in sheep.
Collapse
Affiliation(s)
- Seyed Mohammad Ghoreishifar
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj 31587-11167, Iran;
| | - Christina Marie Rochus
- Animal Breeding and Genomics, Wageningen University and Research, P.O. Box 338, 6700 AH Wageningen, The Netherlands;
| | - Sima Moghaddaszadeh-Ahrabi
- Department of Animal Science, Faculty of Agriculture and Natural Resources, Islamic Azad University, Tabriz Branch, Tabriz 5157944533, Iran;
| | - Pourya Davoudi
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS B2N5E3, Canada; (P.D.); (S.S.A.)
| | - Siavash Salek Ardestani
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS B2N5E3, Canada; (P.D.); (S.S.A.)
| | - Natalia A. Zinovieva
- L.K. Ernst Federal Research Center for Animal Husbandry, 142132 Podolsk, Russia; (N.A.Z.); (T.E.D.)
| | - Tatiana E. Deniskova
- L.K. Ernst Federal Research Center for Animal Husbandry, 142132 Podolsk, Russia; (N.A.Z.); (T.E.D.)
| | - Anna M. Johansson
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, SE-75007 Uppsala, Sweden
- Correspondence:
| |
Collapse
|
23
|
Baazaoui I, Bedhiaf-Romdhani S, Mastrangelo S, Ciani E. Genome-wide analyses reveal population structure and identify candidate genes associated with tail fatness in local sheep from a semi-arid area. Animal 2021; 15:100193. [PMID: 33715983 DOI: 10.1016/j.animal.2021.100193] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2020] [Revised: 01/20/2021] [Accepted: 01/25/2021] [Indexed: 12/21/2022] Open
Abstract
Under a climate change perspective, the genetic make-up of local livestock breeds showing adaptive traits should be explored and preserved as a priority. We used genotype data from the ovine 50 k Illumina BeadChip for assessing breed autozygosity based on runs of homozygosity (ROH) and fine-scale genetic structure and for detecting genomic regions under selection in 63 Tunisian sheep samples. The average genomic inbreeding coefficients based on ROH were estimated at 0.017, 0.021, and 0.024 for Barbarine (BAR, n = 26), Noire de Thibar (NDT, n = 23), and Queue fine de l'Ouest (QFO, n = 14) breeds, respectively. The genomic relationships among individuals based on identity by state (IBS) distance matrix highlighted a recent introgression of QFO into the BAR and a genetic differentiation of NDT samples, possibly explained by past introgression of European gene pools. Genome-wide scan for ROH across breeds and within the BAR sample set identified an outstanding signal on chromosome 13 (46.58-49.61 Mbp). These results were confirmed using FST index, differentiating fat vs. thin-tailed individuals. Candidate genes under selection pressure (CDS2, PROKR1, and BMP2) were associated to lipid storage and probably preferentially selected in fat-tailed BAR animals. Our findings suggest paying more attention to preserve the genetic integrity and adaptive alleles of local sheep breeds.
Collapse
Affiliation(s)
- I Baazaoui
- Faculty of Sciences of Bizerte, University of Carthage, 7021 Zarzouna, Bizerte, Tunisia
| | - S Bedhiaf-Romdhani
- Institut National de la Recherche Agronomique de Tunisie, Laboratoire des Productions Animales et Fourragères, Université de Carthage, 2049 Ariana, Tunisie..
| | - S Mastrangelo
- Department of Agricultural, Food and Forest Sciences, University of Palermo, 90128 Palermo, Italy
| | - E Ciani
- Department of Biosciences, Biotechnologies and Biopharmaceutics, University of Bari, 70121 Bari, Italy
| |
Collapse
|
24
|
Macciotta NPP, Colli L, Cesarani A, Ajmone-Marsan P, Low WY, Tearle R, Williams JL. The distribution of runs of homozygosity in the genome of river and swamp buffaloes reveals a history of adaptation, migration and crossbred events. Genet Sel Evol 2021; 53:20. [PMID: 33639853 PMCID: PMC7912491 DOI: 10.1186/s12711-021-00616-3] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2020] [Accepted: 02/17/2021] [Indexed: 01/03/2023] Open
Abstract
Background Water buffalo is one of the most important livestock species in the world. Two types of water buffalo exist: river buffalo (Bubalus bubalis bubalis) and swamp buffalo (Bubalus bubalis carabanensis). The buffalo genome has been recently sequenced, and thus a new 90 K single nucleotide polymorphism (SNP) bead chip has been developed. In this study, we investigated the genomic population structure and the level of inbreeding of 185 river and 153 swamp buffaloes using runs of homozygosity (ROH). Analyses were carried out jointly and separately for the two buffalo types. Results The SNP bead chip detected in swamp about one-third of the SNPs identified in the river type. In total, 18,116 ROH were detected in the combined data set (17,784 SNPs), and 16,251 of these were unique. ROH were present in both buffalo types mostly detected (~ 59%) in swamp buffalo. The number of ROH per animal was larger and genomic inbreeding was higher in swamp than river buffalo. In the separated datasets (46,891 and 17,690 SNPs for river and swamp type, respectively), 19,760 and 10,581 ROH were found in river and swamp, respectively. The genes that map to the ROH islands are associated with the adaptation to the environment, fitness traits and reproduction. Conclusions Analysis of ROH features in the genome of the two water buffalo types allowed their genomic characterization and highlighted differences between buffalo types and between breeds. A large ROH island on chromosome 2 was shared between river and swamp buffaloes and contained genes that are involved in environmental adaptation and reproduction. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-021-00616-3.
Collapse
Affiliation(s)
| | - Licia Colli
- Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti-DIANA, Università Cattolica del Sacro Cuore, Piacenza, Italia.,Centro di Ricerca sulla Biodiversità e sul DNA Antico-BioDNA, Università Cattolica del Sacro Cuore, Piacenza, Italia
| | - Alberto Cesarani
- Dipartimento di Agraria, Università degli Studi di Sassari, Sassari, Italia. .,Department of Animal and Dairy Science, University of Georgia, Athens, GA, USA.
| | - Paolo Ajmone-Marsan
- Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti-DIANA, Università Cattolica del Sacro Cuore, Piacenza, Italia.,Centro di Ricerca Nutrigenomica e Proteomica-PRONUTRIGEN, Università Cattolica del Sacro Cuore, Piacenza, Italia
| | - Wai Y Low
- The Davies Research Centre, School of Animal and Veterinary Sciences, University of Adelaide, Roseworthy, SA, 5371, Australia
| | - Rick Tearle
- The Davies Research Centre, School of Animal and Veterinary Sciences, University of Adelaide, Roseworthy, SA, 5371, Australia
| | - John L Williams
- Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti-DIANA, Università Cattolica del Sacro Cuore, Piacenza, Italia.,The Davies Research Centre, School of Animal and Veterinary Sciences, University of Adelaide, Roseworthy, SA, 5371, Australia
| |
Collapse
|
25
|
Gorssen W, Meyermans R, Janssens S, Buys N. A publicly available repository of ROH islands reveals signatures of selection in different livestock and pet species. Genet Sel Evol 2021; 53:2. [PMID: 33397285 PMCID: PMC7784028 DOI: 10.1186/s12711-020-00599-7] [Citation(s) in RCA: 45] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2020] [Accepted: 12/11/2020] [Indexed: 01/18/2023] Open
Abstract
BACKGROUND Runs of homozygosity (ROH) have become the state-of-the-art method for analysis of inbreeding in animal populations. Moreover, ROH are suited to detect signatures of selection via ROH islands and are used in other applications, such as genomic prediction and genome-wide association studies (GWAS). Currently, a vast amount of single nucleotide polymorphism (SNP) data is available online, but most of these data have never been used for ROH analysis. Therefore, we performed a ROH analysis on large medium-density SNP datasets in eight animal species (cat, cattle, dog, goat, horse, pig, sheep and water buffalo; 442 different populations) and make these results publicly available. RESULTS The results include an overview of ROH islands per population and a comparison of the incidence of these ROH islands among populations from the same species, which can assist researchers when studying other (livestock) populations or when looking for similar signatures of selection. We were able to confirm many known ROH islands, for example signatures of selection for the myostatin (MSTN) gene in sheep and horses. However, our results also included multiple other ROH islands, which are common to many populations and not identified to date (e.g. on chromosomes D4 and E2 in cats and on chromosome 6 in sheep). CONCLUSIONS We are confident that our repository of ROH islands is a valuable reference for future studies. The discovered ROH island regions represent a unique starting point for new studies or can be used as a reference for future studies. Furthermore, we encourage authors to add their population-specific ROH findings to our repository.
Collapse
Affiliation(s)
- Wim Gorssen
- Livestock Genetics, Department of Biosystems, KU Leuven, Kasteelpark Arenberg 30, Box 2472, 3001, Leuven, Belgium
| | - Roel Meyermans
- Livestock Genetics, Department of Biosystems, KU Leuven, Kasteelpark Arenberg 30, Box 2472, 3001, Leuven, Belgium
| | - Steven Janssens
- Livestock Genetics, Department of Biosystems, KU Leuven, Kasteelpark Arenberg 30, Box 2472, 3001, Leuven, Belgium
| | - Nadine Buys
- Livestock Genetics, Department of Biosystems, KU Leuven, Kasteelpark Arenberg 30, Box 2472, 3001, Leuven, Belgium.
| |
Collapse
|
26
|
Rezaei Topraggaleh T, Numan Bucak M, Shahverdi M, Koohestani Y, Furkan Batur A, Rahimizadeh P, Ili P, Gul M, Ashrafzade AM, Kazem-Allilo A, Garip M, Shahverdi A. Functional and Flow Cytometric Analysis of Buffalo Cryopreserved Spermatozoa: Comparison of Different Breeds and Incubation Times. INTERNATIONAL JOURNAL OF FERTILITY & STERILITY 2021; 15:252-257. [PMID: 34913292 PMCID: PMC8530219 DOI: 10.22074/ijfs.2021.521116.1057] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Subscribe] [Scholar Register] [Received: 12/14/2020] [Accepted: 04/15/2021] [Indexed: 11/18/2022]
Abstract
BACKGROUND The purpose of this research was to compare the functional parameters of frozen-thawed Iranian Azari buffalo spermatozoa with imported semen samples of Italian Mediterranean buffalo (IMB) after the thawing process and 4 hours of incubation. MATERIALS AND METHODS In this experimental study, a total of twenty-four ejaculates from four Iranian Azari buffalo bulls were collected. Semen samples were diluted in AndroMed extender at a concentration of 50×106 spermatozoa/ ml. The diluted samples were filled in 0.5 ml straws and were frozen in a programmable freezer. For imported semen samples, twenty-four samples of four IMB were used, which were diluted in AndroMed extender and frozen by the same procedure. Frozen-thawed sperm motion patterns, mitochondrial activity, membrane integrity, DNA integrity, reactive oxygen species (ROS), and apoptosis status were evaluated immediately after thawing and 4 hours of incubation. RESULTS Post-thawed sperm motility, progressive motility (PM), mitochondrial activity, membrane integrity were significantly higher in imported semen samples in compare with Iranian Azari buffalo. After 4 hours of incubation, sperm velocity patterns were superior in Iranian Azari semen samples. Moreover, the percentage of sperm cells with un-damaged DNA was higher in Iranian semen samples compared to imported samples at the time 0 of incubation. Following 4 hours of incubation, a significant increase in intracellular ROS level leads to reduced membrane integrity, mitochondrial activity, and DNA integrity in both buffalo breeds. At time 4, Iranian samples showed significantly lower apoptosis and higher dead spermatozoa compared to imported semen samples. CONCLUSION Our study showed that the post-thawed quality of Iranian Azari buffalo semen was comparable with imported samples after 4 hours of incubation. Further investigations are recommended to assess the in vitro and in vivo fertility rate of both buffalo breeds.
Collapse
Affiliation(s)
- Tohid Rezaei Topraggaleh
- Department of Anatomical Sciences, Faculty of Medicine, Urmia University of Medical Sciences, Urmia, Iran
| | - Mustafa Numan Bucak
- Department of Reproduction and Artificial Insemination, Faculty of Veterinary Medicine, Selcuk University, Konya, Turkey
| | - Maryam Shahverdi
- Department of Embryology, Reproductive Biomedicine Research Center, Royan Institute for Reproductive Biomedicine, ACECR, Tehran, Iran
| | - Yegane Koohestani
- Department of Anatomical Sciences, Faculty of Medicine, Guilan University of Medical Sciences, Rasht, Iran
| | - Ali Furkan Batur
- Department of Urology, Selcuk University School of Medicine, Konya, Turkey
| | - Pegah Rahimizadeh
- Department of Embryology, Reproductive Biomedicine Research Center, Royan Institute for Reproductive Biomedicine, ACECR, Tehran, Iran
| | - Pinar Ili
- Department of Medical Services and Techniques, Denizli Vocational School of Health Services, Pamukkale University, Denizli, Turkey
| | - Murat Gul
- Department of Urology, Selcuk University School of Medicine, Konya, Turkey
| | - Amir Mahdi Ashrafzade
- Department of Genetics, Reproductive Biomedicine Research Center, Royan Institute for Reproductive Biomedicine, ACECR, Tehran, Iran
| | | | - Mustafa Garip
- Department of Animal Science, Faculty of Veterinary Medicine, Selcuk University, Konya, Turkey
| | - Abdolhossein Shahverdi
- Department of Embryology, Reproductive Biomedicine Research Center, Royan Institute for Reproductive Biomedicine, ACECR, Tehran, Iran,Reproductive Epidemiology Research Center, Royan Institute for Reproductive Biomedicine, ACECR, Tehran, Iran,P.O.Box: 16635-148Department of EmbryologyReproductive Biomedicine Research CenterRoyan Institute for Reproductive BiomedicineACECRTehranIran
| |
Collapse
|
27
|
Genome-Wide Analysis Revealed Homozygosity and Demographic History of Five Chinese Sheep Breeds Adapted to Different Environments. Genes (Basel) 2020; 11:genes11121480. [PMID: 33317115 PMCID: PMC7764688 DOI: 10.3390/genes11121480] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2020] [Revised: 12/06/2020] [Accepted: 12/07/2020] [Indexed: 01/18/2023] Open
Abstract
Homozygosity of long sequence genotypes are a result of parents transmitting identical haplotypes, which can be used to estimate their auto-zygosity. Therefore, we used high-density SNP Chip data to characterize the auto-zygosity of each breed according to the occurrence and distribution of runs of homozygosity (ROH). Subsequently, we identified the genomic regions with high runs of homozygosity frequencies within individuals of each breed. We selected 96 sheep samples from five local Chinese sheep breeds belonging to different geographical locations. We identified 3046 ROHs within the study breed individuals, among which the longer segments (>1–5 Mb) were dominant. On average, ROH segments covered about 12% of the genomes; the coverage rate of OAR20 was the lowest and that of OAR2 was the highest. The distribution analysis of runs of homozygosity showed that the detected ROH mainly distributed between >26 and 28 Mb. The Hetian and Hu sheep showed the lowest ROH distribution. The estimation of homozygosity level reflects the history of modern and ancient inbreeding, which may affect the genomes of Chinese indigenous sheep breeds and indicate that some animals have experienced recent self-pollination events (Yabuyi, Karakul and Wadi). In these sheep breeds, the genomic regions were assumed to be under selection signatures frequently in line with long ROH. These regions included candidate genes associated with disease resistance traits (5S_rRNA), the innate and adaptive immune response (HERC2 and CYFIP1), digestion and metabolism (CENPJ), growth (SPP1), body size and developments (GJB2 and GJA3). This study highlighted new insights into the ROH patterns and provides a basis for future breeding and conservation strategies of Chinese sheep breeds.
Collapse
|
28
|
Minervino AHH, Zava M, Vecchio D, Borghese A. Bubalus bubalis: A Short Story. Front Vet Sci 2020; 7:570413. [PMID: 33335917 PMCID: PMC7736047 DOI: 10.3389/fvets.2020.570413] [Citation(s) in RCA: 53] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2020] [Accepted: 10/27/2020] [Indexed: 12/14/2022] Open
Abstract
The domestic buffalo (Bubalus bubalis), also known as water buffalo or Asian buffalo to prevent confusion with the American bison (Bison bison), wrongly named buffalo in North America, comprises two subspecies: the river buffalo (B. bubalis bubalis) and the swamp buffalo (B. bubalis kerebau). The swamp buffalo has a consistent phenotype and is considered as one type, even if many breeds are recognized within it; conversely, the river buffalo subspecies has many breeds. We found limited information available regarding the worldwide distribution of buffaloes. The best estimate is that 208,098,759 buffalo head are distributed in 77 countries in five continents. In this review, we presented the basic aspects of the water buffalo and unraveled the buffalo path followed from the origin of the species to its current global distribution. We reviewed several data sources to provide a better estimate of the world buffalo count and distribution.
Collapse
Affiliation(s)
| | - Marco Zava
- Argentine Buffalo Breeders Association, Buenos Aires, Argentina
| | - Domenico Vecchio
- Italian National Reference Centre on Water Buffalo Farming and Productions Hygiene and Technologies (CReNBuf), Istituto Zooprofilattico Sperimentale del Mezzogiorno, Portici, Italy
| | - Antonio Borghese
- International Buffalo Federation - IBF, FAO/ESCORENA Buffalo Network, Animal Production Research Institute, Rome, Italy
| |
Collapse
|
29
|
Nascimento AV, Cardoso DF, Santos DJA, Romero ARS, Scalez DCB, Borquis RRA, Neto FRA, Gondro C, Tonhati H. Inbreeding coefficients and runs of homozygosity islands in Brazilian water buffalo. J Dairy Sci 2020; 104:1917-1927. [PMID: 33272579 DOI: 10.3168/jds.2020-18397] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2020] [Accepted: 09/10/2020] [Indexed: 01/03/2023]
Abstract
Characterization of autozygosity is relevant to monitor genetic diversity and manage inbreeding levels in breeding programs. Identification of autozygosity hotspots can unravel genomic regions targeted by selection for economically important traits and can help identify candidate genes for selection. In this study, we estimated the inbreeding levels of a Brazilian population of Murrah buffalo undergoing selection for milk production traits, particularly milk yield. We also studied the distribution of runs of homozygosity (ROH) islands and identified putative genes and quantitative trait loci (QTL) under selection. We genotyped 422 Murrah buffalo for 51,611 SNP; 350 of these had ROH longer than 10 Mb, indicating the occurrence of inbreeding in the last 5 generations. The mean length of the ROH per animal was 4.28 ± 1.85 Mb. Inbreeding coefficients were calculated from the genomic relationship matrix, the pedigree, and the ROH, with estimates varying between 0.242 and 0.035. Inbreeding estimates from the pedigree had a low correlation with the genomic estimates, and estimates from the genomic relationship matrix were much higher than those from the pedigree or the ROH. Signatures of selection were identified in 6 genomic regions, located on chromosomes 1, 2, 3, 5, 16, and 18, encompassing a total of 190 genes and 174 QTL. Many of the genes (e.g., APRT and ACSF3) and QTL identified are related to milk production traits, such as milk yield, milk fat yield and percentage, and milk protein yield and percentage. Other genes are associated with reproduction and immune response traits as well as morphological aspects of the buffalo species. Inbreeding levels in this population are still low but are increasing due to selection and should be managed to avoid future losses due to inbreeding depression. The proximity of genes linked to milk production traits with genes associated with reproduction and immune system traits suggests the need to include these latter genes in the breeding program to avoid negatively affecting them due to selection for production traits.
Collapse
Affiliation(s)
- A V Nascimento
- Department of Animal Science, São Paulo State University (UNESP), Jaboticabal, 14884900, Brazil
| | - D F Cardoso
- Department of Animal Science, São Paulo State University (UNESP), Jaboticabal, 14884900, Brazil
| | - D J A Santos
- Department of Animal Science, University of Maryland, College Park 20742
| | - A R S Romero
- Department of Animal Science, São Paulo State University (UNESP), Jaboticabal, 14884900, Brazil
| | - D C B Scalez
- Department of Animal Science, São Paulo State University (UNESP), Jaboticabal, 14884900, Brazil
| | - R R A Borquis
- College of Agricultural Sciences, Federal University of Grande Dourados (UFGD), Dourados, 79804970, Brazil
| | - F R A Neto
- Goiano Federal Institute, Campus Rio Verde, Rio Verde, 75909120, Brazil
| | - C Gondro
- Department of Animal Science, Michigan State University, East Lansing 48824
| | - H Tonhati
- Department of Animal Science, São Paulo State University (UNESP), Jaboticabal, 14884900, Brazil.
| |
Collapse
|
30
|
Fallahi MH, Shahrbabak HM, Shahrbabak MM, Arpanahi RA, Gholami S. Assessment of Genetic Diversity in Azerbaijani Buffalo Population in Iran Based on Runs of Homozygosity Stretches. RUSS J GENET+ 2020. [DOI: 10.1134/s102279542010004x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
31
|
Sun T, Huang GY, Wang ZH, Teng SH, Cao YH, Sun JL, Hanif Q, Chen NB, Lei CZ, Liao YY. Selection signatures of Fuzhong Buffalo based on whole-genome sequences. BMC Genomics 2020; 21:674. [PMID: 32993537 PMCID: PMC7526191 DOI: 10.1186/s12864-020-07095-8] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2020] [Accepted: 09/23/2020] [Indexed: 12/13/2022] Open
Abstract
BACKGROUND Fuzhong buffalo, a native breed of Guangxi Zhuang Autonomous Region, is traditionally used as a draft animal to provide farm power in the rice cultivation. In addition, the Fuzhong buffalo also prepared for the bullfighting festival organized by the locals. The detection of the selective signatures in its genome can help in elucidating the selection mechanisms in its stamina and muscle development of a draft animal. RESULTS In this study, we analyzed 27 whole genomes of buffalo (including 15 Fuzhong buffalo genomes and 12 published buffalo genomes from Upper Yangtze region). The ZHp, ZFst, π-Ratio, and XP-EHH statistics were used to identify the candidate signatures of positive selection in Fuzhong buffalo. Our results detected a set of candidate genes involving in the pathways and GO terms associated with the response to exercise (e.g., ALDOA, STAT3, AKT2, EIF4E2, CACNA2D2, TCF4, CDH2), immunity (e.g., PTPN22, NKX2-3, PIK3R1, ITK, TMEM173), nervous system (e.g., PTPN21, ROBO1, HOMER1, MAGI2, SLC1A3, NRG3, SNAP47, CTNNA2, ADGRL3). In addition, we also identified several genes related to production and growth traits (e.g., PHLPP1, PRKN, MACF1, UCN3, RALGAPA1, PHKB, PKD1L). Our results depicted several pathways, GO terms, and candidate genes to be associated with response to exercise, immunity, nervous system, and growth traits. CONCLUSIONS The selective sweep analysis of the Fuzhong buffalo demonstrated positive selection pressure on potential target genes involved in behavior, immunity, and growth traits, etc. Our findings provided a valuable resource for future research on buffalo breeding and an insight into the mechanisms of artificial selection.
Collapse
Affiliation(s)
- Ting Sun
- Animal Husbandry Institute of Guangxi Zhuang Autonomous Region, Guangxi Key Laboratory of Livestock Genetic Improvement, Nanning, 530001, China.,College of Animal Science and Technology, Northwest A&F University, Yangling, 712100, Shaanxi, China
| | - Guang-Yun Huang
- Animal Husbandry Institute of Guangxi Zhuang Autonomous Region, Guangxi Key Laboratory of Livestock Genetic Improvement, Nanning, 530001, China
| | - Zi-Hao Wang
- Animal Husbandry Institute of Guangxi Zhuang Autonomous Region, Guangxi Key Laboratory of Livestock Genetic Improvement, Nanning, 530001, China
| | - Shao-Hua Teng
- Animal Husbandry Institute of Guangxi Zhuang Autonomous Region, Guangxi Key Laboratory of Livestock Genetic Improvement, Nanning, 530001, China
| | - Yan-Hong Cao
- Animal Husbandry Institute of Guangxi Zhuang Autonomous Region, Guangxi Key Laboratory of Livestock Genetic Improvement, Nanning, 530001, China
| | - Jun-Li Sun
- Animal Husbandry Institute of Guangxi Zhuang Autonomous Region, Guangxi Key Laboratory of Livestock Genetic Improvement, Nanning, 530001, China
| | - Quratulain Hanif
- Computational Biology Laboratory, Agricultural Biotechnology Division, National Institute for Biotechnology and Genetic Engineering, Faisalabad, Pakistan.,Department of Biotechnology, Pakistan Institute of Engineering and Applied Sciences, Nilore, Islamabad, Pakistan
| | - Ning-Bo Chen
- College of Animal Science and Technology, Northwest A&F University, Yangling, 712100, Shaanxi, China
| | - Chu-Zhao Lei
- College of Animal Science and Technology, Northwest A&F University, Yangling, 712100, Shaanxi, China.
| | - Yu-Ying Liao
- Animal Husbandry Institute of Guangxi Zhuang Autonomous Region, Guangxi Key Laboratory of Livestock Genetic Improvement, Nanning, 530001, China.
| |
Collapse
|
32
|
Ghoreishifar SM, Eriksson S, Johansson AM, Khansefid M, Moghaddaszadeh-Ahrabi S, Parna N, Davoudi P, Javanmard A. Signatures of selection reveal candidate genes involved in economic traits and cold acclimation in five Swedish cattle breeds. Genet Sel Evol 2020; 52:52. [PMID: 32887549 PMCID: PMC7487911 DOI: 10.1186/s12711-020-00571-5] [Citation(s) in RCA: 36] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Accepted: 08/21/2020] [Indexed: 02/01/2023] Open
Abstract
Background Thousands of years of natural and artificial selection have resulted in indigenous cattle breeds that are well-adapted to the environmental challenges of their local habitat and thereby are considered as valuable genetic resources. Understanding the genetic background of such adaptation processes can help us design effective breeding objectives to preserve local breeds and improve commercial cattle. To identify regions under putative selection, GGP HD 150 K single nucleotide polymorphism (SNP) arrays were used to genotype 106 individuals representing five Swedish breeds i.e. native to different regions and covering areas with a subarctic cold climate in the north and mountainous west, to those with a continental climate in the more densely populated south regions. Results Five statistics were incorporated within a framework, known as de-correlated composite of multiple signals (DCMS) to detect signatures of selection. The obtained p-values were adjusted for multiple testing (FDR < 5%), and significant genomic regions were identified. Annotation of genes in these regions revealed various verified and novel candidate genes that are associated with a diverse range of traits, including e.g. high altitude adaptation and response to hypoxia (DCAF8, PPP1R12A, SLC16A3, UCP2, UCP3, TIGAR), cold acclimation (AQP3, AQP7, HSPB8), body size and stature (PLAG1, KCNA6, NDUFA9, AKAP3, C5H12orf4, RAD51AP1, FGF6, TIGAR, CCND2, CSMD3), resistance to disease and bacterial infection (CHI3L2, GBP6, PPFIBP1, REP15, CYP4F2, TIGD2, PYURF, SLC10A2, FCHSD2, ARHGEF17, RELT, PRDM2, KDM5B), reproduction (PPP1R12A, ZFP36L2, CSPP1), milk yield and components (NPC1L1, NUDCD3, ACSS1, FCHSD2), growth and feed efficiency (TMEM68, TGS1, LYN, XKR4, FOXA2, GBP2, GBP5, FGD6), and polled phenotype (URB1, EVA1C). Conclusions We identified genomic regions that may provide background knowledge to understand the mechanisms that are involved in economic traits and adaptation to cold climate in cattle. Incorporating p-values of different statistics in a single DCMS framework may help select and prioritize candidate genes for further analyses.
Collapse
Affiliation(s)
- Seyed Mohammad Ghoreishifar
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Susanne Eriksson
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, SE-75007, Uppsala, Sweden.
| | - Anna M Johansson
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, SE-75007, Uppsala, Sweden
| | - Majid Khansefid
- AgriBio Centre for AgriBioscience, Agriculture Victoria, Bundoora, VIC, 3083, Australia
| | - Sima Moghaddaszadeh-Ahrabi
- Department of Animal Science, Faculty of Agriculture and Natural Resources, Islamic Azad University, Tabriz Branch, Tabriz, Iran
| | - Nahid Parna
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Pourya Davoudi
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, B2N5E3, Canada
| | - Arash Javanmard
- Department of Animal Science, Faculty of Agriculture, University of Tabriz, Tabriz, Iran
| |
Collapse
|
33
|
Estimates of Autozygosity Through Runs of Homozygosity in Farmed Coho Salmon. Genes (Basel) 2020; 11:genes11050490. [PMID: 32365758 PMCID: PMC7290985 DOI: 10.3390/genes11050490] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2020] [Revised: 04/23/2020] [Accepted: 04/24/2020] [Indexed: 11/17/2022] Open
Abstract
The characterization of runs of homozygosity (ROH), using high-density single nucleotide polymorphisms (SNPs) allows inferences to be made about the past demographic history of animal populations and the genomic ROH has become a common approach to characterize the inbreeding. We aimed to analyze and characterize ROH patterns and compare different genomic and pedigree-based methods to estimate the inbreeding coefficient in two pure lines (POP A and B) and one recently admixed line (POP C) of coho salmon (Oncorhynchus kisutch) breeding nuclei, genotyped using a 200 K Affymetrix Axiom® myDesign Custom SNP Array. A large number and greater mean length of ROH were found for the two “pure” lines and the recently admixed line (POP C) showed the lowest number and smaller mean length of ROH. The ROH analysis for different length classes suggests that all three coho salmon lines the genome is largely composed of a high number of short segments (<4 Mb), and for POP C no segment >16 Mb was found. A high variable number of ROH, mean length and inbreeding values across chromosomes; positively the consequence of artificial selection. Pedigree-based inbreeding values tended to underestimate genomic-based inbreeding levels, which in turn varied depending on the method used for estimation. The high positive correlations between different genomic-based inbreeding coefficients suggest that they are consistent and may be more accurate than pedigree-based methods, given that they capture information from past and more recent demographic events, even when there are no pedigree records available.
Collapse
|
34
|
Correction to: Genomic measures of inbreeding coefficients and genome-wide scan for runs of homozygosity islands in Iranian river buffalo, Bubalus bubalis. BMC Genet 2020; 21:42. [PMID: 32268878 PMCID: PMC7140384 DOI: 10.1186/s12863-020-00841-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|