1
|
Deng TX, Ma XY, Duan A, Lu XR, Abdel-Shafy H. Genome-wide copy number variant analysis reveals candidate genes associated with milk production traits in water buffalo (Bubalus bubalis). J Dairy Sci 2024; 107:7022-7037. [PMID: 38762109 DOI: 10.3168/jds.2023-24614] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/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 variants (CNV) spanning 133.29 Mb of the buffalo genome, resulting in 2,100 CNV regions (CNVR), with 1,993 shared CNVR being found within the studied buffalo types. Analyzing CNVR 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 CNVR that exhibited significant divergence between different buffalo breeds, including genes linked to milk production traits. Second, expression quantitative loci analysis revealed differentially expressed CNVR-derived genes (DECG) associated with milk production traits. Notably, known milk production-related genes were among these DECG, validating their relevance. Last, a GWAS identified 3 CNVR 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
|
Guo L, Dao L, Liu B, Wang J, Liu Z, Ma F, Morigen B, Chang C, Bai Y, Guo Y, Shi C, Cao J, Zhang W. Development and application of a 1K functional liquid chip for lactation performance in Bactrian camels. Front Vet Sci 2024; 11:1359923. [PMID: 39021409 PMCID: PMC11253134 DOI: 10.3389/fvets.2024.1359923] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Accepted: 06/12/2024] [Indexed: 07/20/2024] Open
Abstract
Introduction The advancement of high-throughput, high-quality, flexible, and cost-effective genotyping platforms is crucial for the progress of dairy breeding in Bactrian camels. This study focuses on developing and evaluating a 1K functional liquid single nucleotide polymorphism (SNP) array specifically designed for milk performance in Bactrian camels. Methods We utilized RNA sequencing data from 125 lactating camels to identify and select 1,002 loci associated with milk production traits for inclusion in the SNP array. The array's performance was then assessed using 24 randomly selected camels. Additionally, the array was employed to genotype 398 individuals, which allowed for population validation to assess the polymorphism of SNP sites. Results The SNP array demonstrated high overall SNP call rates (> 99%) and a remarkable 100% consistency in genotyping. Population validation results indicate that camels from six breeding areas in Northwest China share a similar genetic background regarding lactation functionality. Discussion This study highlights the potential of the SNP array to accelerate the breeding process of lactating Bactrian camels and provides a robust technical foundation for improving lactation performance.
Collapse
Affiliation(s)
- Lili Guo
- College of Life Science, Inner Mongolia Agricultural University, Hohhot, China
- Inner Mongolia Engineering Research Center of Genomic Big Data for Agriculture, Inner Mongolia Agricultural University, Hohhot, China
- Inner Mongolia Autonomous Region Key Laboratory of Biomanufacturing, Hohhot, China
| | - Lema Dao
- Bactrian Camel Institute of Alsha, Bayanhot, China
| | - Bin Liu
- Inner Mongolia Bionew Technology Co., Ltd., Hohhot, China
| | - Jingyu Wang
- Bactrian Camel Institute of Alsha, Bayanhot, China
| | - Zaixia Liu
- Inner Mongolia Engineering Research Center of Genomic Big Data for Agriculture, Inner Mongolia Agricultural University, Hohhot, China
- College of Animal Science, Inner Mongolia Agricultural University, Hohhot, China
| | - Fengying Ma
- Inner Mongolia Engineering Research Center of Genomic Big Data for Agriculture, Inner Mongolia Agricultural University, Hohhot, China
- College of Animal Science, Inner Mongolia Agricultural University, Hohhot, China
| | - Bielige Morigen
- Animal Disease Prevention and Control Center of Alsha, Bayanhot, China
| | - Chencheng Chang
- Inner Mongolia Engineering Research Center of Genomic Big Data for Agriculture, Inner Mongolia Agricultural University, Hohhot, China
- College of Animal Science, Inner Mongolia Agricultural University, Hohhot, China
| | - Yinbatu Bai
- Inner Mongolia Engineering Research Center of Genomic Big Data for Agriculture, Inner Mongolia Agricultural University, Hohhot, China
- College of Animal Science, Inner Mongolia Agricultural University, Hohhot, China
| | - Yaqiang Guo
- Inner Mongolia Engineering Research Center of Genomic Big Data for Agriculture, Inner Mongolia Agricultural University, Hohhot, China
- College of Animal Science, Inner Mongolia Agricultural University, Hohhot, China
| | - Caixia Shi
- Inner Mongolia Engineering Research Center of Genomic Big Data for Agriculture, Inner Mongolia Agricultural University, Hohhot, China
- College of Animal Science, Inner Mongolia Agricultural University, Hohhot, China
| | - Junwei Cao
- College of Life Science, Inner Mongolia Agricultural University, Hohhot, China
- Inner Mongolia Autonomous Region Key Laboratory of Biomanufacturing, Hohhot, China
| | - Wenguang Zhang
- College of Life Science, Inner Mongolia Agricultural University, Hohhot, China
- Inner Mongolia Engineering Research Center of Genomic Big Data for Agriculture, Inner Mongolia Agricultural University, Hohhot, China
- College of Animal Science, Inner Mongolia Agricultural University, Hohhot, China
| |
Collapse
|
3
|
Khazaei-Koohpar H, Gholizadeh M, Hafezian SH, Esmaeili-Fard SM. Weighted single-step genome-wide association study for direct and maternal genetic effects associated with birth and weaning weights in sheep. Sci Rep 2024; 14:13120. [PMID: 38849438 PMCID: PMC11161479 DOI: 10.1038/s41598-024-63974-0] [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/23/2023] [Accepted: 06/04/2024] [Indexed: 06/09/2024] Open
Abstract
Body weight is an important economic trait for sheep meat production, and its genetic improvement is considered one of the main goals in the sheep breeding program. Identifying genomic regions that are associated with growth-related traits accelerates the process of animal breeding through marker-assisted selection, which leads to increased response to selection. In this study, we conducted a weighted single-step genome-wide association study (WssGWAS) to identify potential candidate genes for direct and maternal genetic effects associated with birth weight (BW) and weaning weight (WW) in Baluchi sheep. The data used in this research included 13,408 birth and 13,170 weaning records collected at Abbas-Abad Baluchi Sheep Breeding Station, Mashhad-Iran. Genotypic data of 94 lambs genotyped by Illumina 50K SNP BeadChip for 54,241 markers were used. The proportion of variance explained by genomic windows was calculated by summing the variance of SNPs within 1 megabase (Mb). The top 10 window genomic regions explaining the highest percentages of additive and maternal genetic variances were selected as candidate window genomic regions associated with body weights. Our findings showed that for BW, the top-ranked genomic regions (1 Mb windows) explained 4.30 and 4.92% of the direct additive and maternal genetic variances, respectively. The direct additive genetic variance explained by the genomic window regions varied from 0.31 on chromosome 1 to 0.59 on chromosome 8. The highest (0.84%) and lowest (0.32%) maternal genetic variances were explained by genomic windows on chromosome 10 and 17, respectively. For WW, the top 10 genomic regions explained 6.38 and 5.76% of the direct additive and maternal genetic variances, respectively. The highest and lowest contribution of direct additive genetic variances were 1.37% and 0.42%, respectively, both explained by genomic regions on chromosome 2. For maternal effects on WW, the highest (1.38%) and lowest (0.41%) genetic variances were explained by genomic windows on chromosome 2. Further investigation of these regions identified several possible candidate genes associated with body weight. Gene ontology analysis using the DAVID database identified several functional terms, such as translation repressor activity, nucleic acid binding, dehydroascorbic acid transporter activity, growth factor activity and SH2 domain binding.
Collapse
Affiliation(s)
- Hava Khazaei-Koohpar
- Department of Animal Science and Fisheries, Sari Agricultural Sciences and Natural Resources University (SANRU), Sari, Iran
| | - Mohsen Gholizadeh
- Department of Animal Science and Fisheries, Sari Agricultural Sciences and Natural Resources University (SANRU), Sari, Iran.
| | - Seyed Hasan Hafezian
- Department of Animal Science and Fisheries, Sari Agricultural Sciences and Natural Resources University (SANRU), Sari, Iran
| | | |
Collapse
|
4
|
Igoshin AV, Mishakova TM, Aitnazarov RB, Ilina AV, Larkin DM, Yudin NS. Association of three single nucleotide polymorphisms in the LPIN1 gene with milk production traits in cows of the Yaroslavl breed. Vavilovskii Zhurnal Genet Selektsii 2024; 28:117-125. [PMID: 38465251 PMCID: PMC10917680 DOI: 10.18699/vjgb-24-14] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2023] [Revised: 01/15/2024] [Accepted: 01/15/2024] [Indexed: 03/12/2024] Open
Abstract
Lipin-1 is a member of the evolutionarily conserved family of proteins and is expressed predominantly in adipose tissue and skeletal muscle. On the one hand, lipin-1 is an enzyme that catalyzes the dephosphorylation of phosphatidic acid to diacylglycerol (DAG) and thus participates in the metabolic pathways of biosynthesis of storage lipids in the cell, membrane phospholipids, and intracellular signaling molecules. On the other hand, lipin-1 is able to be transported from the cytoplasm to the nucleus and is a coactivator of lipid metabolism gene transcription. It was shown, using the analysis of single nucleotide polymorphism (SNP) associations, that the lipin-1 coding gene (LPIN1) is a promising candidate gene for milk production traits in Holstein and Brown Swiss cows. However, it is unclear how much of its effect depends on the breed. The Yaroslavl dairy cattle breed was created in the 18-19 centuries in Russia by breeding northern Great Russian cattle, which were short and poor productive, but well adapted to local climatic conditions and bad food base. It was shown by whole genome genotyping and sequencing that the Yaroslavl breed has unique genetics compared to Russian and other cattle breeds. The aim of the study was to assess the frequency of alleles and genotypes of three SNPs in the LPIN1 gene and to study the association of these SNPs with milk production traits in Yaroslavl cows. Blood samples from 142 cows of the Yaroslavl breed were obtained from two farms in the Yaroslavl region. Genotyping of SNPs was carried out by polymerase chain reaction-restriction fragment length polymorphism method. Associations of SNPs with 305-day milk yield, fat yield, fat percentages, protein yield, and protein percentages were studied from the first to the fourth lactation. Statistical tests were carried out using a mixed linear model, taking into account the relationship between individuals. We identified three SNPs - rs110871255, rs207681322 and rs109039955 with a frequency of a rare allele of 0.042-0.261 in Yaroslavl cows. SNP rs110871255 was associated with fat yield during the third and fourth lactations. SNP rs207681322 was associated with milk yield for the second, third and fourth lactations, as well as protein yield for the third lactation. Thus, we identified significant associations of SNPs rs207681322 and rs110871255 in the LPIN1 gene with a number of milk production traits during several lactations in Yaroslavl cows.
Collapse
Affiliation(s)
- A V Igoshin
- Institute of Cytology and Genetics of the Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia
| | - T M Mishakova
- Institute of Cytology and Genetics of the Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia
| | - R B Aitnazarov
- Institute of Cytology and Genetics of the Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia
| | - A V Ilina
- Federal Williams Research Center for Forage Production and Agroecology, Scientific Research Institute of Livestock Breeding and Forage Production, Yaroslavl Region, Russia
| | - D M Larkin
- Royal Veterinary College, University of London, London, United Kingdom
| | - N S Yudin
- Institute of Cytology and Genetics of the Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia
| |
Collapse
|
5
|
Pallotti S, Picciolini M, Antonini M, Renieri C, Napolioni V. Genome-wide scan for runs of homozygosity in South American Camelids. BMC Genomics 2023; 24:470. [PMID: 37605116 PMCID: PMC10440933 DOI: 10.1186/s12864-023-09547-3] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2023] [Accepted: 07/31/2023] [Indexed: 08/23/2023] Open
Abstract
BACKGROUND Alpaca (Vicugna pacos), llama (Lama glama), vicugna (Vicugna vicugna) and guanaco (Lama guanicoe), are the camelid species distributed over the Andean high-altitude grasslands, the Altiplano, and the Patagonian arid steppes. Despite the wide interest on these animals, most of the loci under selection are still unknown. Using whole-genome sequencing (WGS) data we investigated the occurrence and the distribution of Runs Of Homozygosity (ROHs) across the South American Camelids (SACs) genome to identify the genetic relationship between the four species and the potential signatures of selection. RESULTS A total of 37 WGS samples covering the four species was included in the final analysis. The multi-dimensional scaling approach showed a clear separation between the four species; however, admixture analysis suggested a strong genetic introgression from vicugna and llama to alpaca. Conversely, very low genetic admixture of the guanaco with the other SACs was found. The four species did not show significant differences in the number, length of ROHs (100-500 kb) and genomic inbreeding values. Longer ROHs (> 500 kb) were found almost exclusively in alpaca. Seven overlapping ROHs were shared by alpacas, encompassing nine loci (FGF5, LOC107034918, PRDM8, ANTXR2, LOC102534792, BSN, LOC116284892, DAG1 and RIC8B) while nine overlapping ROHs were found in llama with twenty-five loci annotated (ERC2, FZD9, BAZ1B, BCL7B, LOC116284208, TBL2, MLXIPL, PHF20, TRNAD-AUC, LOC116284365, RBM39, ARFGEF2, DCAF5, EXD2, HSPB11, LRRC42, LDLRAD1, TMEM59, LOC107033213, TCEANC2, LOC102545169, LOC116278408, SMIM15, NDUFAF2 and RCOR1). Four overlapping ROHs, with three annotated loci (DLG1, KAT6B and PDE4D) and three overlapping ROHs, with seven annotated genes (ATP6V1E1, BCL2L13, LOC116276952, BID, KAT6B, LOC116282667 and LOC107034552), were detected for vicugna and guanaco, respectively. CONCLUSIONS The signatures of selection revealed genomic areas potentially selected for production traits as well as for natural adaptation to harsh environment. Alpaca and llama hint a selection driven by environment as well as by farming purpose while vicugna and guanaco showed selection signals for adaptation to harsh environment. Interesting, signatures of selection on KAT6B gene were identified for both vicugna and guanaco, suggesting a positive effect on wild populations fitness. Such information may be of interest to further ecological and animal production studies.
Collapse
Affiliation(s)
- Stefano Pallotti
- Genomic And Molecular Epidemiology (GAME) Lab, School of Biosciences and Veterinary Medicine, University of Camerino, Camerino, Italy.
| | | | - Marco Antonini
- Italian National Agency for New Technologies, Energy and Sustainable Development (ENEA), Roma, Italy
| | - Carlo Renieri
- School of Pharmacy and Health Products, University of Camerino, Camerino, Italy
| | - Valerio Napolioni
- Genomic And Molecular Epidemiology (GAME) Lab, School of Biosciences and Veterinary Medicine, University of Camerino, Camerino, Italy
| |
Collapse
|
6
|
Ahmad SM, Bhat B, Manzoor Z, Dar MA, Taban Q, Ibeagha-Awemu EM, Shabir N, Hussain MI, Shah RA, Ganai NA. Genome wide expression analysis of circular RNAs in mammary epithelial cells of cattle revealed difference in milk synthesis. PeerJ 2022; 10:e13029. [PMID: 35251787 PMCID: PMC8896013 DOI: 10.7717/peerj.13029] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2021] [Accepted: 02/08/2022] [Indexed: 01/11/2023] Open
Abstract
Milk is an excellent source of nutrients for humans. Therefore, in order to enhance the quality and production of milk in cattle, it is interesting to examine the underlying mechanisms. A number of new investigations and research have found that, circRNA; a specific class of non-coding RNAs, is linked with the development of mammary gland and lactation. In the present study, genome wide identification and expression of the circRNAs in mammary epithelial cells of two distinct cattle breeds viz Jersey and Kashmiri at peak lactation was conducted. We reported 1554 and 1286 circRNA in Jersey and Kashmiri cattle, respectively, with 21 circRNAs being differentially expressed in the two breeds. The developmental genes of the established differentially expressed circRNAs were found to be largely enriched in antioxidant activity, progesterone, estradiol, lipid, growth hormone, and drug response. Certain pathways like MAPK, IP3K and immune response pathways were found significantly enriched in KEGG analysis. These results add to our understanding of the controlling mechanisms connected with the lactation process, as well as the function of circRNAs in bovine milk synthesis. Additionally, the comparative analysis of differentially expressed circRNAs showed significant conservation across different species.
Collapse
Affiliation(s)
- Syed Mudasir Ahmad
- Division of Animal Biotechnology, Faculty of Veterinary Sciences, Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir, India, Srinagar, Jammu and Kashmir, India
| | - Basharat Bhat
- Division of Animal Biotechnology, Faculty of Veterinary Sciences, Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir, India, Srinagar, Jammu and Kashmir, India
| | - Zainab Manzoor
- Department of Clinical Biochemistry, University of Kashmir, Srinagar, Jammu and Kashmir, India
| | - Mashooq Ahmad Dar
- Division of Animal Biotechnology, Faculty of Veterinary Sciences, Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir, India, Srinagar, Jammu and Kashmir, India,Department of Clinical Biochemistry, University of Kashmir, Srinagar, Jammu and Kashmir, India
| | - Qamar Taban
- Division of Animal Biotechnology, Faculty of Veterinary Sciences, Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir, India, Srinagar, Jammu and Kashmir, India
| | - Eveline M. Ibeagha-Awemu
- Sherbrooke Research and Development Centre, Agriculture and Agri-Food Canada, Sherbrooke, Canada
| | - Nadeem Shabir
- Division of Animal Biotechnology, Faculty of Veterinary Sciences, Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir, India, Srinagar, Jammu and Kashmir, India
| | - Mohd Isfaqul Hussain
- Division of Veterinary Microbiology, SKUAST-Kashmir, Srinagar, Jammu and Kashmir, India
| | - Riaz A. Shah
- Division of Animal Biotechnology, Faculty of Veterinary Sciences, Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir, India, Srinagar, Jammu and Kashmir, India
| | - Nazir A. Ganai
- Division of Animal Biotechnology, Faculty of Veterinary Sciences, Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir, India, Srinagar, Jammu and Kashmir, India
| |
Collapse
|