1
|
Wang Y, Ma J, Wang J, Zhang L, Xu L, Chen Y, Zhu B, Wang Z, Gao H, Li J, Gao X. Genome-Wide Detection of Copy Number Variations and Their Potential Association with Carcass and Meat Quality Traits in Pingliang Red Cattle. Int J Mol Sci 2024; 25:5626. [PMID: 38891814 PMCID: PMC11172001 DOI: 10.3390/ijms25115626] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2024] [Revised: 05/16/2024] [Accepted: 05/17/2024] [Indexed: 06/21/2024] Open
Abstract
Copy number variation (CNV) serves as a significant source of genetic diversity in mammals and exerts substantial effects on various complex traits. Pingliang red cattle, an outstanding indigenous resource in China, possess remarkable breeding value attributed to their tender meat and superior marbling quality. However, the genetic mechanisms influencing carcass and meat quality traits in Pingliang red cattle are not well understood. We generated a comprehensive genome-wide CNV map for Pingliang red cattle using the GGP Bovine 100K SNP chip. A total of 755 copy number variable regions (CNVRs) spanning 81.03 Mb were identified, accounting for approximately 3.24% of the bovine autosomal genome. Among these, we discovered 270 potentially breed-specific CNVRs in Pingliang red cattle, including 143 gains, 73 losses, and 54 mixed events. Functional annotation analysis revealed significant associations between these specific CNVRs and important traits such as carcass and meat quality, reproduction, exterior traits, growth traits, and health traits. Additionally, our network and transcriptome analysis highlighted CACNA2D1, CYLD, UBXN2B, TG, NADK, and ITGA9 as promising candidate genes associated with carcass weight and intramuscular fat deposition. The current study presents a genome-wide CNV map in Pingliang red cattle, highlighting breed-specific CNVRs, and transcriptome findings provide valuable insights into the underlying genetic characteristics of Pingliang red cattle. These results offer potential avenues for enhancing meat quality through a targeted breeding program.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | | | | | - Junya Li
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China
| | - Xue Gao
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China
| |
Collapse
|
2
|
Zhu H, Lu X, Jiang H, Yang Z, Xu T. Descriptive Statistics and Genome-Wide Copy Number Analysis of Milk Production Traits of Jiangsu Chinese Holstein Cows. Animals (Basel) 2023; 14:17. [PMID: 38200748 PMCID: PMC10778490 DOI: 10.3390/ani14010017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2023] [Revised: 12/05/2023] [Accepted: 12/18/2023] [Indexed: 01/12/2024] Open
Abstract
Milk production traits are the most important quantitative economic traits in dairy cow production; improving the yield and quality of milk is an important way to ensure the production efficiency of the dairy industry. This study carried out a series of in-depth statistical genetics studies and molecular analyses on the Chinese Holstein cows in the Jiangsu Province, such as descriptive statistics and copy number variation analysis. A genetic correlation, phenotypic correlation, and descriptive statistical analysis of five milk production traits (milk yield, milk fat percentage, milk fat yield, milk protein percentage, and milk protein yield) of the dairy cows were analyzed using the SPSS and DMU software. Through quality control, 4173 cows and their genomes were used for genomic study. Then, SNPs were detected using DNA chips, and a copy number variation (CNV) analysis was carried out to locate the quantitative trait loci (QTL) of the milk production traits by Perl program software Penn CNV and hidden Markov model (HMM). The phenotypic means of the milk yield, milk fat percentage, milk fat mass, milk protein percentage, and milk protein mass at the first trimester were lower than those at the other trimesters by 8.821%, 1.031%, 0.930%, 0.003%, and 0.826%, respectively. The five milk production traits showed a significant phenotypic positive correlation (p < 0.01) and a high genetic positive correlation among the three parities. Based on the GGPBovine 100 K SNP data, QTL-detecting research on the fist-parity milk performance of dairy cows was carried out via the CNV. We identified 1731 CNVs and 236 CNVRs in the 29 autosomes of 984 Holstein dairy cows, and 19 CNVRs were significantly associated with the milk production traits (p < 0.05). These CNVRs were analyzed via a bioinformatics analysis; a total of 13 gene ontology (GO) terms and 20 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were significantly enriched (p < 0.05), and these terms and pathways are mainly related to lipid metabolism, amino acid metabolism, and cellular catabolic processes. This study provided a theoretical basis for the molecular-marker-assisted selection of dairy cows by developing descriptive statistics on the milk production traits of dairy cows and by locating the QTL and functional genes that affect the milk production traits of first-born dairy cows. The results describe the basic status of the milk production traits of the Chinese Holstein cows in Jiangsu and locate the QTL and functional genes that affect the milk production traits of the first-born cows, providing a theoretical basis for the molecular-marker-assisted selection of dairy cows.
Collapse
Affiliation(s)
- Hao Zhu
- Joint International Research Laboratory of Agriculture and Agri-Product Safety, Yangzhou University, Yangzhou 225009, China; (H.Z.); (Z.Y.)
- College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China;
| | - Xubin Lu
- College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China;
| | - Hui Jiang
- Center for Quantitative Genetics and Genomics, Aarhus University, 8000 Aarhus C, Denmark;
| | - Zhangping Yang
- Joint International Research Laboratory of Agriculture and Agri-Product Safety, Yangzhou University, Yangzhou 225009, China; (H.Z.); (Z.Y.)
- College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China;
| | - Tianle Xu
- Joint International Research Laboratory of Agriculture and Agri-Product Safety, Yangzhou University, Yangzhou 225009, China; (H.Z.); (Z.Y.)
- College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China;
- International Joint Research Laboratory, Universities of Jiangsu Province of China for Domestic Animal Germplasm Resources and Genetic Improvement, Yangzhou 225009, China
| |
Collapse
|
3
|
Braga LG, Chud TCS, Watanabe RN, Savegnago RP, Sena TM, do Carmo AS, Machado MA, Panetto JCDC, da Silva MVGB, Munari DP. Identification of copy number variations in the genome of Dairy Gir cattle. PLoS One 2023; 18:e0284085. [PMID: 37036840 PMCID: PMC10085049 DOI: 10.1371/journal.pone.0284085] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2022] [Accepted: 03/23/2023] [Indexed: 04/11/2023] Open
Abstract
Studying structural variants that can control complex traits is relevant for dairy cattle production, especially for animals that are tolerant to breeding conditions in the tropics, such as the Dairy Gir cattle. This study identified and characterized high confidence copy number variation regions (CNVR) in the Gir breed genome. A total of 38 animals were whole-genome sequenced, and 566 individuals were genotyped with a high-density SNP panel, among which 36 animals had both sequencing and SNP genotyping data available. Two sets of high confidence CNVR were established: one based on common CNV identified in the studied population (CNVR_POP), and another with CNV identified in sires with both sequence and SNP genotyping data available (CNVR_ANI). We found 10 CNVR_POP and 45 CNVR_ANI, which covered 1.05 Mb and 4.4 Mb of the bovine genome, respectively. Merging these CNV sets for functional analysis resulted in 48 unique high confidence CNVR. The overlapping genes were previously related to embryonic mortality, environmental adaptation, evolutionary process, immune response, longevity, mammary gland, resistance to gastrointestinal parasites, and stimuli recognition, among others. Our results contribute to a better understanding of the Gir breed genome. Moreover, the CNV identified in this study can potentially affect genes related to complex traits, such as production, health, and reproduction.
Collapse
Affiliation(s)
- Larissa G Braga
- Departamento de Engenharia e Ciências Exatas, Universidade Estadual Paulista, Jaboticabal, São Paulo, Brazil
| | - Tatiane C S Chud
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, Ontario, Canada
| | - Rafael N Watanabe
- Departamento de Engenharia e Ciências Exatas, Universidade Estadual Paulista, Jaboticabal, São Paulo, Brazil
| | - Rodrigo P Savegnago
- Department of Animal Science, Michigan State University, East Lansing, Michigan, United States of America
| | - Thomaz M Sena
- Departamento de Engenharia e Ciências Exatas, Universidade Estadual Paulista, Jaboticabal, São Paulo, Brazil
| | - Adriana S do Carmo
- Departamento de Zootecnia, Universidade Federal de Goiás, Goiânia, Goiás, Brazil
| | | | | | | | - Danísio P Munari
- Departamento de Engenharia e Ciências Exatas, Universidade Estadual Paulista, Jaboticabal, São Paulo, Brazil
| |
Collapse
|
4
|
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
|
5
|
Cumer T, Boyer F, Pompanon F. Genome-Wide Detection of Structural Variations Reveals New Regions Associated with Domestication in Small Ruminants. Genome Biol Evol 2021; 13:evab165. [PMID: 34264322 PMCID: PMC8350358 DOI: 10.1093/gbe/evab165] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/08/2021] [Indexed: 11/28/2022] Open
Abstract
During domestication processes, changes in selective pressures induce multiple phenotypical, physiological, and behavioral changes in target species. The rise of next-generation sequencing has provided a chance to study the genetics bases of these changes, most of the time based on single nucleotide polymorphisms (SNPs). However, several studies have highlighted the impact of structural variations (SVs) on individual fitness, particularly in domestic species. We aimed at unraveling the role of SVs during the domestication and later improvement of small ruminants by analyzing whole-genome sequences of 40 domestic sheep and 11 of their close wild relatives (Ovis orientalis), and 40 goats and 18 of their close wild relatives (Capra aegagrus). Using a combination of detection tools, we called 45,796 SVs in Ovis and 15,047 SVs in Capra genomes, including insertions, deletions, inversions, copy number variations, and chromosomal translocations. Most of these SVs were previously unreported in small ruminants. 69 and 45 SVs in sheep and goats, respectively, were in genomic regions with neighboring SNPs highly differentiated between wilds and domestics (i.e., putatively related to domestication). Among them, 25 and 20 SVs were close to or overlapping with genes related to physiological and morpho-anatomical traits linked with productivity (e.g., size, meat or milk quality, wool color), reproduction, or immunity. Finally, several of the SVs differentiated between wilds and domestics would not have been detected by screening only the differentiation of SNPs surrounding them, highlighting the complementarity of SVs and SNPs based approaches to detect signatures of selection.
Collapse
Affiliation(s)
- Tristan Cumer
- Université Grenoble Alpes, Université Savoie Mont-Blanc, CNRS, LECA, Grenoble, France
| | - Frédéric Boyer
- Université Grenoble Alpes, Université Savoie Mont-Blanc, CNRS, LECA, Grenoble, France
| | - François Pompanon
- Université Grenoble Alpes, Université Savoie Mont-Blanc, CNRS, LECA, Grenoble, France
| |
Collapse
|
6
|
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
|
7
|
Chen L, Pryce JE, Hayes BJ, Daetwyler HD. Investigating the Effect of Imputed Structural Variants from Whole-Genome Sequence on Genome-Wide Association and Genomic Prediction in Dairy Cattle. Animals (Basel) 2021; 11:ani11020541. [PMID: 33669735 PMCID: PMC7922624 DOI: 10.3390/ani11020541] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2021] [Revised: 02/09/2021] [Accepted: 02/12/2021] [Indexed: 02/06/2023] Open
Abstract
Simple Summary Structural variants are large changes to the DNA sequences that differ from individual to individual. We discovered and quality-controlled a set of 24,908 structural variants and used a technique called imputation to infer them into 35,588 Holstein and Jersey cattle. We then investigated whether the structural variants affected key dairy cattle traits such as milk production, fertility and overall conformation. Structural variants explained generally less than 10 percent of the phenotypic variation in these traits. Four of the structural variants were significantly associated with dairy cattle production traits. However, the inclusion of the structural variants in the genomic prediction model did not increase genomic prediction accuracy. Abstract Structural variations (SVs) are large DNA segments of deletions, duplications, copy number variations, inversions and translocations in a re-sequenced genome compared to a reference genome. They have been found to be associated with several complex traits in dairy cattle and could potentially help to improve genomic prediction accuracy of dairy traits. Imputation of SVs was performed in individuals genotyped with single-nucleotide polymorphism (SNP) panels without the expense of sequencing them. In this study, we generated 24,908 high-quality SVs in a total of 478 whole-genome sequenced Holstein and Jersey cattle. We imputed 4489 SVs with R2 > 0.5 into 35,568 Holstein and Jersey dairy cattle with 578,999 SNPs with two pipelines, FImpute and Eagle2.3-Minimac3. Genome-wide association studies for production, fertility and overall type with these 4489 SVs revealed four significant SVs, of which two were highly linked to significant SNP. We also estimated the variance components for SNP and SV models for these traits using genomic best linear unbiased prediction (GBLUP). Furthermore, we assessed the effect on genomic prediction accuracy of adding SVs to GBLUP models. The estimated percentage of genetic variance captured by SVs for production traits was up to 4.57% for milk yield in bulls and 3.53% for protein yield in cows. Finally, no consistent increase in genomic prediction accuracy was observed when including SVs in GBLUP.
Collapse
Affiliation(s)
- Long Chen
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC 3083, Australia; (L.C.); (J.E.P.); (B.J.H.)
- School of Applied Systems Biology, La Trobe University, Bundoora, VIC 3083, Australia
| | - Jennie E. Pryce
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC 3083, Australia; (L.C.); (J.E.P.); (B.J.H.)
- School of Applied Systems Biology, La Trobe University, Bundoora, VIC 3083, Australia
| | - Ben J. Hayes
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC 3083, Australia; (L.C.); (J.E.P.); (B.J.H.)
- Queensland Alliance for Agriculture and Food Innovation, Centre for Animal Science, The University of Queensland, St. Lucia, QLD 4067, Australia
| | - Hans D. Daetwyler
- Agriculture Victoria, AgriBio, Centre for AgriBioscience, Bundoora, VIC 3083, Australia; (L.C.); (J.E.P.); (B.J.H.)
- School of Applied Systems Biology, La Trobe University, Bundoora, VIC 3083, Australia
- Correspondence:
| |
Collapse
|
8
|
Berton MP, de Antunes Lemos MV, Seleguim Chud TC, Bonvino Stafuzza N, Kluska S, Amorim ST, Silva Ferlin Lopes L, Cravo Pereira AS, Bickhart D, Liu G, Galvão de Albuquerque L, Baldi F. Genome-wide association study between copy number variation regions and carcass- and meat-quality traits in Nellore cattle. ANIMAL PRODUCTION SCIENCE 2021. [DOI: 10.1071/an20275] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Context
Indicine breeds are the main source of beef products in tropical and subtropical regions. However, genetic improvement for carcass- and meat-quality traits in zebu cattle have been limited and genomics studies concerning structural variations that influence these traits are essential.
Aim
The aim of this study was to perform a genome-wide association study between copy number variation regions (CNVRs) and carcass- and meat quality-traits in Nellore cattle.
Methods
In total, 3794 animals, males and females included, were genotyped using a 777962 single-nucleotide polymorphism platform of BovineHD BeadChip (777k; Illumina Inc.). Of these, 1751 Nellore bulls were slaughtered at 24 months of age for further carcass beef analysis. The following traits were studied: beef tenderness, marbling, rib-eye area, backfat thickness and meat colour (lightness, redness and yellowness). The CNV detection was conducted through PennCNV software. The association analyses were performed using CNVRuler software.
Key results
Several identified genomic regions were linked to quantitative trait loci associated with fat deposition (FABP7) and lipid metabolism (PPARA; PLA2 family; BCHE), extracellular matrix (INS; COL10A1), contraction (SLC34A3; TRDN) and muscle development (CAPZP). The gene-enrichment analyses highlighted biological mechanisms directly related to the metabolism and synthesis of lipids and fatty acids.
Conclusions
The large number of potential candidate genes identified within the CNVRs, as well as the functions and pathways identified, should help better elucidate the genetic mechanisms involved in the expression of beef and carcass traits in Nellore cattle. Several CNVRs harboured genes that might have a functional impact to improve the beef and carcass traits.
Implications
The results obtained contribute to upgrade the sensorial and organoleptic attributes of Nellore cattle and make feasible the genetic improvement of carcass- and meat-quality traits.
Collapse
|
9
|
Li J, Fan Z, Shen F, Pendleton AL, Song Y, Xing J, Yue B, Kidd JM, Li J. Genomic Copy Number Variation Study of Nine Macaca Species Provides New Insights into Their Genetic Divergence, Adaptation, and Biomedical Application. Genome Biol Evol 2020; 12:2211-2230. [PMID: 32970804 PMCID: PMC7846157 DOI: 10.1093/gbe/evaa200] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/19/2020] [Indexed: 02/06/2023] Open
Abstract
Copy number variation (CNV) can promote phenotypic diversification and adaptive evolution. However, the genomic architecture of CNVs among Macaca species remains scarcely reported, and the roles of CNVs in adaptation and evolution of macaques have not been well addressed. Here, we identified and characterized 1,479 genome-wide hetero-specific CNVs across nine Macaca species with bioinformatic methods, along with 26 CNV-dense regions and dozens of lineage-specific CNVs. The genes intersecting CNVs were overrepresented in nutritional metabolism, xenobiotics/drug metabolism, and immune-related pathways. Population-level transcriptome data showed that nearly 46% of CNV genes were differentially expressed across populations and also mainly consisted of metabolic and immune-related genes, which implied the role of CNVs in environmental adaptation of Macaca. Several CNVs overlapping drug metabolism genes were verified with genomic quantitative polymerase chain reaction, suggesting that these macaques may have different drug metabolism features. The CNV-dense regions, including 15 first reported here, represent unstable genomic segments in macaques where biological innovation may evolve. Twelve gains and 40 losses specific to the Barbary macaque contain genes with essential roles in energy homeostasis and immunity defense, inferring the genetic basis of its unique distribution in North Africa. Our study not only elucidated the genetic diversity across Macaca species from the perspective of structural variation but also provided suggestive evidence for the role of CNVs in adaptation and genome evolution. Additionally, our findings provide new insights into the application of diverse macaques to drug study.
Collapse
Affiliation(s)
- Jing Li
- Key Laboratory of Bio-Resources and Eco-Environment (Ministry of Education), College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
- Institute of Animal Genetics and Breeding, College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, Sichuan, China
| | - Zhenxin Fan
- Key Laboratory of Bio-Resources and Eco-Environment (Ministry of Education), College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
- Sichuan Key Laboratory of Conservation Biology on Endangered Wildlife, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Feichen Shen
- Department of Human Genetics, Medical School, University of Michigan
| | | | - Yang Song
- Key Laboratory of Bio-Resources and Eco-Environment (Ministry of Education), College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Jinchuan Xing
- Department of Genetics and the Human Genetics Institute of New Jersey, Rutgers, The State University of New Jersey, Piscataway
| | - Bisong Yue
- Key Laboratory of Bio-Resources and Eco-Environment (Ministry of Education), College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| | - Jeffrey M Kidd
- Department of Human Genetics, Medical School, University of Michigan
| | - Jing Li
- Key Laboratory of Bio-Resources and Eco-Environment (Ministry of Education), College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
- Sichuan Key Laboratory of Conservation Biology on Endangered Wildlife, College of Life Sciences, Sichuan University, Chengdu, Sichuan, China
| |
Collapse
|
10
|
Biscarini F, Mastrangelo S, Catillo G, Senczuk G, Ciampolini R. Insights into Genetic Diversity, Runs of Homozygosity and Heterozygosity-Rich Regions in Maremmana Semi-Feral Cattle Using Pedigree and Genomic Data. Animals (Basel) 2020; 10:E2285. [PMID: 33287320 PMCID: PMC7761732 DOI: 10.3390/ani10122285] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2020] [Revised: 11/27/2020] [Accepted: 11/30/2020] [Indexed: 12/24/2022] Open
Abstract
Semi-feral local livestock populations, like Maremmana cattle, are the object of renewed interest for the conservation of biological diversity and the preservation and exploitation of unique and potentially relevant genetic material. The aim of this study was to estimate genetic diversity parameters in semi-feral Maremmana cattle using both pedigree- and genomic-based approaches (FIS and FROH), and to detect regions of homozygosity (ROH) and heterozygosity (ROHet) in the genome. The average heterozygosity estimates were in the range reported for other cattle breeds (HE=0.261, HO=0.274). Pedigree-based average inbreeding (F) was estimated at 4.9%. The correlation was low between F and genomic-based approaches (r=0.03 with FIS, r=0.21 with FROH), while it was higher between FIS and FROH (r=0.78). The low correlation between F and FROH coefficients may be the result of the limited pedigree depth available for the animals involved in this study. The ROH islands identified in Maremmana cattle included candidate genes associated with climate adaptation, carcass traits or the regulation of body weight, fat and energy metabolism. The ROHet islands contained candidate genes associated with nematode resistance and reproduction traits in livestock. The results of this study confirm that genome-based measures like FROH may be useful estimators of individual autozygosity, and may provide insights on pedigree-based inbreeding estimates in cases when animals' pedigree data are unavailable, thus providing a more detailed picture of the genetic diversity.
Collapse
Affiliation(s)
- Filippo Biscarini
- CNR-IBBA (National Research Council, Institute of Agricultural Biology and Biotechnology), 20133 Milan, Italy
| | - Salvatore Mastrangelo
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, 90128 Palermo, Italy;
| | - Gennaro Catillo
- CREA Research Centre for Animal Production and Acquaculture, CREA, 00015 Monterotondo, Italy;
| | - Gabriele Senczuk
- Dipartimento di Agricoltura, Ambiente e Alimenti, University of Molise, 86100 Campobasso, Italy;
| | - Roberta Ciampolini
- Dipartimento di Scienze Veterinarie—Università di Pisa, 56124 Pisa, Italy;
| |
Collapse
|
11
|
Yaman Y, Aymaz R, Keleş M, Bay V, Özüiçli M, Şenlik B. Association between ovine Toll-like receptor 4 ( TLR4) gene coding variants and presence of Eimeria spp. in naturally infected adult Turkish native sheep. Anim Biotechnol 2020; 32:375-380. [PMID: 32967565 DOI: 10.1080/10495398.2020.1824921] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Abstract
Coccidiosis caused by Eimeria spp. is a protozoan disease prevalent in farm animals, and it is responsible for serious economic losses especially in young animals. It has been popular to breed disease-resistant animals due to the concern about food safety, animal welfare, and public health. Toll-like receptor (TLR) gene family plays a key role in the innate immune system participating in host-antigen interaction, therefore, they are candidate genes for breeding disease-resistant animals. In the present study, possible genetic associations between TLR4 gene coding variants and the presence of Eimeria spp. in adult Turkish sheep were investigated. For this purpose, the presence of Eimeria spp. in fecal samples from six native Turkish sheep were determined, and approximately 1450 bp region in the 3rd exon of the ovine TLR4 gene was sequenced. Ten nonsynonymous and four synonymous single nucleotide polymorphisms (SNPs) were detected in the targeted region. Statistical analyses revealed that the SNP at the codon at 356th position encoding Leucine instead of Phenylalanine (F356L) was significantly associated with the presence of Eimeria spp. It was found that the individuals carrying at least one Leucine amino acid sequence at this position have 2.3-fold more risk for the presence of Eimeria spp.
Collapse
Affiliation(s)
- Yalçın Yaman
- Department of Biometry and Genetics, Sheep Breeding and Research Institute, Bandirma, Turkey
| | - Ramazan Aymaz
- Department of Biometry and Genetics, Sheep Breeding and Research Institute, Bandirma, Turkey
| | - Murat Keleş
- Department of Biometry and Genetics, Sheep Breeding and Research Institute, Bandirma, Turkey
| | - Veysel Bay
- Department of Biometry and Genetics, Sheep Breeding and Research Institute, Bandirma, Turkey
| | - Mehmet Özüiçli
- Veterinary Faculty, Department of Parasitology, Uludag University, Bursa, Turkey
| | - Bayram Şenlik
- Veterinary Faculty, Department of Parasitology, Uludag University, Bursa, Turkey
| |
Collapse
|
12
|
Butty AM, Chud TCS, Miglior F, Schenkel FS, Kommadath A, Krivushin K, Grant JR, Häfliger IM, Drögemüller C, Cánovas A, Stothard P, Baes CF. High confidence copy number variants identified in Holstein dairy cattle from whole genome sequence and genotype array data. Sci Rep 2020; 10:8044. [PMID: 32415111 PMCID: PMC7229195 DOI: 10.1038/s41598-020-64680-3] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2019] [Accepted: 04/15/2020] [Indexed: 12/15/2022] Open
Abstract
Multiple methods to detect copy number variants (CNV) relying on different types of data have been developed and CNV have been shown to have an impact on phenotypes of numerous traits of economic importance in cattle, such as reproduction and immunity. Further improvements in CNV detection are still needed in regard to the trade-off between high-true and low-false positive variant identification rates. Instead of improving single CNV detection methods, variants can be identified in silico with high confidence when multiple methods and datasets are combined. Here, CNV were identified from whole-genome sequences (WGS) and genotype array (GEN) data on 96 Holstein animals. After CNV detection, two sets of high confidence CNV regions (CNVR) were created that contained variants found in both WGS and GEN data following an animal-based (n = 52) and a population-based (n = 36) pipeline. Furthermore, the change in false positive CNV identification rates using different GEN marker densities was evaluated. The population-based approach characterized CNVR, which were more often shared among animals (average 40% more samples per CNVR) and were more often linked to putative functions (48 vs 56% of CNVR) than CNV identified with the animal-based approach. Moreover, false positive identification rates up to 22% were estimated on GEN information. Further research using larger datasets should use a population-wide approach to identify high confidence CNVR.
Collapse
Affiliation(s)
- Adrien M Butty
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | - Tatiane C S Chud
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | - Filippo Miglior
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | - Flavio S Schenkel
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | - Arun Kommadath
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada.,Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB, Canada
| | - Kirill Krivushin
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Jason R Grant
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Irene M Häfliger
- Institute of Genetics, Vetsuisse Faculty, University of Bern, Bern, BE, Switzerland
| | - Cord Drögemüller
- Institute of Genetics, Vetsuisse Faculty, University of Bern, Bern, BE, Switzerland
| | - Angela Cánovas
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada
| | - Paul Stothard
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Christine F Baes
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, ON, Canada. .,Institute of Genetics, Vetsuisse Faculty, University of Bern, Bern, BE, Switzerland.
| |
Collapse
|
13
|
Kommadath A, Grant JR, Krivushin K, Butty AM, Baes CF, Carthy TR, Berry DP, Stothard P. A large interactive visual database of copy number variants discovered in taurine cattle. Gigascience 2020; 8:5523204. [PMID: 31241156 PMCID: PMC6593363 DOI: 10.1093/gigascience/giz073] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2018] [Revised: 02/27/2019] [Accepted: 05/28/2019] [Indexed: 12/19/2022] Open
Abstract
BACKGROUND Copy number variants (CNVs) contribute to genetic diversity and phenotypic variation. We aimed to discover CNVs in taurine cattle using a large collection of whole-genome sequences and to provide an interactive database of the identified CNV regions (CNVRs) that includes visualizations of sequence read alignments, CNV boundaries, and genome annotations. RESULTS CNVs were identified in each of 4 whole-genome sequencing datasets, which together represent >500 bulls from 17 breeds, using a popular multi-sample read-depth-based algorithm, cn.MOPS. Quality control and CNVR construction, performed dataset-wise to avoid batch effects, resulted in 26,223 CNVRs covering 107.75 unique Mb (4.05%) of the bovine genome. Hierarchical clustering of samples by CNVR genotypes indicated clear separation by breeds. An interactive HTML database was created that allows data filtering options, provides graphical and tabular data summaries including Hardy-Weinberg equilibrium tests on genotype proportions, and displays genes and quantitative trait loci at each CNVR. Notably, the database provides sequence read alignments at each CNVR genotype and the boundaries of constituent CNVs in individual samples. Besides numerous novel discoveries, we corroborated the genotypes reported for a CNVR at the KIT locus known to be associated with the piebald coat colour phenotype in Hereford and some Simmental cattle. CONCLUSIONS We present a large comprehensive collection of taurine cattle CNVs in a novel interactive visual database that displays CNV boundaries, read depths, and genome features for individual CNVRs, thus providing users with a powerful means to explore and scrutinize CNVRs of interest more thoroughly.
Collapse
Affiliation(s)
- Arun Kommadath
- Department of Agricultural, Food and Nutritional Science (AFNS), University of Alberta, Edmonton, AB, Canada.,Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, Alberta, Canada
| | - Jason R Grant
- Department of Agricultural, Food and Nutritional Science (AFNS), University of Alberta, Edmonton, AB, Canada
| | - Kirill Krivushin
- Department of Agricultural, Food and Nutritional Science (AFNS), University of Alberta, Edmonton, AB, Canada
| | - Adrien M Butty
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada
| | - Christine F Baes
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada.,Institute of Genetics, Vetsuisse Faculty, University of Bern, Bern, Switzerland
| | - Tara R Carthy
- Teagasc, Animal & Grassland Research and Innovation Centre, Moorepark, Fermoy, Ireland
| | - Donagh P Berry
- Teagasc, Animal & Grassland Research and Innovation Centre, Moorepark, Fermoy, Ireland
| | - Paul Stothard
- Department of Agricultural, Food and Nutritional Science (AFNS), University of Alberta, Edmonton, AB, Canada
| |
Collapse
|
14
|
Sun T, Hanif Q, Chen H, Lei C, Dang R. Copy Number Variations of Four Y-Linked Genes in Swamp Buffaloes. Animals (Basel) 2019; 10:ani10010031. [PMID: 31877875 PMCID: PMC7023270 DOI: 10.3390/ani10010031] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2019] [Revised: 12/12/2019] [Accepted: 12/18/2019] [Indexed: 12/25/2022] Open
Abstract
Simple Summary The amplification of the male-specific region of the Y chromosome was a unique phenomenon during mammalian sex chromosome evolution. The Y-linked copy number variations of many species have been confirmed. However, the Y-linked copy number variations (CNVs) in water buffalo are still unknown. In this study, we investigated the copy number variations of four Y-linked genes (SRY, UTY, DBY, and OFD1Y) in buffalo. Our results showed that UTY was a single-copy gene in buffalo, while DBY, OFD1Y, and SRY exhibited copy number variations in buffalo. Abstract Copy number variation (CNV), a significant source of genetic diversity in the mammalian Y chromosome, is associated with the development of many complex phenotypes, such as spermatogenesis and male fertility. The contribution of Y-linked CNVs has been studied in various species, however, water buffalo has not been explored in this area and the genetic information still remains unknown. The aim of the current study was to investigate the CNVs of four Y-linked genes, including, sex determining Region of Y-Chromosome (SRY), ubiquitously transcribed tetratricopeptide repeat gene protein on the chromosome Y (UTY), DEAD-box helicase 3 Y-linked (DDX3Y, also known as DBY), and oral-facial-digital syndrome 1 Y-linked (OFD1Y) in 254 swamp buffaloes from 15 populations distributed across China, Vietnam, and Laos using quantitative real-time PCR (qPCR). Our results revealed the prevalence of a single-copy UTY gene in buffaloes. The DBY and OFD1Y represented CNVs among and within different buffalo breeds. The SRY showed CNVs only in Vietnamese and Laotian buffaloes. In conclusion, this study indicated that DBY, OFD1Y, and SRY showed CNVs, while the UTY was a single-copy gene in swamp buffaloes.
Collapse
Affiliation(s)
- Ting Sun
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, Shaanxi, China
| | - Quratulain Hanif
- National Institute for Biotechnology and Genetic Engineering, Pakistan Institute of Engineering and Applied Sciences, Faisalabad 577, Pakistan
| | - Hong Chen
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, Shaanxi, China
| | - Chuzhao Lei
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, Shaanxi, China
| | - Ruihua Dang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling 712100, Shaanxi, China
- Correspondence: ; Tel.: +86-153-8862-7637
| |
Collapse
|
15
|
Guo J, Zhong J, Li L, Zhong T, Wang L, Song T, Zhang H. Comparative genome analyses reveal the unique genetic composition and selection signals underlying the phenotypic characteristics of three Chinese domestic goat breeds. Genet Sel Evol 2019; 51:70. [PMID: 31771503 PMCID: PMC6880376 DOI: 10.1186/s12711-019-0512-4] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2019] [Accepted: 11/15/2019] [Indexed: 12/31/2022] Open
Abstract
BACKGROUND As one of the important livestock species around the world, goats provide abundant meat, milk, and fiber to fulfill basic human needs. However, the genetic loci that underlie phenotypic variations in domestic goats are largely unknown, particularly for economically important traits. In this study, we sequenced the whole genome of 38 goats from three Chinese breeds (Chengdu Brown, Jintang Black, and Tibetan Cashmere) and downloaded the genome sequence data of 30 goats from five other breeds (four non-Chinese and one Chinese breed) and 21 Bezoar ibexes to investigate the genetic composition and selection signatures of the Chinese goat breeds after domestication. RESULTS Based on population structure analysis and FST values (average FST = 0.22), the genetic composition of Chengdu Brown goats differs considerably from that of Bezoar ibexes as a result of geographic isolation. Strikingly, the genes under selection that we identified in Tibetan Cashmere goats were significantly enriched in the categories hair growth and bone and nervous system development, possibly because they are involved in adaptation to high-altitude. In particular, we found a large difference in allele frequency of one novel SNP (c.-253G>A) in the 5'-UTR of FGF5 between Cashmere goats and goat breeds with short hair. The mutation at this site introduces a start codon that results in the occurrence of a premature FGF5 protein and is likely a natural causal variant that is involved in the long hair phenotype of cashmere goats. The haplotype tagged with the AGG-allele in exon 12 of DSG3, which encodes a cell adhesion molecule that is expressed mainly in the skin, was almost fixed in Tibetan Cashmere goats, whereas this locus still segregates in the lowland goat breeds. The pigmentation gene KITLG showed a strong signature of selection in Tibetan Cashmere goats. The genes ASIP and LCORL were identified as being under positive selection in Jintang Black goats. CONCLUSIONS After domestication, geographic isolation of some goat breeds has resulted in distinct genetic structures. Furthermore, our work highlights several positively selected genes that likely contributed to breed-related traits in domestic goats.
Collapse
Affiliation(s)
- Jiazhong Guo
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, 611130 China
| | - Jie Zhong
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, 611130 China
| | - Li Li
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, 611130 China
| | - Tao Zhong
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, 611130 China
| | - Linjie Wang
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, 611130 China
| | - Tianzeng Song
- Institute of Animal Science, Tibet Academy of Agricultural and Animal Husbandry Sciences, Lhasa, 850009 China
| | - Hongping Zhang
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu, 611130 China
| |
Collapse
|
16
|
Copy Number Variation of the CADM2 Gene and Its Association with Growth Traits in Yak. Animals (Basel) 2019; 9:ani9121008. [PMID: 31766342 PMCID: PMC6940794 DOI: 10.3390/ani9121008] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2019] [Revised: 11/16/2019] [Accepted: 11/18/2019] [Indexed: 12/16/2022] Open
Abstract
Simple Summary Cell adhesion molecule 2 (CADM2), also known as synaptic cell adhesion molecule 2 (SYNCAM2), is the mediator of synaptic signals enriched in the brain. Overlaps between copy number variation (CNV) regions in CADM2 and quantitative trait loci (QTL) related to body weight have been clarified in a previous study. In this study, two loci were amplified in the CADM2 gene (CNV1: 235,915 bp, exon 1 and partial intron 1; CNV2: 60,430 bp, intron 9) to explore the relationship between CNV types in the CADM2 gene and growth traits in 350 Ashidan yaks. Association analysis illustrated that no significant effect was found on growth traits in CNV1. However, the CNV2 mutation had a significant effect on body weight at the sixth month (p < 0.05). Individuals with the gain-type copy number variation CNV2 were significantly superior to those with loss- or normal-type in terms of body weight (p < 0.05). In summary, this study confirmed that CADM2-CNVs affect growth traits in yaks, and may be candidate genes for successful yak breeding and genetics projects. Abstract Copy number variation (CNV) is currently accepted as a common source of genetic variation. It is reported that CNVs may influence the resistance to disease and complex economic traits, such as residual feed intake, muscle formation, and fat deposition in livestock. Cell adhesion molecule 2 (CADM2) is expressed widely in the brain and adipose tissue and can regulate body weight through the central nervous system. Growth traits are important economic traits for animal selection. In this study, we aimed to explore the effect of CADM2 gene copy number variants on yak growth traits. Here, two CNVs in the CADM2 gene were investigated using the quantitative polymerase chain reaction (qPCR), and the association of the CNVs with growth traits in yak was analyzed using statistical methods by SPSS software. Differences were considered significant if the p value was < 0.05. Statistical analysis indicated significant association of CADM2-CNV2 with the body weight of the Chinese Ashidan yak. A significant effect of CNV2 (p < 0.05) was found on body weight at 6 months. In CNV2, the gain-type copy number variation exhibited greater performance than the other variants, with greater body weight observed at 6 months (p < 0.05). To the best of our knowledge, this is the first attempt to investigate the function of CADM2-CNVs and their association with growth traits in animals. This may be a useful candidate marker in marker-assisted selection of yaks.
Collapse
|
17
|
Keel BN, Nonneman DJ, Lindholm-Perry AK, Oliver WT, Rohrer GA. A Survey of Copy Number Variation in the Porcine Genome Detected From Whole-Genome Sequence. Front Genet 2019; 10:737. [PMID: 31475038 PMCID: PMC6707380 DOI: 10.3389/fgene.2019.00737] [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: 02/19/2019] [Accepted: 07/12/2019] [Indexed: 12/11/2022] Open
Abstract
Copy number variations (CNVs) are gains and losses of large regions of genomic sequence between individuals of a species. Although CNVs have been associated with various phenotypic traits in humans and other species, the extent to which CNVs impact phenotypic variation remains unclear. In swine, as well as many other species, relatively little is understood about the frequency of CNV in the genome, sizes, locations, and other chromosomal properties. In this work, we identified and characterized CNV by utilizing whole-genome sequence from 240 members of an intensely phenotyped experimental swine herd at the U.S. Meat Animal Research Center (USMARC). These animals included all 24 of the purebred founding boars (12 Duroc and 12 Landrace), 48 of the founding Yorkshire-Landrace composite sows, 109 composite animals from generations 4 through 9, 29 composite animals from generation 15, and 30 purebred industry boars (15 Landrace and 15 Yorkshire) used as sires in generations 10 through 15. Using a combination of split reads, paired-end mapping, and read depth approaches, we identified a total of 3,538 copy number variable regions (CNVRs), including 1,820 novel CNVRs not reported in previous studies. The CNVRs covered 0.94% of the porcine genome and overlapped 1,401 genes. Gene ontology analysis identified that CNV-overlapped genes were enriched for functions related to organism development. Additionally, CNVRs overlapped with many known quantitative trait loci (QTL). In particular, analysis of QTL previously identified in the USMARC herd showed that CNVRs were most overlapped with reproductive traits, such as age of puberty and ovulation rate, and CNVRs were significantly enriched for reproductive QTL.
Collapse
Affiliation(s)
- Brittney N Keel
- USDA, ARS, U.S. Meat Animal Research Center, Clay Center, NE, United States
| | - Dan J Nonneman
- USDA, ARS, U.S. Meat Animal Research Center, Clay Center, NE, United States
| | | | - William T Oliver
- USDA, ARS, U.S. Meat Animal Research Center, Clay Center, NE, United States
| | - Gary A Rohrer
- USDA, ARS, U.S. Meat Animal Research Center, Clay Center, NE, United States
| |
Collapse
|
18
|
Pei SW, Qin F, Li WH, Li FD, Yue XP. Copy number variation of ZNF280AY across 21 cattle breeds and its association with the reproductive traits of Holstein and Simmental bulls. J Dairy Sci 2019; 102:7226-7236. [PMID: 31202648 DOI: 10.3168/jds.2018-16063] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2018] [Accepted: 04/23/2019] [Indexed: 11/19/2022]
Abstract
The mammalian Y chromosome gene families in the ampliconic region are expressed predominantly or exclusively in the testis, and their copy number variations (CNV) are significantly associated with male reproductive traits, suggesting they have important roles in spermatogenesis and testicular development. ZNF280AY (zinc finger protein 280A, Y-linked) is a member of the zinc finger protein family and has been identified as a bovid-specific Y-chromosome gene. The current study applied a reliable quantitative real-time PCR method to estimate the CNV of ZNF280AY in 715 bulls across 21 cattle breeds and to further investigate the association of the CNV of ZNF280AY with bull reproductive traits and ZNF280AY mRNA expression levels in adult testis. The results revealed that the median copy number of ZNF280AY was 47, and the copy number varied from 11 to 154, showing significant CNV between and within the investigated cattle breeds. In addition, all 715 bulls were classified into Y1, Y2, and Y3 lineage groups based on a rapid genotyping method described previously. Pairwise comparisons indicated that bulls belonging to the Y1 lineage had a significantly lower median copy number (40) than bulls belonging to the Y2 (52) and Y3 lineages (57). Association analysis revealed that the CNV of ZNF280AY was correlated negatively with the percentage of normal sperm and sperm concentration in Holstein bulls, whereas no significant correlation was observed with ejaculation volume, total sperm count, sperm motility, postthaw motility (PTM), and scrotal circumference in Holstein and Simmental bulls. Furthermore, no correlation was observed between ZNF280AY copy number and ZNF280AY mRNA expression levels in the testis. The current study suggests that the CNV of the ZNF280AY gene family is associated with male reproductive traits and may serve as a valuable marker for early bull fertility selection in Holstein breeding programs.
Collapse
Affiliation(s)
- S W Pei
- State Key Laboratory of Grassland Agro-Ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs; Engineering Research Center of Grassland Industry, Ministry of Education; College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, 730020, P. R. China
| | - F Qin
- School of Pharmacy, Lanzhou University, Lanzhou, 730000, P. R. China
| | - W H Li
- State Key Laboratory of Grassland Agro-Ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs; Engineering Research Center of Grassland Industry, Ministry of Education; College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, 730020, P. R. China
| | - F D Li
- State Key Laboratory of Grassland Agro-Ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs; Engineering Research Center of Grassland Industry, Ministry of Education; College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, 730020, P. R. China; Engineering Laboratory of Sheep Breeding and Reproduction Biotechnology in Gansu Province, Minqin 733300, P. R. China
| | - X P Yue
- State Key Laboratory of Grassland Agro-Ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs; Engineering Research Center of Grassland Industry, Ministry of Education; College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, 730020, P. R. China.
| |
Collapse
|
19
|
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
|
20
|
May K, Scheper C, Brügemann K, Yin T, Strube C, Korkuć P, Brockmann GA, König S. Genome-wide associations and functional gene analyses for endoparasite resistance in an endangered population of native German Black Pied cattle. BMC Genomics 2019; 20:277. [PMID: 30961534 PMCID: PMC6454736 DOI: 10.1186/s12864-019-5659-4] [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: 09/03/2018] [Accepted: 03/29/2019] [Indexed: 12/14/2022] Open
Abstract
Background Gastrointestinal nematodes (GIN), liver flukes (Fasciola hepatica) and bovine lungworms (Dictyocaulus viviparus) are the most important parasitic agents in pastured dairy cattle. Endoparasite infections are associated with reduced milk production and detrimental impacts on female fertility, contributing to economic losses in affected farms. In quantitative-genetic studies, the heritabilities for GIN and F. hepatica were moderate, encouraging studies on genomic scales. Genome-wide association studies (GWAS) based on dense single nucleotide polymorphism (SNP) marker panels allow exploration of the underlying genomic architecture of complex disease traits. The current GWAS combined the identification of potential candidate genes with pathway analyses to obtain deeper insights into bovine immune response and the mechanisms of resistance against endoparasite infections. Results A 2-step approach was applied to infer genome-wide associations in an endangered dual-purpose cattle subpopulation [Deutsches Schwarzbuntes Niederungsrind (DSN)] with a limited number of phenotypic records. First, endoparasite traits from a population of 1166 Black and White dairy cows [including Holstein Friesian (HF) and DSN] naturally infected with GIN, F. hepatica and D. viviparus were precorrected for fixed effects using linear mixed models. Afterwards, the precorrected phenotypes were the dependent traits (rFEC-GIN, rFEC-FH, and rFLC-DV) in GWAS based on 423,654 SNPs from 148 DSN cows. We identified 44 SNPs above the genome-wide significance threshold (pBonf = 4.47 × 10− 7), and 145 associations surpassed the chromosome-wide significance threshold (range: 7.47 × 10− 6 on BTA 1 to 2.18 × 10− 5 on BTA 28). The associated SNPs identified were annotated to 23 candidate genes. The DAVID analysis inferred four pathways as being related to immune response mechanisms or involved in host-parasite interactions. SNP effect correlations considering specific chromosome segments indicate that breeding for resistance to GIN or F. hepatica as measured by fecal egg counts is genetically associated with a higher risk for udder infections. Conclusions We detected a large number of loci with small to moderate effects for endoparasite resistance. The potential candidate genes regulating resistance identified were pathogen-specific. Genetic antagonistic associations between disease resistance and productivity were specific for specific chromosome segments. The 2-step approach was a valid methodological approach to infer genetic mechanisms in an endangered breed with a limited number of phenotypic records. Electronic supplementary material The online version of this article (10.1186/s12864-019-5659-4) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Katharina May
- Institute of Animal Breeding and Genetics, Justus-Liebig-University of Gießen, 35390, Gießen, Germany.,Institute for Parasitology, Center for Infection Medicine, University of Veterinary Medicine Hanover, 30559, Hannover, Germany
| | - Carsten Scheper
- Institute of Animal Breeding and Genetics, Justus-Liebig-University of Gießen, 35390, Gießen, Germany
| | - Kerstin Brügemann
- Institute of Animal Breeding and Genetics, Justus-Liebig-University of Gießen, 35390, Gießen, Germany
| | - Tong Yin
- Institute of Animal Breeding and Genetics, Justus-Liebig-University of Gießen, 35390, Gießen, Germany
| | - Christina Strube
- Institute for Parasitology, Center for Infection Medicine, University of Veterinary Medicine Hanover, 30559, Hannover, Germany
| | - Paula Korkuć
- Department for Crop and Animal Sciences, Breeding Biology and Molecular Genetics, Faculty of Live Science, Humboldt-Universität of Berlin, 10115, Berlin, Germany
| | - Gudrun A Brockmann
- Department for Crop and Animal Sciences, Breeding Biology and Molecular Genetics, Faculty of Live Science, Humboldt-Universität of Berlin, 10115, Berlin, Germany
| | - Sven König
- Institute of Animal Breeding and Genetics, Justus-Liebig-University of Gießen, 35390, Gießen, Germany.
| |
Collapse
|
21
|
Goshu HA, Chu M, Xiaoyun W, Pengjia B, Zhi DX, Yan P. Genomic copy number variation of the CHKB gene alters gene expression and affects growth traits of Chinese domestic yak (Bos grunniens) breeds. Mol Genet Genomics 2019; 294:549-561. [DOI: 10.1007/s00438-018-01530-y] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2018] [Accepted: 12/29/2018] [Indexed: 12/22/2022]
|
22
|
Xu L, Yang L, Wang L, Zhu B, Chen Y, Gao H, Gao X, Zhang L, Liu GE, Li J. Probe-based association analysis identifies several deletions associated with average daily gain in beef cattle. BMC Genomics 2019; 20:31. [PMID: 30630414 PMCID: PMC6327516 DOI: 10.1186/s12864-018-5403-5] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2018] [Accepted: 12/20/2018] [Indexed: 12/20/2022] Open
Abstract
BACKGROUND Average daily gain (ADG) is an important trait that contributes to the production efficiency and economic benefits in the beef cattle industry. The molecular mechanisms of ADG have not yet been fully explored because most recent association studies for ADG are based on SNPs or haplotypes. We reported a systematic CNV discovery and association analysis for ADG in Chinese Simmental beef cattle. RESULTS Our study identified 4912 nonredundant CNVRs with a total length of ~ 248.7 Mb, corresponding to ~ 8.9% of the cattle genome. Using probe-based CNV association, we identified 24 and 12 significant SNP probes within five deletions and two duplications for ADG, respectively. Among them, we found one common deletion with 89 kb imbedded in LHFPL Tetraspan Subfamily Member 6 (LHFPL6) at 22.9 Mb on BTA12, which has high frequency (12.9%) dispersing across population. CNV selection test using VST statistic suggested this common deletion may be under positive selection in Chinese Simmental cattle. Moreover, this deletion was not overlapped with any candidate SNP for ADG compared with previous SNPs-based association studies, suggesting its important role for ADG. In addition, we identified one rare deletion near gene Growth Factor Receptor-bound Protein 10 (GRB10) at 5.1 Mb on BTA4 for ADG using both probe-based association and region-based approaches. CONCLUSIONS Our results provided some valuable insights to elucidate the genetic basis of ADG in beef cattle, and these findings offer an alternative perspective to understand the genetic mechanism of complex traits in terms of copy number variations in farm animals.
Collapse
Affiliation(s)
- Lingyang Xu
- Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China.
| | - Liu Yang
- Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China.,Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, 611130, China
| | - Lei Wang
- Beijing Genecast Biotechnology Co., Beijing, 100191, China
| | - Bo Zhu
- Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Yan Chen
- Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Huijiang Gao
- Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Xue Gao
- Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Lupei Zhang
- Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - George E Liu
- U.S. Department of Agriculture-Agricultural Research Services, Animal Genomics and Improvement Laboratory, Beltsville, MD, 20705, USA.
| | - Junya Li
- Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China.
| |
Collapse
|
23
|
Di Gerlando R, Sardina MT, Tolone M, Sutera AM, Mastrangelo S, Portolano B. Genome-wide detection of copy-number variations in local cattle breeds. ANIMAL PRODUCTION SCIENCE 2019. [DOI: 10.1071/an17603] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
The aim of the present study was to identify copy-number variations (CNVs) in Cinisara (CIN) and Modicana (MOD) cattle breeds on the basis of signal intensity (logR ratio) and B allele frequency of each marker, using Illumina’s BovineSNP50K Genotyping BeadChip. The CNVs were detected with the PennCNV and SVS 8.7.0 software and were aggregated into CNV regions (CNVRs). PennCNV identified 487 CNVs in CIN that aggregated into 86 CNVRs, and 424 CNVs in MOD that aggregated into 81 CNVRs. SVS identified a total of 207 CNVs in CIN that aggregated into 39 CNVRs, and 181 CNVs in MOD that aggregated into 41 CNVRs. The CNVRs identified with the two softwares contained 29 common CNVRs in CIN and 17 common CNVRs in MOD. Only a small number of CNVRs identified in the present study have been identified elsewhere, probably because of the limitations of the array used. In total, 178 and 208 genes were found within the CNVRs of CIN and MOD respectively. Gene Ontology and KEGG pathway analyses showed that several of these genes are involved in milk production, reproduction and behaviour, the immune response, and resistance/susceptibility to infectious diseases. Our results have provided significant information for the construction of more-complete CNV maps of the bovine genome and offer an important resource for the investigation of genomic changes and traits of interest in the CIN and MOD cattle breeds. Our results will also be valuable for future studies and constitute a preliminary report of the CNV distribution resources in local cattle genomes.
Collapse
|
24
|
Bai W, Zhang C, Chen H. Transcriptomic analysis of Momordica charantia polysaccharide on streptozotocin-induced diabetic rats. Gene 2018; 675:208-216. [DOI: 10.1016/j.gene.2018.06.106] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2017] [Revised: 06/11/2018] [Accepted: 06/29/2018] [Indexed: 10/28/2022]
|
25
|
Immunity to gastrointestinal nematodes in ruminants: effector cell mechanisms and cytokines. J Parasit Dis 2018; 42:471-482. [PMID: 30538343 DOI: 10.1007/s12639-018-1023-x] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2018] [Accepted: 08/02/2018] [Indexed: 12/18/2022] Open
Abstract
Gastrointestinal nematodes (GINs) of ruminants are prevalent and have major economic impacts worldwide. The insight studies of immune responses triggered against GINs are of great concern to understand interaction between host's immune system and parasite. T-helper 2 cytokines drive the effector cell mechanisms which include eosinophils and mast cells. The immune responses are controlled by Th2 secreted interleukins (IL); IL3, IL-4, IL-5, IL-9, IL-10 and IL-13. B-Cell immune response is incorporated in defense mechanisms developed against GINs specially immunoglobulins (Ig); IgA, IgE and IgG. The immune resistance of the infected host is presented by failure of larval establishment or hypobiosis, low worm burden and decreased female fecundity. The host-parasite interaction is a complex series that affected by host's genetic constitution, nutrition, age and physiological status. The GINs have different immune evasion mechanisms to improve their survival within the host. Also, management of the host influences GINs parasitism. Thus, the aim of this review is to highlight the hallmarks of immune responses that endorse GINs parasitism. The insights studies of the triggered immune responses developed against GINs will improve the appropriate protective immune strategy.
Collapse
|
26
|
Old origin of a protective endogenous retrovirus (enJSRV) in the Ovis genus. Heredity (Edinb) 2018; 122:187-194. [PMID: 29976957 DOI: 10.1038/s41437-018-0112-z] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2018] [Revised: 06/07/2018] [Accepted: 06/15/2018] [Indexed: 11/08/2022] Open
Abstract
Sheep, the Jaagsiekte sheep retrovirus (JSRV) and its endogenous forms (enJSRVs) are a good model to study long-time relationships between retroviruses and their hosts. Taking advantage of 76 whole genome resequencing data of wild and domestic Ovis, we investigated the evolution of this relationship. An innovative analysis of re-sequencing data allowed characterizing 462 enJSRVs insertion sites (including 435 newly described insertions) in the Ovis genus. We focused our study on endogenous copies inserted in the q13 locus of chromosome 6 (6q13). Those copies are known to confer resistance against exogenous JSRV thanks to alleles bearing a mutation in the gag gene. We characterized (i) the distribution of protective and non-protective alleles across Ovis species and (ii) the copy number variation of the 6q13 locus. Our results challenged the previous hypothesis of fixation and amplification of the protective copies in relation with domestication, and allowed building a new model for the evolution of the 6q13 locus. JSRV would have integrated the 6q13 locus after the Ovis-Capra divergence (5-11 MYA) and before the Ovis diversification (2.4-5 MYA). The protective mutation in the enJSRV 6q13 copy appeared shortly after its insertion and was followed by genomic amplifications, after the divergence between Pachyform lineage on one side and the Argaliform and moufloniform lineages on the other (2.4-5 MYA). Considering the potential selective advantage of the protective mutation, its fixation in both sheep and its closest wild relative Ovis orientalis may be due to natural selection before domestication from O. orientalis populations.
Collapse
|
27
|
Genomic predictions combining SNP markers and copy number variations in Nellore cattle. BMC Genomics 2018; 19:441. [PMID: 29871610 PMCID: PMC5989480 DOI: 10.1186/s12864-018-4787-6] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2017] [Accepted: 05/14/2018] [Indexed: 11/26/2022] Open
Abstract
Background Due to the advancement in high throughput technology, single nucleotide polymorphism (SNP) is routinely being incorporated along with phenotypic information into genetic evaluation. However, this approach often cannot achieve high accuracy for some complex traits. It is possible that SNP markers are not sufficient to predict these traits due to the missing heritability caused by other genetic variations such as microsatellite and copy number variation (CNV), which have been shown to affect disease and complex traits in humans and other species. Results In this study, CNVs were included in a SNP based genomic selection framework. A Nellore cattle dataset consisting of 2230 animals genotyped on BovineHD SNP array was used, and 9 weight and carcass traits were analyzed. A total of six models were implemented and compared based on their prediction accuracy. For comparison, three models including only SNPs were implemented: 1) BayesA model, 2) Bayesian mixture model (BayesB), and 3) a GBLUP model without polygenic effects. The other three models incorporating both SNP and CNV included 4) a Bayesian model similar to BayesA (BayesA+CNV), 5) a Bayesian mixture model (BayesB+CNV), and 6) GBLUP with CNVs modeled as a covariable (GBLUP+CNV). Prediction accuracies were assessed based on Pearson’s correlation between de-regressed EBVs (dEBVs) and direct genomic values (DGVs) in the validation dataset. For BayesA, BayesB and GBLUP, accuracy ranged from 0.12 to 0.62 across the nine traits. A minimal increase in prediction accuracy for some traits was noticed when including CNVs in the model (BayesA+CNV, BayesB+CNV, GBLUP+CNV). Conclusions This study presents the first genomic prediction study integrating CNVs and SNPs in livestock. Combining CNV and SNP marker information proved to be beneficial for genomic prediction of some traits in Nellore cattle. Electronic supplementary material The online version of this article (10.1186/s12864-018-4787-6) contains supplementary material, which is available to authorized users.
Collapse
|
28
|
Mielczarek M, Frąszczak M, Nicolazzi E, Williams JL, Szyda J. Landscape of copy number variations in Bos taurus: individual - and inter-breed variability. BMC Genomics 2018; 19:410. [PMID: 29843606 PMCID: PMC5975385 DOI: 10.1186/s12864-018-4815-6] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2018] [Accepted: 05/22/2018] [Indexed: 11/24/2022] Open
Abstract
Background The number of studies of Copy Number Variation in cattle has increased in recent years. This has been prompted by the increased availability of data on polymorphisms and their relationship with phenotypes. In addition, livestock species are good models for some human phenotypes. In the present study, we described the landscape of CNV driven genetic variation in a large population of 146 individuals representing 13 cattle breeds, using whole genome DNA sequence. Results A highly significant variation among all individuals and within each breed was observed in the number of duplications (P < 10−15) and in the number of deletions (P < 10−15). We also observed significant differences between breeds for duplication (P = 0.01932) and deletion (P = 0.01006) counts. The same variation CNV length - inter-individual and inter-breed differences were significant for duplications (P < 10−15) and deletions (P < 10−15). Moreover, breed-specific variants were identified, with the largest proportion of breed-specific duplications (9.57%) found for Fleckvieh and breed-specific deletions found for Brown Swiss (5.00%). Such breed-specific CNVs were predominantly located in intragenic regions, however in Simmental, one deletion present in five individuals was found in the coding sequence of a novel gene ENSBTAG00000000688 on chromosome 18. In Brown Swiss, Norwegian Red and Simmental breed-specific deletions were located within KIT and MC1R genes, which are responsible for a coat colour. The functional annotation of coding regions underlying the breed-specific CNVs showed that in Norwegian Red, Guernsey, and Simmental significantly under- and overrepresented GO terms were related to chemical stimulus involved in sensory perception of smell and the KEGG pathways for olfactory transduction. In addition, specifically for the Norwegian Red breed, the dopaminergic synapse KEGG pathway was significantly enriched within deleted parts of the genome. Conclusions The CNV landscape in Bos taurus genome revealed by this study was highly complex, with inter-breed differences, but also a significant variation within breeds. The former, may explain some of the phenotypic differences among analysed breeds, and the latter contributes to within-breed variation available for selection. Electronic supplementary material The online version of this article (10.1186/s12864-018-4815-6) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- M Mielczarek
- Biostatistics group, Department of Genetics, Wroclaw University of Environmental and Life Sciences, Kozuchowska 7, 51-631, Wroclaw, Poland. .,National Research Institute of Animal Production, Krakowska 1, 32-083, Balice, Poland.
| | - M Frąszczak
- Biostatistics group, Department of Genetics, Wroclaw University of Environmental and Life Sciences, Kozuchowska 7, 51-631, Wroclaw, Poland
| | - E Nicolazzi
- Council on Dairy Cattle Breeding (CDCB), 4201 Northview Dr, Bowie, MD, 20716, USA
| | - J L Williams
- Davies Research Centre, University of Adelaide, School of Animal and Veterinary Sciences, Roseworthy, SA, 5371, Australia
| | - J Szyda
- Biostatistics group, Department of Genetics, Wroclaw University of Environmental and Life Sciences, Kozuchowska 7, 51-631, Wroclaw, Poland.,National Research Institute of Animal Production, Krakowska 1, 32-083, Balice, Poland
| |
Collapse
|
29
|
Zhou Y, Connor EE, Wiggans GR, Lu Y, Tempelman RJ, Schroeder SG, Chen H, Liu GE. Genome-wide copy number variant analysis reveals variants associated with 10 diverse production traits in Holstein cattle. BMC Genomics 2018; 19:314. [PMID: 29716533 PMCID: PMC5930521 DOI: 10.1186/s12864-018-4699-5] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2017] [Accepted: 04/18/2018] [Indexed: 01/08/2023] Open
Abstract
BACKGROUND Copy number variation (CNV) is an important type of genetic variation contributing to phenotypic differences among mammals and may serve as an alternative molecular marker to single nucleotide polymorphism (SNP) for genome-wide association study (GWAS). Recently, GWAS analysis using CNV has been applied in livestock, although few studies have focused on Holstein cattle. RESULTS We describe 191 CNV detected using intensity data from over 700,000 SNP genotypes generated with the BovineHD Genotyping BeadChip (Illumina, San Diego, CA) in 528 Holstein cows. The CNV were used for GWAS analysis of 10 important production traits of 473 cattle related to feed intake, milk quality, and female fertility, as well as 2 composite traits of net merit and productive life. In total, we detected 57 CNV associated (P < 0.05 after false discovery rate correction) with at least one of the 10 phenotypes. Focusing on feed efficiency and intake-related phenotypes of residual feed intake and dry matter intake, we detected a single CNV associated with both traits which overlaps a predicted olfactory receptor gene OR2A2 (LOC787786). Additionally, 2 CNV within the RXFP4 (relaxin/insulin like family peptide receptor 4) and 2 additional olfactory receptor gene regions, respectively, were associated with residual feed intake. The RXFP4 gene encodes a receptor for an orexigenic peptide, insulin-like peptide 5 produced by intestinal L cells, which is expressed by enteric neurons. Olfactory receptors are critical for transmitting the effects of odorants, contributing to the sense of smell, and have been implicated in participating in appetite regulation. CONCLUSIONS Our results identify CNV for genomic evaluation in Holstein cattle, and provide candidate genes, such as RXFP4, contributing to variation in feed efficiency and feed intake-related traits. These results indicate potential novel targets for manipulating feed intake-related traits of livestock.
Collapse
Affiliation(s)
- Yang Zhou
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, 10300 Baltimore Avenue, Bldg. 306, BARC-East, Beltsville, MD, 20705, USA.,Shaanxi Key Laboratory of Agricultural Molecular Biology, College of Animal Science and Technology, Northwest A&F University, Yangling, 712100, Shaanxi, China.,Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction, Education Ministry of China, Huazhong Agricultural University, Wuhan, 430070, Hubei, China
| | - Erin E Connor
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, 10300 Baltimore Avenue, Bldg. 306, BARC-East, Beltsville, MD, 20705, USA
| | - George R Wiggans
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, 10300 Baltimore Avenue, Bldg. 306, BARC-East, Beltsville, MD, 20705, USA
| | - Yongfang Lu
- Department of Animal Science, Michigan State University, East Lansing, MI, 48824, USA
| | - Robert J Tempelman
- Department of Animal Science, Michigan State University, East Lansing, MI, 48824, USA
| | - Steven G Schroeder
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, 10300 Baltimore Avenue, Bldg. 306, BARC-East, Beltsville, MD, 20705, USA
| | - Hong Chen
- Shaanxi Key Laboratory of Agricultural Molecular Biology, College of Animal Science and Technology, Northwest A&F University, Yangling, 712100, Shaanxi, China
| | - George E Liu
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, 10300 Baltimore Avenue, Bldg. 306, BARC-East, Beltsville, MD, 20705, USA.
| |
Collapse
|
30
|
Tsairidou S, Allen AR, Pong‐Wong R, McBride SH, Wright DM, Matika O, Pooley CM, McDowell SWJ, Glass EJ, Skuce RA, Bishop SC, Woolliams JA. An analysis of effects of heterozygosity in dairy cattle for bovine tuberculosis resistance. Anim Genet 2018; 49:103-109. [PMID: 29368428 PMCID: PMC5888165 DOI: 10.1111/age.12637] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/23/2017] [Indexed: 12/13/2022]
Abstract
Genetic selection of cattle more resistant to bovine tuberculosis (bTB) may offer a complementary control strategy. Hypothesising underlying non-additive genetic variation, we present an approach using genome-wide high density markers to identify genomic loci with dominance effects on bTB resistance and to test previously published regions with heterozygote advantage in bTB. Our data comprised 1151 Holstein-Friesian cows from Northern Ireland, confirmed bTB cases and controls, genotyped with the 700K Illumina BeadChip. Genome-wide markers were tested for associations between heterozygosity and bTB status using marker-based relationships. Results were tested for robustness against genetic structure, and the genotypic frequencies of a significant locus were tested for departures from Hardy-Weinberg equilibrium. Genomic regions identified in our study and in previous publications were tested for dominance effects. Genotypic effects were estimated through ASReml mixed models. A SNP (rs43032684) on chromosome 6 was significant at the chromosome-wide level, explaining 1.7% of the phenotypic variance. In the controls, there were fewer heterozygotes for rs43032684 (P < 0.01) with the genotypic values suggesting that heterozygosity confers a heterozygote disadvantage. The region surrounding rs43032684 had a significant dominance effect (P < 0.01). SNP rs43032684 resides within a pseudogene with a parental gene involved in macrophage response to infection and within a copy-number-variation region previously associated with nematode resistance. No dominance effect was found for the region on chromosome 11, as indicated by a previous candidate region bTB study. These findings require further validation with large-scale data.
Collapse
Affiliation(s)
- S. Tsairidou
- The Roslin Institute and R(D)SVSUniversity of EdinburghEdinburghEH259RGUK
| | - A. R. Allen
- Veterinary Sciences DivisionAgri‐Food and Biosciences InstituteBelfastBT95PXUK
| | - R. Pong‐Wong
- The Roslin Institute and R(D)SVSUniversity of EdinburghEdinburghEH259RGUK
| | - S. H. McBride
- Veterinary Sciences DivisionAgri‐Food and Biosciences InstituteBelfastBT95PXUK
| | - D. M. Wright
- School of Biological SciencesQueen's University of BelfastBelfastBT71NNUK
| | - O. Matika
- The Roslin Institute and R(D)SVSUniversity of EdinburghEdinburghEH259RGUK
| | - C. M. Pooley
- The Roslin Institute and R(D)SVSUniversity of EdinburghEdinburghEH259RGUK
| | - S. W. J. McDowell
- Veterinary Sciences DivisionAgri‐Food and Biosciences InstituteBelfastBT95PXUK
| | - E. J. Glass
- The Roslin Institute and R(D)SVSUniversity of EdinburghEdinburghEH259RGUK
| | - R. A. Skuce
- Veterinary Sciences DivisionAgri‐Food and Biosciences InstituteBelfastBT95PXUK
- School of Biological SciencesQueen's University of BelfastBelfastBT71NNUK
| | - S. C. Bishop
- The Roslin Institute and R(D)SVSUniversity of EdinburghEdinburghEH259RGUK
| | - J. A. Woolliams
- The Roslin Institute and R(D)SVSUniversity of EdinburghEdinburghEH259RGUK
| |
Collapse
|
31
|
Karimi K, Esmailizadeh A, Wu DD, Gondro C. Mapping of genome-wide copy number variations in the Iranian indigenous cattle using a dense SNP data set. ANIMAL PRODUCTION SCIENCE 2018. [DOI: 10.1071/an16384] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
The objective of this study was to present the first map of the copy number variations (CNVs) in Iranian indigenous cattle based on a high-density single nucleotide polymorphism (SNP) dataset. A total of 90 individuals were genotyped using the Illumina BovineHD BeadChip containing 777 962 SNPs. The QuantiSNP algorithm was used to perform a genome-wide CNV detection across autosomal genome. After merging the overlapping CNV, a total of 221 CNV regions were identified encompassing 36.4 Mb or 1.44% of the bovine autosomal genome. The length of the CNV regions ranged from 3.5 to 2252.8 Kb with an average of 163.8 Kb. These regions included 147 loss (66.52%) and 74 gain (33.48%) events containing a total of 637 annotated Ensembl genes. Gene ontology analysis revealed that most of genes in the CNV regions were involved in environmental responses, disease susceptibility and immune system functions. Furthermore, 543 of these genes corresponded to the human orthologous genes, which involved in a wide range of biological functions. Altogether, 73% of the 221 CNV regions overlapped either completely or partially with those previously reported in other cattle studies. Moreover, novel CNV regions involved several quantitative trait loci (QTL)-related to adaptative traits of Iranian indigenous cattle. These results provided a basis to conduct future studies on association between CNV regions and phenotypic variations in the Iranian indigenous cattle.
Collapse
|
32
|
Yang L, Xu L, Zhu B, Niu H, Zhang W, Miao J, Shi X, Zhang M, Chen Y, Zhang L, Gao X, Gao H, Li L, Liu GE, Li J. Genome-wide analysis reveals differential selection involved with copy number variation in diverse Chinese Cattle. Sci Rep 2017; 7:14299. [PMID: 29085051 PMCID: PMC5662686 DOI: 10.1038/s41598-017-14768-0] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2017] [Accepted: 10/12/2017] [Indexed: 12/20/2022] Open
Abstract
Copy number variations (CNVs) are defined as deletions, insertions, and duplications between two individuals of a species. To investigate the diversity and population-genetic properties of CNVs and their diverse selection patterns, we performed a genome-wide CNV analysis using high density SNP array in Chinese native cattle. In this study, we detected a total of 13,225 CNV events and 3,356 CNV regions (CNVRs), overlapping with 1,522 annotated genes. Among them, approximately 71.43 Mb of novel CNVRs were detected in the Chinese cattle population for the first time, representing the unique genomic resources in cattle. A new V i statistic was proposed to estimate the region-specific divergence in CNVR for each group based on unbiased estimates of pairwise V ST . We obtained 12 and 62 candidate CNVRs at the top 1% and top 5% of genome-wide V i value thresholds for each of four groups (North, Northwest, Southwest and South). Moreover, we identified many lineage-differentiated CNV genes across four groups, which were associated with several important molecular functions and biological processes, including metabolic process, response to stimulus, immune system, and others. Our findings provide some insights into understanding lineage-differentiated CNVs under divergent selection in the Chinese native cattle.
Collapse
Affiliation(s)
- Liu Yang
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China.,Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan, 611130, China
| | - Lingyang Xu
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China.
| | - Bo Zhu
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Hong Niu
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Wengang Zhang
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Jian Miao
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China.,College of Animal Sciences, Fujian Agriculture and Forestry University, Fuzhou, Fujian, 350002, China
| | - Xinping Shi
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China.,College of Animal Science and Technology, Agricultural University of Hebei, Baoding, Hebei, 071001, China
| | - Ming Zhang
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan, 611130, China
| | - Yan Chen
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Lupei Zhang
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Xue Gao
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Huijiang Gao
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Li Li
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan, 611130, China
| | - George E Liu
- Animal Genomics and Improvement Laboratory, Agricultural Research Service, USDA, Beltsville, Maryland, 20705, USA
| | - Junya Li
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China.
| |
Collapse
|
33
|
Letaief R, Rebours E, Grohs C, Meersseman C, Fritz S, Trouilh L, Esquerré D, Barbieri J, Klopp C, Philippe R, Blanquet V, Boichard D, Rocha D, Boussaha M. Identification of copy number variation in French dairy and beef breeds using next-generation sequencing. Genet Sel Evol 2017; 49:77. [PMID: 29065859 PMCID: PMC5655909 DOI: 10.1186/s12711-017-0352-z] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2017] [Accepted: 10/17/2017] [Indexed: 11/15/2022] Open
Abstract
Background Copy number variations (CNV) are known to play a major role in genetic variability and disease pathogenesis in several species including cattle. In this study, we report the identification and characterization of CNV in eight French beef and dairy breeds using whole-genome sequence data from 200 animals. Bioinformatics analyses to search for CNV were carried out using four different but complementary tools and we validated a subset of the CNV by both in silico and experimental approaches.
Results We report the identification and localization of 4178 putative deletion-only, duplication-only and CNV regions, which cover 6% of the bovine autosomal genome; they were validated by two in silico approaches and/or experimentally validated using array-based comparative genomic hybridization and single nucleotide polymorphism genotyping arrays. The size of these variants ranged from 334 bp to 7.7 Mb, with an average size of ~ 54 kb. Of these 4178 variants, 3940 were deletions, 67 were duplications and 171 corresponded to both deletions and duplications, which were defined as potential CNV regions. Gene content analysis revealed that, among these variants, 1100 deletions and duplications encompassed 1803 known genes, which affect a wide spectrum of molecular functions, and 1095 overlapped with known QTL regions. Conclusions Our study is a large-scale survey of CNV in eight French dairy and beef breeds. These CNV will be useful to study the link between genetic variability and economically important traits, and to improve our knowledge on the genomic architecture of cattle. Electronic supplementary material The online version of this article (doi:10.1186/s12711-017-0352-z) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Rabia Letaief
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78352, Jouy-en-Josas, France.
| | - Emmanuelle Rebours
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78352, Jouy-en-Josas, France
| | - Cécile Grohs
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78352, Jouy-en-Josas, France
| | - Cédric Meersseman
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78352, Jouy-en-Josas, France.,GMA, INRA, Université de Limoges, UMR1061, Unité Génétique Moléculaire Animale, 123 avenue Albert Thomas, 87060, Limoges Cedex, France
| | - Sébastien Fritz
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78352, Jouy-en-Josas, France.,Allice, Maison Nationale des Eleveurs, 75012, Paris, France
| | - Lidwine Trouilh
- LISBP, CNRS, INRA, INSA, Université de Toulouse, Toulouse, France
| | - Diane Esquerré
- GenPhySE, INRA, Université de Toulouse INPT ENSAT, Université de Toulouse INPT ENVT, 52627, Castanet-Tolosan, France
| | - Johanna Barbieri
- GenPhySE, INRA, Université de Toulouse INPT ENSAT, Université de Toulouse INPT ENVT, 52627, Castanet-Tolosan, France
| | | | - Romain Philippe
- GMA, INRA, Université de Limoges, UMR1061, Unité Génétique Moléculaire Animale, 123 avenue Albert Thomas, 87060, Limoges Cedex, France
| | - Véronique Blanquet
- GMA, INRA, Université de Limoges, UMR1061, Unité Génétique Moléculaire Animale, 123 avenue Albert Thomas, 87060, Limoges Cedex, France
| | - Didier Boichard
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78352, Jouy-en-Josas, France
| | - Dominique Rocha
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78352, Jouy-en-Josas, France
| | - Mekki Boussaha
- GABI, INRA, AgroParisTech, Université Paris-Saclay, 78352, Jouy-en-Josas, France
| |
Collapse
|
34
|
Yan J, Blair HT, Liu M, Li W, He S, Chen L, Dittmer KE, Garrick DJ, Biggs PJ, Dukkipati VS. Genome-wide detection of autosomal copy number variants in several sheep breeds using Illumina OvineSNP50 BeadChips. Small Rumin Res 2017. [DOI: 10.1016/j.smallrumres.2017.08.022] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
35
|
Huang J, Li R, Zhang X, Huang Y, Dang R, Lan X, Chen H, Lei C. Copy number veriation regions detection in Qinchuan cattle. Livest Sci 2017. [DOI: 10.1016/j.livsci.2017.08.016] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
36
|
Genome-Wide SNP Signal Intensity Scanning Revealed Genes Differentiating Cows with Ovarian Pathologies from Healthy Cows. SENSORS 2017; 17:s17081920. [PMID: 28825658 PMCID: PMC5580201 DOI: 10.3390/s17081920] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/19/2017] [Revised: 07/04/2017] [Accepted: 07/06/2017] [Indexed: 01/09/2023]
Abstract
Hypoplasia and ovarian cysts are the most common ovarian pathologies in cattle. In this genome-wide study we analyzed the signal intensity of 648,315 Single Nucleotide Polymorphisms (SNPs) and identified 1338 genes differentiating cows with ovarian pathologies from healthy cows. The sample consisted of six cows presenting an ovarian pathology and six healthy cows. SNP signal intensities were measured with a genotyping process using the Axiom Genome-Wide BOS 1 SNPchip. Statistical tests for equality of variance and mean were applied to SNP intensities, and significance p-values were obtained. A Benjamini-Hochberg multiple testing correction reveled significant SNPs. Corresponding genes were identified using the Bovine Genome UMD 3.1 annotation. Principal Components Analysis (PCA) confirmed differentiation. An analysis of Copy Number Variations (CNVs), obtained from signal intensities, revealed no evidence of association between ovarian pathologies and CNVs. In addition, a haplotype frequency analysis showed no association with ovarian pathologies. Results show that SNP signal intensity, which captures not only information for base-pair genotypes elucidation, but the amount of fluorescence nucleotide synthetization produced in an enzymatic reaction, is a rich source of information that, by itself or in combination with base-pair genotypes, might be used to implement differentiation, prediction and diagnostic procedures, increasing the scope of applications for Genotyping Microarrays.
Collapse
|
37
|
Zhang B, Chamba Y, Shang P, Wang Z, Ma J, Wang L, Zhang H. Comparative transcriptomic and proteomic analyses provide insights into the key genes involved in high-altitude adaptation in the Tibetan pig. Sci Rep 2017. [PMID: 28623314 PMCID: PMC5473931 DOI: 10.1038/s41598-017-03976-3] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Tibetan pigs that inhabit the Tibetan Plateau exhibit striking phenotypic and physiological differences from lowland pigs, and have adapted well to extreme conditions. However, the mechanisms involved in regulating gene expression at high altitude in these animals are not fully understood. In this study, we obtained transcriptomic and proteomic data from the heart tissues of Tibetan and Yorkshire pigs raised in the highlands (TH and YH) and lowlands (TL and YL) via RNA-seq and iTRAQ (isobaric tags for relative and absolute quantitation) analyses, respectively. Comparative analyses of TH vs. YH, TH vs.TL, TL vs. YL, and YH vs. YL yielded 299, 169, 242, and 368 differentially expressed genes (DEGs), and 473, 297, 394, and 297 differentially expressed proteins (DEPs), respectively. By functional annotation of these DEGs and DEPs, genes that were enriched in the HIF-1 signaling pathway (NPPA, ERK2, ENO3, and EGLN3), VEGF signaling pathway (ERK2, A2M, FGF1, CTGF, and DPP4), and hypoxia-related processes (CRYAB, EGLN3, TGFB2, DPP4, and ACE) were identified as important candidate genes for high-altitude adaptation in the Tibetan pig. This study enhances our understanding of the molecular mechanisms involved in hypoxic adaptation in pigs, and furthers our understanding of human hypoxic diseases.
Collapse
Affiliation(s)
- Bo Zhang
- National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Rd., Beijing, 100193, China
| | - Yangzom Chamba
- Tibet Agriculture and Animal Husbandry College, Linzhi, Tibet, 860000, China
| | - Peng Shang
- National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Rd., Beijing, 100193, China.,Tibet Agriculture and Animal Husbandry College, Linzhi, Tibet, 860000, China
| | - Zhixiu Wang
- National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Rd., Beijing, 100193, China
| | - Jun Ma
- National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Rd., Beijing, 100193, China
| | - Liyuang Wang
- National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Rd., Beijing, 100193, China
| | - Hao Zhang
- National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Rd., Beijing, 100193, China.
| |
Collapse
|
38
|
In vivo and in vitro ageing results in accumulation of de novo copy number variations in bulls. Sci Rep 2017; 7:1631. [PMID: 28487564 PMCID: PMC5431667 DOI: 10.1038/s41598-017-01793-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2016] [Accepted: 04/03/2017] [Indexed: 11/24/2022] Open
Abstract
We have identified de novo copy number variations (CNVs) generated in bulls as they age. Blood samples from eight bulls were collected and SNP arrayed in a prospective design over 30 months allowing us to differentiate de novo CNVs from constant CNVs that are present throughout the sampling period. Quite remarkably, the total number of CNVs doubled over the 30-month period, as we observed an almost equal number of de novo and constant CNVs (107 and 111, respectively, i.e. 49% and 51%). Twice as many de novo CNVs emerged during the second half of the sampling schedule as in the first part. It suggests a dynamic generation of de novo CNVs in the bovine genome that becomes more frequent as the age of the animal progresses. In a second experiment de novo CNVs were detected through in vitro ageing of bovine fibroblasts by sampling passage #5, #15 and #25. De novo CNVs also became more frequent, but the proportion of them was only ~25% of the total number of CNVs (21 out of 85). Temporal generation of de novo CNVs resulted in increasing genome coverage. Genes and quantitative trait loci overlapping de novo CNVs were further investigated for ageing related functions.
Collapse
|
39
|
Gao Y, Jiang J, Yang S, Hou Y, Liu GE, Zhang S, Zhang Q, Sun D. CNV discovery for milk composition traits in dairy cattle using whole genome resequencing. BMC Genomics 2017; 18:265. [PMID: 28356085 PMCID: PMC5371188 DOI: 10.1186/s12864-017-3636-3] [Citation(s) in RCA: 55] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2016] [Accepted: 03/17/2017] [Indexed: 01/08/2023] Open
Abstract
Background Copy number variations (CNVs) are important and widely distributed in the genome. CNV detection opens a new avenue for exploring genes associated with complex traits in humans, animals and plants. Herein, we present a genome-wide assessment of CNVs that are potentially associated with milk composition traits in dairy cattle. Results In this study, CNVs were detected based on whole genome re-sequencing data of eight Holstein bulls from four half- and/or full-sib families, with extremely high and low estimated breeding values (EBVs) of milk protein percentage and fat percentage. The range of coverage depth per individual was 8.2–11.9×. Using CNVnator, we identified a total of 14,821 CNVs, including 5025 duplications and 9796 deletions. Among them, 487 differential CNV regions (CNVRs) comprising ~8.23 Mb of the cattle genome were observed between the high and low groups. Annotation of these differential CNVRs were performed based on the cattle genome reference assembly (UMD3.1) and totally 235 functional genes were found within the CNVRs. By Gene Ontology and KEGG pathway analyses, we found that genes were significantly enriched for specific biological functions related to protein and lipid metabolism, insulin/IGF pathway-protein kinase B signaling cascade, prolactin signaling pathway and AMPK signaling pathways. These genes included INS, IGF2, FOXO3, TH, SCD5, GALNT18, GALNT16, ART3, SNCA and WNT7A, implying their potential association with milk protein and fat traits. In addition, 95 CNVRs were overlapped with 75 known QTLs that are associated with milk protein and fat traits of dairy cattle (Cattle QTLdb). Conclusions In conclusion, based on NGS of 8 Holstein bulls with extremely high and low EBVs for milk PP and FP, we identified a total of 14,821 CNVs, 487 differential CNVRs between groups, and 10 genes, which were suggested as promising candidate genes for milk protein and fat traits. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3636-3) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Yahui Gao
- Key Laboratory of Animal Genetics and Breeding of Ministry of Agriculture, National Engineering Laboratory of Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Jianping Jiang
- Key Laboratory of Animal Genetics and Breeding of Ministry of Agriculture, National Engineering Laboratory of Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Shaohua Yang
- Key Laboratory of Animal Genetics and Breeding of Ministry of Agriculture, National Engineering Laboratory of Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Yali Hou
- CAS Key Laboratory of Genomic and Precision Medicine, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, 100101, China
| | - George E Liu
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, Md, 20705, USA
| | - Shengli Zhang
- Key Laboratory of Animal Genetics and Breeding of Ministry of Agriculture, National Engineering Laboratory of Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Qin Zhang
- Key Laboratory of Animal Genetics and Breeding of Ministry of Agriculture, National Engineering Laboratory of Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Dongxiao Sun
- Key Laboratory of Animal Genetics and Breeding of Ministry of Agriculture, National Engineering Laboratory of Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China.
| |
Collapse
|
40
|
Chen L, Chamberlain AJ, Reich CM, Daetwyler HD, Hayes BJ. Detection and validation of structural variations in bovine whole-genome sequence data. Genet Sel Evol 2017; 49:13. [PMID: 28122487 PMCID: PMC5267451 DOI: 10.1186/s12711-017-0286-5] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2016] [Accepted: 01/09/2017] [Indexed: 02/05/2023] Open
Abstract
BACKGROUND Several examples of structural variation (SV) affecting phenotypic traits have been reported in cattle. Currently the identification of SV from whole-genome sequence data (WGS) suffers from a high false positive rate. Our aim was to construct a high quality set of SV calls in cattle using WGS data. First, we tested two SV detection programs, Breakdancer and Pindel, and the overlap of these methods, on simulated sequence data to determine their precision and sensitivity. We then identified population SV from WGS of 252 Holstein and 64 Jersey bulls based on the overlapping calls from the two programs. In addition, we validated an overlapped SV set in 28 twice-sequenced Holstein individuals, and in another two validated sets (one for each breed) that were transmitted from sire to son. We also tested whether highly conserved gene sets across eukaryotes and recently expanded gene families in bovine were depleted and enriched, respectively, for SV. RESULTS In empirical WGS data, 17,518 SV covering 27.36 Mb were found in the Holstein population and 4285 SV covering 8.74 Mb in the Jersey population, of which 4.62 Mb of SV overlapped between Holsteins and Jerseys. A total of 11,534 candidate SV covering 5.64 Mb were validated in the 28 twice-sequenced individuals, while 3.49 and 0.67 Mb of SV were validated from Holstein and Jersey sire-son transmission, respectively. Only eight of 237 core eukaryotic genes had at least a 50-bp overlap with an SV from our validated sets, suggesting that conserved genes are depleted for SV (p < 0.05). In addition, we observed that recently expanded gene families were significantly more associated with SV than other genes. Long interspersed nuclear elements-1 were enriched for deletions when compared to the rest of the genome (p = 0.0035). CONCLUSIONS We reported SV from 252 Holstein and 64 Jersey individuals. A considerable proportion of Jersey population SV (53.5%) were also found in Holstein. In contrast, about 76.90% sire-son transmission validated SV were present in Jerseys and Holsteins. The enrichment of SV in expanding gene families suggests that SV can be a source of genetic variation for evolution.
Collapse
Affiliation(s)
- Long Chen
- AgriBio, Centre for AgriBioscience, Biosciences Research, Department of Economic Development, Jobs, Transport and Resources, Bundoora, VIC, Australia. .,School of Applied Systems Biology, La Trobe University, Bundoora, VIC, Australia.
| | - Amanda J Chamberlain
- AgriBio, Centre for AgriBioscience, Biosciences Research, Department of Economic Development, Jobs, Transport and Resources, Bundoora, VIC, Australia
| | - Coralie M Reich
- AgriBio, Centre for AgriBioscience, Biosciences Research, Department of Economic Development, Jobs, Transport and Resources, Bundoora, VIC, Australia
| | - Hans D Daetwyler
- AgriBio, Centre for AgriBioscience, Biosciences Research, Department of Economic Development, Jobs, Transport and Resources, Bundoora, VIC, Australia.,School of Applied Systems Biology, La Trobe University, Bundoora, VIC, Australia
| | - Ben J Hayes
- AgriBio, Centre for AgriBioscience, Biosciences Research, Department of Economic Development, Jobs, Transport and Resources, Bundoora, VIC, Australia.,School of Applied Systems Biology, La Trobe University, Bundoora, VIC, Australia
| |
Collapse
|
41
|
Keel BN, Lindholm-Perry AK, Snelling WM. Evolutionary and Functional Features of Copy Number Variation in the Cattle Genome. Front Genet 2016; 7:207. [PMID: 27920798 PMCID: PMC5118444 DOI: 10.3389/fgene.2016.00207] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2016] [Accepted: 11/08/2016] [Indexed: 01/18/2023] Open
Abstract
Genomic structural variations are an important source of genetic diversity. Copy number variations (CNVs), gains and losses of large regions of genomic sequence between individuals of a species, have been associated with a wide variety of phenotypic traits. However, in cattle, as well as many other species, relatively little is understood about CNV, including frequency of CNVs in the genome, sizes, and locations, chromosomal properties, and evolutionary processes acting to shape CNV. In this work, we focused on copy number variation in the bovine genome, with the aim to detect CNVs in Bos taurus coding sequence and explore potential evolutionary mechanisms shaping these CNV. We identified and characterized CNV regions by utilizing exome sequence from 175 influential sires used in the Germplasm Evaluation project, representing 10 breeds. We examined various evolutionary and functional aspects of these CNVs, including selective constraint on CNV-overlapped genes, centrality of CNV genes in protein-protein interaction networks, and tissue-specific expression of CNV genes. Patterns of CNV in the Bos taurus genome reveal that reduced functional constraint and mutational bias may play a prominent role in shaping this type of structural variation.
Collapse
Affiliation(s)
- Brittney N Keel
- Agricultural Research Service (USDA), Meat Animal Research Center Clay Center, NE, USA
| | | | - Warren M Snelling
- Agricultural Research Service (USDA), Meat Animal Research Center Clay Center, NE, USA
| |
Collapse
|
42
|
Keel BN, Keele JW, Snelling WM. Genome-wide copy number variation in the bovine genome detected using low coverage sequence of popular beef breeds,. Anim Genet 2016; 48:141-150. [DOI: 10.1111/age.12519] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/27/2016] [Indexed: 12/19/2022]
Affiliation(s)
- B. N. Keel
- USDA; ARS; U.S. Meat Animal Research Center; Clay Center NE 68933 USA
| | - J. W. Keele
- USDA; ARS; U.S. Meat Animal Research Center; Clay Center NE 68933 USA
| | - W. M. Snelling
- USDA; ARS; U.S. Meat Animal Research Center; Clay Center NE 68933 USA
| |
Collapse
|
43
|
Durán Aguilar M, Román Ponce SI, Ruiz López FJ, González Padilla E, Vásquez Peláez CG, Bagnato A, Strillacci MG. Genome-wide association study for milk somatic cell score in holstein cattle using copy number variation as markers. J Anim Breed Genet 2016; 134:49-59. [PMID: 27578198 DOI: 10.1111/jbg.12238] [Citation(s) in RCA: 43] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2016] [Accepted: 08/03/2016] [Indexed: 12/21/2022]
Abstract
Mastitis, the most common and expensive disease in dairy cows, implies significant losses in the dairy industry worldwide. Many efforts have been made to improve genetic mastitis resistance in dairy populations, but low heritability of this trait made this process not as effective as desired. The purpose of this study was to identify genomic regions explaining genetic variation of somatic cell count using copy number variations (CNVs) as markers in the Holstein population, genotyped with the Illumina BovineHD BeadChip. We found 24 and 47 copy number variation regions significantly associated with estimated breeding values for somatic cell score (SCS_EBVs) using SVS 8.3.1 and PennCNV-CNVRuler software, respectively. The association analysis performed with these two software allowed the identification of 18 candidate genes (TERT, NOTCH1, SLC6A3, CLPTM1L, PPARα, BCL-2, ABO, VAV2, CACNA1S, TRAF2, RELA, ELF3, DBH, CDK5, NF2, FASN, EWSR1 and MAP3K11) that result classified in the same functional cluster. These genes are also part of two gene networks, whose genes share the 'stress', 'cell death', 'inflammation' and 'immune response' GO terms. Combining CNV detection/association analysis based on two different algorithms helps towards a more complete identification of genes linked to phenotypic variation of the somatic cell count.
Collapse
Affiliation(s)
- M Durán Aguilar
- Facultad de Estudios Superiores Cuautitlán, UNAM, Cuautitlán Izcalli, México
| | - S I Román Ponce
- Centro Nacional de Investigación en Fisiología y Mejoramiento Animal, INIFAP, Auchitlán, Querétaro, México
| | - F J Ruiz López
- Centro Nacional de Investigación en Fisiología y Mejoramiento Animal, INIFAP, Auchitlán, Querétaro, México
| | - E González Padilla
- Departamento de Genética y Bioestadística, Facultad Medicina Veterinaria y Zootecnia, Universidad Nacional Autónoma de México, México DF, México
| | - C G Vásquez Peláez
- Departamento de Genética y Bioestadística, Facultad Medicina Veterinaria y Zootecnia, Universidad Nacional Autónoma de México, México DF, México
| | - A Bagnato
- Department of Veterinary Medicine (DiMeVet), University of Milan, Milan, Italy
| | - M G Strillacci
- Department of Veterinary Medicine (DiMeVet), University of Milan, Milan, Italy
| |
Collapse
|
44
|
Zvinorova PI, Halimani TE, Muchadeyi FC, Matika O, Riggio V, Dzama K. Breeding for resistance to gastrointestinal nematodes - the potential in low-input/output small ruminant production systems. Vet Parasitol 2016; 225:19-28. [PMID: 27369571 PMCID: PMC4938797 DOI: 10.1016/j.vetpar.2016.05.015] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2015] [Revised: 04/28/2016] [Accepted: 05/12/2016] [Indexed: 02/07/2023]
Abstract
The control of gastrointestinal nematodes (GIN) is mainly based on the use of drugs, grazing management, use of copper oxide wire particles and bioactive forages. Resistance to anthelmintic drugs in small ruminants is documented worldwide. Host genetic resistance to parasites, has been increasingly used as a complementary control strategy, along with the conventional intervention methods mentioned above. Genetic diversity in resistance to GIN has been well studied in experimental and commercial flocks in temperate climates and more developed economies. However, there are very few report outputs from the more extensive low-input/output smallholder systems in developing and emerging countries. Furthermore, results on quantitative trait loci (QTL) associated with nematode resistance from various studies have not always been consistent, mainly due to the different nematodes studied, different host breeds, ages, climates, natural infections versus artificial challenges, infection level at sampling periods, among others. The increasing use of genetic markers (Single Nucleotide Polymorphisms, SNPs) in GWAS or the use of whole genome sequence data and a plethora of analytic methods offer the potential to identify loci or regions associated nematode resistance. Genomic selection as a genome-wide level method overcomes the need to identify candidate genes. Benefits in genomic selection are now being realised in dairy cattle and sheep under commercial settings in the more advanced countries. However, despite the commercial benefits of using these tools, there are practical problems associated with incorporating the use of marker-assisted selection or genomic selection in low-input/output smallholder farming systems breeding schemes. Unlike anthelmintic resistance, there is no empirical evidence suggesting that nematodes will evolve rapidly in response to resistant hosts. The strategy of nematode control has evolved to a more practical manipulation of host-parasite equilibrium in grazing systems by implementation of various strategies, in which improvement of genetic resistance of small ruminant should be included. Therefore, selection for resistant hosts can be considered as one of the sustainable control strategy, although it will be most effective when used to complement other control strategies such as grazing management and improving efficiency of anthelmintics currently.
Collapse
Affiliation(s)
- P I Zvinorova
- Department of Animal Sciences, University of Stellenbosch, Private Bag X1, Matieland, 7602, South Africa; Department of Para-clinical Veterinary Studies, University of Zimbabwe, P.O. MP167, Mt. Pleasant, Harare, Zimbabwe.
| | - T E Halimani
- Department of Animal Science, University of Zimbabwe, P.O. MP167, Mt. Pleasant, Harare, Zimbabwe.
| | - F C Muchadeyi
- Biotechnology Platform, Agriculture Research Council Private Bag X5, Onderstepoort, 0110, South Africa.
| | - O Matika
- The Roslin Institute and R(D)SVS, University of Edinburgh, Easter Bush, MidlothianEH25 9RG, UK.
| | - V Riggio
- The Roslin Institute and R(D)SVS, University of Edinburgh, Easter Bush, MidlothianEH25 9RG, UK.
| | - K Dzama
- Department of Animal Sciences, University of Stellenbosch, Private Bag X1, Matieland, 7602, South Africa.
| |
Collapse
|
45
|
da Silva VH, Regitano LCDA, Geistlinger L, Pértille F, Giachetto PF, Brassaloti RA, Morosini NS, Zimmer R, Coutinho LL. Genome-Wide Detection of CNVs and Their Association with Meat Tenderness in Nelore Cattle. PLoS One 2016; 11:e0157711. [PMID: 27348523 PMCID: PMC4922624 DOI: 10.1371/journal.pone.0157711] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2015] [Accepted: 06/03/2016] [Indexed: 12/20/2022] Open
Abstract
Brazil is one of the largest beef producers and exporters in the world with the Nelore breed representing the vast majority of Brazilian cattle (Bos taurus indicus). Despite the great adaptability of the Nelore breed to tropical climate, meat tenderness (MT) remains to be improved. Several factors including genetic composition can influence MT. In this article, we report a genome-wide analysis of copy number variation (CNV) inferred from Illumina® High Density SNP-chip data for a Nelore population of 723 males. We detected >2,600 CNV regions (CNVRs) representing ≈6.5% of the genome. Comparing our results with previous studies revealed an overlap in ≈1400 CNVRs (>50%). A total of 1,155 CNVRs (43.6%) overlapped 2,750 genes. They were enriched for processes involving guanosine triphosphate (GTP), previously reported to influence skeletal muscle physiology and morphology. Nelore CNVRs also overlapped QTLs for MT reported in other breeds (8.9%, 236 CNVRs) and from a previous study with this population (4.1%, 109 CNVRs). Two CNVRs were also proximal to glutathione metabolism genes that were previously associated with MT. Genome-wide association study of CN state with estimated breeding values derived from meat shear force identified 6 regions, including a region on BTA3 that contains genes of the cAMP and cGMP pathway. Ten CNVRs that overlapped regions associated with MT were successfully validated by qPCR. Our results represent the first comprehensive CNV study in Bos taurus indicus cattle and identify regions in which copy number changes are potentially of importance for the MT phenotype.
Collapse
Affiliation(s)
- Vinicius Henrique da Silva
- Animal Biotechnology Laboratory, Animal Science Department, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo, Brazil
- * E-mail: (LLC); (VHS)
| | | | - Ludwig Geistlinger
- Institute of Bioinformatics, Department of Informatics, Ludwig-Maximilians-Universität München (LMU), Amalienstrasse 17, 80333, München, Germany
| | - Fábio Pértille
- Animal Biotechnology Laboratory, Animal Science Department, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo, Brazil
| | | | - Ricardo Augusto Brassaloti
- Animal Biotechnology Laboratory, Animal Science Department, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo, Brazil
| | - Natália Silva Morosini
- Animal Biotechnology Laboratory, Animal Science Department, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo, Brazil
| | - Ralf Zimmer
- Institute of Bioinformatics, Department of Informatics, Ludwig-Maximilians-Universität München (LMU), Amalienstrasse 17, 80333, München, Germany
| | - Luiz Lehmann Coutinho
- Animal Biotechnology Laboratory, Animal Science Department, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo, Brazil
- * E-mail: (LLC); (VHS)
| |
Collapse
|
46
|
da Silva JM, Giachetto PF, da Silva LO, Cintra LC, Paiva SR, Yamagishi MEB, Caetano AR. Genome-wide copy number variation (CNV) detection in Nelore cattle reveals highly frequent variants in genome regions harboring QTLs affecting production traits. BMC Genomics 2016; 17:454. [PMID: 27297173 PMCID: PMC4907077 DOI: 10.1186/s12864-016-2752-9] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2016] [Accepted: 05/19/2016] [Indexed: 11/10/2022] Open
Abstract
Background Copy number variations (CNVs) have been shown to account for substantial portions of observed genomic variation and have been associated with qualitative and quantitative traits and the onset of disease in a number of species. Information from high-resolution studies to detect, characterize and estimate population-specific variant frequencies will facilitate the incorporation of CNVs in genomic studies to identify genes affecting traits of importance. Results Genome-wide CNVs were detected in high-density single nucleotide polymorphism (SNP) genotyping data from 1,717 Nelore (Bos indicus) cattle, and in NGS data from eight key ancestral bulls. A total of 68,007 and 12,786 distinct CNVs were observed, respectively. Cross-comparisons of results obtained for the eight resequenced animals revealed that 92 % of the CNVs were observed in both datasets, while 62 % of all detected CNVs were observed to overlap with previously validated cattle copy number variant regions (CNVRs). Observed CNVs were used for obtaining breed-specific CNV frequencies and identification of CNVRs, which were subsequently used for gene annotation. A total of 688 of the detected CNVRs were observed to overlap with 286 non-redundant QTLs associated with important production traits in cattle. All of 34 CNVs previously reported to be associated with milk production traits in Holsteins were also observed in Nelore cattle. Comparisons of estimated frequencies of these CNVs in the two breeds revealed 14, 13, 6 and 14 regions in high (>20 %), low (<20 %) and divergent (NEL > HOL, NEL < HOL) frequencies, respectively. Conclusions Obtained results significantly enriched the bovine CNV map and enabled the identification of variants that are potentially associated with traits under selection in Nelore cattle, particularly in genome regions harboring QTLs affecting production traits. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2752-9) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Joaquim Manoel da Silva
- Faculdade de Ciências Agrárias, Biológicas e Sociais Aplicadas, Universidade do Estado de Mato Grosso (UNEMAT), Av. Prof Dr. Renato Figueiro Varella, CEP 78.690-000, Nova Xavantina, Mato Grosso, Brazil.,Programa de Pós-Graduação em Genética e Biologia Molecular-Instituto de Biologia, Universidade Estadual de Campinas (UNICAMP), Campinas, São Paulo, Brazil
| | - Poliana Fernanda Giachetto
- Embrapa Informática Agropecuária - Laboratório Multiusuário de Bioinformática (LMB), Campinas, São Paulo, Brazil
| | | | - Leandro Carrijo Cintra
- Embrapa Informática Agropecuária - Laboratório Multiusuário de Bioinformática (LMB), Campinas, São Paulo, Brazil
| | - Samuel Rezende Paiva
- Embrapa - Secretaria de Relações Internacionais, Brasília, Distrito Federal, Brazil.,Embrapa Recursos Genéticos e Biotecnologia, Brasília, Distrito Federal, Brazil.,CNPq Fellow, ᅟ, ᅟ
| | | | | |
Collapse
|
47
|
Zhou Y, Utsunomiya YT, Xu L, Hay EHA, Bickhart DM, Alexandre PA, Rosen BD, Schroeder SG, Carvalheiro R, de Rezende Neves HH, Sonstegard TS, Van Tassell CP, Ferraz JBS, Fukumasu H, Garcia JF, Liu GE. Genome-wide CNV analysis reveals variants associated with growth traits in Bos indicus. BMC Genomics 2016; 17:419. [PMID: 27245577 PMCID: PMC4888316 DOI: 10.1186/s12864-016-2461-4] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2015] [Accepted: 02/11/2016] [Indexed: 12/19/2022] Open
Abstract
BACKGROUND Apart from single nucleotide polymorphism (SNP), copy number variation (CNV) is another important type of genetic variation, which may affect growth traits and play key roles for the production of beef cattle. To date, no genome-wide association study (GWAS) for CNV and body traits in beef cattle has been reported, so the present study aimed to investigate this type of association in one of the most important cattle subspecies: Bos indicus (Nellore breed). RESULTS We have used intensity data from over 700,000 SNP probes across the bovine genome to detect common CNVs in a sample of 2230 Nellore cattle, and performed GWAS between the detected CNVs and nine growth traits. After filtering for frequency and length, a total of 231 CNVs ranging from 894 bp to 4,855,088 bp were kept and tested as predictors for each growth trait using linear regression analysis with principal components correction. There were 49 significant associations identified among 17 CNVs and seven body traits after false discovery rate correction (P < 0.05). Among the 17 CNVs, three were significant or marginally significant for all the traits. We have compared the locations of associated CNVs with quantitative trait locus and the RefGene database, and found two sets of 9 CNVs overlapping with either known QTLs or genes, respectively. The gene overlapping with CNV100, KCNJ12, is a functional candidate for muscle development and plays critical roles in muscling traits. CONCLUSION This study presents the first CNV-based GWAS of growth traits using high density SNP microarray data in cattle. We detected 17 CNVs significantly associated with seven growth traits and one of them (CNV100) may be involved in growth traits through KCNJ12.
Collapse
Affiliation(s)
- Yang Zhou
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Room 111, BARC-East, Beltsville, Maryland, 20705, USA.,College of Animal Science and Technology, Northwest A&F University, Shaanxi Key Laboratory of Agricultural Molecular Biology, Yangling, Shaanxi, 712100, China
| | - Yuri T Utsunomiya
- Departamento de Medicina Veterinária Preventiva e Reprodução Animal, Faculdade de Ciências Agrárias e Veterinárias, UNESP - Univ Estadual Paulista, Jaboticabal, São Paulo, 14884-900, Brazil
| | - Lingyang Xu
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Room 111, BARC-East, Beltsville, Maryland, 20705, USA.,Institute of Animal Science, Chinese Academy of Agricultural Science, Beijing, 100193, China
| | - El Hamidi Abdel Hay
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Room 111, BARC-East, Beltsville, Maryland, 20705, USA
| | - Derek M Bickhart
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Room 111, BARC-East, Beltsville, Maryland, 20705, USA
| | - Pamela Almeida Alexandre
- Department of Veterinary Medicine, College of Animal Science and Food Engineering, University of São Paulo, Pirassununga, SP, 13635, Brazil
| | - Benjamin D Rosen
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Room 111, BARC-East, Beltsville, Maryland, 20705, USA
| | - Steven G Schroeder
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Room 111, BARC-East, Beltsville, Maryland, 20705, USA
| | - Roberto Carvalheiro
- Departamento de Zootecnia, Faculdade de Ciências Agrárias e Veterinárias, UNESP - Univ Estadual Paulista, Jaboticabal, São Paulo, 14884-900, Brazil
| | - Haroldo Henrique de Rezende Neves
- Departamento de Zootecnia, Faculdade de Ciências Agrárias e Veterinárias, UNESP - Univ Estadual Paulista, Jaboticabal, São Paulo, 14884-900, Brazil
| | - Tad S Sonstegard
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Room 111, BARC-East, Beltsville, Maryland, 20705, USA.,Present address: Recombinetics, Inc., St Paul, MN, 55104, USA
| | - Curtis P Van Tassell
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Room 111, BARC-East, Beltsville, Maryland, 20705, USA
| | - José Bento Sterman Ferraz
- Department of Veterinary Medicine, College of Animal Science and Food Engineering, University of São Paulo, Pirassununga, SP, 13635, Brazil
| | - Heidge Fukumasu
- Department of Veterinary Medicine, College of Animal Science and Food Engineering, University of São Paulo, Pirassununga, SP, 13635, Brazil
| | - Jose Fernando Garcia
- Departamento de Medicina Veterinária Preventiva e Reprodução Animal, Faculdade de Ciências Agrárias e Veterinárias, UNESP - Univ Estadual Paulista, Jaboticabal, São Paulo, 14884-900, Brazil. .,Departamento de Apoio, Produção e Saúde Animal, Faculdade de Medicina Veterinária de Araçatuba, UNESP - Univ Estadual Paulista, Araçatuba, São Paulo, 16050-680, Brazil. .,International Atomic Energy Agency (IAEA) Collaborating Centre on Animal Genomics and Bioinformatics, Araçatuba, SP, Brazil.
| | - George E Liu
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Room 111, BARC-East, Beltsville, Maryland, 20705, USA.
| |
Collapse
|
48
|
Sasaki S, Watanabe T, Nishimura S, Sugimoto Y. Genome-wide identification of copy number variation using high-density single-nucleotide polymorphism array in Japanese Black cattle. BMC Genet 2016; 17:26. [PMID: 26809925 PMCID: PMC4727303 DOI: 10.1186/s12863-016-0335-z] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2015] [Accepted: 01/14/2016] [Indexed: 12/12/2022] Open
Abstract
Background Copy number variation (CNV) is an important source of genetic variability associated with phenotypic variation and disease susceptibility. Comprehensive genome-wide CNV maps provide valuable information for genetic and functional studies. To identify CNV in Japanese Black cattle, we performed a genome-wide autosomal screen using genomic data from 1,481 animals analyzed with the Illumina Bovine High-Density (HD) BeadChip Array (735,293 single-nucleotide polymorphisms (SNPs) with an average marker interval of 3.4 kb on the autosomes). Results We identified a total of 861 CNV regions (CNVRs) across all autosomes, which covered 43.65 Mb of the UMD3.1 genome assembly and corresponded to 1.74 % of the 29 bovine autosomes. Overall, 35 % of the CNVRs were present at a frequency of > 1 % in 1,481 animals. The estimated lengths of CNVRs ranged from 1.1 kb to 1.4 Mb, with an average of 50.7 kb. The average number of CNVR events per animal was 35. Comparisons with previously reported cattle CNV showed that 72 % of the CNVR calls detected in this study were within or overlapped with known CNVRs. Experimentally, three CNVRs were validated using quantitative PCR, and one CNVR was validated using PCR with flanking primers for the deleted region. Out of the 861 CNVRs, 390 contained 717 Ensembl-annotated genes significantly enriched for stimulus response, cellular defense response, and immune response in the Gene Ontology (GO) database. To associate genes contained in CNVRs with phenotypes, we converted 560 bovine Ensembl gene IDs to their 438 orthologous associated mouse gene IDs, and 195 of these mouse orthologous genes were categorized into 1,627 phenotypes in the Mouse Genome Informatics (MGI) database. Conclusions We identified 861 CNVRs in 1,481 Japanese Black cattle using the Illumina BovineHD BeadChip Array. The genes contained in CNVRs were characterized using GO analysis and the mouse orthologous genes were characterized using the MGI database. The comprehensive genome-wide CNVRs map will facilitate identification of genetic variation and disease-susceptibility alleles in Japanese Black cattle. Electronic supplementary material The online version of this article (doi:10.1186/s12863-016-0335-z) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Shinji Sasaki
- National Livestock Breeding Center, Odakura, Nishigo, Fukushima, 961-8511, Japan.
| | - Toshio Watanabe
- National Livestock Breeding Center, Odakura, Nishigo, Fukushima, 961-8511, Japan.
| | - Shota Nishimura
- Shirakawa Institute of Animal Genetics, Japan Livestock Technology Association, Odakura, Nishigo, Fukushima, 961-8061, Japan.
| | - Yoshikazu Sugimoto
- Shirakawa Institute of Animal Genetics, Japan Livestock Technology Association, Odakura, Nishigo, Fukushima, 961-8061, Japan.
| |
Collapse
|
49
|
Wang MD, Dzama K, Rees DJG, Muchadeyi FC. Tropically adapted cattle of Africa: perspectives on potential role of copy number variations. Anim Genet 2015; 47:154-64. [DOI: 10.1111/age.12391] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/01/2015] [Indexed: 12/12/2022]
Affiliation(s)
- M. D. Wang
- Department of Animal Sciences; University of Stellenbosch; Private Bag X1 Matieland 7602 South Africa
- Biotechnology Platform; Agricultural Research Council; Private Bag X5 Onderstepoort 0110 South Africa
| | - K. Dzama
- Department of Animal Sciences; University of Stellenbosch; Private Bag X1 Matieland 7602 South Africa
| | - D. J. G. Rees
- Biotechnology Platform; Agricultural Research Council; Private Bag X5 Onderstepoort 0110 South Africa
| | - F. C. Muchadeyi
- Biotechnology Platform; Agricultural Research Council; Private Bag X5 Onderstepoort 0110 South Africa
| |
Collapse
|
50
|
Wang MD, Dzama K, Hefer CA, Muchadeyi FC. Genomic population structure and prevalence of copy number variations in South African Nguni cattle. BMC Genomics 2015; 16:894. [PMID: 26531252 PMCID: PMC4632335 DOI: 10.1186/s12864-015-2122-z] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2015] [Accepted: 10/22/2015] [Indexed: 12/21/2022] Open
Abstract
Background Copy number variations (CNVs) are modifications in DNA structure comprising of deletions, duplications, insertions and complex multi-site variants. Although CNVs are proven to be involved in a variety of phenotypic discrepancies, the full extent and consequence of CNVs is yet to be understood. To date, no such genomic characterization has been performed in indigenous South African Nguni cattle. Nguni cattle are recognized for their ability to sustain harsh environmental conditions while exhibiting enhanced resistance to disease and parasites and are thought to comprise of up to nine different ecotypes. Methods Illumina BovineSNP50 Beadchip data was utilized to investigate genomic population structure and the prevalence of CNVs in 492 South African Nguni cattle. PLINK, ADMIXTURE, R, gPLINK and Haploview software was utilized for quality control, population structure and haplotype block determination. PennCNV hidden Markov model identified CNVs and genes contained within and 10 Mb downstream from reported CNVs. PANTHER and Ensembl databases were subsequently utilized for gene annotation analyses. Results Population structure analyses on Nguni cattle revealed 5 sub-populations with a possible sub-structure evident at K equal to 8. Four hundred and thirty three CNVs that formed 334 CNVRs ranging from 30 kb to 1 Mb in size are reported. Only 231 of the 492 animals demonstrated CNVRs. Two hundred and eighty nine genes were observed within CNVRs identified. Of these 149, 28, 44, 2 and 14 genes were unique to sub-populations A, B, C, D and E respectively. Gene ontology analyses demonstrated a number of pathways to be represented by respective genes, including immune response, response to abiotic stress and biological regulation processess. Conclusions CNVs may explain part of the phenotypic diversity and the enhanced adaptation evident in Nguni cattle. Genes involved in a number of cellular components, biological processes and molecular functions are reported within CNVRs identified. The significance of such CNVRs and the possible effect thereof needs to be ascertained and may hold interesting insight into the functional and adaptive consequence of CNVs in cattle. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2122-z) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Magretha Diane Wang
- Department of Animal Sciences, University of Stellenbosch, Private Bag X1, Matieland, Stellenbosch, 7602, South Africa. .,Biotechnology Platform, Agricultural Research Council, Private Bag X5, Onderstepoort, 0110, South Africa.
| | - Kennedy Dzama
- Department of Animal Sciences, University of Stellenbosch, Private Bag X1, Matieland, Stellenbosch, 7602, South Africa.
| | - Charles A Hefer
- Biotechnology Platform, Agricultural Research Council, Private Bag X5, Onderstepoort, 0110, South Africa.
| | - Farai C Muchadeyi
- Biotechnology Platform, Agricultural Research Council, Private Bag X5, Onderstepoort, 0110, South Africa.
| |
Collapse
|