1
|
de Oliveira LF, Veroneze R, Sousa KRS, Mulim HA, Araujo AC, Huang Y, Johnson JS, Brito LF. Genomic regions, candidate genes, and pleiotropic variants associated with physiological and anatomical indicators of heat stress response in lactating sows. BMC Genomics 2024; 25:467. [PMID: 38741036 DOI: 10.1186/s12864-024-10365-4] [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: 02/12/2024] [Accepted: 04/29/2024] [Indexed: 05/16/2024] Open
Abstract
BACKGROUND Heat stress (HS) poses significant threats to the sustainability of livestock production. Genetically improving heat tolerance could enhance animal welfare and minimize production losses during HS events. Measuring phenotypic indicators of HS response and understanding their genetic background are crucial steps to optimize breeding schemes for improved climatic resilience. The identification of genomic regions and candidate genes influencing the traits of interest, including variants with pleiotropic effects, enables the refinement of genotyping panels used to perform genomic prediction of breeding values and contributes to unraveling the biological mechanisms influencing heat stress response. Therefore, the main objectives of this study were to identify genomic regions, candidate genes, and potential pleiotropic variants significantly associated with indicators of HS response in lactating sows using imputed whole-genome sequence (WGS) data. Phenotypic records for 18 traits and genomic information from 1,645 lactating sows were available for the study. The genotypes from the PorcineSNP50K panel containing 50,703 single nucleotide polymorphisms (SNPs) were imputed to WGS and after quality control, 1,622 animals and 7,065,922 SNPs were included in the analyses. RESULTS A total of 1,388 unique SNPs located on sixteen chromosomes were found to be associated with 11 traits. Twenty gene ontology terms and 11 biological pathways were shown to be associated with variability in ear skin temperature, shoulder skin temperature, rump skin temperature, tail skin temperature, respiration rate, panting score, vaginal temperature automatically measured every 10 min, vaginal temperature measured at 0800 h, hair density score, body condition score, and ear area. Seven, five, six, two, seven, 15, and 14 genes with potential pleiotropic effects were identified for indicators of skin temperature, vaginal temperature, animal temperature, respiration rate, thermoregulatory traits, anatomical traits, and all traits, respectively. CONCLUSIONS Physiological and anatomical indicators of HS response in lactating sows are heritable but highly polygenic. The candidate genes found are associated with important gene ontology terms and biological pathways related to heat shock protein activities, immune response, and cellular oxidative stress. Many of the candidate genes with pleiotropic effects are involved in catalytic activities to reduce cell damage from oxidative stress and cellular mechanisms related to immune response.
Collapse
Affiliation(s)
- Letícia Fernanda de Oliveira
- Department of Animal Science, Federal University of Viçosa, Viçosa, MG, Brazil
- Department of Animal Sciences, Purdue University, West Lafayette, IN, USA
| | - Renata Veroneze
- Department of Animal Science, Federal University of Viçosa, Viçosa, MG, Brazil
| | - Katiene Régia Silva Sousa
- Department of Animal Sciences, Purdue University, West Lafayette, IN, USA
- Department of Oceanography and Limnology, Federal University of Maranhão, São Luís, MA, Brazil
| | - Henrique A Mulim
- Department of Animal Sciences, Purdue University, West Lafayette, IN, USA
| | | | | | - Jay S Johnson
- USDA-ARS Livestock Behavior Research Unit, West Lafayette, IN, USA
| | - Luiz F Brito
- Department of Animal Sciences, Purdue University, West Lafayette, IN, USA.
| |
Collapse
|
2
|
Sousa LPB, Pinto LFB, Cruz VAR, Oliveira GA, Rojas de Oliveira H, Chud TS, Pedrosa VB, Miglior F, Schenkel FS, Brito LF. Genome-wide association and functional genomic analyses for various hoof health traits in North American Holstein cattle. J Dairy Sci 2024; 107:2207-2230. [PMID: 37939841 DOI: 10.3168/jds.2023-23806] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2023] [Accepted: 10/19/2023] [Indexed: 11/10/2023]
Abstract
Hoof diseases are a major welfare and economic issue in the global dairy cattle production industry, which can be minimized through improved management and breeding practices. Optimal genetic improvement of hoof health could benefit from a deep understanding of the genetic background and biological underpinning of indicators of hoof health. Therefore, the primary objectives of this study were to perform genome-wide association studies, using imputed high-density genetic markers data from North American Holstein cattle, for 8 hoof-related traits: digital dermatitis, sole ulcer, sole hemorrhage, white line lesion, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, and toe ulcer, and a hoof health index. De-regressed estimated breeding values from 25,580 Holstein animals were used as pseudo-phenotypes for the association analyses. The genomic quality control, genotype phasing, and genotype imputation were performed using the PLINK (version 1.9), Eagle (version 2.4.1), and Minimac4 software, respectively. The functional genomic analyses were performed using the GALLO R package and the DAVID platform. We identified 22, 34, 14, 22, 28, 33, 24, 43, and 15 significant markers for digital dermatitis, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, sole hemorrhage, sole ulcer, toe ulcer, white line lesion disease, and the hoof health index, respectively. The significant markers were located across all autosomes, except BTA10, BTA12, BTA20, BTA26, BTA27, and BTA28. Moreover, the genomic regions identified overlap with various previously reported quantitative trait loci for exterior, health, meat and carcass, milk, production, and reproduction traits. The enrichment analyses identified 44 significant gene ontology terms. These enriched genomic regions harbor various candidate genes previously associated with bone development, metabolism, and infectious and immunological diseases. These findings indicate that hoof health traits are highly polygenic and influenced by a wide range of biological processes.
Collapse
Affiliation(s)
- Luis Paulo B Sousa
- Department of Animal Sciences, Federal University of Bahia, Salvador, BA, 40170-110, Brazil
| | - Luis Fernando B Pinto
- Department of Animal Sciences, Federal University of Bahia, Salvador, BA, 40170-110, Brazil
| | - Valdecy A R Cruz
- Department of Animal Sciences, Federal University of Bahia, Salvador, BA, 40170-110, Brazil
| | - Gerson A Oliveira
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Hinayah Rojas de Oliveira
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada; Department of Animal Sciences, Purdue University, West Lafayette, IN 47907
| | - Tatiane S Chud
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada; PEAK, Madison, WI 53718
| | - Victor B Pedrosa
- Department of Animal Sciences, Purdue University, West Lafayette, IN 47907
| | - Filippo Miglior
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada; Lactanet Canada, Guelph, ON, N1K 1E5, Canada
| | - Flávio S Schenkel
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Luiz F Brito
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada; Department of Animal Sciences, Purdue University, West Lafayette, IN 47907.
| |
Collapse
|
3
|
Chen SY, Gloria LS, Pedrosa VB, Doucette J, Boerman JP, Brito LF. Unraveling the genomic background of resilience based on variability in milk yield and milk production levels in North American Holstein cattle through genome-wide association study and Mendelian randomization analyses. J Dairy Sci 2024; 107:1035-1053. [PMID: 37776995 DOI: 10.3168/jds.2023-23650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2023] [Accepted: 09/04/2023] [Indexed: 10/02/2023]
Abstract
Breeding more resilient animals will benefit the dairy cattle industry in the long term, especially as global climate changes become more severe. Previous studies have reported genetic parameters for various milk yield-based resilience indicators, but the underlying genomic background of these traits remain unknown. In this study, we conducted GWAS of 62,029 SNPs with 4 milk yield-based resilience indicators, including the weighted occurrence frequency (wfPert) and accumulated milk losses (dPert) of milk yield perturbations, and log-transformed variance (LnVar) and lag-1 autocorrelation (rauto) of daily yield residuals. These variables were previously derived from 5.6 million daily milk yield records from 21,350 lactations (parities 1-3) of 11,787 North American Holstein cows. The average daily milk yield (ADMY) throughout lactation was also included to compare the shared genetic background of resilience indicators with milk yield. The differential genetic background of these indicators was first revealed by the significant genomic regions identified and significantly enriched biological pathways of positional candidate genes, which confirmed the genetic difference among resilience indicators. Interestingly, the functional analyses of candidate genes suggested that the regulation of intestinal homeostasis is most likely affecting resilience derived based on variability in milk yield. Based on Mendelian randomization analyses of multiple instrumental SNPs, we further found an unfavorable causal association of ADMY with LnVar. In conclusion, the resilience indicators evaluated are genetically different traits, and there are causal associations of milk yield with some of the resilience indicators evaluated. In addition to providing biological insights into the molecular regulation mechanisms of resilience derived based on variability in milk yield, this study also indicates the need for developing selection indexes combining multiple indicator traits and taking into account their genetic relationship for breeding more resilient dairy cattle.
Collapse
Affiliation(s)
- Shi-Yi Chen
- Department of Animal Sciences, Purdue University, West Lafayette, IN 47907; Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan, 611130, China
| | - Leonardo S Gloria
- Department of Animal Sciences, Purdue University, West Lafayette, IN 47907
| | - Victor B Pedrosa
- Department of Animal Sciences, Purdue University, West Lafayette, IN 47907
| | - Jarrod Doucette
- Agriculture Information Technology (AgIT), Purdue University, West Lafayette, IN 47907
| | | | - Luiz F Brito
- Department of Animal Sciences, Purdue University, West Lafayette, IN 47907.
| |
Collapse
|
4
|
Singh A, Malla WA, Kumar A, Jain A, Thakur MS, Khare V, Tiwari SP. Review: genetic background of milk fatty acid synthesis in bovines. Trop Anim Health Prod 2023; 55:328. [PMID: 37749432 DOI: 10.1007/s11250-023-03754-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Accepted: 09/12/2023] [Indexed: 09/27/2023]
Abstract
Milk fat composition is an important trait for the dairy industry as it directly influences the nutritional and technological properties of milk and other dairy products. The synthesis of milk fat is a complex process regulated by a network of genes. Thus, understanding the genetic variation and molecular mechanisms regulating milk fat synthesis will help to improve the nutritional quality of dairy products. In this review, we provide an overview of milk fat synthesis in bovines along with the candidate genes involved in the pathway. We also discuss de novo synthesis of fatty acids (ACSS, ACACA, FASN), uptake of FAs (FATP, FAT, LPL), intracellular activation and channelling of FAs (ACSL, FABP), elongation (EVOLV6), desaturation (SCD, FADS), formation of triglycerides (GPAM, AGPAT, LIPIN, DGAT), and milk lipid secretion (BTN1A1, XDH, PLIN2). The genetic variability of individual fatty acids will help to develop selection strategies for obtaining a healthier milk fat profile in bovines. Thus, this review will offer a potential understanding of the molecular mechanisms that regulate milk fat synthesis in bovines.
Collapse
Affiliation(s)
- Akansha Singh
- College of Veterinary Science and Animal Husbandry, NDVSU, Jabalpur, M.P, 482001, India.
| | - Waseem Akram Malla
- ICMR-National Institute of Malaria Research, Field Unit Guwahati, Assam, 781022, India
| | - Amit Kumar
- ICAR- Indian Veterinary Research Institute, Izatnagar, Bareilly, U.P, 243122, India
| | - Asit Jain
- College of Veterinary Science and Animal Husbandry, NDVSU, Jabalpur, M.P, 482001, India
| | - Mohan Singh Thakur
- College of Veterinary Science and Animal Husbandry, NDVSU, Jabalpur, M.P, 482001, India
| | - Vaishali Khare
- College of Veterinary Science and Animal Husbandry, NDVSU, Jabalpur, M.P, 482001, India
| | - Sita Prasad Tiwari
- College of Veterinary Science and Animal Husbandry, NDVSU, Jabalpur, M.P, 482001, India
| |
Collapse
|
5
|
Chen Y, Atashi H, Grelet C, Mota RR, Vanderick S, Hu H, Gengler N. Genome-wide association study and functional annotation analyses for nitrogen efficiency index and its composition traits in dairy cattle. J Dairy Sci 2023; 106:3397-3410. [PMID: 36894424 DOI: 10.3168/jds.2022-22351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2022] [Accepted: 10/24/2022] [Indexed: 03/09/2023]
Abstract
The aims of this study were (1) to identify genomic regions associated with a N efficiency index (NEI) and its composition traits and (2) to analyze the functional annotation of identified genomic regions. The NEI included N intake (NINT1), milk true protein N (MTPN1), milk urea N yield (MUNY1) in primiparous cattle, and N intake (NINT2+), milk true protein N (MTPN2+), and milk urea N yield (MUNY2+) in multiparous cattle (2 to 5 parities). The edited data included 1,043,171 records on 342,847 cows distributed in 1,931 herds. The pedigree consisted of 505,125 animals (17,797 males). Data of 565,049 SNPs were available for 6,998 animals included in the pedigree (5,251 females and 1,747 males). The SNP effects were estimated using a single-step genomic BLUP approach. The proportion of the total additive genetic variance explained by windows of 50 consecutive SNPs (with an average size of about 240 kb) was calculated. The top 3 genomic regions explaining the largest rate of the total additive genetic variance of the NEI and its composition traits were selected for candidate gene identification and quantitative trait loci (QTL) annotation. The selected genomic regions explained from 0.17% (MTPN2+) to 0.58% (NEI) of the total additive genetic variance. The largest explanatory genomic regions of NEI, NINT1, NINT2+, MTPN1, MTPN2+, MUNY1, and MUNY2+ were Bos taurus autosome 14 (1.52-2.09 Mb), 26 (9.24-9.66 Mb), 16 (75.41-75.51 Mb), 6 (8.73-88.92 Mb), 6 (8.73-88.92 Mb), 11 (103.26-103.41 Mb), 11 (103.26-103.41 Mb). Based on the literature, gene ontology, Kyoto Encyclopedia of Genes and Genomes, and protein-protein interaction, 16 key candidate genes were identified for NEI and its composition traits, which are mainly expressed in the milk cell, mammary, and liver tissues. The number of enriched QTL related to NEI, NINT1, NINT2+, MTPN1, and MTPN2+ were 41, 6, 4, 11, 36, 32, and 32, respectively, and most of them were related to the milk, health, and production classes. In conclusion, this study identified genomic regions associated with NEI and its composition traits, and identified key candidate genes describing the genetic mechanisms of N use efficiency-related traits. Furthermore, the NEI reflects not only its composition traits but also the interactions among them.
Collapse
Affiliation(s)
- Y Chen
- TERRA Teaching and Research Center, University of Liège, Gembloux Agro-Bio Tech (ULiège-GxABT), 5030 Gembloux, Belgium.
| | - H Atashi
- TERRA Teaching and Research Center, University of Liège, Gembloux Agro-Bio Tech (ULiège-GxABT), 5030 Gembloux, Belgium; Department of Animal Science, Shiraz University, 71441-65186 Shiraz, Iran
| | - C Grelet
- Walloon Agricultural Research Center (CRA-W), 5030 Gembloux, Belgium
| | - R R Mota
- Council on Dairy Cattle Breeding, Bowie, MD 20716
| | - S Vanderick
- TERRA Teaching and Research Center, University of Liège, Gembloux Agro-Bio Tech (ULiège-GxABT), 5030 Gembloux, Belgium
| | - H Hu
- TERRA Teaching and Research Center, University of Liège, Gembloux Agro-Bio Tech (ULiège-GxABT), 5030 Gembloux, Belgium
| | | | - N Gengler
- TERRA Teaching and Research Center, University of Liège, Gembloux Agro-Bio Tech (ULiège-GxABT), 5030 Gembloux, Belgium
| |
Collapse
|
6
|
Nazar M, Abdalla IM, Chen Z, Ullah N, Liang Y, Chu S, Xu T, Mao Y, Yang Z, Lu X. Genome-Wide Association Study for Udder Conformation Traits in Chinese Holstein Cattle. Animals (Basel) 2022; 12:2542. [PMID: 36230283 PMCID: PMC9559277 DOI: 10.3390/ani12192542] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2022] [Revised: 09/14/2022] [Accepted: 09/20/2022] [Indexed: 11/16/2022] Open
Abstract
Udder conformation traits are one of the most economic traits in dairy cows, greatly affecting animal health, milk production, and producer profitability in the dairy industry. Genetic analysis of udder structure and scores have been developed in Holstein cattle. In our research, we conducted a genome-wide association study for five udder traits, including anterior udder attachment (AUA), central suspensory ligament (CSL), posterior udder attachment height (PUAH), posterior udder attachment width (PUAW), and udder depth (UD), in which the fixed and random model circulating probability unification (FarmCPU) model was applied for the association analysis. The heritability and the standard errors of these five udder traits ranged from 0.04 ± 0.00 to 0.49 ± 0.03. Phenotype data were measured from 1000 Holstein cows, and the GeneSeek Genomic Profiler (GGP) Bovine 100 K SNP chip was used to analyze genotypic data in Holstein cattle. For GWAS analysis, 984 individual cows and 84,407 single-nucleotide polymorphisms (SNPs) remained after quality control; a total of 18 SNPs were found at the GW significant threshold (p < 5.90 × 10−7). Many candidate genes were identified within 200kb upstream or downstream of the significant SNPs, which include MGST1, MGST2, MTUS1, PRKN, STXBP6, GRID2, E2F8, CDH11, FOXP1, SLF1, TMEM117, SBF2, GC, ADGRB3, and GCLC. Pathway analysis revealed that 58 Gene Ontology (GO) terms and 18 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched with adjusted p values, and these GO terms and the KEGG pathway analysis were associated with biological information, metabolism, hormonal growth, and development processes. These results could give valuable biological information for the genetic architecture of udder conformation traits in dairy Holstein cattle.
Collapse
Affiliation(s)
- Mudasir Nazar
- College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China
| | | | - Zhi Chen
- College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China
| | - Numan Ullah
- College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China
| | - Yan Liang
- College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China
| | - Shuangfeng Chu
- 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
| | - Yongjiang Mao
- College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China
| | - Zhangping Yang
- College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China
| | - Xubin Lu
- College of Animal Science and Technology, Yangzhou University, Yangzhou 225009, China
| |
Collapse
|
7
|
Chen SY, Schenkel FS, Melo ALP, Oliveira HR, Pedrosa VB, Araujo AC, Melka MG, Brito LF. Identifying pleiotropic variants and candidate genes for fertility and reproduction traits in Holstein cattle via association studies based on imputed whole-genome sequence genotypes. BMC Genomics 2022; 23:331. [PMID: 35484513 PMCID: PMC9052698 DOI: 10.1186/s12864-022-08555-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Accepted: 04/12/2022] [Indexed: 02/06/2023] Open
Abstract
Background Genetic progress for fertility and reproduction traits in dairy cattle has been limited due to the low heritability of most indicator traits. Moreover, most of the quantitative trait loci (QTL) and candidate genes associated with these traits remain unknown. In this study, we used 5.6 million imputed DNA sequence variants (single nucleotide polymorphisms, SNPs) for genome-wide association studies (GWAS) of 18 fertility and reproduction traits in Holstein cattle. Aiming to identify pleiotropic variants and increase detection power, multiple-trait analyses were performed using a method to efficiently combine the estimated SNP effects of single-trait GWAS based on a chi-square statistic. Results There were 87, 72, and 84 significant SNPs identified for heifer, cow, and sire traits, respectively, which showed a wide and distinct distribution across the genome, suggesting that they have relatively distinct polygenic nature. The biological functions of immune response and fatty acid metabolism were significantly enriched for the 184 and 124 positional candidate genes identified for heifer and cow traits, respectively. No known biological function was significantly enriched for the 147 positional candidate genes found for sire traits. The most important chromosomes that had three or more significant QTL identified are BTA22 and BTA23 for heifer traits, BTA8 and BTA17 for cow traits, and BTA4, BTA7, BTA17, BTA22, BTA25, and BTA28 for sire traits. Several novel and biologically important positional candidate genes were strongly suggested for heifer (SOD2, WTAP, DLEC1, PFKFB4, TRIM27, HECW1, DNAH17, and ADAM3A), cow (ANXA1, PCSK5, SPESP1, and JMJD1C), and sire (ELMO1, CFAP70, SOX30, DGCR8, SEPTIN14, PAPOLB, JMJD1C, and NELL2) traits. Conclusions These findings contribute to better understand the underlying biological mechanisms of fertility and reproduction traits measured in heifers, cows, and sires, which may contribute to improve genomic evaluation for these traits in dairy cattle. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08555-z.
Collapse
Affiliation(s)
- Shi-Yi Chen
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA.,Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, 611130, Sichuan, China
| | - Flavio S Schenkel
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Ana L P Melo
- Department of Reproduction and Animal Evaluation, Rural Federal University of Rio de Janeiro, Seropédica, RJ, 23897-000, Brazil
| | - Hinayah R Oliveira
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA.,Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Victor B Pedrosa
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA.,Department of Animal Sciences, State University of Ponta Grossa, Ponta Grossa, PR, 84030-900, Brazil
| | - Andre C Araujo
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA
| | - Melkaye G Melka
- Department of Animal and Food Science, University of Wisconsin River Falls, River Falls, WI, 54022, USA
| | - Luiz F Brito
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907-2041, USA. .,Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada.
| |
Collapse
|
8
|
Raza SHA, Liang C, Alshammari AM, Aloufi BH, Gui L, Khan R, Zan L. Genetic variants in the serum amyloid A2 (SAA2) gene as a potential marker for milk production traits in Chinese Holstein cows. Vet Med Sci 2022; 8:1835-1840. [PMID: 35470977 PMCID: PMC9297766 DOI: 10.1002/vms3.796] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022] Open
Abstract
Background This study was conducted to detect potential polymorphisms of the serum amyloid A2 (SAA2) gene and explore their relationships with milk production traits in Chinese Holstein cows. Objectives: This study used sequencing technology conducted in 532 Chinese Holstein cows. Methods Three single nucleotide polymorphisms (SNPs) were identified within intron 1, named g.14061A>G, g.14072G>C and g.14819C>T. Eight estimated haplotypes were identified, of which three major haplotypes had a frequency of Hap3 (‐ACC‐), Hap5 (‐GCC‐) and Hap2 (‐AGT‐), with 17.9%, 12.30% and 8.10%, respectively. Results The association analysis of single markers (g.14061A>G and g.14819C>T) and combined genotypes (Hap1/5) revealed prominent effects on milk production traits in Chinese Holstein cows (p < 0.05). Conclusions Our results suggest that the SAA2 gene is associated with economic traits in Chinese Holstein cows and may be used as candidate gene for marker‐assisted selection and management in breeding programs.
Collapse
Affiliation(s)
- Sayed Haidar Abbas Raza
- State Key Laboratory of Animal Genetics Breeding & Reproduction, College of Animal Science and TechnologyNorthwest A&F UniversityYanglingP. R. China
- National Beef Cattle Improvement CenterNorthwest A&F UniversityYanglingP. R. China
| | - Chengcheng Liang
- State Key Laboratory of Animal Genetics Breeding & Reproduction, College of Animal Science and TechnologyNorthwest A&F UniversityYanglingP. R. China
- National Beef Cattle Improvement CenterNorthwest A&F UniversityYanglingP. R. China
| | | | - Bandar H. Aloufi
- Faculty of Science, Biology DepartmentUniversity of HailHailSaudi Arabia
| | - Linsheng Gui
- State Key Laboratory of Plateau Ecology and AgricultureQinghai UniversityXiningP. R. China
| | - Rajwali Khan
- Department of Livestock Management, Breeding and GeneticsThe University of AgriculturePeshawarPakistan
| | - Linsen Zan
- State Key Laboratory of Animal Genetics Breeding & Reproduction, College of Animal Science and TechnologyNorthwest A&F UniversityYanglingP. R. China
- National Beef Cattle Improvement CenterNorthwest A&F UniversityYanglingP. R. China
| |
Collapse
|
9
|
Pedrosa VB, Schenkel FS, Chen SY, Oliveira HR, Casey TM, Melka MG, Brito LF. Genomewide Association Analyses of Lactation Persistency and Milk Production Traits in Holstein Cattle Based on Imputed Whole-Genome Sequence Data. Genes (Basel) 2021; 12:genes12111830. [PMID: 34828436 PMCID: PMC8624223 DOI: 10.3390/genes12111830] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2021] [Revised: 11/13/2021] [Accepted: 11/17/2021] [Indexed: 12/22/2022] Open
Abstract
Lactation persistency and milk production are among the most economically important traits in the dairy industry. In this study, we explored the association of over 6.1 million imputed whole-genome sequence variants with lactation persistency (LP), milk yield (MILK), fat yield (FAT), fat percentage (FAT%), protein yield (PROT), and protein percentage (PROT%) in North American Holstein cattle. We identified 49, 3991, 2607, 4459, 805, and 5519 SNPs significantly associated with LP, MILK, FAT, FAT%, PROT, and PROT%, respectively. Various known associations were confirmed while several novel candidate genes were also revealed, including ARHGAP35, NPAS1, TMEM160, ZC3H4, SAE1, ZMIZ1, PPIF, LDB2, ABI3, SERPINB6, and SERPINB9 for LP; NIM1K, ZNF131, GABRG1, GABRA2, DCHS1, and SPIDR for MILK; NR6A1, OLFML2A, EXT2, POLD1, GOT1, and ETV6 for FAT; DPP6, LRRC26, and the KCN gene family for FAT%; CDC14A, RTCA, HSTN, and ODAM for PROT; and HERC3, HERC5, LALBA, CCL28, and NEURL1 for PROT%. Most of these genes are involved in relevant gene ontology (GO) terms such as fatty acid homeostasis, transporter regulator activity, response to progesterone and estradiol, response to steroid hormones, and lactation. The significant genomic regions found contribute to a better understanding of the molecular mechanisms related to LP and milk production in North American Holstein cattle.
Collapse
Affiliation(s)
- Victor B. Pedrosa
- Department of Animal Sciences, Purdue University, West Lafayette, IN 47907, USA; (V.B.P.); (S.-Y.C.); (H.R.O.); (T.M.C.)
- Department of Animal Sciences, State University of Ponta Grossa, Ponta Grossa 84030-900, Brazil
| | - Flavio S. Schenkel
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON N1G2W1, Canada;
| | - Shi-Yi Chen
- Department of Animal Sciences, Purdue University, West Lafayette, IN 47907, USA; (V.B.P.); (S.-Y.C.); (H.R.O.); (T.M.C.)
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, College of Animal Science & Technology, Sichuan Agricultural University, Chengdu 611130, China
| | - Hinayah R. Oliveira
- Department of Animal Sciences, Purdue University, West Lafayette, IN 47907, USA; (V.B.P.); (S.-Y.C.); (H.R.O.); (T.M.C.)
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON N1G2W1, Canada;
| | - Theresa M. Casey
- Department of Animal Sciences, Purdue University, West Lafayette, IN 47907, USA; (V.B.P.); (S.-Y.C.); (H.R.O.); (T.M.C.)
| | - Melkaye G. Melka
- Department of Animal and Food Science, University of Wisconsin River Falls, River Falls, WI 54022, USA;
| | - Luiz F. Brito
- Department of Animal Sciences, Purdue University, West Lafayette, IN 47907, USA; (V.B.P.); (S.-Y.C.); (H.R.O.); (T.M.C.)
- Correspondence:
| |
Collapse
|
10
|
Genome-Wide Association Study Identifies Candidate Genes Associated with Feet and Leg Conformation Traits in Chinese Holstein Cattle. Animals (Basel) 2021; 11:ani11082259. [PMID: 34438715 PMCID: PMC8388412 DOI: 10.3390/ani11082259] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2021] [Revised: 07/24/2021] [Accepted: 07/28/2021] [Indexed: 12/12/2022] Open
Abstract
Simple Summary Feet and leg problems are among the major reasons for dairy cows leaving the herd, as well as having direct association with production and reproduction efficiency, health (e.g., claw disorders and lameness) and welfare. Hence, understanding the genetic architecture underlying feet and conformation traits in dairy cattle offers new opportunities toward the genetic improvement and long-term selection. Through a genome-wide association study on Chinese Holstein cattle, we identified several candidate genes associated with feet and leg conformation traits. These results could provide useful information about the molecular breeding basis of feet and leg traits, thus improving the longevity and productivity of dairy cattle. Abstract Feet and leg conformation traits are considered one of the most important economical traits in dairy cattle and have a great impact on the profitability of milk production. Therefore, identifying the single nucleotide polymorphisms (SNPs), genes and pathways analysis associated with these traits might contribute to the genomic selection and long-term plan selection for dairy cattle. We conducted genome-wide association studies (GWASs) using the fixed and random model circulating probability unification (FarmCPU) method to identify SNPs associated with bone quality, heel depth, rear leg side view and rear leg rear view of Chinese Holstein cows. Phenotypic measurements were collected from 1000 individuals of Chinese Holstein cattle and the GeneSeek Genomic Profiler Bovine 100 K SNP chip was utilized for individual genotyping. After quality control, 984 individual cows and 84,906 SNPs remained for GWAS work; as a result, we identified 20 significant SNPs after Bonferroni correction. Several candidate genes were identified within distances of 200 kb upstream or downstream to the significant SNPs, including ADIPOR2, INPP4A, DNMT3A, ALDH1A2, PCDH7, XKR4 and CADPS. Further bioinformatics analyses showed 34 gene ontology terms and two signaling pathways were significantly enriched (p ≤ 0.05). Many terms and pathways are related to biological quality, metabolism and development processes; these identified SNPs and genes could provide useful information about the genetic architecture of feet and leg traits, thus improving the longevity and productivity of Chinese Holstein dairy cattle.
Collapse
|
11
|
AGPAT3 Gene polymorphisms are associated with milk production traits in Chinese Holstein cows. J DAIRY RES 2021; 88:247-252. [PMID: 34261571 DOI: 10.1017/s0022029921000546] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
The current study reports the identification of previously undiscovered single-nucleotide polymorphisms (SNPs) in the bovine AGPAT3 gene and further investigates their associations with milk production traits. Our results demonstrate that the major allele C of the SNP g.12264 C > T is positively correlated with test-day milk yield, protein percentage and 305-day milk yield. Importantly, in silico analysis showed that the C/T transition at this locus gives rise to two new transcription factor binding sites (TFBS), E2F1 and Nkx3-2. Polymorphism g.18658 G > A was the only SNP associated with milk urea nitrogen (MUN) with the G allele related to an increase in milk urea nitrogen as well as fat percentage. The GG genotype of SNP g.28731 A > G was associated with the highest fat and protein percentage and lowest 305-day milk yield and somatic cell score (SCS). The association between AGPAT3 locus and milk production traits could be utilized in marker-assisted selection for the genetic improvement of milk production traits and, probably in conjunction with other traits, for selection to improve fitness of dairy cattle.
Collapse
|
12
|
Shi L, Wu X, Yang Y, Ma Z, Lv X, Liu L, Li Y, Zhao F, Han B, Sun D. A post-GWAS confirming the genetic effects and functional polymorphisms of AGPAT3 gene on milk fatty acids in dairy cattle. J Anim Sci Biotechnol 2021; 12:24. [PMID: 33522959 PMCID: PMC7849138 DOI: 10.1186/s40104-020-00540-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2020] [Accepted: 12/14/2020] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND People are paying more attention to the healthy and balanced diet with the improvement of their living standards. Milk fatty acids (FAs) have been reported that they were related to some atherosclerosis and coronary heart diseases in human. In our previous genome-wide association study (GWAS) on milk FAs in dairy cattle, 83 genome-wide significant single nucleotide polymorphisms (SNPs) were detected. Among them, two SNPs, ARS-BFGL-NGS-109493 and BTA-56389-no-rs associated with C18index (P = 0.0459), were located in the upstream of 1-acylglycerol-3-phosphate O-acyltransferase 3 (AGPAT3) gene. AGPAT3 is involved in glycerol-lipid, glycerol-phospholipid metabolism and phospholipase D signaling pathways. Hence, it was inferred as a candidate gene for milk FAs. The aim of this study was to further confirm the genetic effects of the AGPAT3 gene on milk FA traits in dairy cattle. RESULTS Through re-sequencing the complete coding region, and 3000 bp of 5' and 3' regulatory regions of the AGPAT3 gene, a total of 17 SNPs were identified, including four in 5' regulatory region, one in 5' untranslated region (UTR), three in introns, one in 3' UTR, and eight in 3' regulatory region. By the linkage disequilibrium (LD) analysis with Haploview4.1 software, two haplotype blocks were observed that were formed by four and 12 identified SNPs, respectively. Using SAS9.2, we performed single locus-based and haplotype-based association analysis on 24 milk FAs in 1065 Chinese Holstein cows, and discovered that all the SNPs and the haplotype blocks were significantly associated with C6:0, C8:0 and C10:0 (P < 0.0001-0.0384). Further, with Genomatix, we predicted that four SNPs in 5' regulatory region (g.146702957G > A, g.146704373A > G, g.146704618A > G and g.146704699G > A) changed the transcription factor binding sites (TFBSs) for transcription factors SMARCA3, REX1, VMYB, BRACH, NKX26, ZBED4, SP1, USF1, ARNT and FOXA1. Out of them, two SNPs were validated to impact transcriptional activity by performing luciferase assay that the alleles A of both SNPs, g.146704373A > G and g.146704618A > G, increased the transcriptional activities of AGPAT3 promoter compared with alleles G (P = 0.0004). CONCLUSIONS In conclusion, our findings first demonstrated the significant genetic associations of the AGPAT3 gene with milk FAs in dairy cattle, and two potential causal mutations were detected.
Collapse
Affiliation(s)
- Lijun Shi
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Road, Haidian District, Beijing, 100193, China
- Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Xin Wu
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Road, Haidian District, Beijing, 100193, China
| | - Yuze Yang
- Beijing General Station of Animal Husbandry, Beijing, 100101, China
| | - Zhu Ma
- Beijing Dairy Cattle Center, Beijing, 100192, China
| | - Xiaoqing Lv
- Beijing Dairy Cattle Center, Beijing, 100192, China
| | - Lin Liu
- Beijing Dairy Cattle Center, Beijing, 100192, China
| | - Yanhua Li
- Beijing Dairy Cattle Center, Beijing, 100192, China
| | - Feng Zhao
- Beijing Dairy Cattle Center, Beijing, 100192, China
| | - Bo Han
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Road, Haidian District, Beijing, 100193, China
| | - Dongxiao Sun
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Road, Haidian District, Beijing, 100193, China.
| |
Collapse
|
13
|
Chen SY, Oliveira HR, Schenkel FS, Pedrosa VB, Melka MG, Brito LF. Using imputed whole-genome sequence variants to uncover candidate mutations and genes affecting milking speed and temperament in Holstein cattle. J Dairy Sci 2020; 103:10383-10398. [PMID: 32952011 DOI: 10.3168/jds.2020-18897] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2020] [Accepted: 07/10/2020] [Indexed: 12/12/2022]
Abstract
Milking speed (MS) and temperament (MT) are 2 workability traits of great importance in dairy cattle production and breeding. This is mainly due to an increased intensification of the worldwide production systems and greater adoption of precision technologies with less human-cattle interaction. Both MS and MT are heritable traits and thus, genomic selection is a promising tool to expedite their genetic progress. However, the genetic architecture and biological mechanisms underlying the phenotypic expression of these traits remain underexplored. In this study, we investigated the association of >5.7 million imputed whole-genome sequence variants with MT and MS in 4,381 and 4,219 North American Holstein cattle, respectively. The statistical analyses were performed using a mixed linear model fitting a polygenic effect. We detected 40 and 35 significant SNPs independently associated with MT and MS, respectively, which were distributed across 26 chromosomes. Eight candidate genes (GRIN3A, KCNJ3, BOSTAUV1R417, BOSTAUV1R419, MAP2K5, KCTD3, GAP43, and LSAMP) were suggested to play an important role in MT as they are involved in biologically relevant pathways, such as glutamatergic synapse, vomeronasal receptor and oxytocin signaling. Within their coding and upstream sequences, we used an independent data set to further detect or validate significantly differentiated SNP between cattle breeds with known differences in MT. There were fewer candidate genes potentially implicated in MS, but immunity-related genes (e.g., BOLA-NC1 and LOC512672), also identified in other populations, were validated in this study. The significant SNP and novel candidate genes identified contribute to a better understanding of the biological mechanisms underlying both traits in dairy cattle. This information will also be useful for the optimization of prediction of genomic breeding values by giving greater weights to SNP located in the genomic regions identified.
Collapse
Affiliation(s)
- Shi-Yi Chen
- Department of Animal Sciences, Purdue University, West Lafayette, IN, 47907; Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan, 611130, China
| | - Hinayah R Oliveira
- Department of Animal Sciences, Purdue University, West Lafayette, IN, 47907; Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Flavio S Schenkel
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Victor B Pedrosa
- Department of Animal Sciences, State University of Ponta Grossa, Ponta Grossa, PR, 84030-900, Brazil
| | - Melkaye G Melka
- Department of Animal and Food Science, University of Wisconsin River Falls, 54022
| | - Luiz F Brito
- Department of Animal Sciences, Purdue University, West Lafayette, IN, 47907.
| |
Collapse
|
14
|
Shi L, Liu L, Lv X, Ma Z, Li C, Li Y, Zhao F, Sun D, Han B. Identification of genetic effects and potential causal polymorphisms of CPM gene impacting milk fatty acid traits in Chinese Holstein. Anim Genet 2020; 51:491-501. [PMID: 32301146 DOI: 10.1111/age.12936] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2019] [Revised: 02/03/2020] [Accepted: 03/15/2020] [Indexed: 11/27/2022]
Abstract
Our previous GWAS revealed 83 significant SNPs and 20 promising candidate genes associated with milk fatty acid traits in dairy cattle. Out of them, the carboxypeptidase M (CPM) gene contains a genome-wide significant SNP, Hapmap49848-BTA-106779, which is strongly associated with myristic acid (C14:0; P = 0.0064). Herein, we aimed to confirm the genetic effects of CPM on milk fatty acids in Chinese Holstein. Seven SNPs were detected by re-sequencing the sequences of entire exons and 3000 bp of up-/downstream flanking regions of the CPM gene, of which three were in 5' flanking region, one in the 3' UTR and three were in the 3' flanking region. Using the Haploview 4.1, we estimated the LD among the identified SNPs and found two haplotype blocks. With the animal model, we performed the SNP- and haplotype-based association analyses, and observed that these SNPs and haplotype blocks mainly had strong genetic associations with medium-chain saturated fatty acids (caproic acid, C6:0; caprylic acid, C8:0; capric acid, C10:0; and lauric acid, C12:0) (P < 0.0001-0.0257). In addition, using the Genomatix software, we predicted that three SNPs in the 5' flanking region of CPM (g.45079507A>G, g.45080228C>A and g.45080335C>G) changed the transcription factor binding sites for PREF (progesterone receptor biding site), ZBRK1 (transcription factor with eight central zinc fingers and an N-terminal KRAB domain), SOX9 (sex-determining region Y-box 9, dimeric binding sites), SOX6 (sex-determining region Y-box 6) and FOXP1-ES (alternative splicing variant of FOXP1, activated in ESCs). Further, the dual-luciferase reporter assay showed these three SNPs altered the transcriptional activity of CPM gene (P ≤ 0.0006). In summary, using the post-GWAS strategy, we first confirmed the significant genetic effects of CPM with milk fatty acids in dairy cattle, and identified three potential causal mutations.
Collapse
Affiliation(s)
- L Shi
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing, 100193, China.,Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - L Liu
- Beijing Dairy Cattle Center, Beijing, 100192, China
| | - X Lv
- Beijing Dairy Cattle Center, Beijing, 100192, China
| | - Z Ma
- Beijing Dairy Cattle Center, Beijing, 100192, China
| | - C Li
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing, 100193, China
| | - Y Li
- Beijing Dairy Cattle Center, Beijing, 100192, China
| | - F Zhao
- Beijing Dairy Cattle Center, Beijing, 100192, China
| | - D Sun
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing, 100193, China
| | - B Han
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing, 100193, China
| |
Collapse
|
15
|
Chen Q, Qu K, Ma Z, Zhan J, Zhang F, Shen J, Ning Q, Jia P, Zhang J, Chen N, Chen H, Huang B, Lei C. Genome-Wide Association Study Identifies Genomic Loci Associated With Neurotransmitter Concentration in Cattle. Front Genet 2020; 11:139. [PMID: 32292413 PMCID: PMC7119242 DOI: 10.3389/fgene.2020.00139] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2019] [Accepted: 02/06/2020] [Indexed: 11/30/2022] Open
Abstract
Abnormal neurotransmitter concentration is one of the factors that affect the health status, behavioral personality, and welfare level of animals, but the genetic basis of the abnormality is still largely unknown. The objective of this study is to identify putative genomic loci associated with neurotransmitter concentration in cattle. We measured serotonin (5HT), dopamine (DA), cortisol, glutamate (Glu), and ACTH concentrations in blood serum using double-antibody sandwich ELISA in 30 Brahman cattle and 127 Yunling cattle. Interestingly, we found that ACTH concentration was positively correlated with body weight, cannon circumference, and hip width (P < 0.05). Genome-wide association study (GWAS) was performed with mixed linear models using autosomal SNPs derived from the whole-genome sequence. We identified five, five, two, three, and five suggestive loci associated with 5HT, DA, cortisol, Glu, and ACTH concentration, respectively. These 20 associated loci implicated 18 candidate genes. For Glu concentration, the most significant association locus was assigned to MCHR1, a G-coupled receptor that could modulate glutamate release. For dopamine concentration, a very strong association locus was located in the intron of SLC18A2, which is a critical mediator of dopamine dynamics. However, for ACTH concentration, a very strong association locus was assigned to HTR1F, a G protein-coupled receptor that can influence the release of ACTH. Other candidate genes of interest identified for neurotransmitter concentration were PRMT6, GADD45A, PCCA, ANGPT1, ACCS, LOC100336971, TNR, GSDMA, CNTN3, CARMIL1, CDKAL1, RBFOX1, PCDH15, and LGALS12. Our findings will provide targets for the genetic improvement of neurotransmitter-related traits in domestic cattle and basic materials for studying the mechanism of neurotransmitter synthesis, release, and transport in human and animals.
Collapse
Affiliation(s)
- Qiuming Chen
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Kaixing Qu
- Yunnan Academy of Grassland and Animal Science, Kunming, China
| | - Zhijie Ma
- Academy of Animal Science and Veterinary Medicine, Qinghai University, Xining, China
| | - Jingxi Zhan
- Yunnan Academy of Grassland and Animal Science, Kunming, China
| | - Fengwei Zhang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Jiafei Shen
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Qingqing Ning
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Peng Jia
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Jicai Zhang
- Yunnan Academy of Grassland and Animal Science, Kunming, China
| | - Ningbo Chen
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Hong Chen
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Bizhi Huang
- Yunnan Academy of Grassland and Animal Science, Kunming, 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, China
| |
Collapse
|
16
|
Wang T, Li J, Gao X, Song W, Chen C, Yao D, Ma J, Xu L, Ma Y. Genome-wide association study of milk components in Chinese Holstein cows using single nucleotide polymorphism. Livest Sci 2020. [DOI: 10.1016/j.livsci.2020.103951] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
17
|
Shi L, Liu L, Lv X, Ma Z, Yang Y, Li Y, Zhao F, Sun D, Han B. Polymorphisms and genetic effects of PRLR, MOGAT1, MINPP1 and CHUK genes on milk fatty acid traits in Chinese Holstein. BMC Genet 2019; 20:69. [PMID: 31419940 PMCID: PMC6698030 DOI: 10.1186/s12863-019-0769-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2018] [Accepted: 08/06/2019] [Indexed: 01/28/2023] Open
Abstract
BACKGROUND Our initial genome-wide association study (GWAS) identified 20 promising candidate genes for milk fatty acid (FA) traits in a Chinese Holstein population, including PRLR, MOGAT1, MINPP1 and CHUK genes. In this study, we performed whether they had significant genetic effects on milk FA traits in Chinese Holstein. RESULTS We re-sequenced the entire exons and 3000 bp of the 5' and 3' flanking regions, and identified 11 single nucleotide polymorphisms (SNPs), containing four in PRLR, two in MOGAT1, two in MINPP1, and three in CHUK. The SNP-based association analyses showed that all the 11 SNPs were significantly associated with at least one milk FA trait (P = 0.0456 ~ < 0.0001), and none of them had association with C11:0, C13:0, C15:0 and C16:0 (P > 0.05). By the linkage disequilibrium (LD) analyses, we found two, one, one, and one haplotype blocks in PRLR, MOGAT1, MINPP1, and CHUK, respectively, and each haplotype block was significantly associated with at least one milk FA trait (P = 0.0456 ~ < 0.0001). Further, g.38949011G > A in PRLR, and g.111599360A > G and g.111601747 T > A in MOGAT1 were predicted to alter the transcription factor binding sites (TFBSs). A missense mutation, g.39115344G > A, could change the PRLR protein structure. The g.20966385C > G of CHUK varied the binding sequences for microRNAs. Therefore, we deduced the five SNPs as the potential functional mutations. CONCLUSION In summary, we first detected the genetic effects of PRLR, MOGAT1, MINPP1 and CHUK genes on milk FA traits, and researched the potential functional mutations. These data provided the basis for further investigation on function validation of the four genes in Chinese Holstein.
Collapse
Affiliation(s)
- Lijun Shi
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Road, Haidian District, Beijing, 100193 China
| | - Lin Liu
- Beijing Dairy Cattle Center, Beijing, 100192 China
| | - Xiaoqing Lv
- Beijing Dairy Cattle Center, Beijing, 100192 China
| | - Zhu Ma
- Beijing Dairy Cattle Center, Beijing, 100192 China
| | - Yuze Yang
- Beijing General Station of Animal Husbandry, Beijing, 100101 China
| | - Yanhua Li
- Beijing Dairy Cattle Center, Beijing, 100192 China
| | - Feng Zhao
- Beijing Dairy Cattle Center, Beijing, 100192 China
| | - Dongxiao Sun
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Road, Haidian District, Beijing, 100193 China
| | - Bo Han
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Road, Haidian District, Beijing, 100193 China
| |
Collapse
|
18
|
Shi L, Liu L, Ma Z, Lv X, Li C, Xu L, Han B, Li Y, Zhao F, Yang Y, Sun D. Identification of genetic associations of ECHS1 gene with milk fatty acid traits in dairy cattle. Anim Genet 2019; 50:430-438. [PMID: 31392738 DOI: 10.1111/age.12833] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/15/2019] [Indexed: 12/20/2022]
Abstract
Our previous genome-wide association study identified 83 genome-wide significant SNPs and 20 novel promising candidate genes for milk fatty acids in Chinese Holstein. Among them, the enoyl-CoA hydratase, short chain 1 (ECHS1) and enoyl-CoA hydratase and 3-hydroxyacyl CoA dehydrogenase (EHHADH) genes were located near two SNPs and one SNP respectively, and they play important roles in fatty acid metabolism pathways. We herein validated whether the two genes have genetic effects on milk fatty acid traits in dairy cattle. By re-sequencing the full-length coding region, partially adjacent introns and 3000 bp up/downstream flanking sequences, we identified 12 SNPs in ECHS1: two in exons, four in the 3' flanking region and six in introns. The g.25858322C>T SNP results in an amino acid replacement from leucine to phenylalanine and changes the secondary structure of the ECHS1 protein, and single-locus association analysis showed that it was significantly associated with three milk fatty acids (P = 0.0002-0.0013). The remaining 11 SNPs were found to be significantly associated with at least one milk fatty acid (P = <0.0001-0.0040). Also, we found that two haplotype blocks, consisting of nine and two SNPs respectively, were significantly associated with eight milk fatty acids (P = <0.0001-0.0125). However, none of polymorphisms was observed in the EHHADH gene. In conclusion, our findings are the first to indicate that the ECHS1 gene has a significant genetic impact on long-chain unsaturated and medium-chain saturated fatty acid traits in dairy cattle, although the biological mechanism is still undetermined and requires further in-depth validation.
Collapse
Affiliation(s)
- L Shi
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing, 100193, China
| | - L Liu
- Beijing Dairy Cattle Center, Beijing, 100192, China
| | - Z Ma
- Beijing Dairy Cattle Center, Beijing, 100192, China
| | - X Lv
- Beijing Dairy Cattle Center, Beijing, 100192, China
| | - C Li
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing, 100193, China
| | - L Xu
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing, 100193, China
| | - B Han
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing, 100193, China
| | - Y Li
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing, 100193, China.,Beijing Dairy Cattle Center, Beijing, 100192, China
| | - F Zhao
- Beijing Dairy Cattle Center, Beijing, 100192, China
| | - Y Yang
- Beijing Municipal Bureau of Agriculture, Beijing, 100101, China
| | - D Sun
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing, 100193, China
| |
Collapse
|
19
|
Shi L, Lv X, Liu L, Yang Y, Ma Z, Han B, Sun D. A post-GWAS confirming effects of PRKG1 gene on milk fatty acids in a Chinese Holstein dairy population. BMC Genet 2019; 20:53. [PMID: 31269900 PMCID: PMC6610796 DOI: 10.1186/s12863-019-0755-7] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2019] [Accepted: 06/20/2019] [Indexed: 01/03/2023] Open
Abstract
BACKGROUND We previously conducted a genome-wide association study (GWAS) strategy for milk fatty acids in Chinese Holstein, and identified 83 genome-wide significant single nucleotide polymorphisms (SNPs) and 314 suggestive significant SNPs. Among them, two SNPs, BTB-01077939 and BTA-11275-no-rs associated with C10:0, C12:0, and C14 index (P = 0.000014 ~ 0.000024), were within and close to (0.85 Mb) protein kinase, cGMP-dependent, type І (PRKG1) gene on BTA26, respectively. PRKG1 gene plays a key role in lipolysis to release fatty acids and glycerol through the hydrolysis of triacyglycerol in adipocytes. We herein considered it as a promising candidate for milk fatty acids. The purpose of this study was to investigate whether PRKG1 had effects on milk fatty acids. RESULTS By direct sequencing the PCR products of pooled DNA, we identified a total of six SNPs, including one in 5' flanking region, four in 3' untranslated region (UTR), and one in 3' flanking region. The single-locus association analysis was carried out, and showed that the six SNPs mainly had significant associations with C6:0, C8:0 and C17:1 (P < 0.0001 ~ 0.0035). In addition, we observed a haplotype block formed by g.6903810G > A and g.6904047G > T with Haploview 4.1, and it was strongly associated with C8:0, C10:0, C16:1, C17:1, C20:0 and C16 index (P = < 0.0001 ~ 0.0123). The SNP, g.8344262A > T, was predicted to alter the binding site (BS) of transcription factor (TF) GAGA box with Genomatix software, and the subsequent luciferase assay verified that it really changed the transcriptional activity of PRKG1 gene (P = 0.0009). CONCLUSION In conclusion, to our best of knowledge, we are the first who identified the significant effects of PRKG1 on milk fatty acids in dairy cattle.
Collapse
Affiliation(s)
- Lijun Shi
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Road, Haidian District, Beijing, 100193 China
| | - Xiaoqing Lv
- Beijing Dairy Cattle Center, Beijing, 100192 China
| | - Lin Liu
- Beijing Dairy Cattle Center, Beijing, 100192 China
| | - Yuze Yang
- Beijing Municipal Bureau of Agriculture, Beijing, 100101 China
| | - Zhu Ma
- Beijing Dairy Cattle Center, Beijing, 100192 China
| | - Bo Han
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Road, Haidian District, Beijing, 100193 China
| | - Dongxiao Sun
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, No. 2 Yuanmingyuan West Road, Haidian District, Beijing, 100193 China
| |
Collapse
|
20
|
Poulsen NA, Robinson RC, Barile D, Larsen LB, Buitenhuis B. A genome-wide association study reveals specific transferases as candidate loci for bovine milk oligosaccharides synthesis. BMC Genomics 2019; 20:404. [PMID: 31117955 PMCID: PMC6532250 DOI: 10.1186/s12864-019-5786-y] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2018] [Accepted: 05/08/2019] [Indexed: 12/11/2022] Open
Abstract
BACKGROUND Human milk oligosaccharides (OS) play a key role in brain and gut microbiota development of the neonate, but the underlying biosynthetic steps of OS in the mammary gland are still largely unknown. As bovine milk contains OS with somewhat similar structures and functionalities there is increased interest in further understanding the genetic basis underlying the OS content of milk for eventual extraction and generation of value-added ingredients for infant formulas and nutraceuticals. The present study is the first to report on genetic parameter estimation as well as on a genome wide association study (GWAS) from the largest bovine milk OS dataset analyzed to date. RESULTS In total 15 different bovine milk OS were monitored. Heritabilities ranged from 0 to 0.68 in Danish Holstein and from 0 to 0.92 in Danish Jersey. The GWAS identified in total 1770 SNPs (FDR < 0.10) for five different OS in Danish Holstein and 6913 SNPs (FDR < 0.10) for 11 OS in Danish Jersey. In Danish Holstein, a major overlapping QTL was identified on BTA1 for LNH and LNT explaining 24% of the variation in these OS. The most significant SNPs were associated with B3GNT5, a gene encoding a glycosyltransferase involved in glycan synthesis. In Danish Jersey, a very strong QTL was detected for the OS with composition 2 Hex 1 HexNAc (isomer 1) on BTA11. The most significant SNP had -log10(P-value) of 52.88 (BOVINEHD1100030300) and was assigned to ABO, a gene encoding ABO blood group glycosyltransferases. This SNP has been reported to be a missense mutation and explains 56% of the OS variation. Other candidate genes of interest identified for milk OS were ALG3, B3GALNT2, LOC520336, PIGV, MAN1C1, ST6GALNAC6, GLT6D1, GALNT14, GALNT17, COLGALT2, LFNG and SIGLEC. CONCLUSION To our knowledge, this is the first study documenting a solid breeding potential for bovine milk OS and a strong indication of specific candidate genes related to OS synthesis underlying this genetic influence. This new information has the potential to guide breeding strategies to achieve production of milk with higher diversity and concentration of OS and ultimately facilitate large-scale extraction of bovine milk OS.
Collapse
Affiliation(s)
- Nina A. Poulsen
- Department of Food Science, Aarhus University, Blichers Allé 20, P. O. Box 50, DK-8830 Tjele, Denmark
| | - Randall C. Robinson
- Department of Food Science and Technology, University of California, Davis, One Shields Ave., Davis, CA 95616 USA
| | - Daniela Barile
- Department of Food Science and Technology, University of California, Davis, One Shields Ave., Davis, CA 95616 USA
- Foods for Health Institute, University of California, Davis, One Shields Ave., Davis, CA 95616 USA
| | - Lotte B. Larsen
- Department of Food Science, Aarhus University, Blichers Allé 20, P. O. Box 50, DK-8830 Tjele, Denmark
| | - Bart Buitenhuis
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, Blichers Allé 20, P. O. Box 50, DK-8830 Tjele, Denmark
| |
Collapse
|
21
|
Gebreyesus G, Bovenhuis H, Lund MS, Poulsen NA, Sun D, Buitenhuis B. Reliability of genomic prediction for milk fatty acid composition by using a multi-population reference and incorporating GWAS results. Genet Sel Evol 2019; 51:16. [PMID: 31029078 PMCID: PMC6487064 DOI: 10.1186/s12711-019-0460-z] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2018] [Accepted: 04/10/2019] [Indexed: 01/01/2023] Open
Abstract
Background Large-scale phenotyping for detailed milk fatty acid (FA) composition is difficult due to expensive and time-consuming analytical techniques. Reliability of genomic prediction is often low for traits that are expensive/difficult to measure and for breeds with a small reference population size. An effective method to increase reference population size could be to combine datasets from different populations. Prediction models might also benefit from incorporation of information on the biological underpinnings of quantitative traits. Genome-wide association studies (GWAS) show that genomic regions on Bos taurus chromosomes (BTA) 14, 19 and 26 underlie substantial proportions of the genetic variation in milk FA traits. Genomic prediction models that incorporate such results could enable improved prediction accuracy in spite of limited reference population sizes. In this study, we combine gas chromatography quantified FA samples from the Chinese, Danish and Dutch Holstein populations and implement a genomic feature best linear unbiased prediction (GFBLUP) model that incorporates variants on BTA14, 19 and 26 as genomic features for which random genetic effects are estimated separately. Prediction reliabilities were compared to those estimated with traditional GBLUP models. Results Predictions using a multi-population reference and a traditional GBLUP model resulted in average gains in prediction reliability of 10% points in the Dutch, 8% points in the Danish and 1% point in the Chinese populations compared to predictions based on population-specific references. Compared to the traditional GBLUP, implementation of the GFBLUP model with a multi-population reference led to further increases in prediction reliability of up to 38% points in the Dutch, 23% points in the Danish and 13% points in the Chinese populations. Prediction reliabilities from the GFBLUP model were moderate to high across the FA traits analyzed. Conclusions Our study shows that it is possible to predict genetic merits for milk FA traits with reasonable accuracy by combining related populations of a breed and using models that incorporate GWAS results. Our findings indicate that international collaborations that facilitate access to multi-population datasets could be highly beneficial to the implementation of genomic selection for detailed milk composition traits. Electronic supplementary material The online version of this article (10.1186/s12711-019-0460-z) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Grum Gebreyesus
- Department of Molecular Biology and Genetics, Center for Quantitative Genetics and Genomics, Aarhus University, Blichers Allé 20, P.O. Box 50, 8830, Tjele, Denmark. .,Animal Breeding and Genomics Centre, Wageningen University, P.O. Box 338, 6700 AH, Wageningen, The Netherlands.
| | - Henk Bovenhuis
- Animal Breeding and Genomics Centre, Wageningen University, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| | - Mogens S Lund
- Department of Molecular Biology and Genetics, Center for Quantitative Genetics and Genomics, Aarhus University, Blichers Allé 20, P.O. Box 50, 8830, Tjele, Denmark
| | - Nina A Poulsen
- Department of Food Science, Aarhus University, Blichers Allé 20, P.O. Box 50, 8830, Tjele, Denmark
| | - Dongxiao Sun
- Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture of China, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Bart Buitenhuis
- Department of Molecular Biology and Genetics, Center for Quantitative Genetics and Genomics, Aarhus University, Blichers Allé 20, P.O. Box 50, 8830, Tjele, Denmark
| |
Collapse
|
22
|
Iung LHS, Petrini J, Ramírez-Díaz J, Salvian M, Rovadoscki GA, Pilonetto F, Dauria BD, Machado PF, Coutinho LL, Wiggans GR, Mourão GB. Genome-wide association study for milk production traits in a Brazilian Holstein population. J Dairy Sci 2019; 102:5305-5314. [PMID: 30904307 DOI: 10.3168/jds.2018-14811] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2018] [Accepted: 10/19/2018] [Indexed: 12/19/2022]
Abstract
Advances in the molecular area of selection have expanded knowledge of the genetic architecture of complex traits through genome-wide association studies (GWAS). Several GWAS have been performed so far, but confirming these results is not always possible due to several factors, including environmental conditions. Thus, our objective was to identify genomic regions associated with traditional milk production traits, including milk yield, somatic cell score, fat, protein and lactose percentages, and fatty acid composition in a Holstein cattle population producing under tropical conditions. For this, 75,228 phenotypic records from 5,981 cows and genotypic data of 56,256 SNP from 1,067 cows were used in a weighted single-step GWAS. A total of 46 windows of 10 SNP explaining more than 1% of the genetic variance across 10 Bos taurus autosomes (BTA) harbored well-known and novel genes. The MGST1 (BTA5), ABCG2 (BTA6), DGAT1 (BTA14), and PAEP (BTA11) genes were confirmed within some of the regions identified in our study. Potential novel genes involved in tissue damage and repair of the mammary gland (COL18A1), immune response (LTTC19), glucose homeostasis (SLC37A1), synthesis of unsaturated fatty acids (LTBP1), and sugar transport (SLC37A1 and MFSD4A) were found for milk yield, somatic cell score, fat percentage, and fatty acid composition. Our findings may assist genomic selection by using these regions to design a customized SNP array to improve milk production traits on farms with similar environmental conditions.
Collapse
Affiliation(s)
- L H S Iung
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - J Petrini
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - J Ramírez-Díaz
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - M Salvian
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - G A Rovadoscki
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - F Pilonetto
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - B D Dauria
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - P F Machado
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - L L Coutinho
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil
| | - G R Wiggans
- Animal Genomics and Improvement Laboratory, Agricultural Research Service, USDA, Beltsville, MD 20705-2350
| | - G B Mourão
- Department of Animal Science, University of São Paulo (USP)/Luiz de Queiroz College of Agriculture (ESALQ), Piracicaba, São Paulo 13418900, Brazil.
| |
Collapse
|
23
|
Gebreyesus G, Buitenhuis AJ, Poulsen NA, Visker MHPW, Zhang Q, van Valenberg HJF, Sun D, Bovenhuis H. Multi-population GWAS and enrichment analyses reveal novel genomic regions and promising candidate genes underlying bovine milk fatty acid composition. BMC Genomics 2019; 20:178. [PMID: 30841852 PMCID: PMC6404302 DOI: 10.1186/s12864-019-5573-9] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2018] [Accepted: 02/28/2019] [Indexed: 01/23/2023] Open
Abstract
Background The power of genome-wide association studies (GWAS) is often limited by the sample size available for the analysis. Milk fatty acid (FA) traits are scarcely recorded due to expensive and time-consuming analytical techniques. Combining multi-population datasets can enhance the power of GWAS enabling detection of genomic region explaining medium to low proportions of the genetic variation. GWAS often detect broader genomic regions containing several positional candidate genes making it difficult to untangle the causative candidates. Post-GWAS analyses with data on pathways, ontology and tissue-specific gene expression status might allow prioritization among positional candidate genes. Results Multi-population GWAS for 16 FA traits quantified using gas chromatography (GC) in sample populations of the Chinese, Danish and Dutch Holstein with high-density (HD) genotypes detects 56 genomic regions significantly associated to at least one of the studied FAs; some of which have not been previously reported. Pathways and gene ontology (GO) analyses suggest promising candidate genes on the novel regions including OSBPL6 and AGPS on Bos taurus autosome (BTA) 2, PRLH on BTA 3, SLC51B on BTA 10, ABCG5/8 on BTA 11 and ALG5 on BTA 12. Novel genes in previously known regions, such as FABP4 on BTA 14, APOA1/5/7 on BTA 15 and MGST2 on BTA 17, are also linked to important FA metabolic processes. Conclusion Integration of multi-population GWAS and enrichment analyses enabled detection of several novel genomic regions, explaining relatively smaller fractions of the genetic variation, and revealed highly likely candidate genes underlying the effects. Detection of such regions and candidate genes will be crucial in understanding the complex genetic control of FA metabolism. The findings can also be used to augment genomic prediction models with regions collectively capturing most of the genetic variation in the milk FA traits. Electronic supplementary material The online version of this article (10.1186/s12864-019-5573-9) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- G Gebreyesus
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, Blichers Allé 20, P.O. Box 50, DK-8830, Tjele, Denmark. .,Animal Breeding and Genomics, Wageningen University and Research, P.O. Box 338, 6700 AH, Wageningen, the Netherlands.
| | - A J Buitenhuis
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, Blichers Allé 20, P.O. Box 50, DK-8830, Tjele, Denmark
| | - N A Poulsen
- Department of Food Science, Aarhus University, Blichers Allé 20, P.O. Box 50, DK-8830, Tjele, Denmark
| | - M H P W Visker
- Animal Breeding and Genomics, Wageningen University and Research, P.O. Box 338, 6700 AH, Wageningen, the Netherlands
| | - Q Zhang
- Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture of China, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - H J F van Valenberg
- Dairy Science and Technology Group, Wageningen University and Research, P.O. Box 17, 6700 AA, Wageningen, the Netherlands
| | - D Sun
- Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture of China, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - H Bovenhuis
- Animal Breeding and Genomics, Wageningen University and Research, P.O. Box 338, 6700 AH, Wageningen, the Netherlands
| |
Collapse
|
24
|
Wang Z, Zhu B, Niu H, Zhang W, Xu L, Xu L, Chen Y, Zhang L, Gao X, Gao H, Zhang S, Xu L, Li J. Genome wide association study identifies SNPs associated with fatty acid composition in Chinese Wagyu cattle. J Anim Sci Biotechnol 2019; 10:27. [PMID: 30867906 PMCID: PMC6399853 DOI: 10.1186/s40104-019-0322-0] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2018] [Accepted: 01/14/2019] [Indexed: 12/29/2022] Open
Abstract
Background Fatty acids are important traits that affect meat quality and nutritive values in beef cattle. Detection of genetic variants for fatty acid composition can help to elucidate the genetic mechanism underpinning these traits and promote the improvement of fatty acid profiles. In this study, we performed a genome-wide association study (GWAS) on fatty acid composition using high-density single nucleotide polymorphism (SNP) arrays in Chinese Wagyu cattle. Results In total, we detected 15 and 8 significant genome-wide SNPs for individual fatty acids and fatty acid groups in Chinese Wagyu cattle, respectively. Also, we identified nine candidate genes based on 100 kb regions around associated SNPs. Four SNPs significantly associated with C14:1 cis-9 were embedded with stearoyl-CoA desaturase (SCD), while three SNPs in total were identified for C22:6 n-3 within Phospholipid scramblase family member 5 (PLSCR5), Cytoplasmic linker associated protein 1 (CLASP1), and Chymosin (CYM). Notably, we found the top candidate SNP within SCD can explain ~ 7.37% of phenotypic variance for C14:1 cis-9. Moreover, we detected several blocks with high LD in the 100 kb region around SCD. In addition, we found three significant SNPs within a 100 kb region showing pleiotropic effects related to multiple FA groups (PUFA, n-6, and PUFA/SFA), which contains BAI1 associated protein 2 like 2 (BAIAP2L2), MAF bZIP transcription factor F (MAFF), and transmembrane protein 184B (TMEM184B). Conclusions Our study identified several significant SNPs and candidate genes for individual fatty acids and fatty acid groups in Chinese Wagyu cattle, and these findings will further assist the design of breeding programs for meat quality in cattle.
Collapse
Affiliation(s)
- Zezhao Wang
- 1Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193 China.,2National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Bo Zhu
- 1Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193 China
| | - Hong Niu
- 1Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193 China
| | - Wengang Zhang
- 1Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193 China
| | - Ling Xu
- 1Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193 China
| | - Lei Xu
- 1Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193 China.,3Institute of Animal Husbandry and Veterinary Research, Anhui Academy of Agricultural Sciences, Hefei, 230031 China
| | - Yan Chen
- 1Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193 China
| | - Lupei Zhang
- 1Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193 China
| | - Xue Gao
- 1Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193 China
| | - Huijiang Gao
- 1Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193 China
| | - Shengli Zhang
- 2National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Lingyang Xu
- 1Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193 China
| | - Junya Li
- 1Innovation Team of Cattle Genetic Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193 China
| |
Collapse
|
25
|
Shi L, Han B, Liu L, Lv X, Ma Z, Li C, Xu L, Li Y, Zhao F, Yang Y, Sun D. Determination of Genetic Effects of LIPK and LIPJ Genes on Milk Fatty Acids in Dairy Cattle. Genes (Basel) 2019; 10:genes10020086. [PMID: 30696079 PMCID: PMC6409763 DOI: 10.3390/genes10020086] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2018] [Revised: 01/13/2019] [Accepted: 01/14/2019] [Indexed: 12/30/2022] Open
Abstract
In our previous genome-wide association study (GWAS) on milk fatty acids (FAs) in Chinese Holstein, we discovered 83 genome-wide significant single nucleotide polymorphisms (SNPs) associated with milk FAs. Two of them were close to lipase family member K (LIPK) and lipase family member J (LIPJ), respectively. Hence, this study is a follow-up to verify whether the LIPK and LIPJ have significant genetic effects on milk FAs in dairy cattle. By re-sequencing the entire exons, and 3 kb of 5′ and 3′ flanking regions, two and seven SNPs were identified in LIPK and LIPJ, respectively, including a novel SNP, ss158213049726. With the Haploview 4.1 software, we found that five of the SNPs in LIPJ formed a haplotype block (D′ = 0.96 ~ 1.00). Single-locus association analyses revealed that each SNP in LIPK and LIPJ was significantly associated with at least one milk FA (p = < 1.00 × 10−4 ~ 4.88 × 10−2), and the haplotype-based association analyses showed significant genetic effects on nine milk FAs (p = < 1.00 × 10−4 ~ 3.98 × 10−2). Out of these SNPs, the missense mutation in LIPK gene, rs42774527, could change the protein secondary structure and function predicted by SOPMA, SIFT, and PROVEAN softwares. With the Genomatix software, we predicted that two SNPs, rs110322221 in LIPK and rs211373799 in LIPJ, altered the transcription factors binding sites (TFBSs), indicating their potential regulation on promoter activity of the genes. Furthermore, we found that both LIPK and LIPJ had relatively high expressions in the mammary gland. In conclusion, our research is the first to demonstrate that LIPK and LIPJ genes have significant associations with milk FAs, and the identified SNPs might be served as genetic markers to optimize breeding programs for milk FAs in dairy cattle. This research deserves in-depth verification.
Collapse
Affiliation(s)
- Lijun Shi
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing 100193, China.
| | - Bo Han
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing 100193, China.
| | - Lin Liu
- Beijing Dairy Cattle Center, Qinghe'nanzhen Deshengmenwai Street, Chaoyang District, Beijing 100192, China.
| | - Xiaoqing Lv
- Beijing Dairy Cattle Center, Qinghe'nanzhen Deshengmenwai Street, Chaoyang District, Beijing 100192, China.
| | - Zhu Ma
- Beijing Dairy Cattle Center, Qinghe'nanzhen Deshengmenwai Street, Chaoyang District, Beijing 100192, China.
| | - Cong Li
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing 100193, China.
| | - Lingna Xu
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing 100193, China.
| | - Yanhua Li
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing 100193, China.
- Beijing Dairy Cattle Center, Qinghe'nanzhen Deshengmenwai Street, Chaoyang District, Beijing 100192, China.
| | - Feng Zhao
- Beijing Dairy Cattle Center, Qinghe'nanzhen Deshengmenwai Street, Chaoyang District, Beijing 100192, China.
| | - Yuze Yang
- Beijing General Station of Animal Husbandry, N0.96 Huizhongsi, Yayun Village, Chaoyang District, Beijing, 100101, China.
| | - Dongxiao Sun
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science and Technology, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, National Engineering Laboratory for Animal Breeding, China Agricultural University, Beijing 100193, China.
| |
Collapse
|
26
|
Kroezen V, Schenkel F, Miglior F, Baes C, Squires E. Candidate gene association analyses for ketosis resistance in Holsteins. J Dairy Sci 2018; 101:5240-5249. [DOI: 10.3168/jds.2017-13374] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2017] [Accepted: 02/14/2018] [Indexed: 11/19/2022]
|
27
|
SNP co-association and network analyses identify E2F3, KDM5A and BACH2 as key regulators of the bovine milk fatty acid profile. Sci Rep 2017; 7:17317. [PMID: 29230020 PMCID: PMC5725496 DOI: 10.1038/s41598-017-17434-7] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2017] [Accepted: 11/27/2017] [Indexed: 12/19/2022] Open
Abstract
The fatty acid (FA) profile has a considerable impact on the nutritional and technological quality of milk and dairy products. The molecular mechanism underlying the regulation of fat metabolism in bovine mammary gland have been not completely elucidated. We conducted genome-wide association studies (GWAS) across 65 milk FAs and fat percentage in 1,152 Brown Swiss cows. In total, we identified 175 significant single nucleotide polymorphism (SNPs) spanning all chromosomes. Pathway analyses revealed that 12:0 was associated with the greatest number of overrepresented categories/pathways (e.g. mitogen-activated protein kinase (MAPK) activity and protein phosphorylation), suggesting that it might play an important biological role in controlling milk fat composition. An Associated Weight Matrix approach based on SNP co-associations predicted a network of 791 genes related to the milk FA profile, which were involved in several connected molecular pathways (e.g., MAPK, lipid metabolism and hormone signalling) and undetectable through standard GWAS. Analysis of transcription factors and their putative target genes within the network identified BACH2, E2F3 and KDM5A as key regulators of milk FA metabolism. These findings contribute to increasing knowledge of FA metabolism and mammary gland functionality in dairy cows and may be useful in developing targeted breeding practices to improve milk quality.
Collapse
|
28
|
Liu A, Wang Y, Sahana G, Zhang Q, Liu L, Lund MS, Su G. Genome-wide Association Studies for Female Fertility Traits in Chinese and Nordic Holsteins. Sci Rep 2017; 7:8487. [PMID: 28814769 PMCID: PMC5559619 DOI: 10.1038/s41598-017-09170-9] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2017] [Accepted: 07/20/2017] [Indexed: 12/18/2022] Open
Abstract
Reduced female fertility could cause considerable economic loss and has become a worldwide problem in the modern dairy industry. The objective of this study was to detect quantitative trait loci (QTL) for female fertility traits in Chinese and Nordic Holsteins using various strategies. First, single-trait association analyses were performed for female fertility traits in Chinese and Nordic Holsteins. Second, the SNPs with P-value < 0.005 discovered in Chinese Holsteins were validated in Nordic Holsteins. Third, the summary statistics from single-trait association analyses were combined into meta-analyses to: (1) identify common QTL for multiple fertility traits within each Holstein population; (2) detect SNPs which were associated with a female fertility trait across two Holstein populations. A large numbers of QTL were discovered or confirmed for female fertility traits. The QTL segregating at 31.4~34.1 Mb on BTA13, 48.3~51.9 Mb on BTA23 and 34.0~37.6 Mb on BTA28 shared between Chinese and Nordic Holsteins were further ascertained using a validation approach and meta-analyses. Furthermore, multiple novel variants identified in Chinese Holsteins were validated with Nordic data as well as meta-analyses. The genes IL6R, SLC39A12, CACNB2, ZEB1, ZMIZ1 and FAM213A were concluded to be strong candidate genes for female fertility in Holsteins.
Collapse
Affiliation(s)
- Aoxing Liu
- Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture of China, National Engineering Laboratory of Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China.,Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830, Tjele, Denmark
| | - Yachun Wang
- Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture of China, National Engineering Laboratory of Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China.
| | - Goutam Sahana
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830, Tjele, Denmark
| | - Qin Zhang
- Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture of China, National Engineering Laboratory of Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Lin Liu
- Beijing Dairy Cattle Center, Beijing, 100192, China
| | - Mogens Sandø Lund
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830, Tjele, Denmark
| | - Guosheng Su
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830, Tjele, Denmark.
| |
Collapse
|
29
|
Martin P, Palhière I, Maroteau C, Bardou P, Canale-Tabet K, Sarry J, Woloszyn F, Bertrand-Michel J, Racke I, Besir H, Rupp R, Tosser-Klopp G. A genome scan for milk production traits in dairy goats reveals two new mutations in Dgat1 reducing milk fat content. Sci Rep 2017; 7:1872. [PMID: 28500343 PMCID: PMC5431851 DOI: 10.1038/s41598-017-02052-0] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2016] [Accepted: 04/05/2017] [Indexed: 11/21/2022] Open
Abstract
The quantity of milk and milk fat and proteins are particularly important traits in dairy livestock. However, little is known about the regions of the genome that influence these traits in goats. We conducted a genome wide association study in French goats and identified 109 regions associated with dairy traits. For a major region on chromosome 14 closely associated with fat content, the Diacylglycerol O-Acyltransferase 1 (DGAT1) gene turned out to be a functional and positional candidate gene. The caprine reference sequence of this gene was completed and 29 polymorphisms were found in the gene sequence, including two novel exonic mutations: R251L and R396W, leading to substitutions in the protein sequence. The R251L mutation was found in the Saanen breed at a frequency of 3.5% and the R396W mutation both in the Saanen and Alpine breeds at a frequencies of 13% and 7% respectively. The R396W mutation explained 46% of the genetic variance of the trait, and the R251L mutation 6%. Both mutations were associated with a notable decrease in milk fat content. Their causality was then demonstrated by a functional test. These results provide new knowledge on the genetic basis of milk synthesis and will help improve the management of the French dairy goat breeding program.
Collapse
Affiliation(s)
- Pauline Martin
- GenPhySE, Université de Toulouse, INRA, INPT, ENVT, Castanet Tolosan, France
| | - Isabelle Palhière
- GenPhySE, Université de Toulouse, INRA, INPT, ENVT, Castanet Tolosan, France
| | - Cyrielle Maroteau
- GenPhySE, Université de Toulouse, INRA, INPT, ENVT, Castanet Tolosan, France
- Division of Molecular and Clinical Medecine, School of Medecine, University of Dundee, Ninewells Hospital and Medical School, Dundee, UK
| | - Philippe Bardou
- GenPhySE, Université de Toulouse, INRA, INPT, ENVT, Castanet Tolosan, France
- INRA, Sigenae, Castanet-Tolosan, France
| | - Kamila Canale-Tabet
- GenPhySE, Université de Toulouse, INRA, INPT, ENVT, Castanet Tolosan, France
| | - Julien Sarry
- GenPhySE, Université de Toulouse, INRA, INPT, ENVT, Castanet Tolosan, France
| | - Florent Woloszyn
- GenPhySE, Université de Toulouse, INRA, INPT, ENVT, Castanet Tolosan, France
| | | | - Ines Racke
- Protein Expression and Purification Core Facility, EMBL Heidelberg, Heidelberg, Germany
| | - Hüseyin Besir
- Protein Expression and Purification Core Facility, EMBL Heidelberg, Heidelberg, Germany
| | - Rachel Rupp
- GenPhySE, Université de Toulouse, INRA, INPT, ENVT, Castanet Tolosan, France
| | | |
Collapse
|
30
|
Li X, Lund MS, Janss L, Wang C, Ding X, Zhang Q, Su G. The patterns of genomic variances and covariances across genome for milk production traits between Chinese and Nordic Holstein populations. BMC Genet 2017; 18:26. [PMID: 28298201 PMCID: PMC5353867 DOI: 10.1186/s12863-017-0491-9] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2016] [Accepted: 03/07/2017] [Indexed: 02/06/2023] Open
Abstract
BACKGROUND With the development of SNP chips, SNP information provides an efficient approach to further disentangle different patterns of genomic variances and covariances across the genome for traits of interest. Due to the interaction between genotype and environment as well as possible differences in genetic background, it is reasonable to treat the performances of a biological trait in different populations as different but genetic correlated traits. In the present study, we performed an investigation on the patterns of region-specific genomic variances, covariances and correlations between Chinese and Nordic Holstein populations for three milk production traits. RESULTS Variances and covariances between Chinese and Nordic Holstein populations were estimated for genomic regions at three different levels of genome region (all SNP as one region, each chromosome as one region and every 100 SNP as one region) using a novel multi-trait random regression model which uses latent variables to model heterogeneous variance and covariance. In the scenario of the whole genome as one region, the genomic variances, covariances and correlations obtained from the new multi-trait Bayesian method were comparable to those obtained from a multi-trait GBLUP for all the three milk production traits. In the scenario of each chromosome as one region, BTA 14 and BTA 5 accounted for very large genomic variance, covariance and correlation for milk yield and fat yield, whereas no specific chromosome showed very large genomic variance, covariance and correlation for protein yield. In the scenario of every 100 SNP as one region, most regions explained <0.50% of genomic variance and covariance for milk yield and fat yield, and explained <0.30% for protein yield, while some regions could present large variance and covariance. Although overall correlations between two populations for the three traits were positive and high, a few regions still showed weakly positive or highly negative genomic correlations for milk yield and fat yield. CONCLUSIONS The new multi-trait Bayesian method using latent variables to model heterogeneous variance and covariance could work well for estimating the genomic variances and covariances for all genome regions simultaneously. Those estimated genomic parameters could be useful to improve the genomic prediction accuracy for Chinese and Nordic Holstein populations using a joint reference data in the future.
Collapse
Affiliation(s)
- Xiujin Li
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, Tjele, Denmark.,Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture of China, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China.,State Key Laboratory of Biocontrol, School of Life Sciences, Guangzhou Higher Education Mega Center, Sun Yat-sen University, North Third Road, Guangzhou, Guangdong, 510006, People's Republic of China
| | - Mogens Sandø Lund
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, Tjele, Denmark
| | - Luc Janss
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, Tjele, Denmark
| | - Chonglong Wang
- Department of Pig Genetics and Breeding, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China
| | - Xiangdong Ding
- Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture of China, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Qin Zhang
- Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture of China, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China.
| | - Guosheng Su
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, Tjele, Denmark.
| |
Collapse
|
31
|
Olsen HG, Knutsen TM, Kohler A, Svendsen M, Gidskehaug L, Grove H, Nome T, Sodeland M, Sundsaasen KK, Kent MP, Martens H, Lien S. Genome-wide association mapping for milk fat composition and fine mapping of a QTL for de novo synthesis of milk fatty acids on bovine chromosome 13. Genet Sel Evol 2017; 49:20. [PMID: 28193175 PMCID: PMC5307787 DOI: 10.1186/s12711-017-0294-5] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2016] [Accepted: 02/03/2017] [Indexed: 12/02/2022] Open
Abstract
Background Bovine milk is widely regarded as a nutritious food source for humans, although the effects of individual fatty acids on human health is a subject of debate. Based on the assumption that genomic selection offers potential to improve milk fat composition, there is strong interest to understand more about the genetic factors that influence the biosynthesis of bovine milk and the molecular mechanisms that regulate milk fat synthesis and secretion. For this reason, the work reported here aimed at identifying genetic variants that affect milk fatty acid composition in Norwegian Red cattle. Milk fatty acid composition was predicted from the nation-wide recording scheme using Fourier transform infrared spectroscopy data and applied to estimate heritabilities for 36 individual and combined fatty acid traits. The recordings were used to generate daughter yield deviations that were first applied in a genome-wide association (GWAS) study with 17,343 markers to identify quantitative trait loci (QTL) affecting fatty acid composition, and next on high-density and sequence-level datasets to fine-map the most significant QTL on BTA13 (BTA for Bos taurus chromosome). Results The initial GWAS revealed 200 significant associations, with the strongest signals on BTA1, 13 and 15. The BTA13 QTL highlighted a strong functional candidate gene for de novo synthesis of short- and medium-chained saturated fatty acids; acyl-CoA synthetase short-chain family member 2. However, subsequent fine-mapping using single nucleotide polymorphisms (SNPs) from a high-density chip and variants detected by resequencing showed that the effect was more likely caused by a second nearby gene; nuclear receptor coactivator 6 (NCOA6). These findings were confirmed with results from haplotype studies. NCOA6 is a nuclear receptor that interacts with transcription factors such as PPARγ, which is a major regulator of bovine milk fat synthesis. Conclusions An initial GWAS revealed a highly significant QTL for de novo-synthesized fatty acids on BTA13 and was followed by fine-mapping of the QTL within NCOA6. The most significant SNPs were either synonymous or situated in introns; more research is needed to uncover the underlying causal DNA variation(s). Electronic supplementary material The online version of this article (doi:10.1186/s12711-017-0294-5) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Hanne Gro Olsen
- Centre for Integrative Genetics (CIGENE), Department of Animal and Aquacultural Sciences, Norwegian University of Life Sciences, PO Box 5003, 1432, Ås, Norway.
| | - Tim Martin Knutsen
- Centre for Integrative Genetics (CIGENE), Department of Animal and Aquacultural Sciences, Norwegian University of Life Sciences, PO Box 5003, 1432, Ås, Norway
| | - Achim Kohler
- Department of Mathematical Sciences and Technology, Norwegian University of Life Sciences, PO Box 5003, 1432, Ås, Norway.,Centre for Biospectroscopy and Data Modeling, Nofima AS, Osloveien 1, 1430, Ås, Norway
| | | | | | - Harald Grove
- Centre for Integrative Genetics (CIGENE), Department of Animal and Aquacultural Sciences, Norwegian University of Life Sciences, PO Box 5003, 1432, Ås, Norway
| | - Torfinn Nome
- Centre for Integrative Genetics (CIGENE), Department of Animal and Aquacultural Sciences, Norwegian University of Life Sciences, PO Box 5003, 1432, Ås, Norway
| | - Marte Sodeland
- Institute of Marine Research, Flødevigen, 4817, His, Norway.,Department of Natural Sciences, Faculty of Engineering and Science, University of Agder, PO Box 422, 4604, Kristiansand, Norway
| | - Kristil Kindem Sundsaasen
- Centre for Integrative Genetics (CIGENE), Department of Animal and Aquacultural Sciences, Norwegian University of Life Sciences, PO Box 5003, 1432, Ås, Norway
| | - Matthew Peter Kent
- Centre for Integrative Genetics (CIGENE), Department of Animal and Aquacultural Sciences, Norwegian University of Life Sciences, PO Box 5003, 1432, Ås, Norway
| | - Harald Martens
- Department of Engineering Cybernetics, Norwegian University of Science and Technology, 7034, Trondheim, Norway
| | - Sigbjørn Lien
- Centre for Integrative Genetics (CIGENE), Department of Animal and Aquacultural Sciences, Norwegian University of Life Sciences, PO Box 5003, 1432, Ås, Norway
| |
Collapse
|
32
|
Galla SJ, Buckley TR, Elshire R, Hale ML, Knapp M, McCallum J, Moraga R, Santure AW, Wilcox P, Steeves TE. Building strong relationships between conservation genetics and primary industry leads to mutually beneficial genomic advances. Mol Ecol 2016; 25:5267-5281. [PMID: 27641156 DOI: 10.1111/mec.13837] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2016] [Revised: 08/23/2016] [Accepted: 08/24/2016] [Indexed: 02/06/2023]
Abstract
Several reviews in the past decade have heralded the benefits of embracing high-throughput sequencing technologies to inform conservation policy and the management of threatened species, but few have offered practical advice on how to expedite the transition from conservation genetics to conservation genomics. Here, we argue that an effective and efficient way to navigate this transition is to capitalize on emerging synergies between conservation genetics and primary industry (e.g., agriculture, fisheries, forestry and horticulture). Here, we demonstrate how building strong relationships between conservation geneticists and primary industry scientists is leading to mutually-beneficial outcomes for both disciplines. Based on our collective experience as collaborative New Zealand-based scientists, we also provide insight for forging these cross-sector relationships.
Collapse
Affiliation(s)
- Stephanie J Galla
- School of Biological Sciences, University of Canterbury, Private Bag 4800, Christchurch, 8140, New Zealand.
| | - Thomas R Buckley
- Landcare Research, Private Bag 92170, Auckland Mail Centre, Auckland, 1142, New Zealand.,School of Biological Sciences, University of Auckland, Auckland, 1010, New Zealand
| | - Rob Elshire
- The Elshire Group, Ltd., 52 Victoria Avenue, Palmerston North, 4410, New Zealand
| | - Marie L Hale
- School of Biological Sciences, University of Canterbury, Private Bag 4800, Christchurch, 8140, New Zealand
| | - Michael Knapp
- Department of Anatomy, University of Otago, P.O. Box 913, Dunedin, 9054, New Zealand
| | - John McCallum
- Breeding and Genomics, New Zealand Institute for Plant and Food Research, Private Bag 4704, Christchurch, 8140, New Zealand
| | - Roger Moraga
- AgResearch, Ruakura Research Centre, Bisley Road, Private Bag 3115, Hamilton, 3240, New Zealand
| | - Anna W Santure
- School of Biological Sciences, University of Auckland, Auckland, 1010, New Zealand
| | - Phillip Wilcox
- Department of Mathematics and Statistics, University of Otago, P.O. Box 56, 710 Cumberland Street, Dunedin, 9054, New Zealand
| | - Tammy E Steeves
- School of Biological Sciences, University of Canterbury, Private Bag 4800, Christchurch, 8140, New Zealand
| |
Collapse
|
33
|
Sempéré G, Philippe F, Dereeper A, Ruiz M, Sarah G, Larmande P. Gigwa-Genotype investigator for genome-wide analyses. Gigascience 2016; 5:25. [PMID: 27267926 PMCID: PMC4897896 DOI: 10.1186/s13742-016-0131-8] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2016] [Accepted: 05/16/2016] [Indexed: 01/16/2023] Open
Abstract
BACKGROUND Exploring the structure of genomes and analyzing their evolution is essential to understanding the ecological adaptation of organisms. However, with the large amounts of data being produced by next-generation sequencing, computational challenges arise in terms of storage, search, sharing, analysis and visualization. This is particularly true with regards to studies of genomic variation, which are currently lacking scalable and user-friendly data exploration solutions. DESCRIPTION Here we present Gigwa, a web-based tool that provides an easy and intuitive way to explore large amounts of genotyping data by filtering it not only on the basis of variant features, including functional annotations, but also on genotype patterns. The data storage relies on MongoDB, which offers good scalability properties. Gigwa can handle multiple databases and may be deployed in either single- or multi-user mode. In addition, it provides a wide range of popular export formats. CONCLUSIONS The Gigwa application is suitable for managing large amounts of genomic variation data. Its user-friendly web interface makes such processing widely accessible. It can either be simply deployed on a workstation or be used to provide a shared data portal for a given community of researchers.
Collapse
Affiliation(s)
- Guilhem Sempéré
- UMR InterTryp (CIRAD), Campus International de Baillarguet, 34398, Montpellier, Cedex 5, France.
- South Green Bioinformatics Platform, 1000 Avenue Agropolis, 34934, Montpellier, Cedex 5, France.
| | - Florian Philippe
- UMR DIADE (IRD), 911 Avenue Agropolis, 34934, Montpellier, Cedex 5, France
| | - Alexis Dereeper
- South Green Bioinformatics Platform, 1000 Avenue Agropolis, 34934, Montpellier, Cedex 5, France
- UMR IPME (IRD), 911 Avenue Agropolis, 34394, Montpellier, Cedex 5, France
| | - Manuel Ruiz
- South Green Bioinformatics Platform, 1000 Avenue Agropolis, 34934, Montpellier, Cedex 5, France
- UMR AGAP, CIRAD, 34398, Montpellier, Cedex 5, France
- Institut de Biologie Computationnelle, Université de Montpellier, 860 Rue de St Priest, 34095, Montpellier, Cedex 5, France
- Agrobiodiversity Research Area, International Center for Tropical Agriculture (CIAT), 6713, Cali, Colombia
| | - Gautier Sarah
- South Green Bioinformatics Platform, 1000 Avenue Agropolis, 34934, Montpellier, Cedex 5, France
- INRA, UMR AGAP, 34398, Montpellier, Cedex 5, France
| | - Pierre Larmande
- South Green Bioinformatics Platform, 1000 Avenue Agropolis, 34934, Montpellier, Cedex 5, France
- UMR DIADE (IRD), 911 Avenue Agropolis, 34934, Montpellier, Cedex 5, France
- Institut de Biologie Computationnelle, Université de Montpellier, 860 Rue de St Priest, 34095, Montpellier, Cedex 5, France
- INRIA Zenith Team, LIRMM, 161 Rue Ada, 34095, Montpellier, Cedex 5, France
| |
Collapse
|