Identification of Whole‐genome Significant Single Nucleotide Polymorphisms in Candidate Genes Associated with Body Conformation Traits in Chinese Holstein Cattle
Zhengui Yan,Zhonghua Wang,Qin Zhang,Shujian Yue,Bin Yin,Yunliang Jiang,Kerong Shi
DOI: https://doi.org/10.1111/age.12865
2019-01-01
Animal Genetics
Abstract:GWAS has been proven to be a powerful tool for detecting genetic variants associated with economically important traits such as production,1, 2 body conformation3, 4 and reproduction traits,5 and diseases.6 The objective of this study was to identify SNPs with significant association effects on body conformation traits in Chinese Holstein through the GWAS approach. The experimental population consisted of 445 Chinese Holstein cows, with parity no. 2, raised at four farms in Shandong province of China. Hair follicle samples were collected individually within four days in the same season. In total, 29 body conformation traits were analyzed by GWAS. There were 23 phenotypes measured individually on 1–9 scores, which were body height, height at front end, body depth, chest width, stature, angularity, rump angle (RA), pin width, loin strength, foot angle, heel depth, bone quality, rear legs – side view, rear legs – rear view, udder depth, udder texture, median suspensory, fore udder attachment, front teat placement, attachment height, attachment width (AW), rear teat placement and teat length. The remaining six functional body conformation traits, measured on 1–100 scores, were stature score, rump system score, feet and legs score, udder system score, mammary system score and conformation final score (FS). The definitions of these conformation traits for cattle are available on the website of Canadian Dairy Network (https://www.cdn.ca/articles.php). The summary statistics of these phenotypes are listed in Table S1. Each experimental animal was genotyped by the GGP BovineLD V3 SNP chip, containing 26 151 public SNPs. Because all of the animals in these study populations were female, the SNPs from the X chromosome were also counted. The data quality control procedure excluded individuals with more than 10% missing genotypes and SNPs with a call rate of less than 90%, MAF less than 3% or P-value of the Hardy–Weinberg Equilibrium test less than 1.0E − 6. Finally, 421 animals with 20 632 SNP genotypes were retained for the subsequent GWAS analysis. Physical map length, the number of SNPs and SNP density on each chromosome based on reference genome UMD3.1, before and after the data cleaning procedure, are shown in Table S2. Results of pair-wise LD analysis of the 421 animals from four different farms showed high similarity of LD patterns among the whole population, indicating that similar breeding histories were shared by the four subpopulations (data not provided). A fixed-effects linear regression model was used to carry out the GWAS according to a previous method.7, 2 Bonferroni correction for the genome-wide significance and suggestive thresholds were computed to be 2.42E − 06 (= 0.05/20 632) and 4.85E − 05 (= 1/20 632), respectively. The GWAS based on the mixed-effect linear regression model identified 79 SNPs with genome-wide significant (P < 0.05) association effects on 14 body conformation traits (Table 1, Table S3). Of these SNPs, 28 SNPs were found to be associated with RA, and 14 SNPs associated with FS (Table 1, Fig. S1). We also found six single SNPs that were associated with multiple body conformation traits of dairy cow. These SNPs were BovineHD0600000461, BovineHD1900007686, ARS-BFGL-NGS-41612, BovineHD1700010514, BovineHD3000037672 and ARS-BFGL-NGS-109467 (their detailed information is provided in Table S4). In this study, the DNA region from 44.6 to 65.4 Mb of chromosome 18 was identified to be associated with four different body conformation traits (AW, FS, mammary system score and RA) of dairy cows (Table 2, Table S3). This coincided with previously reported QTL or DNA regions associated with calving performance and udder traits of dairy cattle.8-10 Among these SNPs, Hapmap57004-rs29011610, located near MON1B (MON1 secretory trafficking family member B) gene, were also found to be associated with length of productive life in Holstein cows.11 ARS-BFGL-NGS-116541, which was located within the LIG1 (Ligase I) gene region, was reported to be associated with body weight at birth.12 Similarly, DNA regions on chromosomes 21, 25 and X were also identified to be significantly associated with multiple body conformation traits (Table 2, Table S3). In particular, six SNPs in the small region from 57.55 to 57.62 Mb on the chromosome 21 were found to be associated with AW and FS in the study, and these SNPs were all located within or near the gene region of SLC24A4 (Table S5). Interestingly, SLC24A5 and SLC19A2, members of solute carrier family proteins with SLC24A4, were also found to be associated with RA and FS, respectively (Table 1). Additionally, UA-IFASA-6670, located within the region of GABARAPL1 [GABA (A) receptor-associated protein like 1] gene, were found to be significantly associated with median suspensory in this study. It was reported previously to affect udder attachment and height.13, 14 Also, BovineHD2200013812 was identified to have genetic effects on RA. This SNP was located within the CACNA1D (calcium channel alpha 1D subunit) gene region. Coincidentally, this SNP was also reported to be associated with the length of productive life and udder cleft.14 In summary, a GWAS using linear statistical model was conducted on 29 body conformation traits in a Chinese Holstein cattle population, and 79 SNPs were found to have genome-wide-significant (P < 0.05) association effects on 14 body conformation traits. Among these significant SNPs, 74 of them are newly detected in this study, five have been reported in previous literature and 26 are located in genes and are worth further investigation to potentially identify the causative mutations underlying the QTL. The authors would like to thank the owners of the four farms who generously allowed us to sample hair follicles from their cattle. We are also thankful to Dr Chunsheng Ma, Dr Jizhen Wang, Dr Sihao Wu and Dr Enfa Wang at Alta-Agricorp Co. Ltd for their kind support and suggestions on the verification of phenotype data. We are grateful to Dr Changsong Qian at Neogen Bio-Scientific Technology (Shanghai) Co. Ltd for his technical support during the hair follicle collection. This work was financially supported by the Key Project of Agricultural Fine Breeding of Shandong Province (2016LZGC030, 2019LZGC011), the National Natural Science Foundation of China (31402054), the Natural Science Foundation of Shandong (ZR2013CM013), Funds of the Shandong 'Double Tops' Program (SYL2017YSTD08), the Modern Agricultural Industry Technology System (CARS-36) and the Tai Mountain Scholar Program. The authors declare that there is no conflict of interest. Table S1 Descriptive statistics of the 29 conformation traits used in the GWAS. Table S2 Distribution of SNP markers by chromosomes before and after quality control. Table S3 The chromosomal distribution of significant SNPs associated with body conformation traits. Table S4 The SNPs identified associated with multiple body conformation traits of dairy cows. Table S5 Multiple SNPs located in the SLC24A4 gene region were significantly associated with body conformation traits of dairy cows. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.