1
|
Deng TX, Ma XY, Duan A, Lu XR, Abdel-Shafy H. Genome-wide CNV analysis reveals candidate genes associated with milk production traits in water buffalo (Bubalus bubalis). J Dairy Sci 2024:S0022-0302(24)00802-6. [PMID: 38762109 DOI: 10.3168/jds.2023-24614] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2023] [Accepted: 03/28/2024] [Indexed: 05/20/2024]
Abstract
Buffaloes are vital contributors to the global dairy industry. Understanding the genetic basis of milk production traits in buffalo populations is essential for breeding programs and improving productivity. In this study, we conducted whole-genome resequencing on 387 buffalo genomes from 29 diverse Asian breeds, including 132 river buffaloes, 129 swamp buffaloes, and 126 crossbred buffaloes. We identified 36,548 copy number variant (CNVs) spanning 133.29 Mb of the buffalo genome, resulting in 2,100 copy number variant regions (CNVRs), with 1,993 shared CNVRs being found within the studied buffalo types. Analyzing CNVRs highlighted distinct genetic differentiation between river and swamp buffalo subspecies, verified by evolutionary tree and principal component analyses. Admixture analysis grouped buffaloes into river and swamp categories, with crossbred buffaloes displaying mixed ancestry. To identify candidate genes associated with milk production traits, we employed 3 approaches. First, we used Vst-based population differentiation, revealing 11 genes within CNVRs that exhibited significant divergence between different buffalo breeds, including genes linked to milk production traits. Second, expression quantitative loci (eQTL) analysis revealed differential expression of CNVR-driven genes (DECGs) associated with milk production traits. Notably, known milk production-related genes were among these DECGs, validating their relevance. Last, a genome-wide association study (GWAS) identified 3 CNVRs significantly linked to peak milk yield. Our study provides comprehensive genomic insights into buffalo populations and identifies candidate genes associated with milk production traits. These findings facilitate genetic breeding programs aimed at increasing milk yield and improving quality in this economically important livestock species.
Collapse
Affiliation(s)
- 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.
| | - 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
| | - Anqin Duan
- Guangxi Provincial Key Laboratory of Buffalo Genetics, Breeding and Reproduction Technology, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning 530001, China
| | - Xing-Rong Lu
- Guangxi Provincial Key Laboratory of Buffalo Genetics, Breeding and Reproduction Technology, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning 530001, China
| | - Hamdy Abdel-Shafy
- Department of Animal Production, Faculty of Agriculture, Cairo University, 12613, Giza, Egypt
| |
Collapse
|
2
|
Kumar H, Panigrahi M, G Strillacci M, Sonejita Nayak S, Rajawat D, Ghildiyal K, Bhushan B, Dutt T. Detection of genome-wide copy number variation in Murrah buffaloes. Anim Biotechnol 2023; 34:3783-3795. [PMID: 37381739 DOI: 10.1080/10495398.2023.2227670] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/30/2023]
Abstract
Riverine Buffaloes, especially the Murrah breed because of their adaptability to harsh climatic conditions, is farmed in many countries to convert low-quality feed into valuable dairy products and meat. Here, we investigated the copy number variations (CNVs) in 296 Murrah buffalo using the Axiom® Buffalo Genotyping Array 90K (Affymetrix, Santa Clara, CA, USA). The CNVs were detected on the autosomes, using the Copy Number Analysis Module (CNAM) using the univariate analysis. 7937 CNVs were detected in 279 Buffaloes, the average length of the CNVs was 119,048.87 bp that ranged between 7800 and 4,561,030 bp. These CNVs were accounting for 10.33% of the buffalo genome, which was comparable to cattle, sheep, and goat CNV analyses. Further, CNVs were merged and 1541 CNVRs were detected using the Bedtools-mergeBed command. 485 genes were annotated within 196 CNVRs that were identified in at least 10 animals of Murrah population. Out of these, 40 CNVRs contained 59 different genes that were associated with 69 different traits. Overall, the study identified a significant number of CNVs and CNVRs in the Murrah breed of buffalo, with a wide range of lengths and frequencies across the autosomes. The identified CNVRs contained genes associated with important traits related to production and reproduction, making them potentially important targets for future breeding and genetic improvement efforts.
Collapse
Affiliation(s)
- Harshit Kumar
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, India
| | - Manjit Panigrahi
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, India
| | - Maria G Strillacci
- Department of Veterinary Medicine and Animal Sciences, University of Milan, Lodi, Italy
| | | | - Divya Rajawat
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, India
| | - Kanika Ghildiyal
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, India
| | - Bharat Bhushan
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, India
| | - Triveni Dutt
- Livestock Production and Management Section, Indian Veterinary Research Institute, Izatnagar, India
| |
Collapse
|
3
|
Ahmad SF, Chandrababu Shailaja C, Vaishnav S, Kumar A, Gaur GK, Janga SC, Ahmad SM, Malla WA, Dutt T. Read-depth based approach on whole genome resequencing data reveals important insights into the copy number variation (CNV) map of major global buffalo breeds. BMC Genomics 2023; 24:616. [PMID: 37845620 PMCID: PMC10580622 DOI: 10.1186/s12864-023-09720-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2023] [Accepted: 10/05/2023] [Indexed: 10/18/2023] Open
Abstract
BACKGROUND Elucidating genome-wide structural variants including copy number variations (CNVs) have gained increased significance in recent times owing to their contribution to genetic diversity and association with important pathophysiological states. The present study aimed to elucidate the high-resolution CNV map of six different global buffalo breeds using whole genome resequencing data at two coverages (10X and 30X). Post-quality control, the sequence reads were aligned to the latest draft release of the Bubaline genome. The genome-wide CNVs were elucidated using a read-depth approach in CNVnator with different bin sizes. Adjacent CNVs were concatenated into copy number variation regions (CNVRs) in different breeds and their genomic coverage was elucidated. RESULTS Overall, the average size of CNVR was lower at 30X coverage, providing finer details. Most of the CNVRs were either deletion or duplication type while the occurrence of mixed events was lesser in number on a comparative basis in all breeds. The average CNVR size was lower at 30X coverage (0.201 Mb) as compared to 10X (0.013 Mb) with the finest variants in Banni buffaloes. The maximum number of CNVs was observed in Murrah (2627) and Pandharpuri (25,688) at 10X and 30X coverages, respectively. Whereas the minimum number of CNVs were scored in Surti at both coverages (2092 and 17,373). On the other hand, the highest and lowest number of CNVRs were scored in Jaffarabadi (833 and 10,179 events) and Surti (783 and 7553 events) at both coverages. Deletion events overnumbered duplications in all breeds at both coverages. Gene profiling of common overlapped genes and longest CNVRs provided important insights into the evolutionary history of these breeds and indicate the genomic regions under selection in respective breeds. CONCLUSION The present study is the first of its kind to elucidate the high-resolution CNV map in major buffalo populations using a read-depth approach on whole genome resequencing data. The results revealed important insights into the divergence of major global buffalo breeds along the evolutionary timescale.
Collapse
Affiliation(s)
- Sheikh Firdous Ahmad
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, Uttar Pradesh, 243122, India.
| | - Celus Chandrababu Shailaja
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, Uttar Pradesh, 243122, India
| | - Sakshi Vaishnav
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, Uttar Pradesh, 243122, India
| | - Amit Kumar
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, Uttar Pradesh, 243122, India
| | - Gyanendra Kumar Gaur
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, Uttar Pradesh, 243122, India
| | - Sarath Chandra Janga
- Luddy School of Informatics, Computing & Engineering, Indiana University Indianapolis (IUI), Indianapolis, 46202, USA
| | - Syed Mudasir Ahmad
- Division of Animal Biotechnology, Faculty of Veterinary Sciences and AH, Sher-e-Kashmir University of Agricultural Sciences and Technology, Srinagar, Jammu & Kashmir, 190006, India.
| | - Waseem Akram Malla
- Division of Veterinary Biotechnology, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, Uttar Pradesh, 243122, India
| | - Triveni Dutt
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, Uttar Pradesh, 243122, India
| |
Collapse
|