1
|
Wang Y, Ma J, Wang J, Zhang L, Xu L, Chen Y, Zhu B, Wang Z, Gao H, Li J, Gao X. Genome-Wide Detection of Copy Number Variations and Their Potential Association with Carcass and Meat Quality Traits in Pingliang Red Cattle. Int J Mol Sci 2024; 25:5626. [PMID: 38891814 PMCID: PMC11172001 DOI: 10.3390/ijms25115626] [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: 04/30/2024] [Revised: 05/16/2024] [Accepted: 05/17/2024] [Indexed: 06/21/2024] Open
Abstract
Copy number variation (CNV) serves as a significant source of genetic diversity in mammals and exerts substantial effects on various complex traits. Pingliang red cattle, an outstanding indigenous resource in China, possess remarkable breeding value attributed to their tender meat and superior marbling quality. However, the genetic mechanisms influencing carcass and meat quality traits in Pingliang red cattle are not well understood. We generated a comprehensive genome-wide CNV map for Pingliang red cattle using the GGP Bovine 100K SNP chip. A total of 755 copy number variable regions (CNVRs) spanning 81.03 Mb were identified, accounting for approximately 3.24% of the bovine autosomal genome. Among these, we discovered 270 potentially breed-specific CNVRs in Pingliang red cattle, including 143 gains, 73 losses, and 54 mixed events. Functional annotation analysis revealed significant associations between these specific CNVRs and important traits such as carcass and meat quality, reproduction, exterior traits, growth traits, and health traits. Additionally, our network and transcriptome analysis highlighted CACNA2D1, CYLD, UBXN2B, TG, NADK, and ITGA9 as promising candidate genes associated with carcass weight and intramuscular fat deposition. The current study presents a genome-wide CNV map in Pingliang red cattle, highlighting breed-specific CNVRs, and transcriptome findings provide valuable insights into the underlying genetic characteristics of Pingliang red cattle. These results offer potential avenues for enhancing meat quality through a targeted breeding program.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | | | | | - Junya Li
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China
| | - Xue Gao
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China
| |
Collapse
|
2
|
Capra E, Lazzari B, Milanesi M, Nogueira GP, Garcia JF, Utsunomiya YT, Ajmone-Marsan P, Stella A. Comparison between indicine and taurine cattle DNA methylation reveals epigenetic variation associated to differences in morphological adaptive traits. Epigenetics 2023; 18:2163363. [PMID: 36600398 PMCID: PMC9980582 DOI: 10.1080/15592294.2022.2163363] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
Indicine and taurine subspecies present distinct morphological traits as a consequence of environmental adaptation and artificial selection. Although the two subspecies have been characterized and compared at genome-wide level and at specific loci, their epigenetic diversity has not yet been explored. In this work, Reduced Representation Bisulphite Sequencing (RRBS) profiling of the taurine Angus (A) and indicine Nellore (N) cattle breeds was applied to identify methylation differences between the two subspecies. Genotyping by sequencing (GBS) of the same animals was performed to detect single nucleotide polymorphisms (SNPs) at cytosines in CpG dinucleotides and remove them from the differential methylation analysis. A total of 660,845 methylated cytosines were identified within the CpG context (CpGs) across the 10 animals sequenced (5 N and 5 A). A total of 25,765 of these were differentially methylated (DMCs). Most DMCs clustered in CpG stretches nearby genes involved in cellular and anatomical structure morphogenesis. Also, sequences flanking DMC were enriched in SNPs compared to all other CpGs, either methylated or unmethylated in the two subspecies. Our data suggest a contribution of epigenetics to the regulation and divergence of anatomical morphogenesis in the two subspecies relevant for cattle evolution and sub-species differentiation and adaptation.
Collapse
Affiliation(s)
- E Capra
- Institute of Agricultural Biology and Biotechnology, National Research Council IBBA CNR, Lodi, Italy
| | - B Lazzari
- Institute of Agricultural Biology and Biotechnology, National Research Council IBBA CNR, Lodi, Italy
| | - M Milanesi
- School of Veterinary Medicine, Araçatuba, Department of Production and Animal Health, São Paulo State University (Unesp), Araçatuba, Brazil.,International Atomic Energy Agency, Collaborating Centre on Animal Genomics and Bioinformatics, Araçatuba, Brazil.,Department for Innovation in Biological, Agro-Food and Forest Systems (DIBAF), University of Tuscia, Viterbo, Italy
| | - G P Nogueira
- School of Veterinary Medicine, Araçatuba, Department of Production and Animal Health, São Paulo State University (Unesp), Araçatuba, Brazil
| | - J F Garcia
- School of Veterinary Medicine, Araçatuba, Department of Production and Animal Health, São Paulo State University (Unesp), Araçatuba, Brazil.,International Atomic Energy Agency, Collaborating Centre on Animal Genomics and Bioinformatics, Araçatuba, Brazil
| | - Y T Utsunomiya
- School of Veterinary Medicine, Araçatuba, Department of Production and Animal Health, São Paulo State University (Unesp), Araçatuba, Brazil
| | - P Ajmone-Marsan
- Department of Animal Science, Food and Nutrition - DIANA, and Nutrigenomics and Proteomics Research Center - PRONUTRIGEN, Università Cattolica del Sacro Cuore, Piacenza, Italy
| | - A Stella
- Institute of Agricultural Biology and Biotechnology, National Research Council IBBA CNR, Lodi, Italy
| |
Collapse
|
3
|
Singh VK, Singh S, Nandhini PB, Bhatia AK, Dixit SP, Ganguly I. Comparative genomic diversity analysis of copy number variations (CNV) in indicine and taurine cattle thriving in Europe and Indian subcontinent. Anim Biotechnol 2023; 34:3483-3494. [PMID: 36592947 DOI: 10.1080/10495398.2022.2162910] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Copy number variations (CNVs) include deletions, duplications, and insertions that are larger than 50 bp in size causing structural variation responsible for diversity, adaptation, and breed development. Indian cattle breeds are highly diverse from the taurine breeds. The pattern of CNVRs in 191 animals belonging to 39 cattle breeds (four Indicine and 35 Taurine) was studied based on Illumina 777K BovineHD chip data. The Indicine breeds revealed 2590 CNVs and 335 copy number variation regions (CNVRs) in autosomes. Out of the identified CNVs, 50 were found to be novel. Structure analysis revealed admixed nature of Siri. Neighbor joining tree from CNVR data showed that hot (Kankrej and Hallikar) and cold (Ladakhi and Siri) adapted cattle breeds clustered separately. CNVR of Indian and European breeds revealed that Balkan and Italian breeds of Podolian group are admixed with Indian cattle breeds corroborating indicine introgression (6.1-13.5%). CNVRs spanning the regions of olfactory receptors and immune system genes were identified. AMOVA revealed 9% variation among populations which is 2% greater than SNP based studies showing higher inclusion of variation by CNVR. Detailed analysis of CNVs/CNVRs in Indian cattle adapted to hot and cold climate, and their diversity among worldwide cattle is presented in this study.
Collapse
Affiliation(s)
- V K Singh
- Animal Genetics and Breeding Division, ICAR-National Dairy Research Institute, Karnal, India
| | - S Singh
- Animal Genetics Division, ICAR-National Bureau of Animal Genetic Resources, Karnal, India
| | - P B Nandhini
- Animal Genetics and Breeding Division, ICAR-National Dairy Research Institute, Karnal, India
| | - A K Bhatia
- Animal Genetic Resources Division, ICAR-National Bureau of Animal Genetic Resources, Karnal, India
| | - S P Dixit
- Animal Genetics Division, ICAR-National Bureau of Animal Genetic Resources, Karnal, India
| | - I Ganguly
- Animal Genetics Division, ICAR-National Bureau of Animal Genetic Resources, Karnal, India
| |
Collapse
|
4
|
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
|
5
|
Li S, Lei H, Li J, Sun A, Ahmed Z, Duan H, Chen L, Zhang B, Lei C, Yi K. Analysis of genetic diversity and selection signals in Chaling cattle of southern China using whole-genome scan. Anim Genet 2023; 54:284-294. [PMID: 36864643 DOI: 10.1111/age.13305] [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: 10/12/2022] [Revised: 01/12/2023] [Accepted: 01/30/2023] [Indexed: 03/04/2023]
Abstract
China has diversified resources of indigenous cattle, which are classified into Northern, Central, and Southern groups according to their geographical distribution. Chaling cattle belong to Southern group. This breed is famous for the production of good quality meat with elite meat grades. To analyze the genetic diversity of Chaling cattle, 20 samples were sequenced using whole-genome resequencing technology, along with 138 published whole-genome sequencing data of Indian indicine cattle, Chinese indicine cattle, East Asian taurine cattle, Eurasian taurine cattle, and European taurine cattle as control. It was found that Chaling cattle originated from Chinese indicine cattle. The genetic diversity of Chaling cattle is higher than that of Indian indicine cattle, East Asian taurine cattle, Eurasian taurine cattle, and European taurine cattle, but lower than that of Chinese indicine cattle and Xiangxi cattle. Annotating the selection signals obtained by composite likelihood ratio, θπ, FST , π-ratio, and XP-EHH methods, several genes associated with immunity, heat tolerance, reproduction, growth, and meat quality showed strong selection signals. In general, this study provides a theoretical basis for analyzing the genetic mechanism of Chaling cattle with excellent adaptability, rough feeding tolerance, good immune performance, and good meat quality. This work lays a foundation for genetic breeding of Chaling cattle in future.
Collapse
Affiliation(s)
- Shuang Li
- Hunan Institute of Animal and Veterinary Science, Changsha, Hunan, China.,Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
| | - Hong Lei
- Hunan Institute of Animal and Veterinary Science, Changsha, Hunan, China
| | - Jianbo Li
- Hunan Institute of Animal and Veterinary Science, Changsha, Hunan, China
| | - Ao Sun
- Hunan Institute of Animal and Veterinary Science, Changsha, Hunan, China
| | - Zulfiqar Ahmed
- Department of Livestock and Poultry Production, Faculty of Veterinary and Animal Science, University of Poonch Rawalakot, Azad Kashmir, Pakistan
| | - Hongfeng Duan
- Hunan Institute of Animal and Veterinary Science, Changsha, Hunan, China
| | - Lin Chen
- Chaling County Agricultural Development Corporation Ltd, Chaling, Hunan, China
| | - Baizhong Zhang
- Hunan Institute of Animal and Veterinary Science, Changsha, Hunan, China
| | - Chuzhao Lei
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
| | - Kangle Yi
- Hunan Institute of Animal and Veterinary Science, Changsha, Hunan, China
| |
Collapse
|
6
|
Nguyen TV, Vander Jagt CJ, Wang J, Daetwyler HD, Xiang R, Goddard ME, Nguyen LT, Ross EM, Hayes BJ, Chamberlain AJ, MacLeod IM. In it for the long run: perspectives on exploiting long-read sequencing in livestock for population scale studies of structural variants. Genet Sel Evol 2023; 55:9. [PMID: 36721111 PMCID: PMC9887926 DOI: 10.1186/s12711-023-00783-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2022] [Accepted: 01/23/2023] [Indexed: 02/02/2023] Open
Abstract
Studies have demonstrated that structural variants (SV) play a substantial role in the evolution of species and have an impact on Mendelian traits in the genome. However, unlike small variants (< 50 bp), it has been challenging to accurately identify and genotype SV at the population scale using short-read sequencing. Long-read sequencing technologies are becoming competitively priced and can address several of the disadvantages of short-read sequencing for the discovery and genotyping of SV. In livestock species, analysis of SV at the population scale still faces challenges due to the lack of resources, high costs, technological barriers, and computational limitations. In this review, we summarize recent progress in the characterization of SV in the major livestock species, the obstacles that still need to be overcome, as well as the future directions in this growing field. It seems timely that research communities pool resources to build global population-scale long-read sequencing consortiums for the major livestock species for which the application of genomic tools has become cost-effective.
Collapse
Affiliation(s)
- Tuan V. Nguyen
- grid.452283.a0000 0004 0407 2669Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC 3083 Australia
| | - Christy J. Vander Jagt
- grid.452283.a0000 0004 0407 2669Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC 3083 Australia
| | - Jianghui Wang
- grid.452283.a0000 0004 0407 2669Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC 3083 Australia
| | - Hans D. Daetwyler
- grid.452283.a0000 0004 0407 2669Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC 3083 Australia ,grid.1018.80000 0001 2342 0938School of Applied Systems Biology, La Trobe University, Bundoora, VIC 3083 Australia
| | - Ruidong Xiang
- grid.452283.a0000 0004 0407 2669Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC 3083 Australia ,grid.1008.90000 0001 2179 088XFaculty of Veterinary & Agricultural Science, The University of Melbourne, Parkville, VIC 3052 Australia
| | - Michael E. Goddard
- grid.452283.a0000 0004 0407 2669Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC 3083 Australia ,grid.1008.90000 0001 2179 088XFaculty of Veterinary & Agricultural Science, The University of Melbourne, Parkville, VIC 3052 Australia
| | - Loan T. Nguyen
- grid.1003.20000 0000 9320 7537Queensland Alliance for Agriculture and Food Innovation, University of Queensland, St Lucia, QLD 4072 Australia
| | - Elizabeth M. Ross
- grid.1003.20000 0000 9320 7537Queensland Alliance for Agriculture and Food Innovation, University of Queensland, St Lucia, QLD 4072 Australia
| | - Ben J. Hayes
- grid.1003.20000 0000 9320 7537Queensland Alliance for Agriculture and Food Innovation, University of Queensland, St Lucia, QLD 4072 Australia
| | - Amanda J. Chamberlain
- grid.452283.a0000 0004 0407 2669Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC 3083 Australia ,grid.1018.80000 0001 2342 0938School of Applied Systems Biology, La Trobe University, Bundoora, VIC 3083 Australia
| | - Iona M. MacLeod
- grid.452283.a0000 0004 0407 2669Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC 3083 Australia
| |
Collapse
|
7
|
Braga LG, Chud TCS, Watanabe RN, Savegnago RP, Sena TM, do Carmo AS, Machado MA, Panetto JCDC, da Silva MVGB, Munari DP. Identification of copy number variations in the genome of Dairy Gir cattle. PLoS One 2023; 18:e0284085. [PMID: 37036840 PMCID: PMC10085049 DOI: 10.1371/journal.pone.0284085] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2022] [Accepted: 03/23/2023] [Indexed: 04/11/2023] Open
Abstract
Studying structural variants that can control complex traits is relevant for dairy cattle production, especially for animals that are tolerant to breeding conditions in the tropics, such as the Dairy Gir cattle. This study identified and characterized high confidence copy number variation regions (CNVR) in the Gir breed genome. A total of 38 animals were whole-genome sequenced, and 566 individuals were genotyped with a high-density SNP panel, among which 36 animals had both sequencing and SNP genotyping data available. Two sets of high confidence CNVR were established: one based on common CNV identified in the studied population (CNVR_POP), and another with CNV identified in sires with both sequence and SNP genotyping data available (CNVR_ANI). We found 10 CNVR_POP and 45 CNVR_ANI, which covered 1.05 Mb and 4.4 Mb of the bovine genome, respectively. Merging these CNV sets for functional analysis resulted in 48 unique high confidence CNVR. The overlapping genes were previously related to embryonic mortality, environmental adaptation, evolutionary process, immune response, longevity, mammary gland, resistance to gastrointestinal parasites, and stimuli recognition, among others. Our results contribute to a better understanding of the Gir breed genome. Moreover, the CNV identified in this study can potentially affect genes related to complex traits, such as production, health, and reproduction.
Collapse
Affiliation(s)
- Larissa G Braga
- Departamento de Engenharia e Ciências Exatas, Universidade Estadual Paulista, Jaboticabal, São Paulo, Brazil
| | - Tatiane C S Chud
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, Ontario, Canada
| | - Rafael N Watanabe
- Departamento de Engenharia e Ciências Exatas, Universidade Estadual Paulista, Jaboticabal, São Paulo, Brazil
| | - Rodrigo P Savegnago
- Department of Animal Science, Michigan State University, East Lansing, Michigan, United States of America
| | - Thomaz M Sena
- Departamento de Engenharia e Ciências Exatas, Universidade Estadual Paulista, Jaboticabal, São Paulo, Brazil
| | - Adriana S do Carmo
- Departamento de Zootecnia, Universidade Federal de Goiás, Goiânia, Goiás, Brazil
| | | | | | | | - Danísio P Munari
- Departamento de Engenharia e Ciências Exatas, Universidade Estadual Paulista, Jaboticabal, São Paulo, Brazil
| |
Collapse
|
8
|
Identification and Characterization of Copy Number Variations Regions in West African Taurine Cattle. Animals (Basel) 2022; 12:ani12162130. [PMID: 36009719 PMCID: PMC9405125 DOI: 10.3390/ani12162130] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2022] [Revised: 07/29/2022] [Accepted: 08/16/2022] [Indexed: 11/17/2022] Open
Abstract
A total of 106 West African taurine cattle belonging to the Lagunaire breed of Benin (33), the N’Dama population of Burkina Faso (48), and N’Dama cattle sampled in Congo (25) were analyzed for Copy Number Variations (CNVs) using the BovineHDBeadChip of Illumina and two different CNV calling programs: PennCNV and QuantiSNP. Furthermore, 89 West African zebu samples (Bororo cattle of Mali and Zebu Peul sampled in Benin and Burkina Faso) were used as an outgroup to ensure that analyses reflect the taurine cattle genomic background. Analyses identified 307 taurine-specific CNV regions (CNVRs), covering about 56 Mb on all bovine autosomes. Gene annotation enrichment analysis identified a total of 840 candidate genes on 168 taurine-specific CNVRs. Three different statistically significant functional term annotation clusters (from ACt1 to ACt3) involved in the immune function were identified: ACt1 includes genes encoding lipocalins, proteins involved in the modulation of immune response and allergy; ACt2 includes genes encoding coding B-box-type zinc finger proteins and butyrophilins, involved in innate immune processes; and Act3 includes genes encoding lectin receptors, involved in the inflammatory responses to pathogens and B- and T-cell differentiation. The overlap between taurine-specific CNVRs and QTL regions associated with trypanotolerant response and tick-resistance was relatively low, suggesting that the mechanisms underlying such traits may not be determined by CNV alterations. However, four taurine-specific CNVRs overlapped with QTL regions associated with both traits on BTA23, therefore suggesting that CNV alterations in major histocompatibility complex (MHC) genes can partially explain the existence of genetic mechanisms shared between trypanotolerance and tick resistance in cattle. This research contributes to the understanding of the genomic features of West African taurine cattle.
Collapse
|
9
|
Zhou Y, Yang L, Han X, Han J, Hu Y, Li F, Xia H, Peng L, Boschiero C, Rosen BD, Bickhart DM, Zhang S, Guo A, Van Tassell CP, Smith TPL, Yang L, Liu GE. Assembly of a pangenome for global cattle reveals missing sequences and novel structural variations, providing new insights into their diversity and evolutionary history. Genome Res 2022; 32:gr.276550.122. [PMID: 35977842 PMCID: PMC9435747 DOI: 10.1101/gr.276550.122] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2022] [Accepted: 07/21/2022] [Indexed: 02/03/2023]
Abstract
A cattle pangenome representation was created based on the genome sequences of 898 cattle representing 57 breeds. The pangenome identified 83 Mb of sequence not found in the cattle reference genome, representing 3.1% novel sequence compared with the 2.71-Gb reference. A catalog of structural variants developed from this cattle population identified 3.3 million deletions, 0.12 million inversions, and 0.18 million duplications. Estimates of breed ancestry and hybridization between cattle breeds using insertion/deletions as markers were similar to those produced by single nucleotide polymorphism-based analysis. Hundreds of deletions were observed to have stratification based on subspecies and breed. For example, an insertion of a Bov-tA1 repeat element was identified in the first intron of the APPL2 gene and correlated with cattle breed geographic distribution. This insertion falls within a segment overlapping predicted enhancer and promoter regions of the gene, and could affect important traits such as immune response, olfactory functions, cell proliferation, and glucose metabolism in muscle. The results indicate that pangenomes are a valuable resource for studying diversity and evolutionary history, and help to delineate how domestication, trait-based breeding, and adaptive introgression have shaped the cattle genome.
Collapse
Affiliation(s)
- Yang Zhou
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan 430070, China
| | - Lv Yang
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan 430070, China
| | - Xiaotao Han
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan 430070, China
| | - Jiazheng Han
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan 430070, China
| | - Yan Hu
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan 430070, China
| | - Fan Li
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan 430070, China
| | - Han Xia
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan 430070, China
| | - Lingwei Peng
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan 430070, China
| | - Clarissa Boschiero
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, Maryland 20705, USA
| | - Benjamin D Rosen
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, Maryland 20705, USA
| | - Derek M Bickhart
- Dairy Forage Research Center, ARS USDA, Madison, Wisconsin 53706, USA
| | - Shujun Zhang
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan 430070, China
| | - Aizhen Guo
- The State Key Laboratory of Agricultural Microbiology, Huazhong Agricultural University, Wuhan 430070, China
| | - Curtis P Van Tassell
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, Maryland 20705, USA
| | - Timothy P L Smith
- U.S. Meat Animal Research Center, ARS USDA, Clay Center, Nebraska 68933, USA
| | - Liguo Yang
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan 430070, China
| | - George E Liu
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, Maryland 20705, USA
| |
Collapse
|
10
|
Kelly CJ, Chitko-McKown CG, Chuong EB. Ruminant-specific retrotransposons shape regulatory evolution of bovine immunity. Genome Res 2022; 32:gr.276241.121. [PMID: 35948370 PMCID: PMC9435751 DOI: 10.1101/gr.276241.121] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2021] [Accepted: 05/05/2022] [Indexed: 02/03/2023]
Abstract
Cattle are an important livestock species, and mapping the genomic architecture of agriculturally relevant traits such as disease susceptibility is a major challenge in the bovine research community. Lineage-specific transposable elements (TEs) are increasingly recognized to contribute to gene regulatory evolution and variation, but this possibility has been largely unexplored in ruminant genomes. We conducted epigenomic profiling of the type II interferon (IFN) response in bovine cells and found thousands of ruminant-specific TEs including MER41_BT and Bov-A2 elements predicted to act as IFN-inducible enhancer elements. CRISPR knockout experiments in bovine cells established that critical immune factors including IFNAR2 and IL2RB are transcriptionally regulated by TE-derived enhancers. Finally, population genomic analysis of 38 individuals revealed that a subset of polymorphic TE insertions may function as enhancers in modern cattle. Our study reveals that lineage-specific TEs have shaped the evolution of ruminant IFN responses and potentially continue to contribute to immune gene regulatory differences across modern breeds and individuals. Together with previous work in human cells, our findings demonstrate that lineage-specific TEs have been independently co-opted to regulate IFN-inducible gene expression in multiple species, supporting TE co-option as a recurrent mechanism driving the evolution of IFN-inducible transcriptional networks.
Collapse
Affiliation(s)
- Conor J Kelly
- Department of Molecular, Cellular, and Developmental Biology and BioFrontiers Institute, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Carol G Chitko-McKown
- USDA, ARS, Roman L. Hruska US Meat Animal Research Center (MARC), Clay Center, Nebraska 68933, USA
| | - Edward B Chuong
- Department of Molecular, Cellular, and Developmental Biology and BioFrontiers Institute, University of Colorado Boulder, Boulder, Colorado 80309, USA
| |
Collapse
|
11
|
Zhou J, Liu L, Reynolds E, Huang X, Garrick D, Shi Y. Discovering Copy Number Variation in Dual-Purpose XinJiang Brown Cattle. Front Genet 2022; 12:747431. [PMID: 35222511 PMCID: PMC8873982 DOI: 10.3389/fgene.2021.747431] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2021] [Accepted: 11/01/2021] [Indexed: 12/02/2022] Open
Abstract
Copy number variants (CNVs), which are a class of structural variant, can be important in relating genomic variation to phenotype. The primary aims of this study were to discover the common CNV regions (CNVRs) in the dual-purpose XinJiang-Brown cattle population and to detect differences between CNVs inferred using the ARS-UCD 1.2 (ARS) or the UMD 3.1 (UMD) genome assemblies based on the 150K SNP (Single Nucleotide Polymorphisms) Chip. PennCNV and CNVPartition methods were applied to calculate the deviation of the standardized signal intensity of SNPs markers to detect CNV status. Following the discovery of CNVs, we used the R package HandyCNV to generate and visualize CNVRs, compare CNVs and CNVRs between genome assemblies, and identify consensus genes using annotation resources. We identified 38 consensus CNVRs using the ARS assembly with 1.95% whole genome coverage, and 33 consensus CNVRs using the UMD assembly with 1.46% whole genome coverage using PennCNV and CNVPartition. We identified 37 genes that intersected 13 common CNVs (>5% frequency), these included functionally interesting genes such as GBP4 for which an increased copy number has been negatively associated with cattle stature, and the BoLA gene family which has been linked to the immune response and adaption of cattle. The ARS map file of the GGP Bovine 150K Bead Chip maps the genomic position of more SNPs with increased accuracy compared to the UMD map file. Comparison of the CNVRs identified between the two reference assemblies suggests the newly released ARS reference assembly is better for CNV detection. In spite of this, different CNV detection methods can complement each other to generate a larger number of CNVRs than using a single approach and can highlight more genes of interest.
Collapse
Affiliation(s)
- Jinghang Zhou
- School of Agriculture, Ningxia University, Yinchuan, China
- AL Rae Centre for Genetics and Breeding, Massey University, Hamilton, New Zealand
| | - Liyuan Liu
- School of Agriculture, Ningxia University, Yinchuan, China
- AL Rae Centre for Genetics and Breeding, Massey University, Hamilton, New Zealand
| | - Edwardo Reynolds
- AL Rae Centre for Genetics and Breeding, Massey University, Hamilton, New Zealand
| | - Xixia Huang
- College of Animal Science, Xinjiang Agricultural University, Urumqi, China
| | - Dorian Garrick
- AL Rae Centre for Genetics and Breeding, Massey University, Hamilton, New Zealand
- *Correspondence: Yuangang Shi, ; Dorian Garrick, mailto:
| | - Yuangang Shi
- School of Agriculture, Ningxia University, Yinchuan, China
- *Correspondence: Yuangang Shi, ; Dorian Garrick, mailto:
| |
Collapse
|
12
|
Prunier J, Carrier A, Gilbert I, Poisson W, Albert V, Taillon J, Bourret V, Côté SD, Droit A, Robert C. CNVs with adaptive potential in Rangifer tarandus: genome architecture and new annotated assembly. Life Sci Alliance 2021; 5:5/3/e202101207. [PMID: 34911809 PMCID: PMC8711850 DOI: 10.26508/lsa.202101207] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2021] [Revised: 11/29/2021] [Accepted: 11/29/2021] [Indexed: 01/13/2023] Open
Abstract
Rangifer tarandus has experienced recent drastic population size reductions throughout its circumpolar distribution and preserving the species implies genetic diversity conservation. To facilitate genomic studies of the species populations, we improved the genome assembly by combining long read and linked read and obtained a new highly accurate and contiguous genome assembly made of 13,994 scaffolds (L90 = 131 scaffolds). Using de novo transcriptome assembly of RNA-sequencing reads and similarity with annotated human gene sequences, 17,394 robust gene models were identified. As copy number variations (CNVs) likely play a role in adaptation, we additionally investigated these variations among 20 genomes representing three caribou ecotypes (migratory, boreal and mountain). A total of 1,698 large CNVs (length > 1 kb) showing a genome distribution including hotspots were identified. 43 large CNVs were particularly distinctive of the migratory and sedentary ecotypes and included genes annotated for functions likely related to the expected adaptations. This work includes the first publicly available annotation of the caribou genome and the first assembly allowing genome architecture analyses, including the likely adaptive CNVs reported here.
Collapse
Affiliation(s)
- Julien Prunier
- Département de Médecine Moléculaire, Faculté de Médecine, Université Laval, Quebec City, Canada
| | - Alexandra Carrier
- Département des sciences animales, Faculté des Sciences de l'Agriculture et de l'Alimentation, Université Laval, Quebec City, Canada
| | - Isabelle Gilbert
- Département des sciences animales, Faculté des Sciences de l'Agriculture et de l'Alimentation, Université Laval, Quebec City, Canada
| | - William Poisson
- Département des sciences animales, Faculté des Sciences de l'Agriculture et de l'Alimentation, Université Laval, Quebec City, Canada
| | - Vicky Albert
- Ministère des Forêts, de la Faune et des Parcs du Québec, Quebec City, Canada
| | - Joëlle Taillon
- Ministère des Forêts, de la Faune et des Parcs du Québec, Quebec City, Canada
| | - Vincent Bourret
- Ministère des Forêts, de la Faune et des Parcs du Québec, Quebec City, Canada
| | - Steeve D Côté
- Caribou Ungava, département de biologie, Faculté des Sciences et de Génie, Université Laval, Quebec City, Canada
| | - Arnaud Droit
- Département de Médecine Moléculaire, Faculté de Médecine, Université Laval, Quebec City, Canada
| | - Claude Robert
- Département des sciences animales, Faculté des Sciences de l'Agriculture et de l'Alimentation, Université Laval, Quebec City, Canada
| |
Collapse
|
13
|
Zhang X, Liu X, Zhang Z, Zhao Y, Wang C. Atomic-Level Investigation of Reactant Recognition Mechanism and Thermodynamic Property in Glucosamine 6-Phosphate Deaminase Catalysis. Front Chem 2021; 9:737492. [PMID: 34414167 PMCID: PMC8369159 DOI: 10.3389/fchem.2021.737492] [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: 07/07/2021] [Accepted: 07/20/2021] [Indexed: 11/13/2022] Open
Abstract
Glucosamine 6-phosphate deaminase (NagB) influences the direction of N-acetylglucosamine (GlcNAc) metabolism, facilitating the conversion of D-glucosamine 6-phosphate (GlcN6P) to D-fructose 6-phosphate (Fru6P) with the release of ammonia. Here, extensive molecular dynamics simulations combined with various techniques were performed to study the recognition and delivery process of GlcN6P by SmuNagB, due to its guidance of subsequent enzymatic reaction. The key residues Lys194, His130, Arg127, Thr38, and Ser37 stabilize GlcN6P in the active site by hydrogen bond interactions, therein electrostatic and polar solvent effects provide the primary traction. Four delivery channels were identified, with GlcN6P most likely to enter the active site of NagB through a “door” comprising residues 6–10, 122–136, and 222–233. The corresponding mechanism and thermodynamic properties were investigated. An exothermic recognition and delivery process were detected, accompanied by the flipping of GlcN6P and changes in key direct and indirect hydrogen bond interactions, which provide the driving force for the chemical reaction to occur. Furthermore, “the lid motif” was identified that remain open in alkaline condition with different extent of opening at each stage of transfer that induced GlcN6P to move the active site of NagB. The work will assist in the elucidation of the catalytic mechanism of action of NagB, allowing inhibitors to be designed with superior dynamic behavior.
Collapse
Affiliation(s)
- Xiao Zhang
- The Key Laboratory of Natural Medicine and Immuno-Engineering, Henan University, Kaifeng, China
| | - Xiaoyuan Liu
- The Key Laboratory of Natural Medicine and Immuno-Engineering, Henan University, Kaifeng, China
| | - Zhiyang Zhang
- The Key Laboratory of Natural Medicine and Immuno-Engineering, Henan University, Kaifeng, China
| | - Yuan Zhao
- The Key Laboratory of Natural Medicine and Immuno-Engineering, Henan University, Kaifeng, China
| | - Chaojie Wang
- The Key Laboratory of Natural Medicine and Immuno-Engineering, Henan University, Kaifeng, China
| |
Collapse
|
14
|
Jang J, Kim K, Lee YH, Kim H. Population differentiated copy number variation of Bos taurus, Bos indicus and their African hybrids. BMC Genomics 2021; 22:531. [PMID: 34253178 PMCID: PMC8276479 DOI: 10.1186/s12864-021-07808-7] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2021] [Accepted: 06/10/2021] [Indexed: 01/10/2023] Open
Abstract
Background CNV comprises a large proportion in cattle genome and is associated with various traits. However, there were few population-scale comparison studies on cattle CNV. Results Here, autosome-wide CNVs were called by read depth of NGS alignment result and copy number variation regions (CNVRs) defined from 102 Eurasian taurine (EAT) of 14 breeds, 28 Asian indicine (ASI) of 6 breeds, 22 African taurine (AFT) of 2 breeds, and 184 African humped cattle (AFH) of 17 breeds. The copy number of every CNVRs were compared between populations and CNVRs with population differentiated copy numbers were sorted out using the pairwise statistics VST and Kruskal-Wallis test. Three hundred sixty-two of CNVRs were significantly differentiated in both statistics and 313 genes were located on the population differentiated CNVRs. Conclusion For some of these genes, the averages of copy numbers were also different between populations and these may be candidate genes under selection. These include olfactory receptors, pathogen-resistance, parasite-resistance, heat tolerance and productivity related genes. Furthermore, breed- and individual-level comparison was performed using the presence or copy number of the autosomal CNVRs. Our findings were based on identification of CNVs from short Illumina reads of 336 individuals and 39 breeds, which to our knowledge is the largest dataset for this type of analysis and revealed important CNVs that may play a role in cattle adaption to various environments. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07808-7.
Collapse
Affiliation(s)
- Jisung Jang
- Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul, Republic of Korea
| | - Kwondo Kim
- Department of Agricultural Biotechnology and Research Institute of Agriculture and Life Sciences, Seoul National University, Seoul, Republic of Korea
| | - Young Ho Lee
- Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul, Republic of Korea
| | - Heebal Kim
- Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul, Republic of Korea. .,Department of Agricultural Biotechnology and Research Institute of Agriculture and Life Sciences, Seoul National University, Seoul, Republic of Korea. .,eGnome, Inc, Seoul, South Korea.
| |
Collapse
|
15
|
Gutierrez-Reinoso MA, Aponte PM, Garcia-Herreros M. Genomic Analysis, Progress and Future Perspectives in Dairy Cattle Selection: A Review. Animals (Basel) 2021; 11:599. [PMID: 33668747 PMCID: PMC7996307 DOI: 10.3390/ani11030599] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2020] [Revised: 11/23/2020] [Accepted: 11/24/2020] [Indexed: 12/16/2022] Open
Abstract
Genomics comprises a set of current and valuable technologies implemented as selection tools in dairy cattle commercial breeding programs. The intensive progeny testing for production and reproductive traits based on genomic breeding values (GEBVs) has been crucial to increasing dairy cattle productivity. The knowledge of key genes and haplotypes, including their regulation mechanisms, as markers for productivity traits, may improve the strategies on the present and future for dairy cattle selection. Genome-wide association studies (GWAS) such as quantitative trait loci (QTL), single nucleotide polymorphisms (SNPs), or single-step genomic best linear unbiased prediction (ssGBLUP) methods have already been included in global dairy programs for the estimation of marker-assisted selection-derived effects. The increase in genetic progress based on genomic predicting accuracy has also contributed to the understanding of genetic effects in dairy cattle offspring. However, the crossing within inbred-lines critically increased homozygosis with accumulated negative effects of inbreeding like a decline in reproductive performance. Thus, inaccurate-biased estimations based on empirical-conventional models of dairy production systems face an increased risk of providing suboptimal results derived from errors in the selection of candidates of high genetic merit-based just on low-heritability phenotypic traits. This extends the generation intervals and increases costs due to the significant reduction of genetic gains. The remarkable progress of genomic prediction increases the accurate selection of superior candidates. The scope of the present review is to summarize and discuss the advances and challenges of genomic tools for dairy cattle selection for optimizing breeding programs and controlling negative inbreeding depression effects on productivity and consequently, achieving economic-effective advances in food production efficiency. Particular attention is given to the potential genomic selection-derived results to facilitate precision management on modern dairy farms, including an overview of novel genome editing methodologies as perspectives toward the future.
Collapse
Affiliation(s)
- Miguel A. Gutierrez-Reinoso
- Facultad de Ciencias Agropecuarias y Recursos Naturales, Carrera de Medicina Veterinaria, Universidad Técnica de Cotopaxi (UTC), Latacunga 05-0150, Ecuador
- Laboratorio de Biotecnología Animal, Departamento de Ciencia Animal, Facultad de Ciencias Veterinarias, Universidad de Concepción (UdeC), Chillán 3780000, Chile
| | - Pedro M. Aponte
- Colegio de Ciencias Biológicas y Ambientales (COCIBA), Universidad San Francisco de Quito (USFQ), Quito 170157, Ecuador
- Campus Cumbayá, Instituto de Investigaciones en Biomedicina “One-health”, Universidad San Francisco de Quito (USFQ), Quito 170157, Ecuador
| | - Manuel Garcia-Herreros
- Instituto Nacional de Investigação Agrária e Veterinária (INIAV), 2005-048 Santarém, Portugal
| |
Collapse
|