1
|
Nayak SS, Panigrahi M, Rajawat D, Ghildiyal K, Sharma A, Parida S, Bhushan B, Mishra BP, Dutt T. Comprehensive selection signature analyses in dairy cattle exploiting purebred and crossbred genomic data. Mamm Genome 2023; 34:615-631. [PMID: 37843569 DOI: 10.1007/s00335-023-10021-4] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2023] [Accepted: 09/24/2023] [Indexed: 10/17/2023]
Abstract
The main objective of the current research was to locate, annotate, and highlight specific areas of the bovine genome that are undergoing intense positive selection. Here, we are analyzing selection signatures in crossbred (Bos taurus X Bos indicus), taurine (Bos taurus), and indicine (Bos indicus) cattle breeds. Indicine cattle breeds found throughout India are known for their higher heat tolerance and disease resilience. More breeds and more methods can provide a better understanding of the selection signature. So, we have worked on nine distinct cattle breeds utilizing seven different summary statistics, which is a fairly extensive approach. In this study, we carried out a thorough genome-wide investigation of selection signatures using bovine 50K SNP data. We have included the genotyped data of two taurine, two crossbreds, and five indicine cattle breeds, for a total of 320 animals. During the 1950s, these indicine (cebuine) cattle breeds were exported with the aim of enhancing the resilience of taurine breeds in Western countries. For this study, we employed seven summary statistics, including intra-population, i.e., Tajima's D, CLR, iHS, and ROH and inter-population statistics, i.e., FST, XP-EHH, and Rsb. The NCBI database, PANTHER 17.0, and CattleQTL database were used for annotation after finding the important areas under selection. Some genes, including EPHA6, CTNNA2, NPFFR2, HS6ST3, NPR3, KCNIP4, LIPK, SDCBP, CYP7A1, NSMAF, UBXN2B, UGDH, UBE2K, and DAB1, were shown to be shared by three or more different approaches. Therefore, it gives evidence of the most intense selection in these areas. These genes are mostly linked to milk production and adaptability traits. This study also reveals selection regions that contain genes which are crucial to numerous biological functions, including those associated with milk production, coat color, glucose metabolism, oxidative stress response, immunity and circadian rhythms.
Collapse
Affiliation(s)
- Sonali Sonejita Nayak
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, UP, 243122, India
| | - Manjit Panigrahi
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, UP, 243122, India.
| | - Divya Rajawat
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, UP, 243122, India
| | - Kanika Ghildiyal
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, UP, 243122, India
| | - Anurodh Sharma
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, UP, 243122, India
| | - Subhashree Parida
- Division of Pharmacology & Toxicology, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, UP, 243122, India
| | - Bharat Bhushan
- Division of Animal Genetics, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, UP, 243122, India
| | - B P Mishra
- ICAR-National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - Triveni Dutt
- Livestock Production and Management Section, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, UP, 243122, India
| |
Collapse
|
2
|
Kovalchuk SN, Arkhipova AL. Development of TaqMan PCR assay for genotyping SNP rs211250281 of the bovine agpat6 gene. Anim Biotechnol 2023; 34:3250-3255. [PMID: 35635030 DOI: 10.1080/10495398.2022.2077742] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/01/2022]
Abstract
Milk fat percentage is an important production trait of dairy cattle and is one of the goals of breeding programs. Over 95% of the milk fat accounts for triacylglycerols. AGPAT6 (1-acylglycerol-3-phosphate O-acyltransferase 6) catalyzes an intermediary step of triglyceride synthesis in the mammary cells. Genome-wide association studies identified SNP rs211250281 (g27: 36520069 G/T) in the agpat6 gene associated with milk fat content and fat-to-protein ratio in dairy cattle. The article presents data on the development of TaqMan PCR assay for genotyping SNP rs211250281 of the bovine agpat6 gene. In this method, a primer pair, initiating amplification of 75-bp fragments of the agpat6 gene, and two allele-specific TaqMan probes are used. Identification of the G and T alleles is based on a comparison of the final fluorescence intensity of FAM and VIC dyes, respectively. The developed TaqMan PCR assay was validated by Sanger sequencing method. The results of both methods fully coincided, that confirmed high accuracy of the developed TaqMan PCR assay. This reliable, simple, rapid, and high-throughput method could be a suitable tool for studying the distribution of the SNP rs211250281 among different cattle breeds and its association with milk fat content.
Collapse
Affiliation(s)
- Svetlana N Kovalchuk
- Department of Molecular Biotechnology, Institute of Innovative Biotechnologies in Animal Husbandry - the branch of L.K. Ernst Federal Research Center for Animal Husbandry, Moscow, Russia
| | - Anna L Arkhipova
- Department of Molecular Biotechnology, Institute of Innovative Biotechnologies in Animal Husbandry - the branch of L.K. Ernst Federal Research Center for Animal Husbandry, Moscow, Russia
| |
Collapse
|
3
|
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
|
4
|
Naderi S, Moradi MH, Farhadian M, Yin T, Jaeger M, Scheper C, Korkuc P, Brockmann GA, König S, May K. Assessing selection signatures within and between selected lines of dual-purpose black and white and German Holstein cattle. Anim Genet 2020; 51:391-408. [PMID: 32100321 DOI: 10.1111/age.12925] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/02/2020] [Indexed: 12/29/2022]
Abstract
The aim of this study was to detect selection signatures considering cows from the German Holstein (GH) and the local dual-purpose black and white (DSN) population, as well as from generated sub-populations. The 4654 GH and 261 DSN cows were genotyped with the BovineSNP50 Genotyping BeadChip. The geographical herd location was used as an environmental descriptor to create the East-DSN and West-DSN sub-populations. In addition, two further sub-populations of GH cows were generated, using the extreme values for solutions of residual effects of cows for the claw disorder dermatitis digitalis. These groups represented the most susceptible and most resistant cows. We used cross-population extended haplotype homozygosity methodology (XP-EHH) to identify the most recent selection signatures. Furthermore, we calculated Wright's fixation index (FST ). Chromosomal segments for the top 0.1 percentile of negative or positive XP-EHH scores were studied in detail. For gene annotations, we used the Ensembl database and we considered a window of 250 kbp downstream and upstream of each core SNP corresponding to peaks of XP-EHH. In addition, functional interactions among potential candidate genes were inferred via gene network analyses. The most outstanding XP-EHH score was on chromosome 12 (at 77.34 Mb) for DSN and on chromosome 20 (at 36.29-38.42 Mb) for GH. Selection signature locations harbored QTL for several economically important milk and meat quality traits, reflecting the different breeding goals for GH and DSN. The average FST value between GH and DSN was quite low (0.068), indicating shared founders. For group stratifications according to cow health, several identified potential candidate genes influence disease resistance, especially to dermatitis digitalis.
Collapse
Affiliation(s)
- S Naderi
- Institute of Animal Breeding and Genetics, Justus-Liebig University Giessen, Ludwigstr. 21b, Giessen, Germany
| | - M H Moradi
- Department of Animal Sciences, Arak University, Shahid Beheshti Street, Arak, Iran
| | - M Farhadian
- Department of Animal Science, University of Tabriz, 29 Bahman Boulevard, Tabriz, Iran
| | - T Yin
- Institute of Animal Breeding and Genetics, Justus-Liebig University Giessen, Ludwigstr. 21b, Giessen, Germany
| | - M Jaeger
- Institute of Animal Breeding and Genetics, Justus-Liebig University Giessen, Ludwigstr. 21b, Giessen, Germany
| | - C Scheper
- Institute of Animal Breeding and Genetics, Justus-Liebig University Giessen, Ludwigstr. 21b, Giessen, Germany
| | - P Korkuc
- Albrecht Daniel Thaer Institute for Agricultural and Horticultural Sciences, Humboldt University Berlin, Invalidenstr. 42, Berlin, D-10115, Germany
| | - G A Brockmann
- Albrecht Daniel Thaer Institute for Agricultural and Horticultural Sciences, Humboldt University Berlin, Invalidenstr. 42, Berlin, D-10115, Germany
| | - S König
- Institute of Animal Breeding and Genetics, Justus-Liebig University Giessen, Ludwigstr. 21b, Giessen, Germany
| | - K May
- Institute of Animal Breeding and Genetics, Justus-Liebig University Giessen, Ludwigstr. 21b, Giessen, Germany
| |
Collapse
|
5
|
Ghoreishifar SM, Moradi-Shahrbabak H, Fallahi MH, Jalil Sarghale A, Moradi-Shahrbabak M, Abdollahi-Arpanahi R, Khansefid M. Genomic measures of inbreeding coefficients and genome-wide scan for runs of homozygosity islands in Iranian river buffalo, Bubalus bubalis. BMC Genet 2020; 21:16. [PMID: 32041535 PMCID: PMC7011551 DOI: 10.1186/s12863-020-0824-y] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2019] [Accepted: 02/04/2020] [Indexed: 01/06/2023] Open
Abstract
Background Consecutive homozygous fragments of a genome inherited by offspring from a common ancestor are known as runs of homozygosity (ROH). ROH can be used to calculate genomic inbreeding and to identify genomic regions that are potentially under historical selection pressure. The dataset of our study consisted of 254 Azeri (AZ) and 115 Khuzestani (KHZ) river buffalo genotyped for ~ 65,000 SNPs for the following two purposes: 1) to estimate and compare inbreeding calculated using ROH (FROH), excess of homozygosity (FHOM), correlation between uniting gametes (FUNI), and diagonal elements of the genomic relationship matrix (FGRM); 2) to identify frequently occurring ROH (i.e. ROH islands) for our selection signature and gene enrichment studies. Results In this study, 9102 ROH were identified, with an average number of 21.2 ± 13.1 and 33.2 ± 15.9 segments per animal in AZ and KHZ breeds, respectively. On average in AZ, 4.35% (108.8 ± 120.3 Mb), and in KHZ, 5.96% (149.1 ± 107.7 Mb) of the genome was autozygous. The estimated inbreeding values based on FHOM, FUNI and FGRM were higher in AZ than they were in KHZ, which was in contrast to the FROH estimates. We identified 11 ROH islands (four in AZ and seven in KHZ). In the KHZ breed, the genes located in ROH islands were enriched for multiple Gene Ontology (GO) terms (P ≤ 0.05). The genes located in ROH islands were associated with diverse biological functions and traits such as body size and muscle development (BMP2), immune response (CYP27B1), milk production and components (MARS, ADRA1A, and KCTD16), coat colour and pigmentation (PMEL and MYO1A), reproductive traits (INHBC, INHBE, STAT6 and PCNA), and bone development (SUOX). Conclusion The calculated FROH was in line with expected higher inbreeding in KHZ than in AZ because of the smaller effective population size of KHZ. Thus, we find that FROH can be used as a robust estimate of genomic inbreeding. Further, the majority of ROH peaks were overlapped with or in close proximity to the previously reported genomic regions with signatures of selection. This tells us that it is likely that the genes in the ROH islands have been subject to artificial or natural selection.
Collapse
Affiliation(s)
- Seyed Mohammad Ghoreishifar
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Hossein Moradi-Shahrbabak
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran.
| | - Mohammad Hossein Fallahi
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Ali Jalil Sarghale
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Mohammad Moradi-Shahrbabak
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj, 31587-11167, Iran
| | - Rostam Abdollahi-Arpanahi
- Departments of Animal and Poultry Science, College of Aburaihan, University of Tehran, Pakdasht, 33916-53755, Iran
| | - Majid Khansefid
- AgriBio Centre for AgriBioscience, Agriculture Victoria, Bundoora, VIC, 3083, Australia
| |
Collapse
|
6
|
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
|
7
|
Cai Z, Guldbrandtsen B, Lund MS, Sahana G. Retraction Note: Dissecting closely linked association signals in combination with the mammalian phenotype database can identify candidate genes in dairy cattle. BMC Genet 2018; 19:111. [PMID: 30537928 PMCID: PMC6288910 DOI: 10.1186/s12863-018-0698-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
Affiliation(s)
- Zexi Cai
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830, Tjele, Denmark.
| | - Bernt Guldbrandtsen
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830, Tjele, Denmark
| | - Mogens Sandø Lund
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830, Tjele, Denmark
| | - Goutam Sahana
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830, Tjele, Denmark
| |
Collapse
|
8
|
Cai Z, Guldbrandtsen B, Lund MS, Sahana G. Prioritizing candidate genes post-GWAS using multiple sources of data for mastitis resistance in dairy cattle. BMC Genomics 2018; 19:656. [PMID: 30189836 PMCID: PMC6127918 DOI: 10.1186/s12864-018-5050-x] [Citation(s) in RCA: 45] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2018] [Accepted: 08/31/2018] [Indexed: 12/31/2022] Open
Abstract
Background Improving resistance to mastitis, one of the costliest diseases in dairy production, has become an important objective in dairy cattle breeding. However, mastitis resistance is influenced by many genes involved in multiple processes, including the response to infection, inflammation, and post-infection healing. Low genetic heritability, environmental variations, and farm management differences further complicate the identification of links between genetic variants and mastitis resistance. Consequently, studies of the genetics of variation in mastitis resistance in dairy cattle lack agreement about the responsible genes. Results We associated 15,552,968 imputed whole-genome sequencing markers for 5147 Nordic Holstein cattle with mastitis resistance in a genome-wide association study (GWAS). Next, we augmented P-values for markers in genes in the associated regions using Gene Ontology terms, Kyoto Encyclopedia of Genes and Genomes pathway analysis, and mammalian phenotype database. To confirm results of gene-based analyses, we used gene expression data from E. coli-challenged cow udders. We identified 22 independent quantitative trait loci (QTL) that collectively explained 14% of the variance in breeding values for resistance to clinical mastitis (CM). Using association test statistics with multiple pieces of independent information on gene function and differential expression during bacterial infection, we suggested putative causal genes with biological relevance for 12 QTL affecting resistance to CM in dairy cattle. Conclusion Combining information on the nearest positional genes, gene-based analyses, and differential gene expression data from RNA-seq, we identified putative causal genes (candidate genes with biological evidence) in QTL for mastitis resistance in Nordic Holstein cattle. The same strategy can be applied for other traits. Electronic supplementary material The online version of this article (10.1186/s12864-018-5050-x) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Zexi Cai
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830, Tjele, Denmark.
| | - Bernt Guldbrandtsen
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830, Tjele, Denmark
| | - Mogens Sandø Lund
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830, Tjele, Denmark
| | - Goutam Sahana
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830, Tjele, Denmark
| |
Collapse
|