1
|
Yuan H, Wei W, Zhang Y, Li C, Zhao S, Chao Z, Xia C, Quan J, Gao C. Unveiling the Influence of Copy Number Variations on Genetic Diversity and Adaptive Evolution in China's Native Pig Breeds via Whole-Genome Resequencing. Int J Mol Sci 2024; 25:5843. [PMID: 38892031 PMCID: PMC11172908 DOI: 10.3390/ijms25115843] [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/10/2024] [Revised: 05/22/2024] [Accepted: 05/25/2024] [Indexed: 06/21/2024] Open
Abstract
Copy number variations (CNVs) critically influence individual genetic diversity and phenotypic traits. In this study, we employed whole-genome resequencing technology to conduct an in-depth analysis of 50 pigs from five local swine populations [Rongchang pig (RC), Wuzhishan pig (WZS), Tibetan pig (T), Yorkshire (YL) and Landrace (LR)], aiming to assess their genetic potential and explore their prospects in the field of animal model applications. We identified a total of 96,466 CNVs, which were subsequently integrated into 7112 non-redundant CNVRs, encompassing 1.3% of the swine genome. Functional enrichment analysis of the genes within these CNVRs revealed significant associations with sensory perception, energy metabolism, and neural-related pathways. Further selective scan analyses of the local pig breeds RC, T, WZS, along with YL and LR, uncovered that for the RC variety, the genes PLA2G10 and ABCA8 were found to be closely related to fat metabolism and cardiovascular health. In the T breed, the genes NCF2 and CSGALNACT1 were associated with immune response and connective tissue characteristics. As for the WZS breed, the genes PLIN4 and CPB2 were primarily linked to fat storage and anti-inflammatory responses. In summary, this research underscores the pivotal role of CNVs in fostering the diversity and adaptive evolution of pig breeds while also offering valuable insights for further exploration of the advantageous genetic traits inherent to China's local pig breeds. This facilitates the creation of experimental animal models tailored to the specific characteristics of these breeds, contributing to the advancement of livestock and biomedical research.
Collapse
Affiliation(s)
- Haonan Yuan
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou 730030, China; (H.Y.); (W.W.); (Y.Z.); (S.Z.)
| | - Wenjing Wei
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou 730030, China; (H.Y.); (W.W.); (Y.Z.); (S.Z.)
| | - Yue Zhang
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou 730030, China; (H.Y.); (W.W.); (Y.Z.); (S.Z.)
| | - Changwen Li
- State Key Laboratory for Animal Disease Control and Prevention, Harbin Veterinary Research Institute, Chinese Academy of Agricultural Sciences (CAAS), Heilongjiang Provincial Key Laboratory of Laboratory Animal and Comparative Medicine, National Poultry Laboratory Animal Resource Center, Harbin 150069, China; (C.L.); (C.X.)
| | - Shengguo Zhao
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou 730030, China; (H.Y.); (W.W.); (Y.Z.); (S.Z.)
| | - Zhe Chao
- Institute of Animal Science and Veterinary Medicine, Hainan Academy of Agricultural Sciences, Key Laboratory of Tropical Animal Breeding and Disease Research, Haikou 571100, China;
| | - Changyou Xia
- State Key Laboratory for Animal Disease Control and Prevention, Harbin Veterinary Research Institute, Chinese Academy of Agricultural Sciences (CAAS), Heilongjiang Provincial Key Laboratory of Laboratory Animal and Comparative Medicine, National Poultry Laboratory Animal Resource Center, Harbin 150069, China; (C.L.); (C.X.)
| | - Jinqiang Quan
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou 730030, China; (H.Y.); (W.W.); (Y.Z.); (S.Z.)
| | - Caixia Gao
- State Key Laboratory for Animal Disease Control and Prevention, Harbin Veterinary Research Institute, Chinese Academy of Agricultural Sciences (CAAS), Heilongjiang Provincial Key Laboratory of Laboratory Animal and Comparative Medicine, National Poultry Laboratory Animal Resource Center, Harbin 150069, China; (C.L.); (C.X.)
| |
Collapse
|
2
|
Jin M, Liu G, Lu J, Chen Z, Wang H, Li T, Wei C. Characteristics of genetic basis copy number variation in production and adaptation traits of Chinese indigenous sheep. J Genet Genomics 2024; 51:566-569. [PMID: 38280418 DOI: 10.1016/j.jgg.2024.01.005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2023] [Revised: 01/17/2024] [Accepted: 01/19/2024] [Indexed: 01/29/2024]
Affiliation(s)
- Meilin Jin
- State Key Laboratory of Animal Biotech Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China; College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Gang Liu
- National Animal Husbandry Service, National Center of Preservation & Utilization of Animal Genetic Resources, Beijing 100125, China
| | - Jian Lu
- National Animal Husbandry Service, National Center of Preservation & Utilization of Animal Genetic Resources, Beijing 100125, China
| | | | - Huihua Wang
- State Key Laboratory of Animal Biotech Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China
| | - Taotao Li
- State Key Laboratory of Animal Biotech Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China
| | - Caihong Wei
- State Key Laboratory of Animal Biotech Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China.
| |
Collapse
|
3
|
Benfica LF, Brito LF, do Bem RD, de Oliveira LF, Mulim HA, Braga LG, Cyrillo JNSG, Bonilha SFM, Mercadante MEZ. Detection and characterization of copy number variation in three differentially-selected Nellore cattle populations. Front Genet 2024; 15:1377130. [PMID: 38694873 PMCID: PMC11061390 DOI: 10.3389/fgene.2024.1377130] [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/26/2024] [Accepted: 04/05/2024] [Indexed: 05/04/2024] Open
Abstract
Introduction: Nellore cattle (Bos taurus indicus) is the main beef cattle breed raised in Brazil. This breed is well adapted to tropical conditions and, more recently, has experienced intensive genetic selection for multiple performance traits. Over the past 43 years, an experimental breeding program has been developed in the Institute of Animal Science (IZ, Sertaozinho, SP, Brazil), which resulted in three differentially-selected lines known as Nellore Control (NeC), Nellore Selection (NeS), and Nellore Traditional (NeT). The primary goal of this selection experiment was to determine the response to selection for yearling weight (YW) and residual feed intake (RFI) on Nellore cattle. The main objectives of this study were to: 1) identify copy number variation (CNVs) in Nellore cattle from three selection lines; 2) identify and characterize CNV regions (CNVR) on these three lines; and 3) perform functional enrichment analyses of the CNVR identified. Results: A total of 14,914 unique CNVs and 1,884 CNVRs were identified when considering all lines as a single population. The CNVRs were non-uniformly distributed across the chromosomes of the three selection lines included in the study. The NeT line had the highest number of CNVRs (n = 1,493), followed by the NeS (n = 823) and NeC (n = 482) lines. The CNVRs covered 23,449,890 bp (0.94%), 40,175,556 bp (1.61%), and 63,212,273 bp (2.54%) of the genome of the NeC, NeS, and NeT lines, respectively. Two CNVRs were commonly identified between the three lines, and six, two, and four exclusive regions were identified for NeC, NeS, and NeT, respectively. All the exclusive regions overlap with important genes, such as SMARCD3, SLC15A1, and MAPK1. Key biological processes associated with the candidate genes were identified, including pathways related to growth and metabolism. Conclusion: This study revealed large variability in CNVs and CNVRs across three Nellore lines differentially selected for YW and RFI. Gene annotation and gene ontology analyses of the exclusive CNVRs to each line revealed specific genes and biological processes involved in the expression of growth and feed efficiency traits. These findings contribute to the understanding of the genetic mechanisms underlying the phenotypic differences among the three Nellore selection lines.
Collapse
Affiliation(s)
- Lorena F. Benfica
- Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
- Department of Animal Science, Faculty of Agricultural and Veterinary Sciences, Sao Paulo State University, Jaboticabal, São Paulo, Brazil
| | - Luiz F. Brito
- Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
| | - Ricardo D. do Bem
- Department of Animal Science, Faculty of Agricultural and Veterinary Sciences, Sao Paulo State University, Jaboticabal, São Paulo, Brazil
| | | | - Henrique A. Mulim
- Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
| | - Larissa G. Braga
- Department of Animal Science, Faculty of Agricultural and Veterinary Sciences, Sao Paulo State University, Jaboticabal, São Paulo, Brazil
- Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada
| | | | - Sarah F. M. Bonilha
- Beef Cattle Research Center, Institute of Animal Science, Sertaozinho, São Paulo, Brazil
| | - Maria Eugenia Z. Mercadante
- Department of Animal Science, Faculty of Agricultural and Veterinary Sciences, Sao Paulo State University, Jaboticabal, São Paulo, Brazil
- Beef Cattle Research Center, Institute of Animal Science, Sertaozinho, São Paulo, Brazil
| |
Collapse
|
4
|
Yang Y, Tang J, Yang H, Yang S, Cai M, Qi A, Lan X, Huang B, Su C, Chen H. Copy number variation of bovine S100A7 as a positional candidate affected body measurements. Anim Biotechnol 2023; 34:2141-2149. [PMID: 35815693 DOI: 10.1080/10495398.2022.2077740] [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: 11/01/2022]
Abstract
Beef production is closely related to the national economy and the attention has been paid to the improvement of beef cattle by molecular markers associated. Copy number variations (CNVs) recently have been gained many researches and recognized as an important source of genetic variation. Extensive studies have indicated that CNVs have effects on a large range of economic traits by a wide range of gene copy number alteration. S100A7 is a member of S100 family which is a famous family of Ca2+-binding proteins. S100A7 plays a crucial role in many important phenotypes (progress) including inflammatory diseases, psoriasis, obesity, etc. The aim of our study was to explore the phenotypic effects of CNV located in the S100A7 gene of bovine chromosome 3. We detected S100A7 CNV by qPCR in different cattle breeds, including Qinchuan cattle, Yunling cattle, Xianan cattle and a crossbred group Pinan. The copy number was identified as gain, normal and loss type, our results showed that the gain type was the main type in three types of S100A7 CNV of the whole tested breeds. After CNV detection, association analysis between S100A7 CNV and growth traits was carried out in four cattle breeds. We found significant effects of the CNV on cattle growth traits with differently preferred CNV types such as gain type with better chest depth (p = 0.043) in QC, loss type with better body length (p = 0.008) and rump width (p = 0.014) in YL, normal with better chest girth (p = 0.001), gain with better waist width (p = 0.001) and rump width (p = 0.044) in PN. These results suggested that the S100A7 CNV could affect the phenotypic traits and be used as a promising genetic marker for cattle molecular breeding.
Collapse
Affiliation(s)
- Yu Yang
- Shaanxi Key Laboratory of Molecular Biology for Agriculture, College of Animal Science and Technology, Northwest A & F University, Yangling, China
| | - Jia Tang
- Shaanxi Key Laboratory of Molecular Biology for Agriculture, College of Animal Science and Technology, Northwest A & F University, Yangling, China
| | - Haiyan Yang
- Shaanxi Key Laboratory of Molecular Biology for Agriculture, College of Animal Science and Technology, Northwest A & F University, Yangling, China
| | - Shuling Yang
- Shaanxi Key Laboratory of Molecular Biology for Agriculture, College of Animal Science and Technology, Northwest A & F University, Yangling, China
| | - Ming Cai
- Yunnan Academy of Grassland and Animal Science, Kunming, China
| | - Ao Qi
- Shaanxi Key Laboratory of Molecular Biology for Agriculture, College of Animal Science and Technology, Northwest A & F University, Yangling, China
| | - Xianyong Lan
- Shaanxi Key Laboratory of Molecular Biology for Agriculture, College of Animal Science and Technology, Northwest A & F University, Yangling, China
| | - Bizhi Huang
- Yunnan Academy of Grassland and Animal Science, Kunming, China
| | - Chao Su
- Shaanxi Key Laboratory of Molecular Biology for Agriculture, College of Animal Science and Technology, Northwest A & F University, Yangling, China
| | - Hong Chen
- Shaanxi Key Laboratory of Molecular Biology for Agriculture, College of Animal Science and Technology, Northwest A & F University, Yangling, China
- College of Animal Science, Xinjiang Agricultural University, Urumqi, China
| |
Collapse
|
5
|
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
|
6
|
Kooverjee BB, Soma P, van der Nest MA, Scholtz MM, Neser FWC. Copy Number Variation Discovery in South African Nguni-Sired and Bonsmara-Sired Crossbred Cattle. Animals (Basel) 2023; 13:2513. [PMID: 37570321 PMCID: PMC10417447 DOI: 10.3390/ani13152513] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2023] [Revised: 07/26/2023] [Accepted: 08/01/2023] [Indexed: 08/13/2023] Open
Abstract
Crossbreeding forms part of Climate-Smart beef production and is one of the strategies to mitigate the effects of climate change. Two Nguni-sired and three Bonsmara-sired crossbred animals underwent whole genome sequencing. Following quality control and file preparation, the sequence data were investigated for genome-wide copy number variation (CNV) using the panelcn.MOPS tool. A total of 355 CNVs were identified in the crossbreds, of which 274 were unique in Bonsmara-sired crossbreds and 81 unique in the Nguni-sired crossbreds. Genes that differed in copy number in both crossbreds included genes related to growth (SCRN2, LOC109572916) and fertility-related factors (RPS28, LOC1098562432, LOC109570037). Genes that were present only in the Bonsmara-sired crossbreds included genes relating to lipid metabolism (MAF1), olfaction (LOC109569114), body size (HES7), immunity (LOC10957335, LOC109877039) and disease (DMBT1). Genes that were present only in the Nguni-sired crossbreds included genes relating to ketosis (HMBOX1) and amino acid transport (LOC109572916). Results of this study indicate that Nguni and Bonsmara cattle can be utilized in crossbreeding programs as they may enhance the presence of economically important traits associated with both breeds. This will produce crossbred animals that are good meat producers, grow faster, have high fertility, strong immunity and a better chance of producing in South Africa's harsh climate conditions. Ultimately, this study provides new genetic insights into the adaptability of Nguni and Bonsmara crossbred cattle.
Collapse
Affiliation(s)
| | - Pranisha Soma
- Animal Production, Agricultural Research Council, Pretoria 0062, South Africa;
| | - Magrieta A. van der Nest
- Department of Biochemistry, Genetics and Microbiology, University of Pretoria, Pretoria 0028, South Africa;
| | - Michiel M. Scholtz
- Animal Production, Agricultural Research Council, Pretoria 0062, South Africa;
- Department of Animal Science, University of the Free State, Bloemfontein 9300, South Africa;
| | - Frederick W. C. Neser
- Department of Animal Science, University of the Free State, Bloemfontein 9300, South Africa;
| |
Collapse
|
7
|
Peripolli E, Stafuzza NB, Machado MA, do Carmo Panetto JC, do Egito AA, Baldi F, da Silva MVGB. Assessment of copy number variants in three Brazilian locally adapted cattle breeds using whole-genome re-sequencing data. Anim Genet 2023; 54:254-270. [PMID: 36740987 DOI: 10.1111/age.13298] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Revised: 12/13/2021] [Accepted: 01/13/2023] [Indexed: 02/07/2023]
Abstract
Further characterization of genetic structural variations should strongly focus on small and endangered local breeds given their role in unraveling genes and structural variants underlying selective pressures and phenotype variation. A comprehensive genome-wide assessment of copy number variations (CNVs) based on whole-genome re-sequencing data was performed on three Brazilian locally adapted cattle breeds (Caracu Caldeano, Crioulo Lageano, and Pantaneiro) using the ARS-UCD1.2 genome assembly. Data from 36 individuals with an average coverage depth of 14.07× per individual was used. A total of 24 945 CNVs were identified distributed among the breeds (Caracu Caldeano = 7285, Crioulo Lageano = 7297, and Pantaneiro = 10 363). Deletion events were 1.75-2.07-fold higher than duplications, and the total length of CNVs is composed mostly of a high number of segments between 10 and 30 kb. CNV regions (CNVRs) are not uniformly scattered throughout the genomes (n = 463), and 105 CNVRs were found overlapping among the studied breeds. Functional annotation of the CNVRs revealed variants with high consequence on protein sequence harboring relevant genes, in which we highlighted the BOLA-DQB, BOLA-DQA5, CD1A, β-defensins, PRG3, and ULBP21 genes. Enrichment analysis based on the gene list retrieved from the CNVRs disclosed over-represented terms (p < 0.01) strongly associated with immunity and cattle resilience to harsh environments. Additionally, QTL associated with body conformation and dairy-related traits were also unveiled within the CNVRs. These results provide better understanding of the selective forces shaping the genome of such cattle breeds and identify traces of natural selection pressures by which these populations have been exposed to challenging environmental conditions.
Collapse
Affiliation(s)
- Elisa Peripolli
- School of Agricultural and Veterinarian Sciences, São Paulo State University (Unesp), Jaboticabal, Brazil
| | | | | | | | | | - Fernando Baldi
- School of Agricultural and Veterinarian Sciences, São Paulo State University (Unesp), Jaboticabal, Brazil
| | | |
Collapse
|
8
|
Identification of Copy Number Variations in Four Horse Breed Populations in South Korea. Animals (Basel) 2022; 12:ani12243501. [PMID: 36552421 PMCID: PMC9774267 DOI: 10.3390/ani12243501] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2022] [Revised: 11/21/2022] [Accepted: 12/08/2022] [Indexed: 12/14/2022] Open
Abstract
In this study, genome-wide CNVs were identified using a total of 469 horses from four horse populations (Jeju horses, Thoroughbreds, Jeju riding horses, and Hanla horses). We detected a total of 843 CNVRs throughout all autosomes: 281, 30, 301, and 310 CNVRs for Jeju horses, Thoroughbreds, Jeju riding horses, and Hanla horses, respectively. Of the total CNVRs, copy number losses were found to be the most abundant (48.99%), while gains and mixed CNVRs accounted for 41.04% and 9.96% of the total CNVRs, respectively. The length of the CNVRs ranged from 0.39 kb to 2.8 Mb, while approximately 7.2% of the reference horse genome assembly was covered by the total CNVRs. By comparing the CNVRs among the populations, we found a significant portion of the CNVRs (30.13%) overlapped; the highest number of shared CNVRs was between Hanla horses and Jeju riding horses. When compared with the horse CNVRs of previous studies, 26.8% of CNVRs were found to be uniquely detected in this study. The CNVRs were not randomly distributed throughout the genome; in particular, the Equus caballus autosome (ECA) 7 comprised the largest proportion of its genome (16.3%), while ECA 24 comprised the smallest (0.7%). Furthermore, functional analysis was applied to CNVRs that overlapped with genes (genic-CNVRs); these overlapping areas may be potentially associated with the olfactory pathway and nervous system. A racing performance QTL was detected in a CNVR of Thoroughbreds, Jeju riding horses, and Hanla horses, and the CNVR value was mixed for three breeds.
Collapse
|
9
|
Effects of Genetic Variation of the Sorting Nexin 29 ( SNX29) Gene on Growth Traits of Xiangdong Black Goat. Animals (Basel) 2022; 12:ani12243461. [PMID: 36552381 PMCID: PMC9774745 DOI: 10.3390/ani12243461] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2022] [Revised: 11/20/2022] [Accepted: 11/30/2022] [Indexed: 12/13/2022] Open
Abstract
Previous studies have found that the copy number variation (CNV) and insertion/deletion (indels) located in the sorting nexin 29 (SNX29) gene, which is an important candidate gene related to meat production and quality, are associated with growth traits of African goats and Shaanbei white cashmere goats. However, the genetic effects of SNX29 genetic variation on growth traits of Xiangdong black (XDB) goat (a representative meat goat breed in China) are still unclear. The purpose of this study was to detect the mRNA expression level of SNX29 and to explore the genetic effects of CNV and indel within SNX29 on growth traits and gene expression in XDB goat. The SNX29 mRNA expression profile showed that the SNX29 was highly expressed in adipose tissues, indicating that the SNX29 gene could play a key role in subcutaneous adipose deposition of XDB goat. 17 bp indel (g.10559298-10559314), 21 bp indel (g.10918982-10919002) and CNV were detected in 516 individuals of XDB goat by PCR or qPCR. The association analysis of SNX29 CNV with growth traits in XDB goats showed that SNX29 CNV was significantly correlated with chest circumference and abdominal circumference (p < 0.01), and the normal type of SNX29 CNV goat individuals were more advantageous. For the mRNA expression of SNX29 gene, individuals with SNX29 copy number normal type had a higher trend than that of SNX29 gene with copy number gain type in longissimus dorsi muscle (p = 0.07), whereas individuals with SNX29 copy number gain type had a higher trend in abdominal adipose (p = 0.09). Overall, these results suggested that the SNX29 gene could play an important role in growth and development of XDB goats and could be used for marker-assisted selection (MAS) in XDB goats.
Collapse
|
10
|
Liu M, Huang C, Dai R, Ren W, Li X, Wu X, Ma X, Chu M, Bao P, Guo X, Pei J, Xiong L, Yan P, Liang C. Copy Number Variations in the MICALL2 and MOGAT2 Genes Are Associated with Ashidan Yak Growth Traits. Animals (Basel) 2022; 12:ani12202779. [PMID: 36290165 PMCID: PMC9597734 DOI: 10.3390/ani12202779] [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: 09/05/2022] [Revised: 10/12/2022] [Accepted: 10/13/2022] [Indexed: 11/29/2022] Open
Abstract
Copy number variations (CNVs) are a result of genomic rearrangement affecting DNA regions over 1 kb in length, and can include inversions, translocations, deletions, and duplications. The molecule interacting with CasL-like protein 2 (MICALL2) gene is primarily associated with mitochondrial protein targeting and exhibits predicted stress fiber colocalization. The monoacylglycerol O-acyltransferase 2 (MOGAT2) gene encodes an enzyme responsible for catalyzing diacylglycerol synthesis from 2-monoacylglycerol and fatty acyl-CoA. For this study, blood samples were obtained from 315 yaks, and the body weight, body length, withers height, and chest girth of these animals were measured at 6, 12, 18, and 30 months of age. Genomic DNA was harvested from the collected blood samples, and CNVs in these samples were detected by qPCR. The resultant data were compared using ANOVAs, revealing significant associations between MICALL2 gene CNVs and body weight at 6 months of age (p < 0.05), body length and chest girth at 30 months of age (p < 0.05), and withers height at 18 months of age (p < 0.01) in Ashidan yaks. Similarly, MOGAT2 CNVs were significantly associated with body weight at 6 and 30 months of age (p < 0.05), and with withers height at 18 months of age (p < 0.01) in these Ashidan yaks. MICALL2 and MOGAT2 gene expression was further analyzed in yak tissue samples, revealing that MICALL2 was most highly expressed in the adipose tissue, whereas MOGAT2 was most highly expressed in the lung. These results thus confirmed the relationship between CNVs in the MICALL2 and MOGAT2 genes and Ashidan yak growth traits, providing a valuable gene locus that can be leveraged for future marker-assisted yak breeding efforts.
Collapse
Affiliation(s)
- Modian Liu
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Science, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Chun Huang
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Science, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Rongfeng Dai
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Science, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Wenwen Ren
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Science, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Xinyi Li
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Science, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Xiaoyun Wu
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Science, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Xiaoming Ma
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Science, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Min Chu
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Science, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Pengjia Bao
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Science, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Xian Guo
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Science, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Jie Pei
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Science, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Lin Xiong
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Science, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Ping Yan
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Science, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
- Correspondence: (P.Y.); (C.L.)
| | - Chunnian Liang
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Science, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
- Correspondence: (P.Y.); (C.L.)
| |
Collapse
|
11
|
Davoudi P, Do DN, Rathgeber B, Colombo SM, Sargolzaei M, Plastow G, Wang Z, Karimi K, Hu G, Valipour S, Miar Y. Genome-wide detection of copy number variation in American mink using whole-genome sequencing. BMC Genomics 2022; 23:649. [PMID: 36096727 PMCID: PMC9468235 DOI: 10.1186/s12864-022-08874-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2022] [Accepted: 09/05/2022] [Indexed: 02/06/2023] Open
Abstract
BACKGROUND Copy number variations (CNVs) represent a major source of genetic diversity and contribute to the phenotypic variation of economically important traits in livestock species. In this study, we report the first genome-wide CNV analysis of American mink using whole-genome sequence data from 100 individuals. The analyses were performed by three complementary software programs including CNVpytor, DELLY and Manta. RESULTS A total of 164,733 CNVs (144,517 deletions and 20,216 duplications) were identified representing 5378 CNV regions (CNVR) after merging overlapping CNVs, covering 47.3 Mb (1.9%) of the mink autosomal genome. Gene Ontology and KEGG pathway enrichment analyses of 1391 genes that overlapped CNVR revealed potential role of CNVs in a wide range of biological, molecular and cellular functions, e.g., pathways related to growth (regulation of actin cytoskeleton, and cAMP signaling pathways), behavior (axon guidance, circadian entrainment, and glutamatergic synapse), lipid metabolism (phospholipid binding, sphingolipid metabolism and regulation of lipolysis in adipocytes), and immune response (Wnt signaling, Fc receptor signaling, and GTPase regulator activity pathways). Furthermore, several CNVR-harbored genes associated with fur characteristics and development (MYO5A, RAB27B, FGF12, SLC7A11, EXOC2), and immune system processes (SWAP70, FYN, ORAI1, TRPM2, and FOXO3). CONCLUSIONS This study presents the first genome-wide CNV map of American mink. We identified 5378 CNVR in the mink genome and investigated genes that overlapped with CNVR. The results suggest potential links with mink behaviour as well as their possible impact on fur quality and immune response. Overall, the results provide new resources for mink genome analysis, serving as a guideline for future investigations in which genomic structural variations are present.
Collapse
Affiliation(s)
- Pourya Davoudi
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| | - Duy Ngoc Do
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| | - Bruce Rathgeber
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| | - Stefanie M Colombo
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| | - Mehdi Sargolzaei
- Department of Pathobiology, University of Guelph, Guelph, ON, Canada
- Select Sires Inc., Plain City, OH, USA
| | - Graham Plastow
- Livestock Gentec, Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Zhiquan Wang
- Livestock Gentec, Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Karim Karimi
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| | - Guoyu Hu
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| | - Shafagh Valipour
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| | - Younes Miar
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada.
| |
Collapse
|
12
|
Hu L, Zhang L, Li Q, Liu H, Xu T, Zhao N, Han X, Xu S, Zhao X, Zhang C. Genome-wide analysis of CNVs in three populations of Tibetan sheep using whole-genome resequencing. Front Genet 2022; 13:971464. [PMID: 36160022 PMCID: PMC9490000 DOI: 10.3389/fgene.2022.971464] [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: 06/17/2022] [Accepted: 08/23/2022] [Indexed: 01/29/2023] Open
Abstract
Copy number variation (CNV), an important source of genomic structural variation, can disturb genetic structure, dosage, regulation and expression, and is associated with phenotypic diversity and adaptation to local environments in mammals. In the present study, 24 resequencing datasets were used to characterize CNVs in three ecotypic populations of Tibetan sheep and assess CNVs related to domestication and adaptation in Qinghai-Tibetan Plateau. A total of 87,832 CNV events accounting for 0.3% of the sheep genome were detected. After merging the overlapping CNVs, 2777 CNV regions (CNVRs) were obtained, among which 1098 CNVRs were shared by the three populations. The average length of these CNVRs was more than 3 kb, and duplication events were more frequent than deletions. Functional analysis showed that the shared CNVRs were significantly enriched in 56 GO terms and 18 KEGG pathways that were mainly concerned with ABC transporters, olfactory transduction and oxygen transport. Moreover, 188 CNVRs overlapped with 97 quantitative trait loci (QTLs), such as growth and carcass QTLs, immunoglobulin QTLs, milk yield QTLs and fecal egg counts QTLs. PCDH15, APP and GRID2 overlapped with body weight QTLs. Furthermore, Vst analysis showed that RUNX1, LOC101104348, LOC105604082 and PAG11 were highly divergent between Highland-type Tibetan Sheep (HTS) and Valley-type Tibetan sheep (VTS), and RUNX1 and LOC101111988 were significantly differentiated between VTS and Oura-type Tibetan sheep (OTS). The duplication of RUNX1 may facilitate the hypoxia adaptation of OTS and HTS in Qinghai-Tibetan Plateau, which deserves further research in detail. In conclusion, for the first time, we represented the genome-wide distribution characteristics of CNVs in Tibetan sheep by resequencing, and provided a valuable genetic variation resource, which will facilitate the elucidation of the genetic basis underlying the distinct phenotypic traits and local adaptation of Tibetan sheep.
Collapse
Affiliation(s)
- Linyong Hu
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Liangzhi Zhang
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Qi Li
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Hongjin Liu
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Tianwei Xu
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Na Zhao
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Xueping Han
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
- Technology Extension Service of Animal Husbandry of Qinghai, Xining, China
| | - Shixiao Xu
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Xinquan Zhao
- Key Laboratory of Adaptation and Evolution of Plateau Biota, Northwest Institute of Plateau Biology, Chinese Academy of Sciences, Xining, China
| | - Cunfang Zhang
- State Key Laboratory of Plateau Ecology and Agriculture, Qinghai University, Xining, China
- *Correspondence: Cunfang Zhang,
| |
Collapse
|
13
|
Moradi MH, Mahmodi R, Farahani AHK, Karimi MO. Genome-wide evaluation of copy gain and loss variations in three Afghan sheep breeds. Sci Rep 2022; 12:14286. [PMID: 35996004 PMCID: PMC9395407 DOI: 10.1038/s41598-022-18571-4] [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: 02/09/2022] [Accepted: 08/16/2022] [Indexed: 11/23/2022] Open
Abstract
Copy number variation (CNV) is one of the main sources of variation between different individuals that has recently attracted much researcher interest as a major source for heritable variation in complex traits. The aim of this study was to identify CNVs in Afghan indigenous sheep consisting of three Arab, Baluchi, and Gadik breeds using genomic arrays containing 53,862 single nucleotide polymorphism (SNP) markers. Data were analyzed using the Hidden Markov Model (HMM) of PennCNV software. In this study, out of 45 sheep studied, 97.8% (44 animals) have shown CNVs. In total, 411 CNVs were observed for autosomal chromosomes and the entire sequence length of around 144 Mb was identified across the genome. The average number of CNVs per each sheep was 9.13. The identified CNVs for Arab, Baluchi, and Gadik breeds were 306, 62, and 43, respectively. After merging overlapped regions, a total of 376 copy number variation regions (CNVR) were identified, which are 286, 50, and 40 for Arab, Baluchi, and Gadik breeds, respectively. Bioinformatics analysis was performed to identify the genes and QTLs reported in these regions and the biochemical pathways involved by these genes. The results showed that many of these CNVRs overlapped with the genes or QTLs that are associated with various pathways such as immune system development, growth, reproduction, and environmental adaptions. Furthermore, to determine a genome-wide pattern of selection signatures in Afghan sheep breeds, the unbiased estimates of FST was calculated and the results indicated that 37 of the 376 CNVRs (~ 10%) have been also under selection signature, most of those overlapped with the genes influencing production, reproduction and immune system. Finally, the statistical methods used in this study was applied in an external dataset including 96 individuals of the Iranian sheep breed. The results indicated that 20 of the 114 CNVRs (18%) identified in Iranian sheep breed were also identified in our study, most of those overlapped with the genes influencing production, reproduction and immune system. Overall, this is the first attempts to develop the genomic map of loss and gain variation in the genome of Afghan indigenous sheep breeds, and may be important to shed some light on the genomic regions associated with some economically important traits in these breeds.
Collapse
Affiliation(s)
- Mohammad Hossein Moradi
- Department of Animal Science, Faculty of Agriculture and Natural Resources, Arak University, Arak, 38156-8-8349, Iran.
| | - Roqiah Mahmodi
- Department of Animal Science, Faculty of Agriculture and Natural Resources, Arak University, Arak, 38156-8-8349, Iran
| | | | - Mohammad Osman Karimi
- Department of Animal Science, Faculty of Agriculture and Natural Resources, Herat University, Herat, Afghanistan
| |
Collapse
|
14
|
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
|
15
|
Ahmad SF, Singh A, Panda S, Malla WA, Kumar A, Dutt T. Genome-wide elucidation of CNV regions and their association with production and reproduction traits in composite Vrindavani cattle. Gene 2022; 830:146510. [PMID: 35447249 DOI: 10.1016/j.gene.2022.146510] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2021] [Revised: 03/23/2022] [Accepted: 04/14/2022] [Indexed: 11/17/2022]
Abstract
The present study was aimed to analyze the genome-wide copy number variations (CNVs) in Vrindavani composite cattle and concatenate them into CNV regions (CNVRs), and finally test the association of CNVRs with different production and reproduction traits. Genotypic data, generated on BovineSNP50 Beadchip (v3) array for 96 Vrindavani animals, was used to elucidate the CNVs at the genome level. Intensity data covering over 53,218 SNP genotypes on bovine genome was used. Algorithm based on Hidden Markov Model was employed in PennCNV program to detect, normalize and filter CNVs across the genome. 252 putative CNVs, detected via PennCNV program, in different individuals were concatenated into 71 CNV regions (CNVRs) using CNVRuler program. Association of CNVRs with important (re)production traits in Vrindavani animals was assessed using linear regression. Five CNVRs were found to be significantly associated with ten important (re)production traits. The genes harbored in these regions provided useful insights into the association of CNVRs with genes and ultimately the variation at phenotype level. Important genes that overlapped with CNVRs included WASHC4, HS6ST3, MBNL2, TOLLIP, PIDD1 and TSPAN4. Furthermore, the CNVRs were found to overlap with important QTLs available in AnimalQTL database which affect milk yield and composition along with reproduction and immune function traits. The copy number states of three enes were validated using digital droplet PCR technique. The results from the present study significantly enhance the understanding about CNVs in Vrindavani cattle and should help establish its CNV map. The study will also enable further investigation on association of these variants with important traits of economic interest including disease incidence.
Collapse
Affiliation(s)
- Sheikh Firdous Ahmad
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, Uttar Pradesh, India.
| | - Akansha Singh
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, Uttar Pradesh, India
| | - Snehasmita Panda
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, Uttar Pradesh, India
| | - Waseem Akram Malla
- Division of Veterinary Biotechnology, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, Uttar Pradesh, India
| | - Amit Kumar
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, Uttar Pradesh, India.
| | - Triveni Dutt
- Livestock Production and Management Section, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, Uttar Pradesh, India
| |
Collapse
|
16
|
Guo X, Pei J, Wu X, Bao P, Ding X, Xiong L, Chu M, Lan X, Yan P. Detection of InDel and CNV of SPAG17 gene and their associations with bovine growth traits. Anim Biotechnol 2022; 33:440-447. [PMID: 32820682 DOI: 10.1080/10495398.2020.1803342] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Sperm-associated antigen 17 (SPAG17) gene encodes a central pair protein, which is involved in flagellar motility, male fertility and skeletal growth in ruminants. The insertions/deletions (indels) and copy number variations (CNVs) influence phenotypic traits by altering the sequences and copy numbers of functional genes, respectively. This study identified a novel 8-bp indel of SPAG17 gene in 1520 individuals from eight different cattle breeds, as well as a novel CNV region in 355 animals. The correlation analysis of indel showed that the individuals of ID genotype had superior performance traits such as body height (p = 0.038) and body slanting length (p = 0.041) as compared to other genotypes in Xianan cattle. For the CNV, different copy numbers were closely related to the body height in Qinchuan (p = 0.045) and body weight in Xianan (p = 0.036) breeds. Importantly, significant difference was observed between the 8-bp indel and the copy number loss in Xianan breed (p < 0.01). These findings indicated that the variations within the bovine SPAG17 gene can be considered as an effective DNA molecular marker for beef cattle breeding.
Collapse
Affiliation(s)
- Xian Guo
- Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, Gansu, China
- Key Laboratory of Yak Breeding Engineering of Gansu Province, Lanzhou, Gansu, China
| | - Jie Pei
- Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, Gansu, China
- Key Laboratory of Yak Breeding Engineering of Gansu Province, Lanzhou, Gansu, China
| | - Xiaoyun Wu
- Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, Gansu, China
- Key Laboratory of Yak Breeding Engineering of Gansu Province, Lanzhou, Gansu, China
| | - Pengjia Bao
- Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, Gansu, China
- Key Laboratory of Yak Breeding Engineering of Gansu Province, Lanzhou, Gansu, China
| | - Xuezhi Ding
- Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, Gansu, China
- Key Laboratory of Yak Breeding Engineering of Gansu Province, Lanzhou, Gansu, China
| | - Lin Xiong
- Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, Gansu, China
- Key Laboratory of Yak Breeding Engineering of Gansu Province, Lanzhou, Gansu, China
| | - Min Chu
- Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, Gansu, China
- Key Laboratory of Yak Breeding Engineering of Gansu Province, Lanzhou, Gansu, China
| | - Xianyong Lan
- Shaanxi Key Laboratory of Molecular Biology for Agriculture, College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China
| | - Ping Yan
- Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, Gansu, China
- Key Laboratory of Yak Breeding Engineering of Gansu Province, Lanzhou, Gansu, China
| |
Collapse
|
17
|
Wang M, Liu Y, Bi X, Ma H, Zeng G, Guo J, Guo M, Ling Y, Zhao C. Genome-Wide Detection of Copy Number Variants in Chinese Indigenous Horse Breeds and Verification of CNV-Overlapped Genes Related to Heat Adaptation of the Jinjiang Horse. Genes (Basel) 2022; 13:genes13040603. [PMID: 35456409 PMCID: PMC9033042 DOI: 10.3390/genes13040603] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2022] [Revised: 03/23/2022] [Accepted: 03/25/2022] [Indexed: 02/05/2023] Open
Abstract
In the present study, genome-wide CNVs were detected in a total of 301 samples from 10 Chinese indigenous horse breeds using the Illumina Equine SNP70 Bead Array, and the candidate genes related to adaptability to high temperature and humidity in Jinjiang horses were identified and validated. We determined a total of 577 CNVs ranging in size from 1.06 Kb to 2023.07 Kb on the 31 pairs of autosomes. By aggregating the overlapping CNVs for each breed, a total of 495 CNVRs were detected in the 10 Chinese horse breeds. As many as 211 breed-specific CNVRs were determined, of which 64 were found in the Jinjiang horse population. By removing repetitive CNV regions between breeds, a total of 239 CNVRs were identified in the Chinese indigenous horse breeds including 102 losses, 133 gains and 4 of both events (losses and gains in the same region), in which 131 CNVRs were novel and only detected in the present study compared with previous studies. The total detected CNVR length was 41.74 Mb, accounting for 1.83% of the total length of equine autosomal chromosomes. The coverage of CNVRs on each chromosome varied from 0.47% to 15.68%, with the highest coverage on ECA 12, but the highest number of CNVRs was detected on ECA1 and ECA24. A total of 229 genes overlapping with CNVRs were detected in the Jinjiang horse population, which is an indigenous horse breed unique to the southeastern coast of China exhibiting adaptability to high temperature and humidity. The functional annotation of these genes showed significant relation to cellular heat acclimation and immunity. The expression levels of the candidate genes were validated by heat shock treatment of various durations on fibroblasts of horses. The results show that the expression levels of HSPA1A were significantly increased among the different heat shock durations. The expression level of NFKBIA and SOCS4 declined from the beginning of heat shock to 2 h after heat shock and then showed a gradual increase until it reached the highest value at 6 h and 10 h of heat shock, respectively. Breed-specific CNVRs of Chinese indigenous horse breeds were revealed in the present study, and the results facilitate mapping CNVs on the whole genome and also provide valuable insights into the molecular mechanisms of adaptation to high temperature and humidity in the Jinjiang horse.
Collapse
Affiliation(s)
- Min Wang
- College of Animal Science and Technology, China Agricultural University, Beijing 100193, China; (M.W.); (Y.L.); (X.B.); (Y.L.)
- Equine Center, China Agricultural University, Beijing 100193, China
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, Beijing 100193, China
- National Engineering Laboratory for Animal Breeding, Beijing 100193, China
- Beijing Key Laboratory for Genetic Improvement of Livestock and Poultry, Beijing 100193, China
| | - Yu Liu
- College of Animal Science and Technology, China Agricultural University, Beijing 100193, China; (M.W.); (Y.L.); (X.B.); (Y.L.)
- Equine Center, China Agricultural University, Beijing 100193, China
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, Beijing 100193, China
- National Engineering Laboratory for Animal Breeding, Beijing 100193, China
- Beijing Key Laboratory for Genetic Improvement of Livestock and Poultry, Beijing 100193, China
| | - Xiaokun Bi
- College of Animal Science and Technology, China Agricultural University, Beijing 100193, China; (M.W.); (Y.L.); (X.B.); (Y.L.)
- Equine Center, China Agricultural University, Beijing 100193, China
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, Beijing 100193, China
- National Engineering Laboratory for Animal Breeding, Beijing 100193, China
- Beijing Key Laboratory for Genetic Improvement of Livestock and Poultry, Beijing 100193, China
| | - Hongying Ma
- Shaanxi Key Laboratory for Animal Conservation, Shaanxi Institute of Zoology, Xi’an 710032, China;
| | - Guorong Zeng
- Jinjiang Animal Husbandry and Veterinary Station, Quanzhou 362200, China; (G.Z.); (J.G.); (M.G.)
| | - Jintu Guo
- Jinjiang Animal Husbandry and Veterinary Station, Quanzhou 362200, China; (G.Z.); (J.G.); (M.G.)
| | - Minghao Guo
- Jinjiang Animal Husbandry and Veterinary Station, Quanzhou 362200, China; (G.Z.); (J.G.); (M.G.)
| | - Yao Ling
- College of Animal Science and Technology, China Agricultural University, Beijing 100193, China; (M.W.); (Y.L.); (X.B.); (Y.L.)
| | - Chunjiang Zhao
- College of Animal Science and Technology, China Agricultural University, Beijing 100193, China; (M.W.); (Y.L.); (X.B.); (Y.L.)
- Equine Center, China Agricultural University, Beijing 100193, China
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, Beijing 100193, China
- National Engineering Laboratory for Animal Breeding, Beijing 100193, China
- Beijing Key Laboratory for Genetic Improvement of Livestock and Poultry, Beijing 100193, China
- Correspondence:
| |
Collapse
|
18
|
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
|
19
|
Identification of Copy Number Variations and Genetic Diversity in Italian Insular Sheep Breeds. Animals (Basel) 2022; 12:ani12020217. [PMID: 35049839 PMCID: PMC8773107 DOI: 10.3390/ani12020217] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2021] [Revised: 01/14/2022] [Accepted: 01/16/2022] [Indexed: 02/05/2023] Open
Abstract
Copy number variants (CNVs) are one of the major contributors to genetic diversity and phenotypic variation in livestock. The aim of this work is to identify CNVs and perform, for the first time, a CNV-based population genetics analysis with five Italian sheep breeds (Barbaresca, Comisana, Pinzirita, Sarda, and Valle del Belìce). We identified 10,207 CNVs with an average length of 1.81 Mb. The breeds showed similar mean numbers of CNVs, ranging from 20 (Sarda) to 27 (Comisana). A total of 365 CNV regions (CNVRs) were determined. The length of the CNVRs varied among breeds from 2.4 Mb to 124.1 Mb. The highest number of shared CNVRs was between Comisana and Pinzirita, and only one CNVR was shared among all breeds. Our results indicated that segregating CNVs expresses a certain degree of diversity across all breeds. Despite the low/moderate genetic differentiation among breeds, the different approaches used to disclose the genetic relationship showed that the five breeds tend to cluster in distinct groups, similar to the previous studies based on single-nucleotide polymorphism markers. Gene enrichment was described for the 37 CNVRs selected, considering the top 10%. Out of 181 total genes, 67 were uncharacterized loci. Gene Ontology analysis showed that several of these genes are involved in lipid metabolism, immune response, and the olfactory pathway. Our results corroborated previous studies and showed that CNVs represent valuable molecular resources for providing useful information for separating the population and could be further used to explore the function and evolutionary aspect of sheep genome.
Collapse
|
20
|
Igoshin AV, Deniskova TE, Yurchenko AA, Yudin NS, Dotsev AV, Selionova MI, Zinovieva NA, Larkin DM. Copy number variants in genomes of local sheep breeds from Russia. Anim Genet 2021; 53:119-132. [PMID: 34904242 DOI: 10.1111/age.13163] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/28/2021] [Indexed: 01/21/2023]
Abstract
Copy number variants (CNVs) are genomic structural variations that contribute to many adaptive and economically important traits in livestock. In this study, we detected CNVs in 354 animals from 16 Russian indigenous sheep breeds and analysed their possible functional roles. Our analysis of the entire sample set resulted in 4527 CNVs forming 1450 CNV regions (CNVRs). When constructing CNVRs for individual breeds, a total of 2715 regions ranging from 88 in Groznensk to 337 in Osetin breeds were identified. To make interbreed CNVR frequency comparison possible, we also identified core CNVRs using CNVs with overlapping chromosomal locations found in different breeds. This resulted in 137 interbreed CNVRs with frequency >15% in at least one breed. Functional enrichment analysis of genes affected by CNVRs in individual breeds revealed 12 breeds with significant enrichments in olfactory perception, PRAME family proteins, and immune response. Function of genes affected by interbreed and breed-specific CNVRs revealed candidates related to domestication, adaptation to high altitudes and cold climates, reproduction, parasite resistance, milk and meat qualities, wool traits, fat storage, and fat metabolism. Our work is the first attempt to uncover and characterise the CNV makeup of Russian indigenous sheep breeds. Further experimental and functional validation of CNVRs would help in developing new and improving existing sheep breeds.
Collapse
Affiliation(s)
- A V Igoshin
- The Federal Research Center Institute of Cytology and Genetics SB RAS, Novosibirsk, 630090, Russia
| | - T E Deniskova
- L.K. Ernst Federal Research Center for Animal Husbandry, Podolsk, 142132, Russia
| | - A A Yurchenko
- The Federal Research Center Institute of Cytology and Genetics SB RAS, Novosibirsk, 630090, Russia
| | - N S Yudin
- The Federal Research Center Institute of Cytology and Genetics SB RAS, Novosibirsk, 630090, Russia.,Novosibirsk State University, Novosibirsk, 630090, Russia
| | - A V Dotsev
- L.K. Ernst Federal Research Center for Animal Husbandry, Podolsk, 142132, Russia
| | - M I Selionova
- Russian State Agrarian University, Moscow Timiryazev Agricultural Academy, Moscow, 127550, Russia
| | - N A Zinovieva
- L.K. Ernst Federal Research Center for Animal Husbandry, Podolsk, 142132, Russia
| | - D M Larkin
- The Federal Research Center Institute of Cytology and Genetics SB RAS, Novosibirsk, 630090, Russia.,Royal Veterinary College, University of London, London, NW1 0TU, UK
| |
Collapse
|
21
|
Dixit SP, Bhatia AK, Ganguly I, Singh S, Dash S, Sharma A, Anandkumar N, Dang AK, Jayakumar S. Genome analyses revealed genetic admixture and selection signatures in Bos indicus. Sci Rep 2021; 11:21924. [PMID: 34753978 PMCID: PMC8578574 DOI: 10.1038/s41598-021-01144-2] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2021] [Accepted: 10/21/2021] [Indexed: 11/20/2022] Open
Abstract
The genomic diversity and relationship among seven diverse cattle breeds viz. Sahiwal, Tharparkar, Gir, Vechur, Ongole, Kangayam and Hariana were investigated in 132 random samples based on high density SNP array comprising > 777 K SNPs. A total of 1993 SNPs (0.25% of the total) having greater power (FST ≥ 0.20) to differentiate these cattle populations were identified, and utilized to partition genome of each animal into a predefined number of clusters. The structure of these cattle indicated shared ancestry of dairy breeds viz. Gir, Tharparkar and Sahiwal. Most of the animals (> 76%) of different populations under study except Vechur clustered into their own group of animals called breed. Vechur population retained highest rate of admixture, consistent with its crossing with other breeds. Ongole, Kangayam and Hariana shared comparatively less of their genome (≤ 15%) with other breeds. The study indicated that all seven breeds evolved from their independent ancestry but there was intermixing of these breeds in the recent past. The selection signatures identified between draft (Kangayam) and dairy breeds included several genes like FAM19A2, RAB31P, BEST3, DGKA, AHCY, PIGU and PFKP which are involved in immune response, metabolic pathway, transportation of glucose and sugars, signaling pathways, cellular processes, cell division and glycolysis regulation, respectively. Moreover, these genomic regions also harbour QTLs affecting milk performance traits. The signatures were also identified even between the dairy breeds. In comparison to large-sized cattle, there were significant differences in the number of QTLs affecting production (body weight, growth rate etc.) and morphological traits (height) in short-statured Vechur breed. The presence of HMGA2 gene in the selection signature on chromosome 5 may explain the variations in stature between these cattle.
Collapse
Affiliation(s)
- S P Dixit
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India.
| | - A K Bhatia
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - Indrajit Ganguly
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - Sanjeev Singh
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - Soumya Dash
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - Anurodh Sharma
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - N Anandkumar
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - A K Dang
- ICAR - National Dairy Research Institute, Karnal, Haryana, 132001, India
| | - S Jayakumar
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| |
Collapse
|
22
|
Kava R, Peripolli E, Berton MP, Lemos M, Lobo RB, Stafuzza NB, Pereira AS, Baldi F. Genome-wide structural variations in Brazilian Senepol cattle, a tropically adapted taurine breed. Livest Sci 2021. [DOI: 10.1016/j.livsci.2021.104708] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
23
|
Guo S, Wu X, Pei J, Wang X, Bao P, Xiong L, Chu M, Liang C, Yan P, Guo X. Genome-wide CNV analysis reveals variants associated with high-altitude adaptation and meat traits in Qaidam cattle. ELECTRON J BIOTECHN 2021. [DOI: 10.1016/j.ejbt.2021.07.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022] Open
|
24
|
Upadhyay M, Derks MFL, Andersson G, Medugorac I, Groenen MAM, Crooijmans RPMA. Introgression contributes to distribution of structural variations in cattle. Genomics 2021; 113:3092-3102. [PMID: 34242710 DOI: 10.1016/j.ygeno.2021.07.005] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2021] [Revised: 06/24/2021] [Accepted: 07/03/2021] [Indexed: 11/19/2022]
Abstract
Structural variations (SVs) are an important source of phenotypic diversity in cattle. Here, 72 whole genome sequences representing taurine and zebu cattle were used to identify SVs. Applying multiple approaches, 16,738 SVs were identified. A comparison against the Database of Genomic Variants archives revealed that 1575 SVs were novel in our data. A novel duplication covering the entire GALNT15 gene, was observed only in N'Dama. A duplication, which was previously reported only in zebu and associated with navel length, was also observed in N'Dama. Investigation of a novel deletion located upstream of CAST13 gene and identified only in Italian cattle and zebu, revealed its introgressed origin in the former. Overall, our data highlights how the SVs distribution in cattle is also shaped by forces such as demographical differences and gene flow. The cattle SVs of this study and its meta-data can be visualized on an interactive genome browser at https://tinyurl.com/svCowArs.
Collapse
Affiliation(s)
- Maulik Upadhyay
- Animal Breeding and Genomics, Wageningen University & Research, Droevendaalsesteeg 1, 6708PB, Wageningen, The Netherlands; Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, 75007 Uppsala, Sweden; Population Genomics Group, Department of Veterinary Sciences, Ludwig-Maximilians-University Munich, 80539 Munich, Germany.
| | - Martijn F L Derks
- Animal Breeding and Genomics, Wageningen University & Research, Droevendaalsesteeg 1, 6708PB, Wageningen, The Netherlands.
| | - Göran Andersson
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, 75007 Uppsala, Sweden.
| | - Ivica Medugorac
- Population Genomics Group, Department of Veterinary Sciences, Ludwig-Maximilians-University Munich, 80539 Munich, Germany.
| | - Martien A M Groenen
- Animal Breeding and Genomics, Wageningen University & Research, Droevendaalsesteeg 1, 6708PB, Wageningen, The Netherlands.
| | - Richard P M A Crooijmans
- Animal Breeding and Genomics, Wageningen University & Research, Droevendaalsesteeg 1, 6708PB, Wageningen, The Netherlands.
| |
Collapse
|
25
|
Kumar H, Panigrahi M, Saravanan KA, Rajawat D, Parida S, Bhushan B, Gaur GK, Dutt T, Mishra BP, Singh RK. Genome-wide detection of copy number variations in Tharparkar cattle. Anim Biotechnol 2021; 34:448-455. [PMID: 34191685 DOI: 10.1080/10495398.2021.1942027] [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] [Indexed: 10/21/2022]
Abstract
Copy number variations (CNVs) are major forms of genetic variation with an increasing importance in animal genomics. This study used the Illumina BovineSNP 50 K BeadChip to detect the genome-wide CNVs in the Tharparkar cattle. With the aid of PennCNV software, we noticed a total of 447 copy number variation regions (CNVRs) across the autosomal genome, occupying nearly 2.17% of the bovine genome. The average size of detected CNVRs was found to be 122.2 kb, the smallest CNVR being 50.02 kb in size, to the largest being 1,232.87 Kb. Enrichment analyses of the genes in these CNVRs gave significant associations with molecular adaptation-related Gene Ontology (GO) terms. Most CNVR genes were significantly enriched for specific biological functions; signaling pathways, sensory responses to stimuli, and various cellular processes. In addition, QTL analysis of CNVRs described them to be linked with economically essential traits in cattle. The findings here provide crucial information for constructing a more comprehensive CNVR map for the indigenous cattle genome.
Collapse
Affiliation(s)
- Harshit Kumar
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, India
| | - Manjit Panigrahi
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, India
| | - K A Saravanan
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, India
| | - Divya Rajawat
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, India
| | - Subhashree Parida
- Division of Pharmacology & Toxicology, ICAR-Indian Veterinary Research Institute, Bareilly, India
| | - Bharat Bhushan
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, India
| | - G K Gaur
- Division of Animal Genetics, Indian Veterinary Research Institute, Bareilly, India
| | - Triveni Dutt
- Livestock Production and Management Section, ICAR-Indian Veterinary Research Institute, Bareilly, India
| | - B P Mishra
- Division of Animal Biotechnology, ICAR-Indian Veterinary Research Institute, Bareilly, India
| | - R K Singh
- Division of Animal Biotechnology, ICAR-Indian Veterinary Research Institute, Bareilly, India
| |
Collapse
|
26
|
Qiu Y, Ding R, Zhuang Z, Wu J, Yang M, Zhou S, Ye Y, Geng Q, Xu Z, Huang S, Cai G, Wu Z, Yang J. Genome-wide detection of CNV regions and their potential association with growth and fatness traits in Duroc pigs. BMC Genomics 2021; 22:332. [PMID: 33964879 PMCID: PMC8106131 DOI: 10.1186/s12864-021-07654-7] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2020] [Accepted: 04/27/2021] [Indexed: 11/19/2022] Open
Abstract
Background In the process of pig breeding, the average daily gain (ADG), days to 100 kg (AGE), and backfat thickness (BFT) are directly related to growth rate and fatness. However, the genetic mechanisms involved are not well understood. Copy number variation (CNV), an important source of genetic diversity, can affect a variety of complex traits and diseases and has gradually been thrust into the limelight. In this study, we reported the genome-wide CNVs of Duroc pigs using SNP genotyping data from 6627 animals. We also performed a copy number variation region (CNVR)-based genome-wide association studies (GWAS) for growth and fatness traits in two Duroc populations. Results Our study identified 953 nonredundant CNVRs in U.S. and Canadian Duroc pigs, covering 246.89 Mb (~ 10.90%) of the pig autosomal genome. Of these, 802 CNVRs were in U.S. Duroc pigs with 499 CNVRs were in Canadian Duroc pigs, indicating 348 CNVRs were shared by the two populations. Experimentally, 77.8% of nine randomly selected CNVRs were validated through quantitative PCR (qPCR). We also identified 35 CNVRs with significant association with growth and fatness traits using CNVR-based GWAS. Ten of these CNVRs were associated with both ADG and AGE traits in U.S. Duroc pigs. Notably, four CNVRs showed significant associations with ADG, AGE, and BFT, indicating that these CNVRs may play a pleiotropic role in regulating pig growth and fat deposition. In Canadian Duroc pigs, nine CNVRs were significantly associated with both ADG and AGE traits. Further bioinformatic analysis identified a subset of potential candidate genes, including PDGFA, GPER1, PNPLA2 and BSCL2. Conclusions The present study provides a necessary supplement to the CNV map of the Duroc genome through large-scale population genotyping. In addition, the CNVR-based GWAS results provide a meaningful way to elucidate the genetic mechanisms underlying complex traits. The identified CNVRs can be used as molecular markers for genetic improvement in the molecular-guided breeding of modern commercial pigs. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07654-7.
Collapse
Affiliation(s)
- Yibin Qiu
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, Guangdong, 510642, People's Republic of China
| | - Rongrong Ding
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, Guangdong, 510642, People's Republic of China.,Guangdong Wens Breeding Swine Technology Co., Ltd., Yunfu, Guangdong, 527400, People's Republic of China
| | - Zhanwei Zhuang
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, Guangdong, 510642, People's Republic of China
| | - Jie Wu
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, Guangdong, 510642, People's Republic of China
| | - Ming Yang
- Guangdong Wens Breeding Swine Technology Co., Ltd., Yunfu, Guangdong, 527400, People's Republic of China
| | - Shenping Zhou
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, Guangdong, 510642, People's Republic of China
| | - Yong Ye
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, Guangdong, 510642, People's Republic of China
| | - Qian Geng
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, Guangdong, 510642, People's Republic of China
| | - Zheng Xu
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, Guangdong, 510642, People's Republic of China.,Lingnan Guangdong Laboratory of Modern Agriculture, Guangzhou, 510642, People's Republic of China
| | - Sixiu Huang
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, Guangdong, 510642, People's Republic of China
| | - Gengyuan Cai
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, Guangdong, 510642, People's Republic of China.,Guangdong Wens Breeding Swine Technology Co., Ltd., Yunfu, Guangdong, 527400, People's Republic of China.,Lingnan Guangdong Laboratory of Modern Agriculture, Guangzhou, 510642, People's Republic of China
| | - Zhenfang Wu
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, Guangdong, 510642, People's Republic of China. .,Guangdong Wens Breeding Swine Technology Co., Ltd., Yunfu, Guangdong, 527400, People's Republic of China. .,Lingnan Guangdong Laboratory of Modern Agriculture, Guangzhou, 510642, People's Republic of China.
| | - Jie Yang
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, Guangdong, 510642, People's Republic of China. .,Lingnan Guangdong Laboratory of Modern Agriculture, Guangzhou, 510642, People's Republic of China.
| |
Collapse
|
27
|
Serres-Armero A, Davis BW, Povolotskaya IS, Morcillo-Suarez C, Plassais J, Juan D, Ostrander EA, Marques-Bonet T. Copy number variation underlies complex phenotypes in domestic dog breeds and other canids. Genome Res 2021; 31:762-774. [PMID: 33863806 PMCID: PMC8092016 DOI: 10.1101/gr.266049.120] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2020] [Accepted: 02/26/2021] [Indexed: 01/02/2023]
Abstract
Extreme phenotypic diversity, a history of artificial selection, and socioeconomic value make domestic dog breeds a compelling subject for genomic research. Copy number variation (CNV) is known to account for a significant part of inter-individual genomic diversity in other systems. However, a comprehensive genome-wide study of structural variation as it relates to breed-specific phenotypes is lacking. We have generated whole genome CNV maps for more than 300 canids. Our data set extends the canine structural variation landscape to more than 100 dog breeds, including novel variants that cannot be assessed using microarray technologies. We have taken advantage of this data set to perform the first CNV-based genome-wide association study (GWAS) in canids. We identify 96 loci that display copy number differences across breeds, which are statistically associated with a previously compiled set of breed-specific morphometrics and disease susceptibilities. Among these, we highlight the discovery of a long-range interaction involving a CNV near MED13L and TBX3, which could influence breed standard height. Integration of the CNVs with chromatin interactions, long noncoding RNA expression, and single nucleotide variation highlights a subset of specific loci and genes with potential functional relevance and the prospect to explain trait variation between dog breeds.
Collapse
Affiliation(s)
- Aitor Serres-Armero
- IBE, Institut de Biologia Evolutiva (Universitat Pompeu Fabra/CSIC), Ciencies Experimentals i de la Salut, Barcelona 08003, Spain
| | - Brian W Davis
- Cancer Genetics and Comparative Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, Maryland 20892, USA.,Department of Veterinary Integrative Biosciences, College of Veterinary Medicine, Texas A&M University, College Station, Texas 77843, USA
| | - Inna S Povolotskaya
- Veltischev Research and Clinical Institute for Pediatrics of the Pirogov Russian National Research Medical University, Moscow 117997, Russia
| | - Carlos Morcillo-Suarez
- IBE, Institut de Biologia Evolutiva (Universitat Pompeu Fabra/CSIC), Ciencies Experimentals i de la Salut, Barcelona 08003, Spain
| | - Jocelyn Plassais
- Cancer Genetics and Comparative Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - David Juan
- IBE, Institut de Biologia Evolutiva (Universitat Pompeu Fabra/CSIC), Ciencies Experimentals i de la Salut, Barcelona 08003, Spain
| | - Elaine A Ostrander
- Cancer Genetics and Comparative Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Tomas Marques-Bonet
- IBE, Institut de Biologia Evolutiva (Universitat Pompeu Fabra/CSIC), Ciencies Experimentals i de la Salut, Barcelona 08003, Spain.,CNAG-CRG, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona 08028, Spain.,Institucio Catalana de Recerca i Estudis Avançats (ICREA), Barcelona, Catalonia 08010, Spain.,Institut Català de Paleontologia Miquel Crusafont, Universitat Autònoma de Barcelona, Cerdanyola del Vallès, Catalonia 08201, Spain
| |
Collapse
|
28
|
Yuan C, Lu Z, Guo T, Yue Y, Wang X, Wang T, Zhang Y, Hou F, Niu C, Sun X, Zhao H, Zhu S, Liu J, Yang B. A global analysis of CNVs in Chinese indigenous fine-wool sheep populations using whole-genome resequencing. BMC Genomics 2021; 22:78. [PMID: 33485316 PMCID: PMC7825165 DOI: 10.1186/s12864-021-07387-7] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2019] [Accepted: 01/13/2021] [Indexed: 12/13/2022] Open
Abstract
Background Copy number variation (CNV) is an important source of genetic variation that has a significant influence on phenotypic diversity, economically important traits and the evolution of livestock species. In this study, the genome-wide CNV distribution characteristics of 32 fine-wool sheep from three breeds were analyzed using resequencing. Results A total of 1,747,604 CNVs were detected in this study, and 7228 CNV regions (CNVR) were obtained after merging overlapping CNVs; these regions accounted for 2.17% of the sheep reference genome. The average length of the CNVRs was 4307.17 bp. “Deletion” events took place more frequently than “duplication” or “both” events. The CNVRs obtained overlapped with previously reported sheep CNVRs to variable extents (4.39–55.46%). Functional enrichment analysis showed that the CNVR-harboring genes were mainly involved in sensory perception systems, nutrient metabolism processes, and growth and development processes. Furthermore, 1855 of the CNVRs were associated with 166 quantitative trait loci (QTL), including milk QTLs, carcass QTLs, and health-related QTLs, among others. In addition, the 32 fine-wool sheep were divided into horned and polled groups to analyze for the selective sweep of CNVRs, and it was found that the relaxin family peptide receptor 2 (RXFP2) gene was strongly influenced by selection. Conclusions In summary, we constructed a genomic CNV map for Chinese indigenous fine-wool sheep using resequencing, thereby providing a valuable genetic variation resource for sheep genome research, which will contribute to the study of complex traits in sheep. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07387-7.
Collapse
Affiliation(s)
- Chao Yuan
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences of Chinese Academy of Agricultural Sciences, Sheep Breeding Engineering Technology Research Center, Lanzhou, 730050, China
| | - Zengkui Lu
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences of Chinese Academy of Agricultural Sciences, Sheep Breeding Engineering Technology Research Center, Lanzhou, 730050, China
| | - Tingting Guo
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences of Chinese Academy of Agricultural Sciences, Sheep Breeding Engineering Technology Research Center, Lanzhou, 730050, China
| | - Yaojing Yue
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences of Chinese Academy of Agricultural Sciences, Sheep Breeding Engineering Technology Research Center, Lanzhou, 730050, China
| | - Xijun Wang
- Gansu Provincial Sheep Breeding Technology Extension Station, Sunan, 734031, China
| | - Tianxiang Wang
- Gansu Provincial Sheep Breeding Technology Extension Station, Sunan, 734031, China
| | - Yajun Zhang
- Xinjiang Gongnaisi Breeding Sheep Farm, Xinyuan, 835808, China
| | - Fujun Hou
- Aohan Banner Breeding Sheep Farm, Chifeng, 024300, China
| | - Chune Niu
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences of Chinese Academy of Agricultural Sciences, Sheep Breeding Engineering Technology Research Center, Lanzhou, 730050, China
| | - Xiaopin Sun
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences of Chinese Academy of Agricultural Sciences, Sheep Breeding Engineering Technology Research Center, Lanzhou, 730050, China
| | - Hongchang Zhao
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences of Chinese Academy of Agricultural Sciences, Sheep Breeding Engineering Technology Research Center, Lanzhou, 730050, China
| | - Shaohua Zhu
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences of Chinese Academy of Agricultural Sciences, Sheep Breeding Engineering Technology Research Center, Lanzhou, 730050, China
| | - Jianbin Liu
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences of Chinese Academy of Agricultural Sciences, Sheep Breeding Engineering Technology Research Center, Lanzhou, 730050, China.
| | - Bohui Yang
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences of Chinese Academy of Agricultural Sciences, Sheep Breeding Engineering Technology Research Center, Lanzhou, 730050, China.
| |
Collapse
|
29
|
Orlando L. The Evolutionary and Historical Foundation of the Modern Horse: Lessons from Ancient Genomics. Annu Rev Genet 2020; 54:563-581. [PMID: 32960653 DOI: 10.1146/annurev-genet-021920-011805] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The domestication of the horse some 5,500 years ago followed those of dogs, sheep, goats, cattle, and pigs by ∼2,500-10,000 years. By providing fast transportation and transforming warfare, the horse had an impact on human history with no equivalent in the animal kingdom. Even though the equine sport industry has considerable economic value today, the evolutionary history underlying the emergence of the modern domestic horse remains contentious. In the last decade, novel sequencing technologies have revolutionized our capacity to sequence the complete genome of organisms, including from archaeological remains. Applied to horses, these technologies have provided unprecedented levels of information and have considerably changed models of horse domestication. This review illustrates how ancient DNA, especially ancient genomes, has inspired researchers to rethink the process by which horses were first domesticated and then diversified into a variety of breeds showing a range of traits that are useful to humans.
Collapse
Affiliation(s)
- Ludovic Orlando
- Laboratoire d'Anthropobiologie Moléculaire et Imagerie de Synthèse, Faculté de Médecine Purpan, Université Toulouse III-Paul Sabatier, 31000 Toulouse, France;
| |
Collapse
|
30
|
Hu Y, Xia H, Li M, Xu C, Ye X, Su R, Zhang M, Nash O, Sonstegard TS, Yang L, Liu GE, Zhou Y. Comparative analyses of copy number variations between Bos taurus and Bos indicus. BMC Genomics 2020; 21:682. [PMID: 33004001 PMCID: PMC7528262 DOI: 10.1186/s12864-020-07097-6] [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: 02/21/2020] [Accepted: 09/23/2020] [Indexed: 12/15/2022] Open
Abstract
Background Bos taurus and Bos indicus are two main sub-species of cattle. However, the differential copy number variations (CNVs) between them are not yet well studied. Results Based on the new high-quality cattle reference genome ARS-UCD1.2, we identified 13,234 non-redundant CNV regions (CNVRs) from 73 animals of 10 cattle breeds (4 Bos taurus and 6 Bos indicus), by integrating three detection strategies. While 6990 CNVRs (52.82%) were shared by Bos taurus and Bos indicus, large CNV differences were discovered between them and these differences could be used to successfully separate animals into two subspecies. We found that 2212 and 538 genes uniquely overlapped with either indicine-specific CNVRs and or taurine-specific CNVRs, respectively. Based on FST, we detected 16 candidate lineage-differential CNV segments (top 0.1%) under selection, which overlapped with eight genes (CTNNA1, ENSBTAG00000004415, PKN2, BMPER, PDE1C, DNAJC18, MUSK, and PLCXD3). Moreover, we obtained 1.74 Mbp indicine-specific sequences, which could only be mapped on the Bos indicus reference genome UOA_Brahman_1. We found these sequences and their associated genes were related to heat resistance, lipid and ATP metabolic process, and muscle development under selection. We further analyzed and validated the top significant lineage-differential CNV. This CNV overlapped genes related to muscle cell differentiation, which might be generated from a retropseudogene of CTH but was deleted along Bos indicus lineage. Conclusions This study presents a genome wide CNV comparison between Bos taurus and Bos indicus. It supplied essential genome diversity information for understanding of adaptation and phenotype differences between the Bos taurus and Bos indicus populations.
Collapse
Affiliation(s)
- Yan Hu
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education & College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, 430070, China
| | - Han Xia
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education & College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, 430070, China
| | - Mingxun Li
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Building 306, Room 111, BARC-East, Beltsville, MD, 20705, USA.,College of Animal Science and Technology, Yangzhou University, Yangzhou, 225009, China
| | - Chang Xu
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education & College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, 430070, China
| | - Xiaowei Ye
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education & College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, 430070, China
| | - Ruixue Su
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education & College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, 430070, China
| | - Mai Zhang
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education & College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, 430070, China
| | - Oyekanmi Nash
- Centre for Genomics Research and Innovation, National Biotechnology Development Agency, Abuja, Nigeria
| | | | - Liguo Yang
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education & College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, 430070, China
| | - George E Liu
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Building 306, Room 111, BARC-East, Beltsville, MD, 20705, USA.
| | - Yang Zhou
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education & College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, 430070, China.
| |
Collapse
|
31
|
Guan D, Martínez A, Castelló A, Landi V, Luigi-Sierra MG, Fernández-Álvarez J, Cabrera B, Delgado JV, Such X, Jordana J, Amills M. A genome-wide analysis of copy number variation in Murciano-Granadina goats. Genet Sel Evol 2020; 52:44. [PMID: 32770942 PMCID: PMC7414533 DOI: 10.1186/s12711-020-00564-4] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2020] [Accepted: 07/28/2020] [Indexed: 02/06/2023] Open
Abstract
BACKGROUND In this work, our aim was to generate a map of the copy number variations (CNV) segregating in a population of Murciano-Granadina goats, the most important dairy breed in Spain, and to ascertain the main biological functions of the genes that map to copy number variable regions. RESULTS Using a dataset that comprised 1036 Murciano-Granadina goats genotyped with the Goat SNP50 BeadChip, we were able to detect 4617 and 7750 autosomal CNV with the PennCNV and QuantiSNP software, respectively. By applying the EnsembleCNV algorithm, these CNV were assembled into 1461 CNV regions (CNVR), of which 486 (33.3% of the total CNVR count) were consistently called by PennCNV and QuantiSNP and used in subsequent analyses. In this set of 486 CNVR, we identified 78 gain, 353 loss and 55 gain/loss events. The total length of all the CNVR (95.69 Mb) represented 3.9% of the goat autosomal genome (2466.19 Mb), whereas their size ranged from 2.0 kb to 11.1 Mb, with an average size of 196.89 kb. Functional annotation of the genes that overlapped with the CNVR revealed an enrichment of pathways related with olfactory transduction (fold-enrichment = 2.33, q-value = 1.61 × 10-10), ABC transporters (fold-enrichment = 5.27, q-value = 4.27 × 10-04) and bile secretion (fold-enrichment = 3.90, q-value = 5.70 × 10-03). CONCLUSIONS A previous study reported that the average number of CNVR per goat breed was ~ 20 (978 CNVR/50 breeds), which is much smaller than the number we found here (486 CNVR). We attribute this difference to the fact that the previous study included multiple caprine breeds that were represented by small to moderate numbers of individuals. Given the low frequencies of CNV (in our study, the average frequency of CNV is 1.44%), such a design would probably underestimate the levels of the diversity of CNV at the within-breed level. We also observed that functions related with sensory perception, metabolism and embryo development are overrepresented in the set of genes that overlapped with CNV, and that these loci often belong to large multigene families with tens, hundreds or thousands of paralogous members, a feature that could favor the occurrence of duplications or deletions by non-allelic homologous recombination.
Collapse
Affiliation(s)
- Dailu Guan
- Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Amparo Martínez
- Departamento de Genética, Universidad de Córdoba, 14071, Córdoba, Spain
| | - Anna Castelló
- Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain.,Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Vincenzo Landi
- Departamento de Genética, Universidad de Córdoba, 14071, Córdoba, Spain.,Department of Veterinary Medicine, University of Bari "Aldo Moro", SP. 62 per Casamassima km. 3, 70010, Valenzano, BA, Italy
| | - María Gracia Luigi-Sierra
- Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Javier Fernández-Álvarez
- Asociación Nacional de Criadores de Caprino de Raza Murciano-Granadina (CAPRIGRAN), 18340, Granada, Spain
| | - Betlem Cabrera
- Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain.,Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | | | - Xavier Such
- Group of Research in Ruminants (G2R), Department of Animal and Food Science, Universitat Autònoma de Barcelona (UAB), Bellaterra, Barcelona, Spain
| | - Jordi Jordana
- Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Marcel Amills
- Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain. .,Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain.
| |
Collapse
|
32
|
Hao D, Wang X, Thomsen B, Kadarmideen HN, Wang X, Lan X, Huang Y, Qi X, Chen H. Copy Number Variations and Expression Levels of Guanylate-Binding Protein 6 Gene Associated with Growth Traits of Chinese Cattle. Animals (Basel) 2020; 10:E566. [PMID: 32230930 PMCID: PMC7222342 DOI: 10.3390/ani10040566] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2020] [Revised: 03/18/2020] [Accepted: 03/22/2020] [Indexed: 11/16/2022] Open
Abstract
Association studies have indicated profound effects of copy number variations (CNVs) on various phenotypes in different species. In this study, we identified the CNV distributions and expression levels of guanylate-binding protein 6 (GBP6) associated with the growth traits of Chinese cattle. The results showed that the phenotypic values of body size and weight of Xianan (XN) cattle were higher than those of Nanyang (NY) cattle. The medium CNV types were mostly identified in the XN and NY breeds, but their CNV distributions were significantly different (adjusted p < 0.05). The association analysis revealed that the body weight, cannon circumference and chest circumference of XN cattle had significantly different values in different CNV types (p < 0.05), with CNV gain types (Log22-ΔΔCt > 0.5) displaying superior phenotypic values. We also found that transcription levels varied in different tissues (p < 0.001) and the CNV gain types showed the highest relative gene expression levels in the muscle tissue, consistent with the highest phenotypic values of body weight and cannon circumference among the three CNV types. Consequently, our results suggested that CNV gain types of GBP6 could be used as the candidate markers in the cattle-breeding program for growth traits.
Collapse
Affiliation(s)
- Dan Hao
- College of Animal Science and Technology, Northwest A&F University, Shaanxi Key Laboratory of Animal Genetics, Breeding and Reproduction, Yangling 712100, Shaanxi, China; (D.H.); (X.W.); (X.L.); (Y.H.)
- Department of Molecular Biology and Genetics, Aarhus University, 8000 Aarhus C, Denmark;
| | - Xiao Wang
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, 2800 Kongens Lyngby, Denmark; (X.W.); (H.N.K.)
| | - Bo Thomsen
- Department of Molecular Biology and Genetics, Aarhus University, 8000 Aarhus C, Denmark;
| | - Haja N. Kadarmideen
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, 2800 Kongens Lyngby, Denmark; (X.W.); (H.N.K.)
| | - Xiaogang Wang
- College of Animal Science and Technology, Northwest A&F University, Shaanxi Key Laboratory of Animal Genetics, Breeding and Reproduction, Yangling 712100, Shaanxi, China; (D.H.); (X.W.); (X.L.); (Y.H.)
| | - Xianyong Lan
- College of Animal Science and Technology, Northwest A&F University, Shaanxi Key Laboratory of Animal Genetics, Breeding and Reproduction, Yangling 712100, Shaanxi, China; (D.H.); (X.W.); (X.L.); (Y.H.)
| | - Yongzhen Huang
- College of Animal Science and Technology, Northwest A&F University, Shaanxi Key Laboratory of Animal Genetics, Breeding and Reproduction, Yangling 712100, Shaanxi, China; (D.H.); (X.W.); (X.L.); (Y.H.)
| | - Xinglei Qi
- Bureau of Animal Husbandry of Biyang County, Biyang 463700, Henan, China;
| | - Hong Chen
- College of Animal Science and Technology, Northwest A&F University, Shaanxi Key Laboratory of Animal Genetics, Breeding and Reproduction, Yangling 712100, Shaanxi, China; (D.H.); (X.W.); (X.L.); (Y.H.)
| |
Collapse
|
33
|
Dixit SP, Singh S, Ganguly I, Bhatia AK, Sharma A, Kumar NA, Dang AK, Jayakumar S. Genome-Wide Runs of Homozygosity Revealed Selection Signatures in Bos indicus. Front Genet 2020; 11:92. [PMID: 32153647 PMCID: PMC7046685 DOI: 10.3389/fgene.2020.00092] [Citation(s) in RCA: 39] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2019] [Accepted: 01/28/2020] [Indexed: 02/04/2023] Open
Abstract
Genome-wide runs of homozygosity (ROH) are suitable for understanding population history, calculating genomic inbreeding, deciphering genetic architecture of complex traits and diseases as well as identifying genes linked with agro-economic traits. Autozygosity and ROH islands, genomic regions with elevated ROH frequencies, were characterized in 112 animals of seven Indian native cattle breeds (B. indicus) using BovineHD BeadChip. In total, 4138 ROH were detected. The average number of ROH per animal was maximum in draft breed, Kangayam (63.62 ± 22.71) and minimum in dairy breed, Sahiwal (24.62 ± 11.03). The mean ROH length was maximum in Vechur (6.97 Mb) and minimum in Hariana (4.04 Mb). Kangayam revealed the highest ROH based inbreeding (FROH > 1Mb = 0.113 ± 0.059), whereas Hariana (FROH > 1Mb = 0.042 ± 0.031) and Sahiwal (FROH > 1Mb = 0.043 ± 0.048) showed the lowest. The high standard deviation observed in each breed highlights a considerable variability in autozygosity. Out of the total autozygous segments observed in each breed except Vechur, > 80% were of short length (< 8 Mb) and contributed almost 50% of the genome proportion under ROH. However, in Vechur cattle, long ROH contributed 75% of the genome proportion under ROH. ROH patterns revealed Hariana and Sahiwal breeds as less consanguineous, while recent inbreeding was apparent in Vechur. Maximum autozygosity observed in Kangayam is attributable to both recent and ancient inbreeding. The ROH islands were harbouring higher proportion of QTLs for production traits (20.68% vs. 14.64%; P≤ 0.05) but lower for reproductive traits (11.49% vs. 15.76%; P≤ 0.05) in dairy breeds compared to draft breed. In draft cattle, genes associated with resistant to diseases/higher immunity (LYZL1, SVIL, and GPX4) and stress tolerant (CCT4) were identified in ROH islands; while in dairy breeds, for milk production (PTGFR, CSN1S1, CSN2, CSN1S2, and CSN3). Significant difference in ROH islands among large and short statured breeds was observed at chromosome 3 and 5 involving genes like PTGFR and HMGA2 responsible for milk production and stature, respectively. PCA analysis on consensus ROH regions revealed distinct clustering of dairy, draft and short stature cattle breeds.
Collapse
Affiliation(s)
- S P Dixit
- Animal Genetics Division, ICAR-National Bureau of Animal Genetic Resources, Karnal, India
| | - Sanjeev Singh
- Animal Genetics Division, ICAR-National Bureau of Animal Genetic Resources, Karnal, India
| | - Indrajit Ganguly
- Animal Genetics Division, ICAR-National Bureau of Animal Genetic Resources, Karnal, India
| | - Avnish Kumar Bhatia
- Animal Genetics Division, ICAR-National Bureau of Animal Genetic Resources, Karnal, India
| | - Anurodh Sharma
- Animal Genetics Division, ICAR-National Bureau of Animal Genetic Resources, Karnal, India
| | - N Anand Kumar
- Animal Genetics & Breeding Division, ICAR-National Dairy Research Institute, Karnal, India
| | - Ajay Kumar Dang
- Animal Physiology Division, ICAR-National Dairy Research Institute, Karnal, India
| | - S Jayakumar
- Animal Genetics Division, ICAR-National Bureau of Animal Genetic Resources, Karnal, India
| |
Collapse
|
34
|
Zhang Y, Hu Y, Wang X, Jiang Q, Zhao H, Wang J, Ju Z, Yang L, Gao Y, Wei X, Bai J, Zhou Y, Huang J. Population Structure, and Selection Signatures Underlying High-Altitude Adaptation Inferred From Genome-Wide Copy Number Variations in Chinese Indigenous Cattle. Front Genet 2020; 10:1404. [PMID: 32117428 PMCID: PMC7033542 DOI: 10.3389/fgene.2019.01404] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2019] [Accepted: 12/23/2019] [Indexed: 12/12/2022] Open
Abstract
Copy number variations (CNVs) have been demonstrated as crucial substrates for evolution, adaptation and breed formation. Chinese indigenous cattle breeds exhibit a broad geographical distribution and diverse environmental adaptability. Here, we analyzed the population structure and adaptation to high altitude of Chinese indigenous cattle based on genome-wide CNVs derived from the high-density BovineHD SNP array. We successfully detected the genome-wide CNVs of 318 individuals from 24 Chinese indigenous cattle breeds and 37 yaks as outgroups. A total of 5,818 autosomal CNV regions (683 bp-4,477,860 bp in size), covering ~14.34% of the bovine genome (UMD3.1), were identified, showing abundant CNV resources. Neighbor-joining clustering, principal component analysis (PCA), and population admixture analysis based on these CNVs support that most Chinese cattle breeds are hybrids of Bos taurus taurus (hereinafter to be referred as Bos taurus) and Bos taurus indicus (Bos indicus). The distribution patterns of the CNVs could to some extent be related to the geographical backgrounds of the habitat of the breeds, and admixture among cattle breeds from different districts. We analyzed the selective signatures of CNVs positively involved in high-altitude adaptation using pairwise Fst analysis within breeds with a strong Bos taurus background (taurine-type breeds) and within Bos taurus×Bos indicus hybrids, respectively. CNV-overlapping genes with strong selection signatures (at top 0.5% of Fst value), including LETM1 (Fst = 0.490), TXNRD2 (Fst = 0.440), and STUB1 (Fst = 0.420) within taurine-type breeds, and NOXA1 (Fst = 0.233), RUVBL1 (Fst = 0.222), and SLC4A3 (Fst=0.154) within hybrids, were potentially involved in the adaptation to hypoxia. Thus, we provide a new profile of population structure from the CNV aspects of Chinese indigenous cattle and new insights into high-altitude adaptation in cattle.
Collapse
Affiliation(s)
- Yaran Zhang
- Dairy Cattle Research Center, Shandong Academy of Agricultural Sciences, Jinan, China
| | - Yan Hu
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education & College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, China
| | - Xiuge Wang
- Dairy Cattle Research Center, Shandong Academy of Agricultural Sciences, Jinan, China
| | - Qiang Jiang
- Dairy Cattle Research Center, Shandong Academy of Agricultural Sciences, Jinan, China
| | - Han Zhao
- Dairy Cattle Research Center, Shandong Academy of Agricultural Sciences, Jinan, China
| | - Jinpeng Wang
- Dairy Cattle Research Center, Shandong Academy of Agricultural Sciences, Jinan, China
| | - Zhihua Ju
- Dairy Cattle Research Center, Shandong Academy of Agricultural Sciences, Jinan, China
| | - Liguo Yang
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education & College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, China
| | - Yaping Gao
- Dairy Cattle Research Center, Shandong Academy of Agricultural Sciences, Jinan, China
| | - Xiaochao Wei
- Dairy Cattle Research Center, Shandong Academy of Agricultural Sciences, Jinan, China
| | - Jiachen Bai
- Dairy Cattle Research Center, Shandong Academy of Agricultural Sciences, Jinan, China
| | - Yang Zhou
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education & College of Animal Science and Technology, Huazhong Agricultural University, Wuhan, China
| | - Jinming Huang
- Dairy Cattle Research Center, Shandong Academy of Agricultural Sciences, Jinan, China.,Engineering Center of Animal Breeding and Reproduction, Jinan, China
| |
Collapse
|
35
|
Lee YL, Bosse M, Mullaart E, Groenen MAM, Veerkamp RF, Bouwman AC. Functional and population genetic features of copy number variations in two dairy cattle populations. BMC Genomics 2020; 21:89. [PMID: 31992181 PMCID: PMC6988284 DOI: 10.1186/s12864-020-6496-1] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2019] [Accepted: 01/14/2020] [Indexed: 12/27/2022] Open
Abstract
BACKGROUND Copy Number Variations (CNVs) are gain or loss of DNA segments that are known to play a role in shaping a wide range of phenotypes. In this study, we used two dairy cattle populations, Holstein Friesian and Jersey, to discover CNVs using the Illumina BovineHD Genotyping BeadChip aligned to the ARS-UCD1.2 assembly. The discovered CNVs were investigated for their functional impact and their population genetics features. RESULTS We discovered 14,272 autosomal CNVs, which were aggregated into 1755 CNV regions (CNVR) from 451 animals. These CNVRs together cover 2.8% of the bovine autosomes. The assessment of the functional impact of CNVRs showed that rare CNVRs (MAF < 0.01) are more likely to overlap with genes, than common CNVRs (MAF ≥ 0.05). The Population differentiation index (Fst) based on CNVRs revealed multiple highly diverged CNVRs between the two breeds. Some of these CNVRs overlapped with candidate genes such as MGAM and ADAMTS17 genes, which are related to starch digestion and body size, respectively. Lastly, linkage disequilibrium (LD) between CNVRs and BovineHD BeadChip SNPs was generally low, close to 0, although common deletions (MAF ≥ 0.05) showed slightly higher LD (r2 = ~ 0.1 at 10 kb distance) than the rest. Nevertheless, this LD is still lower than SNP-SNP LD (r2 = ~ 0.5 at 10 kb distance). CONCLUSIONS Our analyses showed that CNVRs detected using BovineHD BeadChip arrays are likely to be functional. This finding indicates that CNVs can potentially disrupt the function of genes and thus might alter phenotypes. Also, the population differentiation index revealed two candidate genes, MGAM and ADAMTS17, which hint at adaptive evolution between the two populations. Lastly, low CNVR-SNP LD implies that genetic variation from CNVs might not be fully captured in routine animal genetic evaluation, which relies solely on SNP markers.
Collapse
Affiliation(s)
- Young-Lim Lee
- Wageningen University & Research, Animal Breeding and Genomics, P.O. Box 338, Wageningen, AH, 6700, the Netherlands.
| | - Mirte Bosse
- Wageningen University & Research, Animal Breeding and Genomics, P.O. Box 338, Wageningen, AH, 6700, the Netherlands
| | | | - Martien A M Groenen
- Wageningen University & Research, Animal Breeding and Genomics, P.O. Box 338, Wageningen, AH, 6700, the Netherlands
| | - Roel F Veerkamp
- Wageningen University & Research, Animal Breeding and Genomics, P.O. Box 338, Wageningen, AH, 6700, the Netherlands
| | - Aniek C Bouwman
- Wageningen University & Research, Animal Breeding and Genomics, P.O. Box 338, Wageningen, AH, 6700, the Netherlands
| |
Collapse
|
36
|
|
37
|
Solé M, Ablondi M, Binzer-Panchal A, Velie BD, Hollfelder N, Buys N, Ducro BJ, François L, Janssens S, Schurink A, Viklund Å, Eriksson S, Isaksson A, Kultima H, Mikko S, Lindgren G. Inter- and intra-breed genome-wide copy number diversity in a large cohort of European equine breeds. BMC Genomics 2019; 20:759. [PMID: 31640551 PMCID: PMC6805398 DOI: 10.1186/s12864-019-6141-z] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2019] [Accepted: 09/25/2019] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Copy Number Variation (CNV) is a common form of genetic variation underlying animal evolution and phenotypic diversity across a wide range of species. In the mammalian genome, high frequency of CNV differentiation between breeds may be candidates for population-specific selection. However, CNV differentiation, selection and its population genetics have been poorly explored in horses. RESULTS We investigated the patterns, population variation and gene annotation of CNV using the Axiom® Equine Genotyping Array (670,796 SNPs) from a large cohort of individuals (N = 1755) belonging to eight European horse breeds, varying from draught horses to several warmblood populations. After quality control, 152,640 SNP CNVs (individual markers), 18,800 segment CNVs (consecutive SNP CNVs of same gain/loss state or both) and 939 CNV regions (CNVRs; overlapping segment CNVs by at least 1 bp) compared to the average signal of the reference (Belgian draught horse) were identified. Our analyses showed that Equus caballus chromosome 12 (ECA12) was the most enriched in segment CNV gains and losses (~ 3% average proportion of the genome covered), but the highest number of segment CNVs were detected on ECA1 and ECA20 (regardless of size). The Friesian horses showed private SNP CNV gains (> 20% of the samples) on ECA1 and Exmoor ponies displayed private SNP CNV losses on ECA25 (> 20% of the samples). The Warmblood cluster showed private SNP CNV gains located in ECA9 and Draught cluster showed private SNP CNV losses located in ECA7. The length of the CNVRs ranged from 1 kb to 21.3 Mb. A total of 10,612 genes were annotated within the CNVRs. The PANTHER annotation of these genes showed significantly under- and overrepresented gene ontology biological terms related to cellular processes and immunity (Bonferroni P-value < 0.05). We identified 80 CNVRs overlapping with known QTL for fertility, coat colour, conformation and temperament. We also report 67 novel CNVRs. CONCLUSIONS This work revealed that CNV patterns, in the genome of some European horse breeds, occurred in specific genomic regions. The results provide support to the hypothesis that high frequency private CNVs residing in genes may potentially be responsible for the diverse phenotypes seen between horse breeds.
Collapse
Affiliation(s)
- Marina Solé
- Department of Animal Breeding & Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden.
| | - Michela Ablondi
- Department of Veterinary Science, Università di Parma, Parma, Italy
| | - Amrei Binzer-Panchal
- Department of Medical Sciences, Array and Analysis Facility, Uppsala University, Uppsala, Sweden
| | - Brandon D Velie
- Faculty of Life and Environmental Science, University of Sydney, Sydney, NSW, Australia
| | - Nina Hollfelder
- Department of Medical Sciences, Array and Analysis Facility, Uppsala University, Uppsala, Sweden
| | - Nadine Buys
- Livestock Genetics, Department of Biosystems, KU Leuven, 3001, Leuven, Belgium
| | - Bart J Ducro
- Animal Breeding and Genomics, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, the Netherlands
| | - Liesbeth François
- Livestock Genetics, Department of Biosystems, KU Leuven, 3001, Leuven, Belgium
| | - Steven Janssens
- Livestock Genetics, Department of Biosystems, KU Leuven, 3001, Leuven, Belgium
| | - Anouk Schurink
- Animal Breeding and Genomics, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, the Netherlands.,Centre for Genetic Resources, the Netherlands (CGN), Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, the Netherlands
| | - Åsa Viklund
- Department of Animal Breeding & Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Susanne Eriksson
- Department of Animal Breeding & Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Anders Isaksson
- Department of Medical Sciences, Array and Analysis Facility, Uppsala University, Uppsala, Sweden
| | - Hanna Kultima
- Department of Medical Sciences, Array and Analysis Facility, Uppsala University, Uppsala, Sweden
| | - Sofia Mikko
- Department of Animal Breeding & Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - Gabriella Lindgren
- Department of Animal Breeding & Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden.,Livestock Genetics, Department of Biosystems, KU Leuven, 3001, Leuven, Belgium
| |
Collapse
|
38
|
Identification of Copy Number Variation in Domestic Chicken Using Whole-Genome Sequencing Reveals Evidence of Selection in the Genome. Animals (Basel) 2019; 9:ani9100809. [PMID: 31618984 PMCID: PMC6826909 DOI: 10.3390/ani9100809] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2019] [Revised: 10/13/2019] [Accepted: 10/14/2019] [Indexed: 12/12/2022] Open
Abstract
Simple Summary Chickens have been bred for meat and egg production as a source of animal protein. With the increase of productivity as the main purpose of domestication, factors such as metabolism and immunity were boosted, which are detectable signs of selection on the genome. This study focused on copy number variation (CNV) to find evidence of domestication on the genome. CNV was detected from whole-genome sequencing of 65 chickens including Red Jungle Fowl, broilers, and layers. After that, CNV region, the overlapping region of CNV between individuals, was made to identify which genomic regions showed copy number differentiation. The 663 domesticated-specific CNV regions were associated with various functions such as metabolism and organ development. Also, by performing population differentiation analyses such as clustering analysis and ANOVA test, we found that there are a lot of genomic regions with different copy number patterns between broilers and layers. This result indicates that different genetic variations can be found, depending on the purpose of artificial selection and provides considerations for future animal breeding. Abstract Copy number variation (CNV) has great significance both functionally and evolutionally. Various CNV studies are in progress to find the cause of human disease and to understand the population structure of livestock. Recent advances in next-generation sequencing (NGS) technology have made CNV detection more reliable and accurate at whole-genome level. However, there is a lack of CNV studies on chickens using NGS. Therefore, we obtained whole-genome sequencing data of 65 chickens including Red Jungle Fowl, Cornish (broiler), Rhode Island Red (hybrid), and White Leghorn (layer) from the public databases for CNV region (CNVR) detection. Using CNVnator, a read-depth based software, a total of 663 domesticated-specific CNVRs were identified across autosomes. Gene ontology analysis of genes annotated in CNVRs showed that mainly enriched terms involved in organ development, metabolism, and immune regulation. Population analysis revealed that CN and RIR are closer to each other than WL, and many genes (LOC772271, OR52R1, RD3, ADH6, TLR2B, PRSS2, TPK1, POPDC3, etc.) with different copy numbers between breeds found. In conclusion, this study has helped to understand the genetic characteristics of domestic chickens at CNV level, which may provide useful information for the development of breeding systems in chickens.
Collapse
|
39
|
Jia C, Wang H, Li C, Wu X, Zan L, Ding X, Guo X, Bao P, Pei J, Chu M, Liang C, Yan P. Genome-wide detection of copy number variations in polled yak using the Illumina BovineHD BeadChip. BMC Genomics 2019; 20:376. [PMID: 31088363 PMCID: PMC6518677 DOI: 10.1186/s12864-019-5759-1] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2019] [Accepted: 05/02/2019] [Indexed: 01/29/2023] Open
Abstract
Background Copy number variations (CNVs), which are genetic variations present throughout mammalian genomes, are a vital source of phenotypic variation that can lead to the development of unique traits. In this study we used the Illunima BovineHD BeadChip to conduct genome-wide detection of CNVs in 215 polled yaks. Results A total of 1066 CNV regions (CNVRs) were detected with a total length of 181.6 Mb, comprising ~ 7.2% of the bovine autosomal genome. The size of these CNVRs ranged from 5.53 kb to 1148.45 kb, with an average size of 170.31 kb. Eight out of nine randomly chosen CNVRs were successfully validated by qPCR. A functional enrichment analysis of the CNVR-associated genes indicated their relationship to a number of molecular adaptations that enable yaks to thrive at high altitudes. One third of the detected CNVRs were mapped to QTLs associated with six classes of economically important traits, indicating that these CNVRs may play an important role in variations of these traits. Conclusions Our genome-wide yak CNV map may thus provide valuable insights into both the molecular mechanisms of high altitude adaptation and the potential genomic basis of economically important traits in yak. Electronic supplementary material The online version of this article (10.1186/s12864-019-5759-1) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Congjun Jia
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China.,College of Animal Science and Technology, Northwest A&F University, Yangling, 712100, China
| | - Hongbo Wang
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China
| | - Chen Li
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China
| | - Xiaoyun Wu
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China
| | - Linsen Zan
- College of Animal Science and Technology, Northwest A&F University, Yangling, 712100, China
| | - Xuezhi Ding
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China
| | - Xian Guo
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China
| | - Pengjia Bao
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China
| | - Jie Pei
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China
| | - Min Chu
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China
| | - Chunnian Liang
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China.
| | - Ping Yan
- Key Laboratory of Yak Breeding Engineering Gansu Province, Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou, 730050, China.
| |
Collapse
|
40
|
Lye ZN, Purugganan MD. Copy Number Variation in Domestication. TRENDS IN PLANT SCIENCE 2019; 24:352-365. [PMID: 30745056 DOI: 10.1016/j.tplants.2019.01.003] [Citation(s) in RCA: 100] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/10/2018] [Revised: 01/08/2019] [Accepted: 01/10/2019] [Indexed: 05/22/2023]
Abstract
Domesticated plants have long served as excellent models for studying evolution. Many genes and mutations underlying important domestication traits have been identified, and most causal mutations appear to be SNPs. Copy number variation (CNV) is an important source of genetic variation that has been largely neglected in studies of domestication. Ongoing work demonstrates the importance of CNVs as a source of genetic variation during domestication, and during the diversification of domesticated taxa. Here, we review how CNVs contribute to evolutionary processes underlying domestication, and review examples of domestication traits caused by CNVs. We draw from examples in plant species, but also highlight cases in animal systems that could illuminate the roles of CNVs in the domestication process.
Collapse
Affiliation(s)
- Zoe N Lye
- Center for Genomics and Systems Biology, 12 Waverly Place, New York University, New York, NY 10003, USA
| | - Michael D Purugganan
- Center for Genomics and Systems Biology, 12 Waverly Place, New York University, New York, NY 10003, USA; Center for Genomics and Systems Biology, New York University Abu Dhabi, Saadiyat Island, Abu Dhabi, United Arab Emirates.
| |
Collapse
|
41
|
Liu M, Fang L, Liu S, Pan MG, Seroussi E, Cole JB, Ma L, Chen H, Liu GE. Array CGH-based detection of CNV regions and their potential association with reproduction and other economic traits in Holsteins. BMC Genomics 2019; 20:181. [PMID: 30845913 PMCID: PMC6407259 DOI: 10.1186/s12864-019-5552-1] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2018] [Accepted: 02/21/2019] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Copy number variations (CNVs) are structural variants consisting of large-scale insertions and deletions of genomic fragments. Exploring CNVs and estimating their effects on phenotypes are useful for genome selection but remain challenging in the livestock. RESULTS We identified 1043 CNV regions (CNVRs) from array comparative genomic hybridization (CGH) data of 47 Holstein bulls. Using a probe-based CNV association approach, we detected 87 CNVRs significantly (Bonferroni-corrected P value < 0.05) associated with at least one out of 41 complex traits. Within them, 39 CNVRs were simultaneously associated with at least 2 complex traits. Notably, 24 CNVRs were markedly related to daughter pregnancy rate (DPR). For example, CNVR661 containing CYP4A11 and CNVR213 containing CTR9, respectively, were associated with DPR and other traits related to reproduction, production, and body conformation. CNVR758 was also significantly related to DPR, with a nearby gene CAPZA3, encoding one of F-actin-capping proteins which play a role in determining sperm architecture and male fertility. We corroborated these CNVRs by examining their overlapped quantitative trait loci and comparing with previously published CNV results. CONCLUSION To our knowledge, this is one of the first genome-wide association studies based on CNVs called by array CGH in Holstein cattle. Our results contribute substantial information about the potential CNV impacts on reproduction, health, production, and body conformation traits, which lay the foundation for incorporating CNV into the future dairy cattle breeding program.
Collapse
Affiliation(s)
- Mei Liu
- College of Animal Science and Technology, Northwest A&F University, Shaanxi Key Laboratory of Agricultural Molecular Biology, Yangling, 712100 Shaanxi China
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
| | - Lingzhao Fang
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD USA
| | - Shuli Liu
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Michael G. Pan
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
| | - Eyal Seroussi
- Agricultural Research Organization (ARO), Volcani Center, Institute of Animal Science, Department of Quantitative and Molecular Genetics, HaMaccabim Road, P.O.B 15159, 7528809 Rishon LeTsiyon, Israel
| | - John B. Cole
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
| | - Li Ma
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD USA
| | - Hong Chen
- College of Animal Science and Technology, Northwest A&F University, Shaanxi Key Laboratory of Agricultural Molecular Biology, Yangling, 712100 Shaanxi China
| | - George E. Liu
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
| |
Collapse
|
42
|
Xu JW, Zheng L, Li LJ, Yao YF, Hua H, Yang SZ, Wen YF, Song CC, Cao XK, Liu KP, Zhang GM, Yang JM, Hao D, Dang RH, Lan XY, Lei CZ, Qi XL, Chen H, Huang YZ. Novel copy number variation of the KLF3 gene is associated with growth traits in beef cattle. Gene 2019; 680:99-104. [DOI: 10.1016/j.gene.2018.08.040] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2017] [Revised: 07/14/2018] [Accepted: 08/08/2018] [Indexed: 01/01/2023]
|
43
|
Zheng L, Xu JW, Li JC, Wang DH, An QM, Xu LN, Ma YL, Wang J, Peng SJ, Lei CZ, Lan XY, Chen H, Huo LJ, Huang YZ. Distribution and association study in copy number variation of KCNJ12 gene across four Chinese cattle populations. Gene 2018; 689:90-96. [PMID: 30572095 DOI: 10.1016/j.gene.2018.12.019] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2018] [Revised: 11/29/2018] [Accepted: 12/06/2018] [Indexed: 10/27/2022]
Abstract
Copy number variation is a large genome variation which usually happens in the noncoding-region, and it may occur at the locus associated with the functional gene to further influence the phenotype. Potassium inwardly-rectifying channel, subfamily J 12 (KCNJ12) gene expressed widely in cardiomyocytes and neurons, plays an important role in tumor therapy and muscle movement regulation. In this study, we detected the distribution of CNVs for KCNJ12 gene in 404 individuals belonging to four Chinese cattle breeds (NY, JX, JA and GF). We also investigated the KCNJ12 gene expression in different tissues of JX cattle. Additionally, we examined the association of two CNV regions (CNV1: 1,600 bp, intron 1; CNV2: 4,800 bp, intergenic) with growth traits. The statistical analyses indicated that the CNV1 is associated with the body length, rump length and weight in JX cattle population (P < 0.05); and there has a significant association with the body length, chest circumference, and body weight in GF cattle (P < 0.05).The CNV2 had a significant effect on the body length and body weight in JX cattle (P < 0.05); the body length, chest circumference, rump length and body weight in GF cattle (P < 0.01 or P < 0.05). The copy numbers of KCNJ12 gene presented the negative correlations with the transcript level of gene in skeletal muscles (P < 0.05). Our results provide evidence that CNV1 and CNV 2 in KCNJ12 are associated with growth traits in two cattle populations and may be used as candidates for marker-assisted selection and breeding management in cattle.
Collapse
Affiliation(s)
- Li Zheng
- Henan University of Animal Husbandry and Economy, Zhengzhou, Henan 450046, People's Republic of China; Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction, Education Ministry of China, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan 430070, People's Republic of China
| | - Jia-Wei Xu
- College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi 712100, People's Republic of China
| | - Ji-Chao Li
- Henan University of Animal Husbandry and Economy, Zhengzhou, Henan 450046, People's Republic of China
| | - Da-Hui Wang
- College of Agriculture and Forestry Engineering, Tongren Unviersity, Tongren, Guizhou 554300, People's Republic of China
| | - Qing-Ming An
- College of Agriculture and Forestry Engineering, Tongren Unviersity, Tongren, Guizhou 554300, People's Republic of China
| | - Lin-Na Xu
- Gansu Animal Husbandry Industry Administration, Lanzhou, Gansu 730000, People's Republic of China
| | - Yi-Lei Ma
- College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi 712100, People's Republic of China
| | - Jian Wang
- College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi 712100, People's Republic of China
| | - Shu-Jun Peng
- College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi 712100, People's Republic of China
| | - Chu-Zhao Lei
- College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi 712100, People's Republic of China
| | - Xian-Yong Lan
- College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi 712100, People's Republic of China
| | - Hong Chen
- College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi 712100, People's Republic of China
| | - Li-Jun Huo
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction, Education Ministry of China, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan 430070, People's Republic of China.
| | - Yong-Zhen Huang
- College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi 712100, People's Republic of China.
| |
Collapse
|
44
|
Copy Number Variations of KLF6 Modulate Gene Transcription and Growth Traits in Chinese Datong Yak (Bos Grunniens). Animals (Basel) 2018; 8:ani8090145. [PMID: 30134528 PMCID: PMC6162419 DOI: 10.3390/ani8090145] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2018] [Revised: 08/06/2018] [Accepted: 08/18/2018] [Indexed: 02/08/2023] Open
Abstract
Copy number variation (CNV) is a significant marker of the genetic and phenotypic diversity among individuals that accounts for complex quantitative traits of phenotype and diseases via modulating gene dosage and disrupting coding regions in the genome. Biochemically, Kruppel-like factor 6 (KLF6) genes plays a significant role in the regulation of cell differentiation and proliferation and muscle development. The aim of this study was to detect the distributions of KLF6 copy number variations (CNVs) in five breeds of domestic yak and to explore their effect on growth traits and gene expression. The data were analyzed by real-time quantitative PCR (qPCR). Our results elucidated that a decreased CNV in the KLF6 gene is more highly associated (p < 0.05) with various growth traits than increased or normal CNVs in six-month-old and five-year-old Datong yak. Nevertheless, negative correlations between the DNA copy number and KLF6 gene expression were observed in the skeletal muscle of adult Datong yak. These results suggest that CNVs of the KLF6 gene could be crucial genomic markers for growth phenotypes of Chinese Datong yak breeds and this finding constitutes the first evidence of the biological role of KLF6 CNVs in Chinese Datong yak breeds.
Collapse
|