Whole-genome Sequencing Reveals Mutational Landscape Underlying Phenotypic Differences Between Two Widespread Chinese Cattle Breeds

Overview

Journal PLoS One

Specialties General Medicine
Science

Date 2017 Aug 26

PMID 28841720

Citations 14

Authors

Yao Xu

Yu Jiang

Tao Shi

Hanfang Cai

Xianyong Lan

Xin Zhao

Martin Plath

Hong Chen

Affiliations

Soon will be listed here.

Abstract

Whole-genome sequencing provides a powerful tool to obtain more genetic variability that could produce a range of benefits for cattle breeding industry. Nanyang (Bos indicus) and Qinchuan (Bos taurus) are two important Chinese indigenous cattle breeds with distinct phenotypes. To identify the genetic characteristics responsible for variation in phenotypes between the two breeds, in the present study, we for the first time sequenced the genomes of four Nanyang and four Qinchuan cattle with 10 to 12 fold on average of 97.86% and 98.98% coverage of genomes, respectively. Comparison with the Bos_taurus_UMD_3.1 reference assembly yielded 9,010,096 SNPs for Nanyang, and 6,965,062 for Qinchuan cattle, 51% and 29% of which were novel SNPs, respectively. A total of 154,934 and 115,032 small indels (1 to 3 bp) were found in the Nanyang and Qinchuan genomes, respectively. The SNP and indel distribution revealed that Nanyang showed a genetically high diversity as compared to Qinchuan cattle. Furthermore, a total of 2,907 putative cases of copy number variation (CNV) were identified by aligning Nanyang to Qinchuan genome, 783 of which (27%) encompassed the coding regions of 495 functional genes. The gene ontology (GO) analysis revealed that many CNV genes were enriched in the immune system and environment adaptability. Among several CNV genes related to lipid transport and fat metabolism, Lepin receptor gene (LEPR) overlapping with CNV_1815 showed remarkably higher copy number in Qinchuan than Nanyang (log2 (ratio) = -2.34988; P value = 1.53E-102). Further qPCR and association analysis investigated that the copy number of the LEPR gene presented positive correlations with transcriptional expression and phenotypic traits, suggesting the LEPR CNV may contribute to the higher fat deposition in muscles of Qinchuan cattle. Our findings provide evidence that the distinct phenotypes of Nanyang and Qinchuan breeds may be due to the different genetic variations including SNPs, indels and CNV.

Citing Articles

Insights into Adaption and Growth Evolution: Genome-Wide Copy Number Variation Analysis in Chinese Hainan Yellow Cattle Using Whole-Genome Re-Sequencing Data.

Zhong Z, Wang Z, Xie X, Pan D, Su Z, Fan J Int J Mol Sci. 2024; 25(22).

PMID: 39595990 PMC: 11594005. DOI: 10.3390/ijms252211919.

Uncovering structural variants in Creole cattle from Guadeloupe and their impact on environmental adaptation through whole genome sequencing.

Ben-Jemaa S, Boussaha M, Mandonnet N, Bardou P, Naves M PLoS One. 2024; 19(8):e0309411.

PMID: 39186744 PMC: 11346954. DOI: 10.1371/journal.pone.0309411.

Genome-wide analysis reveals genomic diversity and signatures of selection in Qinchuan beef cattle.

Yu H, Zhang K, Cheng G, Mei C, Wang H, Zan L BMC Genomics. 2024; 25(1):558.

PMID: 38834950 PMC: 11151602. DOI: 10.1186/s12864-024-10482-0.

Whole-genome sequencing and genomic analysis of Norduz goat (Capra hircus).

Arslan M Mamm Genome. 2023; 34(3):437-448.

PMID: 37004528 DOI: 10.1007/s00335-023-09990-3.

Genome-wide analysis of CNVs in three populations of Tibetan sheep using whole-genome resequencing.

Hu L, Zhang L, Li Q, Liu H, Xu T, Zhao N Front Genet. 2022; 13:971464.

PMID: 36160022 PMC: 9490000. DOI: 10.3389/fgene.2022.971464.

References

Jeon J, Shim S, Nam H, Ryu G, Hong E, Kim H . Copy number variation at leptin receptor gene locus associated with metabolic traits and the risk of type 2 diabetes mellitus. BMC Genomics. 2010; 11:426. PMC: 2996954. DOI: 10.1186/1471-2164-11-426. View

Liu G, Van Tassell C, Sonstegard T, Li R, Alexander L, Keele J . Detection of germline and somatic copy number variations in cattle. Dev Biol (Basel). 2008; 132:231-237. DOI: 10.1159/000317165. View

Du Z, Zhou X, Ling Y, Zhang Z, Su Z . agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010; 38(Web Server issue):W64-70. PMC: 2896167. DOI: 10.1093/nar/gkq310. View

Kawahara-Miki R, Tsuda K, Shiwa Y, Arai-Kichise Y, Matsumoto T, Kanesaki Y . Whole-genome resequencing shows numerous genes with nonsynonymous SNPs in the Japanese native cattle Kuchinoshima-Ushi. BMC Genomics. 2011; 12:103. PMC: 3048544. DOI: 10.1186/1471-2164-12-103. View

Stothard P, Choi J, Basu U, Sumner-Thomson J, Meng Y, Liao X . Whole genome resequencing of black Angus and Holstein cattle for SNP and CNV discovery. BMC Genomics. 2011; 12:559. PMC: 3229636. DOI: 10.1186/1471-2164-12-559. View

Lee K, Chung W, Lee S, Choi J, Kim J, Lim D . Whole-genome resequencing of Hanwoo (Korean cattle) and insight into regions of homozygosity. BMC Genomics. 2013; 14:519. PMC: 3750754. DOI: 10.1186/1471-2164-14-519. View

Xu L, Hou Y, Bickhart D, Zhou Y, Abdel Hay E, Song J . Population-genetic properties of differentiated copy number variations in cattle. Sci Rep. 2016; 6:23161. PMC: 4804293. DOI: 10.1038/srep23161. View

Blattler A, Yao L, Witt H, Guo Y, Nicolet C, Berman B . Global loss of DNA methylation uncovers intronic enhancers in genes showing expression changes. Genome Biol. 2014; 15(9):469. PMC: 4203885. DOI: 10.1186/s13059-014-0469-0. View

Zhou X, Su Z . EasyGO: Gene Ontology-based annotation and functional enrichment analysis tool for agronomical species. BMC Genomics. 2007; 8:246. PMC: 1940007. DOI: 10.1186/1471-2164-8-246. View

10.

Mitchell-Olds T, Willis J, Goldstein D . Which evolutionary processes influence natural genetic variation for phenotypic traits?. Nat Rev Genet. 2007; 8(11):845-56. DOI: 10.1038/nrg2207. View

11.

Hinds D, Stuve L, Nilsen G, Halperin E, Eskin E, Ballinger D . Whole-genome patterns of common DNA variation in three human populations. Science. 2005; 307(5712):1072-9. DOI: 10.1126/science.1105436. View

12.

Xu Y, Zhang L, Shi T, Zhou Y, Cai H, Lan X . Copy number variations of MICAL-L2 shaping gene expression contribute to different phenotypes of cattle. Mamm Genome. 2013; 24(11-12):508-16. DOI: 10.1007/s00335-013-9483-x. View

13.

Bickhart D, Hou Y, Schroeder S, Alkan C, Cardone M, Matukumalli L . Copy number variation of individual cattle genomes using next-generation sequencing. Genome Res. 2012; 22(4):778-90. PMC: 3317159. DOI: 10.1101/gr.133967.111. View

14.

Flicek P, Amode M, Barrell D, Beal K, Brent S, Carvalho-Silva D . Ensembl 2012. Nucleic Acids Res. 2011; 40(Database issue):D84-90. PMC: 3245178. DOI: 10.1093/nar/gkr991. View

15.

Ellegren H, Sheldon B . Genetic basis of fitness differences in natural populations. Nature. 2008; 452(7184):169-75. DOI: 10.1038/nature06737. View

16.

Keel B, Lindholm-Perry A, Snelling W . Evolutionary and Functional Features of Copy Number Variation in the Cattle Genome. Front Genet. 2016; 7:207. PMC: 5118444. DOI: 10.3389/fgene.2016.00207. View

17.

Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N . The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009; 25(16):2078-9. PMC: 2723002. DOI: 10.1093/bioinformatics/btp352. View

18.

McCarthy M, Abecasis G, Cardon L, Goldstein D, Little J, Ioannidis J . Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nat Rev Genet. 2008; 9(5):356-69. DOI: 10.1038/nrg2344. View

19.

Mackay T . The genetic architecture of quantitative traits. Annu Rev Genet. 2001; 35:303-39. DOI: 10.1146/annurev.genet.35.102401.090633. View

20.

Hou Y, Liu G, Bickhart D, Cardone M, Wang K, Kim E . Genomic characteristics of cattle copy number variations. BMC Genomics. 2011; 12:127. PMC: 3053260. DOI: 10.1186/1471-2164-12-127. View