1
|
Zhang Z, Zhao W, Wang Z, Pan Y, Wang Q, Zhang Z. Integration of ssGWAS and ROH analyses for uncovering genetic variants associated with reproduction traits in Large White pigs. Anim Genet 2024. [PMID: 39129705 DOI: 10.1111/age.13465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2024] [Revised: 05/26/2024] [Accepted: 07/05/2024] [Indexed: 08/13/2024]
Abstract
The low heritability of reproduction traits such as total number born (TNB), number born alive (NBA) and adjusted litter weight until 21 days at weaning (ALW) poses a challenge for genetic improvement. In this study, we aimed to identify genetic variants that influence these traits and evaluate the accuracy of genomic selection (GS) using these variants as genomic features. We performed single-step genome-wide association studies (ssGWAS) on 17 823 Large White (LW) pigs, of which 2770 were genotyped by 50K single nucleotide polymorphism (SNP) chips. Additionally, we analyzed runs of homozygosity (ROH) in the population and tested their effects on the traits. The genomic feature best linear unbiased prediction (GFBLUP) was then carried out in an independent population of 350 LW pigs using identified trait-related SNP subsets as genomic features. As a result, our findings identified five, one and four SNP windows that explaining more than 1% of genetic variance for ALW, TNB, and NBA, respectively and discovered 358 hotspots and nine ROH islands. The ROH SSC1:21814570-27186456 and SSC11:7220366-14276394 were found to be significantly associated with ALW and NBA, respectively. We assessed the genomic estimated breeding value accuracy through 20 replicates of five-fold cross-validation. Our findings demonstrate that GFBLUP, incorporating SNPs located in effective ROH (p-value < 0.05) as genomic features, might enhance GS accuracy for ALW compared with GBLUP. Additionally, using SNPs explaining more than 0.1% of the genetic variance in ssGWAS for NBA as genomic features might improve the GS accuracy, too. However, it is important to note that the incorporation of inappropriate genomic features can significantly reduce GS accuracy. In conclusion, our findings provide valuable insights into the genetic mechanisms of reproductive traits in pigs and suggest that the ssGWAS and ROH have the potential to enhance the accuracy of GS for reproductive traits in LW pigs.
Collapse
Affiliation(s)
- Zhenyang Zhang
- Department of Animal Science, College of Animal Science, Zhejiang University, Hangzhou, China
| | - Wei Zhao
- SciGene Biotechnology Co. Ltd, Hefei, China
| | - Zhen Wang
- Department of Animal Science, College of Animal Science, Zhejiang University, Hangzhou, China
| | - Yuchun Pan
- Department of Animal Science, College of Animal Science, Zhejiang University, Hangzhou, China
- Hainan Institute, Zhejiang University, Sanya, China
| | - Qishan Wang
- Department of Animal Science, College of Animal Science, Zhejiang University, Hangzhou, China
- Hainan Institute, Zhejiang University, Sanya, China
| | - Zhe Zhang
- Department of Animal Science, College of Animal Science, Zhejiang University, Hangzhou, China
| |
Collapse
|
2
|
Schneider H, Krizanac AM, Falker-Gieske C, Heise J, Tetens J, Thaller G, Bennewitz J. Genomic dissection of the correlation between milk yield and various health traits using functional and evolutionary information about imputed sequence variants of 34,497 German Holstein cows. BMC Genomics 2024; 25:265. [PMID: 38461236 DOI: 10.1186/s12864-024-10115-6] [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: 08/17/2023] [Accepted: 02/13/2024] [Indexed: 03/11/2024] Open
Abstract
BACKGROUND Over the last decades, it was subject of many studies to investigate the genomic connection of milk production and health traits in dairy cattle. Thereby, incorporating functional information in genomic analyses has been shown to improve the understanding of biological and molecular mechanisms shaping complex traits and the accuracies of genomic prediction, especially in small populations and across-breed settings. Still, little is known about the contribution of different functional and evolutionary genome partitioning subsets to milk production and dairy health. Thus, we performed a uni- and a bivariate analysis of milk yield (MY) and eight health traits using a set of ~34,497 German Holstein cows with 50K chip genotypes and ~17 million imputed sequence variants divided into 27 subsets depending on their functional and evolutionary annotation. In the bivariate analysis, eight trait-combinations were observed that contrasted MY with each health trait. Two genomic relationship matrices (GRM) were included, one consisting of the 50K chip variants and one consisting of each set of subset variants, to obtain subset heritabilities and genetic correlations. In addition, 50K chip heritabilities and genetic correlations were estimated applying merely the 50K GRM. RESULTS In general, 50K chip heritabilities were larger than the subset heritabilities. The largest heritabilities were found for MY, which was 0.4358 for the 50K and 0.2757 for the subset heritabilities. Whereas all 50K genetic correlations were negative, subset genetic correlations were both, positive and negative (ranging from -0.9324 between MY and mastitis to 0.6662 between MY and digital dermatitis). The subsets containing variants which were annotated as noncoding related, splice sites, untranslated regions, metabolic quantitative trait loci, and young variants ranked highest in terms of their contribution to the traits` genetic variance. We were able to show that linkage disequilibrium between subset variants and adjacent variants did not cause these subsets` high effect. CONCLUSION Our results confirm the connection of milk production and health traits in dairy cattle via the animals` metabolic state. In addition, they highlight the potential of including functional information in genomic analyses, which helps to dissect the extent and direction of the observed traits` connection in more detail.
Collapse
Affiliation(s)
- Helen Schneider
- Institute of Animal Science, University of Hohenheim, 70599, Stuttgart, Germany.
| | - Ana-Marija Krizanac
- Department of Animal Sciences, University of Göttingen, 37077, Göttingen, Germany
| | | | - Johannes Heise
- Vereinigte Informationssysteme Tierhaltung w.V. (VIT), 27283, Verden, Germany
| | - Jens Tetens
- Department of Animal Sciences, University of Göttingen, 37077, Göttingen, Germany
| | - Georg Thaller
- Institute of Animal Breeding and Husbandry, Christian-Albrechts University of Kiel, 24098, Kiel, Germany
| | - Jörn Bennewitz
- Institute of Animal Science, University of Hohenheim, 70599, Stuttgart, Germany
| |
Collapse
|
3
|
Zhao X, Luo H, Lu H, Ma L, Li Y, Dou J, Zhang J, Ma Y, Li J, Wang Y. RNA-Seq Analysis of Peripheral Whole Blood from Dairy Bulls with High and Low Antibody-Mediated Immune Responses-A Preliminary Study. Animals (Basel) 2023; 13:2208. [PMID: 37444006 DOI: 10.3390/ani13132208] [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/16/2023] [Revised: 06/28/2023] [Accepted: 07/03/2023] [Indexed: 07/15/2023] Open
Abstract
Enhancing the immune response through breeding is regarded as an effective strategy for improving animal health, as dairy cattle identified as high immune responders are reported to have a decreased prevalence of economically significant diseases. The identification of differentially expressed genes (DEGs) associated with immune responses might be an effective tool for breeding healthy dairy cattle. In this study, antibody-mediated immune responses (AMIRs) were induced by the immunization of hen egg white lysozyme (HEWL) in six Chinese Holstein dairy bulls divided into high- and low-AMIR groups based on their HEWL antibody level. Then, RNA-seq was applied to explore the transcriptome of peripheral whole blood between the two comparison groups. As a result, several major upregulated and downregulated genes were identified and attributed to the regulation of locomotion, tissue development, immune response, and detoxification. In addition, the result of the KEGG pathway analysis revealed that most DEGs were enriched in pathways related to disease, inflammation, and immune response, including antigen processing and presentation, Staphylococcus aureus infection, intestinal immune network for IgA production, cytokine-cytokine receptor interaction, and complement and coagulation cascades. Moreover, six genes (BOLA-DQA5, C5, CXCL2, HBA, LTF, and COL1A1) were validated using RT-qPCR, which may provide information for genomic selection in breeding programs. These results broaden the knowledge of the immune response mechanism in dairy bulls, which has strong implications for breeding cattle with an enhanced AMIR.
Collapse
Affiliation(s)
- Xiuxin Zhao
- Ningxia Key Laboratory of Ruminant Molecular and Cellular Breeding, College of Animal Science and Technology, Ningxia University, Yinchuan 750021, China
- Institute of Animal Science and Veterinary Medicine, Shandong Academy of Agricultural Sciences, No. 23788, Gongyebei Road, Jinan 250100, China
- Shandong Ox Livestock Breeding Co., Ltd., Jinan 250100, China
| | - Hanpeng Luo
- 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
| | - Haibo Lu
- 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
- Beijing Consortium for Innovative Bio-Breeding, Beijing 101206, China
| | - Longgang Ma
- 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
| | - Yanqin Li
- Institute of Animal Science and Veterinary Medicine, Shandong Academy of Agricultural Sciences, No. 23788, Gongyebei Road, Jinan 250100, China
| | - Jinhuan Dou
- 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
| | - Junxing Zhang
- Ningxia Key Laboratory of Ruminant Molecular and Cellular Breeding, College of Animal Science and Technology, Ningxia University, Yinchuan 750021, China
| | - Yun Ma
- Ningxia Key Laboratory of Ruminant Molecular and Cellular Breeding, College of Animal Science and Technology, Ningxia University, Yinchuan 750021, China
| | - Jianbin Li
- Institute of Animal Science and Veterinary Medicine, Shandong Academy of Agricultural Sciences, No. 23788, Gongyebei Road, Jinan 250100, China
| | - Yachun Wang
- Ningxia Key Laboratory of Ruminant Molecular and Cellular Breeding, College of Animal Science and Technology, Ningxia University, Yinchuan 750021, China
- 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
| |
Collapse
|
4
|
Ye H, Xu Z, Bello SF, Zhu Q, Kong S, Zheng M, Fang X, Jia X, Xu H, Zhang X, Nie Q. Haplotype analysis of genomic prediction by incorporating genomic pathway information based on high-density SNP marker in Chinese yellow-feathered chicken. Poult Sci 2023; 102:102549. [PMID: 36907129 PMCID: PMC10024239 DOI: 10.1016/j.psj.2023.102549] [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: 11/12/2022] [Revised: 01/16/2023] [Accepted: 01/27/2023] [Indexed: 02/09/2023] Open
Abstract
Genomic selection using single nucleotide polymorphism (SNP) markers is now intensively investigated in breeding and has been widely utilized for genetic improvement. Currently, several studies have used haplotype (consisting of multiallelic SNPs) for genomic prediction and revealed its performance advantage. In this study, we comprehensively evaluated the performance of haplotype models for genomic prediction in 15 traits, including 6 growth, 5 carcass, and 4 feeding traits in a Chinese yellow-feathered chicken population. We adopted 3 methods to define haplotypes from high-density SNP panels, and our strategy included combining Kyoto Encyclopedia of Genes and Genomes pathway information and considering linkage disequilibrium (LD) information. Our results showed an increase in prediction accuracy due to haplotypes ranging from -0.04∼27.16% in all traits, where the significant improvements were found in 12 traits. The estimates of haplotype epistasis heritability were strongly correlated with the accuracy increase by haplotype models. In addition, incorporating genomic annotation information could further increase the accuracy of the haplotype model, where the further increase in accuracy is significantly relative to the increase of relative haplotype epistasis heritability. The genomic prediction using LD information for constructing haplotypes has the best prediction performance among the 4 traits. These results uncovered that haplotype methods were beneficial for genomic prediction, and the accuracy could be further increased by incorporating genomic annotation information. Moreover, using LD information would potentially improve the performance of genomic prediction.
Collapse
Affiliation(s)
- Haoqiang Ye
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China
| | - Zhenqiang Xu
- Wen's Nanfang Poultry Breeding Co. Ltd, Guangdong Province, Yunfu 527400, China
| | - Semiu Folaniyi Bello
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China
| | - Qianghui Zhu
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China
| | - Shaofen Kong
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China
| | - Ming Zheng
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China
| | - Xiang Fang
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China
| | - Xinzheng Jia
- Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding, Foshan University, Foshan, 528225 China
| | - Haiping Xu
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China
| | - Xiquan Zhang
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China
| | - Qinghua Nie
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou, 510642 China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, Guangzhou, 510642 China.
| |
Collapse
|
5
|
Li W, Mi S, Zhang J, Liu X, Chen S, Liu S, Feng X, Tang Y, Li Y, Liu L, Fang L, Zhang S, Yu Y. Integrating sperm cell transcriptome and seminal plasma metabolome to analyze the molecular regulatory mechanism of sperm motility in Holstein stud bulls. J Anim Sci 2023; 101:skad214. [PMID: 37366074 PMCID: PMC10355371 DOI: 10.1093/jas/skad214] [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: 09/12/2022] [Accepted: 06/26/2023] [Indexed: 06/28/2023] Open
Abstract
Considering that artificial insemination is the most widely used assisted reproductive technique in the dairy industry, the semen quality of bulls is very important for selecting excellent stud bulls. Sperm motility is one of the important traits of semen quality, and related genes may be regulated by environmental factors. Seminal plasma can affect sperm cell transcriptome and further affect sperm motility through exosome or other processes. However, the molecular regulation mechanism of bull sperm motility has not been studied by combining the sperm cell transcriptome with seminal plasma metabolome. The number of motile sperm per ejaculate (NMSPE) is an integrated indicator for assessing sperm motility in stud bulls. In the present study, we selected 7 bulls with higher NMSPE (5,698.55 million +/- 945.40 million) as group H and 7 bulls with lower NMSPE (2,279.76 million +/- 1,305.69 million) as group L from 53 Holstein stud bulls. The differentially expressed genes (DEGs) in sperm cells were evaluated between the two groups (H vs. L). We conducted gene co-expression network analysis (WGCNA) on H and L groups of bulls, as well as two monozygotic twin Holstein bulls with different NMSPE values, to screen candidate genes for NMSPE. The regulatory effect of seminal plasma metabolome on the candidate genes of NMSPE was also investigated. A total of 1,099 DEGs were identified in the sperm cells of H and L groups. These DEGs were primarily concentrated in energy metabolism and sperm cell transcription. The significantly enriched Kyoto encyclopedia of genes and genomes (KEGG) pathways of the 57 differential metabolites were the aminoacyl-tRNA biosynthesis pathway and vitamin B6 metabolism pathway. Our study discovered 14 genes as the potential candidate markers for sperm motility, including FBXO39. We observed a broad correlation between transcriptome of sperm cells and seminal plasma metabolome, such as three metabolites, namely, mesaconic acid, 2-coumaric acid, and 4-formylaminoantipyrine, might regulate FBXO39 expression through potential pathways. The genes related to seminal plasma metabolites expressed in sperm cells are not only located near the quantitative trait loci of reproductive traits, but also enriched in the genome-wide association study signal of sire conception rate. Collectively, this study was the first to investigate the interplays among transcriptome of sperm cells and seminal plasma metabolome from Holstein stud bulls with different sperm motility.
Collapse
Affiliation(s)
- Wenlong Li
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, Department of Animal Breeding and Genetics, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Siyuan Mi
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, Department of Animal Breeding and Genetics, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Jinning Zhang
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, Department of Animal Breeding and Genetics, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Xueqin Liu
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, Department of Animal Breeding and Genetics, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Siqian Chen
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, Department of Animal Breeding and Genetics, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Shuli Liu
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, Department of Animal Breeding and Genetics, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
- School of Life Sciences, Westlake University, Hangzhou, Zhejiang 310024, China
| | - Xia Feng
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, Department of Animal Breeding and Genetics, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Yongjie Tang
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, Department of Animal Breeding and Genetics, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Yanhua Li
- Beijing Dairy Cattle Center, Qinghe’nanzhen Deshengmenwai Road, Beijing 100192, China
| | - Lin Liu
- Beijing Dairy Cattle Center, Qinghe’nanzhen Deshengmenwai Road, Beijing 100192, China
| | - Lingzhao Fang
- MRC Human Genetics Unit, Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh EH4 2XU, UK
- Center for Quantitative Genetics and Genomics (QGG), Aarhus University, Aarhus, Denmark
| | - Shengli Zhang
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, Department of Animal Breeding and Genetics, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Ying Yu
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, Department of Animal Breeding and Genetics, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| |
Collapse
|
6
|
Cai X, Teng J, Ren D, Zhang H, Li J, Zhang Z. Model Comparison of Heritability Enrichment Analysis in Livestock Population. Genes (Basel) 2022; 13:genes13091644. [PMID: 36140810 PMCID: PMC9498849 DOI: 10.3390/genes13091644] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2022] [Revised: 09/03/2022] [Accepted: 09/10/2022] [Indexed: 11/16/2022] Open
Abstract
Heritability enrichment analysis is an important means of exploring the genetic architecture of complex traits in human genetics. Heritability enrichment is typically defined as the proportion of an SNP subset explained heritability, divided by the proportion of SNPs. Heritability enrichment enables better study of underlying complex traits, such as functional variant/gene subsets, biological networks and metabolic pathways detected through integrating explosively increased omics data. This would be beneficial for genomic prediction of disease risk in humans and genetic values estimation of important economical traits in livestock and plant species. However, in livestock, factors affecting the heritability enrichment estimation of complex traits have not been examined. Previous studies on humans reported that the frequencies, effect sizes, and levels of linkage disequilibrium (LD) of underlying causal variants (CVs) would affect the heritability enrichment estimation. Therefore, the distribution of heritability across the genome should be fully considered to obtain the unbiased estimation of heritability enrichment. To explore the performance of different heritability enrichment models in livestock populations, we used the VanRaden, GCTA and α models, assuming different α values, and the LDAK model, considering LD weight. We simulated three types of phenotypes, with CVs from various minor allele frequency (MAF) ranges: genome-wide (0.005 ≤ MAF ≤ 0.5), common (0.05 ≤ MAF ≤ 0.5), and uncommon (0.01 ≤ MAF < 0.05). The performances of the models with two different subsets (one of which contained known CVs and the other consisting of randomly selected markers) were compared to verify the accuracy of heritability enrichment estimation of functional variant sets. Our results showed that models with known CV subsets provided more robust enrichment estimation. Models with different α values tended to provide stable and accurate estimates for common and genome-wide CVs (relative deviation 0.5−2.2%), while tending to underestimate the enrichment of uncommon CVs. As the α value increased, enrichments from 15.73% higher than true value (i.e., 3.00) to 48.93% lower than true value for uncommon CVs were observed. In addition, the long-range LD windows (e.g., 5000 kb) led to large bias of the enrichment estimations for both common and uncommon CVs. Overall, heritability enrichment estimations were sensitive for the α value assumption and LD weight consideration of different models. Accuracy would be greatly improved by using a suitable model. This study would be helpful in understanding the genetic architecture of complex traits and provides a reference for genetic analysis in the livestock population.
Collapse
|
7
|
Yao Y, Liu S, Xia C, Gao Y, Pan Z, Canela-Xandri O, Khamseh A, Rawlik K, Wang S, Li B, Zhang Y, Pairo-Castineira E, D’Mellow K, Li X, Yan Z, Li CJ, Yu Y, Zhang S, Ma L, Cole JB, Ross PJ, Zhou H, Haley C, Liu GE, Fang L, Tenesa A. Comparative transcriptome in large-scale human and cattle populations. Genome Biol 2022; 23:176. [PMID: 35996157 PMCID: PMC9394047 DOI: 10.1186/s13059-022-02745-4] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2020] [Accepted: 08/09/2022] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Cross-species comparison of transcriptomes is important for elucidating evolutionary molecular mechanisms underpinning phenotypic variation between and within species, yet to date it has been essentially limited to model organisms with relatively small sample sizes. RESULTS Here, we systematically analyze and compare 10,830 and 4866 publicly available RNA-seq samples in humans and cattle, respectively, representing 20 common tissues. Focusing on 17,315 orthologous genes, we demonstrate that mean/median gene expression, inter-individual variation of expression, expression quantitative trait loci, and gene co-expression networks are generally conserved between humans and cattle. By examining large-scale genome-wide association studies for 46 human traits (average n = 327,973) and 45 cattle traits (average n = 24,635), we reveal that the heritability of complex traits in both species is significantly more enriched in transcriptionally conserved than diverged genes across tissues. CONCLUSIONS In summary, our study provides a comprehensive comparison of transcriptomes between humans and cattle, which might help decipher the genetic and evolutionary basis of complex traits in both species.
Collapse
Affiliation(s)
- Yuelin Yao
- MRC Human Genetics Unit at the Institute of Genetics and Cancer, The University of Edinburgh, EH4 2XU Edinburgh, UK
- School of Informatics, The University of Edinburgh, Edinburgh, EH8 9AB UK
| | - Shuli Liu
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705 USA
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Charley Xia
- The Roslin Institute, Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Midlothian, EH25 9RG UK
- Department of Psychology, 7 George Square, The University of Edinburgh, Edinburgh, EH8 9JZ UK
| | - Yahui Gao
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705 USA
- Department of Animal and Avian Sciences, University of Maryland, College Park, MA 20742 USA
| | - Zhangyuan Pan
- Department of Animal Science, University of California, Davis, CA 95616 USA
- Present address: Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Oriol Canela-Xandri
- MRC Human Genetics Unit at the Institute of Genetics and Cancer, The University of Edinburgh, EH4 2XU Edinburgh, UK
| | - Ava Khamseh
- MRC Human Genetics Unit at the Institute of Genetics and Cancer, The University of Edinburgh, EH4 2XU Edinburgh, UK
- School of Informatics, The University of Edinburgh, Edinburgh, EH8 9AB UK
| | - Konrad Rawlik
- The Roslin Institute, Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Midlothian, EH25 9RG UK
| | - Sheng Wang
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, 650223 Yunnan China
| | - Bingjie Li
- Scotland’s Rural College (SRUC), Roslin Institute Building, Midlothian, EH25 9RG UK
| | - Yi Zhang
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Erola Pairo-Castineira
- MRC Human Genetics Unit at the Institute of Genetics and Cancer, The University of Edinburgh, EH4 2XU Edinburgh, UK
- The Roslin Institute, Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Midlothian, EH25 9RG UK
| | - Kenton D’Mellow
- MRC Human Genetics Unit at the Institute of Genetics and Cancer, The University of Edinburgh, EH4 2XU Edinburgh, UK
| | - Xiujin Li
- Guangdong Provincial Key Laboratory of Waterfowl Healthy Breeding, College of Animal Science & Technology, Zhongkai University of Agriculture and Engineering, Guangzhou, 510225 Guangdong China
| | - Ze Yan
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Cong-jun Li
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705 USA
| | - Ying Yu
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Shengli Zhang
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Li Ma
- Department of Animal and Avian Sciences, University of Maryland, College Park, MA 20742 USA
| | - John B. Cole
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705 USA
| | - Pablo J. Ross
- Department of Animal Science, University of California, Davis, CA 95616 USA
| | - Huaijun Zhou
- Department of Animal Science, University of California, Davis, CA 95616 USA
| | - Chris Haley
- MRC Human Genetics Unit at the Institute of Genetics and Cancer, The University of Edinburgh, EH4 2XU Edinburgh, UK
- The Roslin Institute, Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Midlothian, EH25 9RG UK
| | - George E. Liu
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705 USA
| | - Lingzhao Fang
- MRC Human Genetics Unit at the Institute of Genetics and Cancer, The University of Edinburgh, EH4 2XU Edinburgh, UK
- Present address: Center for Quantitative Genetics and Genomics, Aarhus University, Aarhus, Denmark
| | - Albert Tenesa
- MRC Human Genetics Unit at the Institute of Genetics and Cancer, The University of Edinburgh, EH4 2XU Edinburgh, UK
- The Roslin Institute, Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Midlothian, EH25 9RG UK
| |
Collapse
|
8
|
Hao X, Liang A, Plastow G, Zhang C, Wang Z, Liu J, Salzano A, Gasparrini B, Campanile G, Zhang S, Yang L. An Integrative Genomic Prediction Approach for Predicting Buffalo Milk Traits by Incorporating Related Cattle QTLs. Genes (Basel) 2022; 13:genes13081430. [PMID: 36011341 PMCID: PMC9408041 DOI: 10.3390/genes13081430] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2022] [Revised: 08/08/2022] [Accepted: 08/09/2022] [Indexed: 11/16/2022] Open
Abstract
Background: The 90K Axiom Buffalo SNP Array is expected to improve and speed up various genomic analyses for the buffalo (Bubalus bubalis). Genomic prediction is an effective approach in animal breeding to improve selection and reduce costs. As buffalo genome research is lagging behind that of the cow and production records are also limited, genomic prediction performance will be relatively poor. To improve the genomic prediction in buffalo, we introduced a new approach (pGBLUP) for genomic prediction of six buffalo milk traits by incorporating QTL information from the cattle milk traits in order to help improve the prediction performance for buffalo. Results: In simulations, the pGBLUP could outperform BayesR and the GBLUP if the prior biological information (i.e., the known causal loci) was appropriate; otherwise, it performed slightly worse than BayesR and equal to or better than the GBLUP. In real data, the heritability of the buffalo genomic region corresponding to the cattle milk trait QTLs was enriched (fold of enrichment > 1) in four buffalo milk traits (FY270, MY270, PY270, and PM) when the EBV was used as the response variable. The DEBV as the response variable yielded more reliable genomic predictions than the traditional EBV, as has been shown by previous research. The performance of the three approaches (GBLUP, BayesR, and pGBLUP) did not vary greatly in this study, probably due to the limited sample size, incomplete prior biological information, and less artificial selection in buffalo. Conclusions: To our knowledge, this study is the first to apply genomic prediction to buffalo by incorporating prior biological information. The genomic prediction of buffalo traits can be further improved with a larger sample size, higher-density SNP chips, and more precise prior biological information.
Collapse
Affiliation(s)
- Xingjie Hao
- Department of Epidemiology and Biostatistics, School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan 430030, China
- Correspondence: (X.H.); (L.Y.)
| | - Aixin Liang
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan 430070, China
| | - Graham Plastow
- Livestock Gentec Center, Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB T6G 2C8, Canada
| | - Chunyan Zhang
- Livestock Gentec Center, Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB T6G 2C8, Canada
| | - Zhiquan Wang
- Livestock Gentec Center, Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB T6G 2C8, Canada
| | - Jiajia Liu
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan 430070, China
| | - Angela Salzano
- Department of Veterinary Medicine and Animal Productions, University of Naples “Federico II”, 80137 Naples, Italy
| | - Bianca Gasparrini
- Department of Veterinary Medicine and Animal Productions, University of Naples “Federico II”, 80137 Naples, Italy
| | - Giuseppe Campanile
- Department of Veterinary Medicine and Animal Productions, University of Naples “Federico II”, 80137 Naples, Italy
| | - Shujun Zhang
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan 430070, China
| | - Liguo Yang
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan 430070, China
- Correspondence: (X.H.); (L.Y.)
| |
Collapse
|
9
|
Gao Y, Li J, Cai G, Wang Y, Yang W, Li Y, Zhao X, Li R, Gao Y, Tuo W, Baldwin RL, Li CJ, Fang L, Liu GE. Single-cell transcriptomic and chromatin accessibility analyses of dairy cattle peripheral blood mononuclear cells and their responses to lipopolysaccharide. BMC Genomics 2022; 23:338. [PMID: 35501711 PMCID: PMC9063233 DOI: 10.1186/s12864-022-08562-0] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2021] [Accepted: 04/13/2022] [Indexed: 11/24/2022] Open
Abstract
Background Gram-negative bacteria are important pathogens in cattle, causing severe infectious diseases, including mastitis. Lipopolysaccharides (LPS) are components of the outer membrane of Gram-negative bacteria and crucial mediators of chronic inflammation in cattle. LPS modulations of bovine immune responses have been studied before. However, the single-cell transcriptomic and chromatin accessibility analyses of bovine peripheral blood mononuclear cells (PBMCs) and their responses to LPS stimulation were never reported. Results We performed single-cell RNA sequencing (scRNA-seq) and single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) in bovine PBMCs before and after LPS treatment and demonstrated that seven major cell types, which included CD4 T cells, CD8 T cells, and B cells, monocytes, natural killer cells, innate lymphoid cells, and dendritic cells. Bioinformatic analyses indicated that LPS could increase PBMC cell cycle progression, cellular differentiation, and chromatin accessibility. Gene analyses further showed significant changes in differential expression, transcription factor binding site, gene ontology, and regulatory interactions during the PBMC responses to LPS. Consistent with the findings of previous studies, LPS induced activation of monocytes and dendritic cells, likely through their upregulated TLR4 receptor. NF-κB was observed to be activated by LPS and an increased transcription of an array of pro-inflammatory cytokines, in agreement that NF-κB is an LPS-responsive regulator of innate immune responses. In addition, by integrating LPS-induced differentially expressed genes (DEGs) with large-scale GWAS of 45 complex traits in Holstein, we detected trait-relevant cell types. We found that selected DEGs were significantly associated with immune-relevant health, milk production, and body conformation traits. Conclusion This study provided the first scRNAseq and scATAC-seq data for cattle PBMCs and their responses to the LPS stimulation to the best of our knowledge. These results should also serve as valuable resources for the future study of the bovine immune system and open the door for discoveries about immune cell roles in complex traits like mastitis at single-cell resolution. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08562-0.
Collapse
Affiliation(s)
- Yahui Gao
- Institute of Animal Science and Veterinary Medicine, Shandong Academy of Agricultural Sciences, No.202, Gongyebei Road, Jinan, 250100, China.,Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, MD, 20705, USA
| | - Jianbin Li
- Institute of Animal Science and Veterinary Medicine, Shandong Academy of Agricultural Sciences, No.202, Gongyebei Road, Jinan, 250100, China.
| | - Gaozhan Cai
- Institute of Animal Science and Veterinary Medicine, Shandong Academy of Agricultural Sciences, No.202, Gongyebei Road, Jinan, 250100, China.,Shandong Ox Livestock Breeding Co., Ltd, Jinan, 250100, China
| | - Yujiao Wang
- Institute of Animal Science and Veterinary Medicine, Shandong Academy of Agricultural Sciences, No.202, Gongyebei Road, Jinan, 250100, China
| | - Wenjing Yang
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Yanqin Li
- Institute of Animal Science and Veterinary Medicine, Shandong Academy of Agricultural Sciences, No.202, Gongyebei Road, Jinan, 250100, China
| | - Xiuxin Zhao
- Shandong Ox Livestock Breeding Co., Ltd, Jinan, 250100, China
| | - Rongling Li
- Institute of Animal Science and Veterinary Medicine, Shandong Academy of Agricultural Sciences, No.202, Gongyebei Road, Jinan, 250100, China
| | - Yundong Gao
- Institute of Animal Science and Veterinary Medicine, Shandong Academy of Agricultural Sciences, No.202, Gongyebei Road, Jinan, 250100, China
| | - Wenbin Tuo
- Animal Parasitic Diseases Laboratory, BARC, USDA-ARS, Beltsville, MD, 20705, USA
| | - Ransom L Baldwin
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, MD, 20705, USA
| | - Cong-Jun Li
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, MD, 20705, USA.
| | - Lingzhao Fang
- MRC Human Genetics Unit at the Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, UK.
| | - George E Liu
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, MD, 20705, USA.
| |
Collapse
|
10
|
Li J, Wang Y, Mukiibi R, Karisa B, Plastow GS, Li C. Integrative analyses of genomic and metabolomic data reveal genetic mechanisms associated with carcass merit traits in beef cattle. Sci Rep 2022; 12:3389. [PMID: 35232965 PMCID: PMC8888742 DOI: 10.1038/s41598-022-06567-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2021] [Accepted: 02/01/2022] [Indexed: 11/09/2022] Open
Abstract
Improvement of carcass merit traits is a priority for the beef industry. Discovering DNA variants and genes associated with variation in these traits and understanding biological functions/processes underlying their associations are of paramount importance for more effective genetic improvement of carcass merit traits in beef cattle. This study integrates 10,488,742 imputed whole genome DNA variants, 31 plasma metabolites, and animal phenotypes to identify genes and biological functions/processes that are associated with carcass merit traits including hot carcass weight (HCW), rib eye area (REA), average backfat thickness (AFAT), lean meat yield (LMY), and carcass marbling score (CMAR) in a population of 493 crossbred beef cattle. Regression analyses were performed to identify plasma metabolites associated with the carcass merit traits, and the results showed that 4 (3-hydroxybutyric acid, acetic acid, citric acid, and choline), 6 (creatinine, L-glutamine, succinic acid, pyruvic acid, L-lactic acid, and 3-hydroxybutyric acid), 4 (fumaric acid, methanol, D-glucose, and glycerol), 2 (L-lactic acid and creatinine), and 5 (succinic acid, fumaric acid, lysine, glycine, and choline) plasma metabolites were significantly associated with HCW, REA, AFAT, LMY, and CMAR (P-value < 0.1), respectively. Combining the results of metabolome-genome wide association studies using the 10,488,742 imputed SNPs, 103, 160, 83, 43, and 109 candidate genes were identified as significantly associated with HCW, REA, AFAT, LMY, and CMAR (P-value < 1 × 10-5), respectively. By applying functional enrichment analyses for candidate genes of each trait, 26, 24, 26, 24, and 28 significant cellular and molecular functions were predicted for HCW, REA, AFAT, LMY, and CMAR, respectively. Among the five topmost significantly enriched biological functions for carcass merit traits, molecular transport and small molecule biochemistry were two top biological functions associated with all carcass merit traits. Lipid metabolism was the most significant biological function for LMY and CMAR and it was also the second and fourth highest biological function for REA and HCW, respectively. Candidate genes and enriched biological functions identified by the integrative analyses of metabolites with phenotypic traits and DNA variants could help interpret the results of previous genome-wide association studies for carcass merit traits. Our integrative study also revealed additional potential novel genes associated with these economically important traits. Therefore, our study improves understanding of the molecular and biological functions/processes that influence carcass merit traits, which could help develop strategies to enhance genomic prediction of carcass merit traits with incorporation of metabolomic data. Similarly, this information could guide management practices, such as nutritional interventions, with the purpose of boosting specific carcass merit traits.
Collapse
Affiliation(s)
- Jiyuan Li
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada
| | - Yining Wang
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada.,Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB, Canada
| | - Robert Mukiibi
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Edinburgh, Scotland, UK
| | - Brian Karisa
- Results Driven Agriculture Research, Edmonton, AB, Canada
| | - Graham S Plastow
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada.
| | - Changxi Li
- Department of Agricultural, Food and Nutritional Science, University of Alberta, Edmonton, AB, Canada. .,Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB, Canada.
| |
Collapse
|
11
|
Functional annotation of regulatory elements in cattle genome reveals the roles of extracellular interaction and dynamic change of chromatin states in rumen development during weaning. Genomics 2022; 114:110296. [PMID: 35143887 DOI: 10.1016/j.ygeno.2022.110296] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Revised: 12/20/2021] [Accepted: 02/01/2022] [Indexed: 12/24/2022]
Abstract
We profiled landscapes of bovine regulatory elements and explored dynamic changes of chromatin states in rumen development during weaning. The regulatory elements (15 chromatin states) and their coordinated activities in cattle were defined through genome-wide profiling of four histone modifications, CTCF-binding, DNA accessibility, DNA methylation, and transcriptome in rumen epithelial tissues. Each chromatin state presented specific enrichment for sequence ontology, methylation, trait-associated variants, transcription, gene expression-associated variants, selection signatures, and evolutionarily conserved elements. During weaning, weak enhancers and flanking active transcriptional start sites (TSS) were the most dynamic chromatin states and occurred in tandem with significant variations in gene expression and DNA methylation, significantly associated with stature, production, and reproduction economic traits. By comparing with in vitro cultured epithelial cells and in vivo rumen tissues, we showed the commonness and uniqueness of these results, especially the roles of cell interactions and mitochondrial activities in tissue development.
Collapse
|
12
|
Shi S, Zhang Z, Li B, Zhang S, Fang L. Incorporation of Trait-Specific Genetic Information into Genomic Prediction Models. Methods Mol Biol 2022; 2467:329-340. [PMID: 35451781 DOI: 10.1007/978-1-0716-2205-6_11] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Due to the rapid development of high-throughput sequencing technology, we can easily obtain not only the genetic variants at the whole-genome sequence level (e.g., from 1000 Genomes project and 1000 Bull Genomes project), but also a wide range of functional annotations (e.g., enhancers and promoters from ENCODE, FAANG, and FarmGTEx projects) across a wide range of tissues, cell types, developmental stages, and environmental conditions. This huge amount of information leads to a revolution in studying genetics and genomics of complex traits in humans, livestock, and plant species. In this chapter, we focused on and reviewed the genomic prediction methods that incorporate external biological information into genomic prediction, such as sequence ontology, linkage disequilibrium (LD) of SNPs, quantitative trait loci (QTL), and multi-layer omics data (e.g., transcriptome, epigenome, and microbiome).
Collapse
Affiliation(s)
- Shaolei Shi
- College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Zhe Zhang
- Department of Animal Breeding and genetics, College of Animal Science, South China Agricultural University (SCAU), Guangzhou, China
| | - Bingjie Li
- The Roslin Institute Building, Scotland's Rural College, Edinburgh, UK
| | - Shengli Zhang
- College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Lingzhao Fang
- MRC Human Genetics Unit at the Institute of Genetics and Cancer, The University of Edinburgh, Edinburgh, UK.
| |
Collapse
|
13
|
Long EM, Bradbury PJ, Romay MC, Buckler ES, Robbins KR. Genome-wide Imputation Using the Practical Haplotype Graph in the Heterozygous Crop Cassava. G3-GENES GENOMES GENETICS 2021; 12:6423990. [PMID: 34751380 PMCID: PMC8728015 DOI: 10.1093/g3journal/jkab383] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/09/2021] [Accepted: 10/14/2021] [Indexed: 11/13/2022]
Abstract
Genomic applications such as genomic selection and genome-wide association have become increasingly common since the advent of genome sequencing. The cost of sequencing has decreased in the past two decades; however, genotyping costs are still prohibitive to gathering large datasets for these genomic applications, especially in nonmodel species where resources are less abundant. Genotype imputation makes it possible to infer whole-genome information from limited input data, making large sampling for genomic applications more feasible. Imputation becomes increasingly difficult in heterozygous species where haplotypes must be phased. The practical haplotype graph (PHG) is a recently developed tool that can accurately impute genotypes, using a reference panel of haplotypes. We showcase the ability of the PHG to impute genomic information in the highly heterozygous crop cassava (Manihot esculenta). Accurately phased haplotypes were sampled from runs of homozygosity across a diverse panel of individuals to populate PHG, which proved more accurate than relying on computational phasing methods. The PHG achieved high imputation accuracy, using sparse skim-sequencing input, which translated to substantial genomic prediction accuracy in cross-validation testing. The PHG showed improved imputation accuracy, compared to a standard imputation tool Beagle, especially in predicting rare alleles.
Collapse
Affiliation(s)
- Evan M Long
- Plant Breeding and Genetics Section, School of Integrative Plant Science, Cornell University, Ithaca, NY 14853, USA
| | - Peter J Bradbury
- Institute for Genomic Diversity, Cornell University, Ithaca, NY 14853, USA.,United States Department of Agriculture-Agricultural Research Service, Robert W. Holley, Center for Agriculture and Health, Ithaca, NY 14853, USA
| | - M Cinta Romay
- Institute for Genomic Diversity, Cornell University, Ithaca, NY 14853, USA
| | - Edward S Buckler
- Plant Breeding and Genetics Section, School of Integrative Plant Science, Cornell University, Ithaca, NY 14853, USA.,Institute for Genomic Diversity, Cornell University, Ithaca, NY 14853, USA.,United States Department of Agriculture-Agricultural Research Service, Robert W. Holley, Center for Agriculture and Health, Ithaca, NY 14853, USA
| | - Kelly R Robbins
- Plant Breeding and Genetics Section, School of Integrative Plant Science, Cornell University, Ithaca, NY 14853, USA
| |
Collapse
|
14
|
Rice BR, Lipka AE. Diversifying maize genomic selection models. MOLECULAR BREEDING : NEW STRATEGIES IN PLANT IMPROVEMENT 2021; 41:33. [PMID: 37309328 PMCID: PMC10236107 DOI: 10.1007/s11032-021-01221-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/16/2020] [Accepted: 03/07/2021] [Indexed: 06/14/2023]
Abstract
Genomic selection (GS) is one of the most powerful tools available for maize breeding. Its use of genome-wide marker data to estimate breeding values translates to increased genetic gains with fewer breeding cycles. In this review, we cover the history of GS and highlight particular milestones during its adaptation to maize breeding. We discuss how GS can be applied to developing superior maize inbreds and hybrids. Additionally, we characterize refinements in GS models that could enable the encapsulation of non-additive genetic effects, genotype by environment interactions, and multiple levels of the biological hierarchy, all of which could ultimately result in more accurate predictions of breeding values. Finally, we suggest the stages in a maize breeding program where it would be beneficial to apply GS. Given the current sophistication of high-throughput phenotypic, genotypic, and other -omic level data currently available to the maize community, now is the time to explore the implications of their incorporation into GS models and thus ensure that genetic gains are being achieved as quickly and efficiently as possible.
Collapse
Affiliation(s)
- Brian R. Rice
- Department of Crop Sciences, University of Illinois, Urbana, IL USA
| | | |
Collapse
|
15
|
Gao Y, Fang L, Baldwin RL, Connor EE, Cole JB, Van Tassell CP, Ma L, Li CJ, Liu GE. Single-cell transcriptomic analyses of dairy cattle ruminal epithelial cells during weaning. Genomics 2021; 113:2045-2055. [PMID: 33933592 DOI: 10.1016/j.ygeno.2021.04.039] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2020] [Revised: 03/20/2021] [Accepted: 04/27/2021] [Indexed: 11/30/2022]
Abstract
Using the 10× Genomics Chromium Controller, we obtained scRNA-seq data of 5064 and 1372 individual cells from two Holstein calf ruminal epithelial tissues before and after weaning, respectively. We detected six distinct cell clusters, designated their cell types, and reported their marker genes. We then examined these clusters' underlining cell types and relationships by performing cell cycle, pseudotime trajectory, regulatory network, weighted gene co-expression network and gene ontology analyses. By integrating these cell marker genes with Holstein GWAS signals, we found they were enriched for animal production and body conformation traits. Finally, we confirmed their cell identities by comparing them with human and mouse stomach epithelial cells. This study presents an initial effort to implement single-cell transcriptomic analysis in cattle, and demonstrates ruminal tissue epithelial cell types and their developments during weaning, opening the door for new discoveries about tissue/cell type roles in complex traits at single-cell resolution.
Collapse
Affiliation(s)
- Yahui Gao
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, MD 20705, USA; Department of Animal and Avian Sciences, University of Maryland, College Park, MD 20742, USA.
| | - Lingzhao Fang
- MRC Human Genetics Unit at the Medical Research Council Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh EH4 2XU, United Kingdom.
| | - Ransom L Baldwin
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, MD 20705, USA.
| | - Erin E Connor
- Department of Animal and Food Sciences, University of Delaware, Newark, DE 19716, USA.
| | - John B Cole
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, MD 20705, USA.
| | - Curtis P Van Tassell
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, MD 20705, USA.
| | - Li Ma
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD 20742, USA.
| | - Cong-Jun Li
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, MD 20705, USA.
| | - George E Liu
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, MD 20705, USA.
| |
Collapse
|
16
|
Milkevych V, Karaman E, Sahana G, Janss L, Cai Z, Lund MS. MeSCoT: The tool for quantitative trait simulation through the mechanistic modelling of genes' regulatory interactions. G3-GENES GENOMES GENETICS 2021; 11:6255744. [PMID: 33905502 PMCID: PMC8496224 DOI: 10.1093/g3journal/jkab133] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Accepted: 04/10/2021] [Indexed: 11/21/2022]
Abstract
This work represents a novel mechanistic approach to simulate and study genomic networks with accompanying regulatory interactions and complex mechanisms of quantitative trait formation. The approach implemented in MeSCoT software is conceptually based on the omnigenic genetic model of quantitative (complex) trait, and closely imitates the basic in vivo mechanisms of quantitative trait realization. The software provides a framework to study molecular mechanisms of gene-by-gene and gene-by-environment interactions underlying quantitative trait’s realization and allows detailed mechanistic studies of impact of genetic and phenotypic variance on gene regulation. MeSCoT performs a detailed simulation of genes’ regulatory interactions for variable genomic architectures and generates complete set of transcriptional and translational data together with simulated quantitative trait values. Such data provide opportunities to study, for example, verification of novel statistical methods aiming to integrate intermediate phenotypes together with final phenotype in quantitative genetic analyses or to investigate novel approaches for exploiting gene-by-gene and gene-by-environment interactions.
Collapse
Affiliation(s)
- Viktor Milkevych
- Center for Quantitative Genetics and Genomics, Aarhus University, Tjele, Denmark
| | - Emre Karaman
- Center for Quantitative Genetics and Genomics, Aarhus University, Tjele, Denmark
| | - Goutam Sahana
- Center for Quantitative Genetics and Genomics, Aarhus University, Tjele, Denmark
| | - Luc Janss
- Center for Quantitative Genetics and Genomics, Aarhus University, Tjele, Denmark
| | - Zexi Cai
- Center for Quantitative Genetics and Genomics, Aarhus University, Tjele, Denmark
| | - Mogens Sandø Lund
- Center for Quantitative Genetics and Genomics, Aarhus University, Tjele, Denmark
| |
Collapse
|
17
|
Farooq M, van Dijk ADJ, Nijveen H, Aarts MGM, Kruijer W, Nguyen TP, Mansoor S, de Ridder D. Prior Biological Knowledge Improves Genomic Prediction of Growth-Related Traits in Arabidopsis thaliana. Front Genet 2021; 11:609117. [PMID: 33552126 PMCID: PMC7855462 DOI: 10.3389/fgene.2020.609117] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2020] [Accepted: 12/21/2020] [Indexed: 01/11/2023] Open
Abstract
Prediction of growth-related complex traits is highly important for crop breeding. Photosynthesis efficiency and biomass are direct indicators of overall plant performance and therefore even minor improvements in these traits can result in significant breeding gains. Crop breeding for complex traits has been revolutionized by technological developments in genomics and phenomics. Capitalizing on the growing availability of genomics data, genome-wide marker-based prediction models allow for efficient selection of the best parents for the next generation without the need for phenotypic information. Until now such models mostly predict the phenotype directly from the genotype and fail to make use of relevant biological knowledge. It is an open question to what extent the use of such biological knowledge is beneficial for improving genomic prediction accuracy and reliability. In this study, we explored the use of publicly available biological information for genomic prediction of photosynthetic light use efficiency (Φ PSII ) and projected leaf area (PLA) in Arabidopsis thaliana. To explore the use of various types of knowledge, we mapped genomic polymorphisms to Gene Ontology (GO) terms and transcriptomics-based gene clusters, and applied these in a Genomic Feature Best Linear Unbiased Predictor (GFBLUP) model, which is an extension to the traditional Genomic BLUP (GBLUP) benchmark. Our results suggest that incorporation of prior biological knowledge can improve genomic prediction accuracy for both Φ PSII and PLA. The improvement achieved depends on the trait, type of knowledge and trait heritability. Moreover, transcriptomics offers complementary evidence to the Gene Ontology for improvement when used to define functional groups of genes. In conclusion, prior knowledge about trait-specific groups of genes can be directly translated into improved genomic prediction.
Collapse
Affiliation(s)
- Muhammad Farooq
- Bioinformatics Group, Wageningen University, Wageningen, Netherlands
- Molecular Virology and Gene Silencing Lab, Agricultural Biotechnology Division, National Institute for Biotechnology and Genetic Engineering (NIBGE), Punjab, Pakistan
| | - Aalt D. J. van Dijk
- Bioinformatics Group, Wageningen University, Wageningen, Netherlands
- Biometris, Wageningen University, Wageningen, Netherlands
| | - Harm Nijveen
- Bioinformatics Group, Wageningen University, Wageningen, Netherlands
| | - Mark G. M. Aarts
- Laboratory of Genetics, Wageningen University, Wageningen, Netherlands
| | - Willem Kruijer
- Biometris, Wageningen University, Wageningen, Netherlands
| | - Thu-Phuong Nguyen
- Laboratory of Genetics, Wageningen University, Wageningen, Netherlands
| | - Shahid Mansoor
- Molecular Virology and Gene Silencing Lab, Agricultural Biotechnology Division, National Institute for Biotechnology and Genetic Engineering (NIBGE), Punjab, Pakistan
| | - Dick de Ridder
- Bioinformatics Group, Wageningen University, Wageningen, Netherlands
| |
Collapse
|
18
|
Xiang R, Breen EJ, Prowse-Wilkins CP, Chamberlain AJ, Goddard ME. Bayesian genome-wide analysis of cattle traits using variants with functional and evolutionary significance. ANIMAL PRODUCTION SCIENCE 2021. [DOI: 10.1071/an21061] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Context
Functional genomics studies have highlighted genomic regions with regulatory and evolutionary significance. Such information independent of association analysis may benefit fine-mapping and genomic selection of economically important traits. However, systematic evaluation of the use of functional information in mapping, and genomic selection of cattle traits, is lacking. Also, single-nucleotide polymorphisms (SNPs) from the high-density (HD) panel are known to tag informative variants, but the performance of genomic prediction using HD SNPs together with variants supported by different functional genomics is unknown.
Aims
We selected six sets of functionally important variants and modelled each set together with HD SNPs in Bayesian models to map and predict protein, fat and milk yield as well as mastitis, somatic cell count and temperament of dairy cattle.
Methods
Two models were used, namely (1) BayesR, which includes priors of four distribution of variant effects, and (2) BayesRC, which includes additional priors of different functional classes of variants. Bayesian models were trained in three breeds of 28 000 cows of Holstein, Jersey and Australian Red and predicted into 2600 independent bulls.
Key results
Adding functionally important variants significantly increased the enrichment of genetic variance explained for mapped variants, suggesting improved genome-wide mapping precision. Such improvement was significantly higher when the same set of variants was modelled by BayesRC than by BayesR. Combining functional variant sets with HD SNPs improves genomic prediction accuracy in the majority of the cases and such improvement was more common and stronger for non-Holstein breeds and traits such as mastitis, somatic cell count and temperament. In contrast, adding a large number of random sequence variants to HD SNPs reduces mapping precision and has a worse or similar prediction accuracy, compared with using HD SNPs alone to map or predict. While BayesRC tended to have better genomic prediction accuracy than did BayesR, the overall difference in prediction accuracy between the two models was insignificant.
Conclusions
Our findings demonstrated the usefulness of functional data in genomic mapping and prediction.
Implications
We have highlighted the need for effective tools exploiting complex functional datasets to improve genomic prediction.
Collapse
|
19
|
Botelho ME, Lopes MS, Mathur PK, Knol EF, Guimarães SEF, Marques DBD, Lopes PS, Silva FF, Veroneze R. Applying an association weight matrix in weighted genomic prediction of boar taint compounds. J Anim Breed Genet 2020; 138:442-453. [PMID: 33285013 DOI: 10.1111/jbg.12528] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2020] [Revised: 10/13/2020] [Accepted: 11/14/2020] [Indexed: 12/14/2022]
Abstract
Biological information regarding markers and gene association may be used to attribute different weights for single nucleotide polymorphism (SNP) in genome-wide selection. Therefore, we aimed to evaluate the predictive ability and the bias of genomic prediction using models that allow SNP weighting in the genomic relationship matrix (G) building, with and without incorporating biological information to obtain the weights. Firstly, we performed a genome-wide association studies (GWAS) in data set containing single- (SL) or a multi-line (ML) pig population for androstenone, skatole and indole levels. Secondly, 1%, 2%, 5%, 10%, 30% and 50% of the markers explaining the highest proportions of the genetic variance for each trait were selected to build gene networks through the association weight matrix (AWM) approach. The number of edges in the network was computed and used to derive weights for G (AWM-WssGBLUP). The single-step GBLUP (ssGBLUP) and weighted ssGBLUP (WssGBLUP) were used as standard scenarios. All scenarios presented predictive abilities different from zero; however, the great overlap in their confidences interval suggests no differences among scenarios. Most of scenarios of based on AWM provide overestimations for skatole in both SL and ML populations. On the other hand, the skatole and indole prediction were no biased in the ssGBLUP (S1) in both SL and ML populations. Most of scenarios based on AWM provide no biased predictions for indole in both SL and ML populations. In summary, using biological information through AWM matrix and gene networks to derive weights for genomic prediction resulted in no increase in predictive ability for boar taint compounds. In addition, this approach increased the number of analyses steps. Thus, we can conclude that ssGBLUP is most appropriate for the analysis of boar taint compounds in comparison with the weighted strategies used in the present work.
Collapse
Affiliation(s)
- Margareth E Botelho
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, Brazil
| | - Marcos S Lopes
- Topigs Norsvin, Curitiba, Brazil.,Topigs Norsvin Research Center, Beuningen, the Netherlands
| | | | - Egbert F Knol
- Topigs Norsvin Research Center, Beuningen, the Netherlands
| | | | - Daniele B D Marques
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, Brazil
| | - Paulo S Lopes
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, Brazil
| | - Fabyano F Silva
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, Brazil
| | - Renata Veroneze
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, Brazil
| |
Collapse
|
20
|
Lee HJ, Chung YJ, Jang S, Seo DW, Lee HK, Yoon D, Lim D, Lee SH. Genome-wide identification of major genes and genomic prediction using high-density and text-mined gene-based SNP panels in Hanwoo (Korean cattle). PLoS One 2020; 15:e0241848. [PMID: 33264312 PMCID: PMC7710051 DOI: 10.1371/journal.pone.0241848] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2020] [Accepted: 10/21/2020] [Indexed: 11/24/2022] Open
Abstract
It was hypothesized that single-nucleotide polymorphisms (SNPs) extracted from text-mined genes could be more tightly related to causal variant for each trait and that differentially weighting of this SNP panel in the GBLUP model could improve the performance of genomic prediction in cattle. Fitting two GRMs constructed by text-mined SNPs and SNPs except text-mined SNPs from 777k SNPs set (exp_777K) as different random effects showed better accuracy than fitting one GRM (Im_777K) for six traits (e.g. backfat thickness: + 0.002, eye muscle area: + 0.014, Warner–Bratzler Shear Force of semimembranosus and longissimus dorsi: + 0.024 and + 0.068, intramuscular fat content of semimembranosus and longissimus dorsi: + 0.008 and + 0.018). These results can suggest that attempts to incorporate text mining into genomic predictions seem valuable, and further study using text mining can be expected to present the significant results.
Collapse
Affiliation(s)
- Hyo Jun Lee
- Division of Animal and Dairy Science, Chungnam National University, Daejeon, Korea
| | - Yoon Ji Chung
- Division of Animal and Dairy Science, Chungnam National University, Daejeon, Korea
| | - Sungbong Jang
- Department of Animal and Dairy Science, University of Georgia, Athens, GA, United States of America
| | - Dong Won Seo
- Division of Animal and Dairy Science, Chungnam National University, Daejeon, Korea
| | - Hak Kyo Lee
- Department of Animal Biotechnology, Chonbuk National University, Jeonju, Korea
| | - Duhak Yoon
- Department of Animal Science, Kyungpook National University, Sangju, Korea
| | - Dajeong Lim
- Animal Genome & Bioinformatics, National Institute of Animal Science, Wanju, Korea
- * E-mail: (DL); (SHL)
| | - Seung Hwan Lee
- Division of Animal and Dairy Science, Chungnam National University, Daejeon, Korea
- * E-mail: (DL); (SHL)
| |
Collapse
|
21
|
Adoption and Optimization of Genomic Selection To Sustain Breeding for Apricot Fruit Quality. G3-GENES GENOMES GENETICS 2020; 10:4513-4529. [PMID: 33067307 PMCID: PMC7718743 DOI: 10.1534/g3.120.401452] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Genomic selection (GS) is a breeding approach which exploits genome-wide information and whose unprecedented success has shaped several animal and plant breeding schemes through delivering their genetic progress. This is the first study assessing the potential of GS in apricot (Prunus armeniaca) to enhance postharvest fruit quality attributes. Genomic predictions were based on a F1 pseudo-testcross population, comprising 153 individuals with contrasting fruit quality traits. They were phenotyped for physical and biochemical fruit metrics in contrasting climatic conditions over two years. Prediction accuracy (PA) varied from 0.31 for glucose content with the Bayesian LASSO (BL) to 0.78 for ethylene production with RR-BLUP, which yielded the most accurate predictions in comparison to Bayesian models and only 10% out of 61,030 SNPs were sufficient to reach accurate predictions. Useful insights were provided on the genetic architecture of apricot fruit quality whose integration in prediction models improved their performance, notably for traits governed by major QTL. Furthermore, multivariate modeling yielded promising outcomes in terms of PA within training partitions partially phenotyped for target traits. This provides a useful framework for the implementation of indirect selection based on easy-to-measure traits. Thus, we highlighted the main levers to take into account for the implementation of GS for fruit quality in apricot, but also to improve the genetic gain in perennial species.
Collapse
|
22
|
Ye S, Li J, Zhang Z. Multi-omics-data-assisted genomic feature markers preselection improves the accuracy of genomic prediction. J Anim Sci Biotechnol 2020; 11:109. [PMID: 33292577 PMCID: PMC7708144 DOI: 10.1186/s40104-020-00515-5] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2020] [Accepted: 09/22/2020] [Indexed: 12/02/2022] Open
Abstract
Background Presently, multi-omics data (e.g., genomics, transcriptomics, proteomics, and metabolomics) are available to improve genomic predictors. Omics data not only offers new data layers for genomic prediction but also provides a bridge between organismal phenotypes and genome variation that cannot be readily captured at the genome sequence level. Therefore, using multi-omics data to select feature markers is a feasible strategy to improve the accuracy of genomic prediction. In this study, simultaneously using whole-genome sequencing (WGS) and gene expression level data, four strategies for single-nucleotide polymorphism (SNP) preselection were investigated for genomic predictions in the Drosophila Genetic Reference Panel. Results Using genomic best linear unbiased prediction (GBLUP) with complete WGS data, the prediction accuracies were 0.208 ± 0.020 (0.181 ± 0.022) for the startle response and 0.272 ± 0.017 (0.307 ± 0.015) for starvation resistance in the female (male) lines. Compared with GBLUP using complete WGS data, both GBLUP and the genomic feature BLUP (GFBLUP) did not improve the prediction accuracy using SNPs preselected from complete WGS data based on the results of genome-wide association studies (GWASs) or transcriptome-wide association studies (TWASs). Furthermore, by using SNPs preselected from the WGS data based on the results of the expression quantitative trait locus (eQTL) mapping of all genes, only the startle response had greater accuracy than GBLUP with the complete WGS data. The best accuracy values in the female and male lines were 0.243 ± 0.020 and 0.220 ± 0.022, respectively. Importantly, by using SNPs preselected based on the results of the eQTL mapping of significant genes from TWAS, both GBLUP and GFBLUP resulted in great accuracy and small bias of genomic prediction. Compared with the GBLUP using complete WGS data, the best accuracy values represented increases of 60.66% and 39.09% for the starvation resistance and 27.40% and 35.36% for startle response in the female and male lines, respectively. Conclusions Overall, multi-omics data can assist genomic feature preselection and improve the performance of genomic prediction. The new knowledge gained from this study will enrich the use of multi-omics in genomic prediction.
Collapse
Affiliation(s)
- Shaopan Ye
- Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, National Engineering Research Centre for Breeding Swine Industry, College of Animal Science, South China Agricultural University, Guangzhou, Guangdong, China
| | - Jiaqi Li
- Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, National Engineering Research Centre for Breeding Swine Industry, College of Animal Science, South China Agricultural University, Guangzhou, Guangdong, China
| | - Zhe Zhang
- Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, National Engineering Research Centre for Breeding Swine Industry, College of Animal Science, South China Agricultural University, Guangzhou, Guangdong, China.
| |
Collapse
|
23
|
Genomic Prediction Informed by Biological Processes Expands Our Understanding of the Genetic Architecture Underlying Free Amino Acid Traits in Dry Arabidopsis Seeds. G3-GENES GENOMES GENETICS 2020; 10:4227-4239. [PMID: 32978264 PMCID: PMC7642941 DOI: 10.1534/g3.120.401240] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Plant growth, development, and nutritional quality depends upon amino acid homeostasis, especially in seeds. However, our understanding of the underlying genetics influencing amino acid content and composition remains limited, with only a few candidate genes and quantitative trait loci identified to date. Improved knowledge of the genetics and biological processes that determine amino acid levels will enable researchers to use this information for plant breeding and biological discovery. Toward this goal, we used genomic prediction to identify biological processes that are associated with, and therefore potentially influence, free amino acid (FAA) composition in seeds of the model plant Arabidopsis thaliana. Markers were split into categories based on metabolic pathway annotations and fit using a genomic partitioning model to evaluate the influence of each pathway on heritability explained, model fit, and predictive ability. Selected pathways included processes known to influence FAA composition, albeit to an unknown degree, and spanned four categories: amino acid, core, specialized, and protein metabolism. Using this approach, we identified associations for pathways containing known variants for FAA traits, in addition to finding new trait-pathway associations. Markers related to amino acid metabolism, which are directly involved in FAA regulation, improved predictive ability for branched chain amino acids and histidine. The use of genomic partitioning also revealed patterns across biochemical families, in which serine-derived FAAs were associated with protein related annotations and aromatic FAAs were associated with specialized metabolic pathways. Taken together, these findings provide evidence that genomic partitioning is a viable strategy to uncover the relative contributions of biological processes to FAA traits in seeds, offering a promising framework to guide hypothesis testing and narrow the search space for candidate genes.
Collapse
|
24
|
Hu G, Do DN, Gray J, Miar Y. Selection for Favorable Health Traits: A Potential Approach to Cope with Diseases in Farm Animals. Animals (Basel) 2020; 10:E1717. [PMID: 32971980 PMCID: PMC7552752 DOI: 10.3390/ani10091717] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2020] [Accepted: 09/21/2020] [Indexed: 12/17/2022] Open
Abstract
Disease is a global problem for animal farming industries causing tremendous economic losses (>USD 220 billion over the last decade) and serious animal welfare issues. The limitations and deficiencies of current non-selection disease control methods (e.g., vaccination, treatment, eradication strategy, genome editing, and probiotics) make it difficult to effectively, economically, and permanently eliminate the adverse influences of disease in the farm animals. These limitations and deficiencies drive animal breeders to be more concerned and committed to dealing with health problems in farm animals by selecting animals with favorable health traits. Both genetic selection and genomic selection contribute to improving the health of farm animals by selecting certain health traits (e.g., disease tolerance, disease resistance, and immune response), although both of them face some challenges. The objective of this review was to comprehensively review the potential of selecting health traits in coping with issues caused by diseases in farm animals. Within this review, we highlighted that selecting health traits can be applied as a method of disease control to help animal agriculture industries to cope with the adverse influences caused by diseases in farm animals. Certainly, the genetic/genomic selection solution cannot solve all the disease problems in farm animals. Therefore, management, vaccination, culling, medical treatment, and other measures must accompany selection solution to reduce the adverse impact of farm animal diseases on profitability and animal welfare.
Collapse
Affiliation(s)
| | | | | | - Younes Miar
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS B2N 5E3, Canada; (G.H.); (D.N.D.); (J.G.)
| |
Collapse
|
25
|
Yan Z, Huang H, Freebern E, Santos DJA, Dai D, Si J, Ma C, Cao J, Guo G, Liu GE, Ma L, Fang L, Zhang Y. Integrating RNA-Seq with GWAS reveals novel insights into the molecular mechanism underpinning ketosis in cattle. BMC Genomics 2020; 21:489. [PMID: 32680461 PMCID: PMC7367229 DOI: 10.1186/s12864-020-06909-z] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2020] [Accepted: 07/13/2020] [Indexed: 01/12/2023] Open
Abstract
Background Ketosis is a common metabolic disease during the transition period in dairy cattle, resulting in long-term economic loss to the dairy industry worldwide. While genetic selection of resistance to ketosis has been adopted by many countries, the genetic and biological basis underlying ketosis is poorly understood. Results We collected a total of 24 blood samples from 12 Holstein cows, including 4 healthy and 8 ketosis-diagnosed ones, before (2 weeks) and after (5 days) calving, respectively. We then generated RNA-Sequencing (RNA-Seq) data and seven blood biochemical indicators (bio-indicators) from leukocytes and plasma in each of these samples, respectively. By employing a weighted gene co-expression network analysis (WGCNA), we detected that 4 out of 16 gene-modules, which were significantly engaged in lipid metabolism and immune responses, were transcriptionally (FDR < 0.05) correlated with postpartum ketosis and several bio-indicators (e.g., high-density lipoprotein and low-density lipoprotein). By conducting genome-wide association signal (GWAS) enrichment analysis among six common health traits (ketosis, mastitis, displaced abomasum, metritis, hypocalcemia and livability), we found that 4 out of 16 modules were genetically (FDR < 0.05) associated with ketosis, among which three were correlated with postpartum ketosis based on WGCNA. We further identified five candidate genes for ketosis, including GRINA, MAF1, MAFA, C14H8orf82 and RECQL4. Our phenome-wide association analysis (Phe-WAS) demonstrated that human orthologues of these candidate genes were also significantly associated with many metabolic, endocrine, and immune traits in humans. For instance, MAFA, which is involved in insulin secretion, glucose response, and transcriptional regulation, showed a significantly higher association with metabolic and endocrine traits compared to other types of traits in humans. Conclusions In summary, our study provides novel insights into the molecular mechanism underlying ketosis in cattle, and highlights that an integrative analysis of omics data and cross-species mapping are promising for illustrating the genetic architecture underpinning complex traits.
Collapse
Affiliation(s)
- Ze Yan
- National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Hetian Huang
- National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China.,College of Animal Science and Veterinary Medicine, Henan Agricultural University, Zhengzhou, 450002, China
| | - Ellen Freebern
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD, 20742, USA
| | - Daniel J A Santos
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD, 20742, USA
| | - Dongmei Dai
- National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Jingfang Si
- National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Chong Ma
- Department of Clinical Veterinary Medicine, College of Veterinary Medicine, China Agricultural University, Beijing, 100193, China
| | - Jie Cao
- Department of Clinical Veterinary Medicine, College of Veterinary Medicine, China Agricultural University, Beijing, 100193, China
| | - Gang Guo
- Beijing Sunlon Livestock Development Co Ltd., Beijing, 100076, China
| | - George E Liu
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, Agricultural Research Service, USDA, Beltsville, MD, 20705, USA
| | - Li Ma
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD, 20742, USA.
| | - Lingzhao Fang
- MRC Human Genetics Unit at the Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, UK.
| | - Yi Zhang
- National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics, Breeding and Reproduction of Ministry of Agriculture and Rural Affairs, College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China.
| |
Collapse
|
26
|
Liu S, Yu Y, Zhang S, Cole JB, Tenesa A, Wang T, McDaneld TG, Ma L, Liu GE, Fang L. Epigenomics and genotype-phenotype association analyses reveal conserved genetic architecture of complex traits in cattle and human. BMC Biol 2020; 18:80. [PMID: 32620158 PMCID: PMC7334855 DOI: 10.1186/s12915-020-00792-6] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2019] [Accepted: 05/12/2020] [Indexed: 02/01/2023] Open
Abstract
Background Lack of comprehensive functional annotations across a wide range of tissues and cell types severely hinders the biological interpretations of phenotypic variation, adaptive evolution, and domestication in livestock. Here we used a combination of comparative epigenomics, genome-wide association study (GWAS), and selection signature analysis, to shed light on potential adaptive evolution in cattle. Results We cross-mapped 8 histone marks of 1300 samples from human to cattle, covering 178 unique tissues/cell types. By uniformly analyzing 723 RNA-seq and 40 whole genome bisulfite sequencing (WGBS) datasets in cattle, we validated that cross-mapped histone marks captured tissue-specific expression and methylation, reflecting tissue-relevant biology. Through integrating cross-mapped tissue-specific histone marks with large-scale GWAS and selection signature results, we for the first time detected relevant tissues and cell types for 45 economically important traits and artificial selection in cattle. For instance, immune tissues are significantly associated with health and reproduction traits, multiple tissues for milk production and body conformation traits (reflecting their highly polygenic architecture), and thyroid for the different selection between beef and dairy cattle. Similarly, we detected relevant tissues for 58 complex traits and diseases in humans and observed that immune and fertility traits in humans significantly correlated with those in cattle in terms of relevant tissues, which facilitated the identification of causal genes for such traits. For instance, PIK3CG, a gene highly specifically expressed in mononuclear cells, was significantly associated with both age-at-menopause in human and daughter-still-birth in cattle. ICAM, a T cell-specific gene, was significantly associated with both allergic diseases in human and metritis in cattle. Conclusion Collectively, our results highlighted that comparative epigenomics in conjunction with GWAS and selection signature analyses could provide biological insights into the phenotypic variation and adaptive evolution. Cattle may serve as a model for human complex traits, by providing additional information beyond laboratory model organisms, particularly when more novel phenotypes become available in the near future.
Collapse
Affiliation(s)
- Shuli Liu
- Animal Genomics and Improvement Laboratory, Agricultural Research Service, USDA, BARC-East, Beltsville, MD, 20705, USA.,College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Ying Yu
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Shengli Zhang
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - John B Cole
- Animal Genomics and Improvement Laboratory, Agricultural Research Service, USDA, BARC-East, Beltsville, MD, 20705, USA
| | - Albert Tenesa
- MRC Human Genetics Unit, Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, UK.,The Roslin Institute, University of Edinburgh, Edinburgh, EH25 9RG, UK
| | - Ting Wang
- Department of Genetics, Washington University School of Medicine, St. Louis, MO, 63110, USA
| | - Tara G McDaneld
- US Meat Animal Research Center, Agricultural Research Service, USDA, Clay Center, NE, 68933, USA
| | - Li Ma
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD, 20742, USA.
| | - George E Liu
- Animal Genomics and Improvement Laboratory, Agricultural Research Service, USDA, BARC-East, Beltsville, MD, 20705, USA.
| | - Lingzhao Fang
- Animal Genomics and Improvement Laboratory, Agricultural Research Service, USDA, BARC-East, Beltsville, MD, 20705, USA. .,MRC Human Genetics Unit, Institute of Genetics and Molecular Medicine, University of Edinburgh, Edinburgh, EH4 2XU, UK. .,Department of Animal and Avian Sciences, University of Maryland, College Park, MD, 20742, USA.
| |
Collapse
|
27
|
Fang L, Cai W, Liu S, Canela-Xandri O, Gao Y, Jiang J, Rawlik K, Li B, Schroeder SG, Rosen BD, Li CJ, Sonstegard TS, Alexander LJ, Van Tassell CP, VanRaden PM, Cole JB, Yu Y, Zhang S, Tenesa A, Ma L, Liu GE. Comprehensive analyses of 723 transcriptomes enhance genetic and biological interpretations for complex traits in cattle. Genome Res 2020; 30:790-801. [PMID: 32424068 PMCID: PMC7263193 DOI: 10.1101/gr.250704.119] [Citation(s) in RCA: 77] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2019] [Accepted: 05/01/2020] [Indexed: 12/13/2022]
Abstract
By uniformly analyzing 723 RNA-seq data from 91 tissues and cell types, we built a comprehensive gene atlas and studied tissue specificity of genes in cattle. We demonstrated that tissue-specific genes significantly reflected the tissue-relevant biology, showing distinct promoter methylation and evolution patterns (e.g., brain-specific genes evolve slowest, whereas testis-specific genes evolve fastest). Through integrative analyses of those tissue-specific genes with large-scale genome-wide association studies, we detected relevant tissues/cell types and candidate genes for 45 economically important traits in cattle, including blood/immune system (e.g., CCDC88C) for male fertility, brain (e.g., TRIM46 and RAB6A) for milk production, and multiple growth-related tissues (e.g., FGF6 and CCND2) for body conformation. We validated these findings by using epigenomic data across major somatic tissues and sperm. Collectively, our findings provided novel insights into the genetic and biological mechanisms underlying complex traits in cattle, and our transcriptome atlas can serve as a primary source for biological interpretation, functional validation, studies of adaptive evolution, and genomic improvement in livestock.
Collapse
Affiliation(s)
- Lingzhao Fang
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705, USA
- Department of Animal and Avian Sciences, University of Maryland, College Park, Maryland 20742, USA
- The Roslin Institute, Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Midlothian EH25 9RG, United Kingdom
- Medical Research Council Human Genetics Unit at the Medical Research Council Institute of Genetics and Molecular Medicine, The University of Edinburgh, Edinburgh EH4 2XU, United Kingdom
| | - Wentao Cai
- Department of Animal and Avian Sciences, University of Maryland, College Park, Maryland 20742, USA
- College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Shuli Liu
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705, USA
- College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Oriol Canela-Xandri
- The Roslin Institute, Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Midlothian EH25 9RG, United Kingdom
- Medical Research Council Human Genetics Unit at the Medical Research Council Institute of Genetics and Molecular Medicine, The University of Edinburgh, Edinburgh EH4 2XU, United Kingdom
| | - Yahui Gao
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705, USA
- Department of Animal and Avian Sciences, University of Maryland, College Park, Maryland 20742, USA
| | - Jicai Jiang
- Department of Animal and Avian Sciences, University of Maryland, College Park, Maryland 20742, USA
| | - Konrad Rawlik
- The Roslin Institute, Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Midlothian EH25 9RG, United Kingdom
| | - Bingjie Li
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705, USA
| | - Steven G Schroeder
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705, USA
| | - Benjamin D Rosen
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705, USA
| | - Cong-Jun Li
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705, USA
| | | | - Leeson J Alexander
- Fort Keogh Livestock and Range Research Laboratory, Agricultural Research Service, USDA, Miles City, Montana 59301, USA
| | - Curtis P Van Tassell
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705, USA
| | - Paul M VanRaden
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705, USA
| | - John B Cole
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705, USA
| | - Ying Yu
- College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Shengli Zhang
- College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Albert Tenesa
- The Roslin Institute, Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Midlothian EH25 9RG, United Kingdom
- Medical Research Council Human Genetics Unit at the Medical Research Council Institute of Genetics and Molecular Medicine, The University of Edinburgh, Edinburgh EH4 2XU, United Kingdom
| | - Li Ma
- Department of Animal and Avian Sciences, University of Maryland, College Park, Maryland 20742, USA
| | - George E Liu
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, Maryland 20705, USA
| |
Collapse
|
28
|
Causes, types, etiological agents, prevalence, diagnosis, treatment, prevention, effects on human health and future aspects of bovine mastitis. Anim Health Res Rev 2020; 21:36-49. [PMID: 32051050 DOI: 10.1017/s1466252319000094] [Citation(s) in RCA: 49] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Mastitis is among the most common and challenging diseases of dairy animals. It is an inflammation of udder tissues due to physical damage, chemical irritation, or infection caused by certain pathogens. Bovine mastitis has been known for ages, but its complex etiology and multi-factorial nature make it difficult to control. Mastitis may have a negative impact on human health by inducing antibiotic-resistant pathogens that may spread, which is threatening. Researchers are continuously struggling to devise suitable methods for mastitis control. Management strategies are mainly focused on disease prevention by farm management which includes proper hygiene, trained staff to monitor minor changes in the udder or milk, and better diagnostic and treatment methods. New technologies which have the potential to unravel this complicated disease include improved diagnostic tools, based on advanced genomics or proteomics, prevention, based on vaccines and immune modulators, and metabolic products of probiotics such as bacteriocins and gene therapy.
Collapse
|
29
|
Grigoletto L, Brito LF, Mattos EC, Eler JP, Bussiman FO, Silva BDCA, da Silva RP, Carvalho FE, Berton MP, Baldi F, Ferraz JBS. Genome-wide associations and detection of candidate genes for direct and maternal genetic effects influencing growth traits in the Montana Tropical® Composite population. Livest Sci 2019. [DOI: 10.1016/j.livsci.2019.09.013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
30
|
Song H, Ye S, Jiang Y, Zhang Z, Zhang Q, Ding X. Using imputation-based whole-genome sequencing data to improve the accuracy of genomic prediction for combined populations in pigs. Genet Sel Evol 2019; 51:58. [PMID: 31638889 PMCID: PMC6805481 DOI: 10.1186/s12711-019-0500-8] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2019] [Accepted: 10/07/2019] [Indexed: 11/17/2022] Open
Abstract
BACKGROUND For genomic selection in populations with a small reference population, combining populations of the same breed or populations of related breeds is an effective way to increase the size of the reference population. However, genomic predictions based on single nucleotide polymorphism (SNP)-chip genotype data using combined populations with different genetic backgrounds or from different breeds have not shown a clear advantage over using within-population or within-breed predictions. The increasing availability of whole-genome sequencing (WGS) data provides new opportunities for combined population genomic prediction. Our objective was to investigate the accuracy of genomic prediction using imputation-based WGS data from combined populations in pigs. Using 80K SNP panel genotypes, WGS genotypes, or genotypes on WGS variants that were pruned based on linkage disequilibrium (LD), three methods [genomic best linear unbiased prediction (GBLUP), single-step (ss)GBLUP, and genomic feature (GF)BLUP] were implemented with different prior information to identify the best method to improve the accuracy of genomic prediction for combined populations in pigs. RESULTS In total, 2089 and 2043 individuals with production and reproduction phenotypes, respectively, from three Yorkshire populations with different genetic backgrounds were genotyped with the PorcineSNP80 panel. Imputation accuracy from 80K to WGS variants reached 92%. The results showed that use of the WGS data compared to the 80K SNP panel did not increase the accuracy of genomic prediction in a single population, but using WGS data with LD pruning and GFBLUP with prior information did yield higher accuracy than the 80K SNP panel. For the 80K SNP panel genotypes, using the combined population resulted in a slight improvement, no change, or even a slight decrease in accuracy in comparison with the single population for GBLUP and ssGBLUP, while accuracy increased by 1 to 2.4% when using WGS data. Notably, the GFBLUP method did not perform well for both the combined population and the single populations. CONCLUSIONS The use of WGS data was beneficial for combined population genomic prediction. Simply increasing the number of SNPs to the WGS level did not increase accuracy for a single population, while using pruned WGS data based on LD and GFBLUP with prior information could yield higher accuracy than the 80K SNP panel.
Collapse
Affiliation(s)
- Hailiang Song
- Key Laboratory of Animal Genetics and Breeding of Ministry of Agriculture, National Engineering Laboratory of Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Shaopan Ye
- Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, National Engineering Research Centre for Breeding Swine Industry, College of Animal Science, South China Agricultural University, Guangzhou, Guangdong China
| | - Yifan Jiang
- Key Laboratory of Animal Genetics and Breeding of Ministry of Agriculture, National Engineering Laboratory of Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Zhe Zhang
- Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, National Engineering Research Centre for Breeding Swine Industry, College of Animal Science, South China Agricultural University, Guangzhou, Guangdong China
| | - Qin Zhang
- Shandong Provincial Key Laboratory of Animal Biotechnology and Disease Control and Prevention, Shandong Agricultural University, Taian, China
| | - Xiangdong Ding
- Key Laboratory of Animal Genetics and Breeding of Ministry of Agriculture, National Engineering Laboratory of Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| |
Collapse
|
31
|
Rohde PD, Jensen IR, Sarup PM, Ørsted M, Demontis D, Sørensen P, Kristensen TN. Genetic Signatures of Drug Response Variability in Drosophila melanogaster. Genetics 2019; 213:633-650. [PMID: 31455722 PMCID: PMC6781897 DOI: 10.1534/genetics.119.302381] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2019] [Accepted: 08/26/2019] [Indexed: 12/27/2022] Open
Abstract
Knowledge of the genetic basis underlying variation in response to environmental exposures or treatments is important in many research areas. For example, knowing the set of causal genetic variants for drug responses could revolutionize personalized medicine. We used Drosophila melanogaster to investigate the genetic signature underlying behavioral variability in response to methylphenidate (MPH), a drug used in the treatment of attention-deficit/hyperactivity disorder. We exposed a wild-type D. melanogaster population to MPH and a control treatment, and observed an increase in locomotor activity in MPH-exposed individuals. Whole-genome transcriptomic analyses revealed that the behavioral response to MPH was associated with abundant gene expression alterations. To confirm these patterns in a different genetic background and to further advance knowledge on the genetic signature of drug response variability, we used a system of inbred lines, the Drosophila Genetic Reference Panel (DGRP). Based on the DGRP, we showed that the behavioral response to MPH was strongly genotype-dependent. Using an integrative genomic approach, we incorporated known gene interactions into the genomic analyses of the DGRP, and identified putative candidate genes for variability in drug response. We successfully validated 71% of the investigated candidate genes by gene expression knockdown. Furthermore, we showed that MPH has cross-generational behavioral and transcriptomic effects. Our findings establish a foundation for understanding the genetic mechanisms driving genotype-specific responses to medical treatment, and highlight the opportunities that integrative genomic approaches have in optimizing medical treatment of complex diseases.
Collapse
Affiliation(s)
- Palle Duun Rohde
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830 Tjele, Denmark
- The Lundbeck Foundation Initiative for Integrative Psychiatric Research, iPSYCH, 8000 Aarhus C, Denmark
- Center for Integrative Sequencing, Aarhus University, 8000, Denmark
| | - Iben Ravnborg Jensen
- Section for Biology and Environmental Science, Department of Chemistry and Bioscience, Aalborg University, 9220, Denmark
| | - Pernille Merete Sarup
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830 Tjele, Denmark
| | - Michael Ørsted
- Section for Biology and Environmental Science, Department of Chemistry and Bioscience, Aalborg University, 9220, Denmark
| | - Ditte Demontis
- The Lundbeck Foundation Initiative for Integrative Psychiatric Research, iPSYCH, 8000 Aarhus C, Denmark
- Center for Integrative Sequencing, Aarhus University, 8000, Denmark
- Department of Biomedicine, Aarhus University, 8000, Denmark
| | - Peter Sørensen
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830 Tjele, Denmark
| | - Torsten Nygaard Kristensen
- Section for Biology and Environmental Science, Department of Chemistry and Bioscience, Aalborg University, 9220, Denmark
- Section for Genetics, Ecology and Evolution, Department of Bioscience, Aarhus University, 8000, Denmark
| |
Collapse
|
32
|
Fang L, Liu S, Liu M, Kang X, Lin S, Li B, Connor EE, Baldwin RL, Tenesa A, Ma L, Liu GE, Li CJ. Functional annotation of the cattle genome through systematic discovery and characterization of chromatin states and butyrate-induced variations. BMC Biol 2019; 17:68. [PMID: 31419979 PMCID: PMC6698049 DOI: 10.1186/s12915-019-0687-8] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2019] [Accepted: 08/05/2019] [Indexed: 12/20/2022] Open
Abstract
BACKGROUND The functional annotation of genomes, including chromatin accessibility and modifications, is important for understanding and effectively utilizing the increased amount of genome sequences reported. However, while such annotation has been well explored in a diverse set of tissues and cell types in human and model organisms, relatively little data are available for livestock genomes, hindering our understanding of complex trait variation, domestication, and adaptive evolution. Here, we present the first complete global landscape of regulatory elements in cattle and explore the dynamics of chromatin states in rumen epithelial cells induced by the rumen developmental regulator-butyrate. RESULTS We established the first global map of regulatory elements (15 chromatin states) and defined their coordinated activities in cattle, through genome-wide profiling for six histone modifications, RNA polymerase II, CTCF-binding sites, DNA accessibility, DNA methylation, and transcriptome in rumen epithelial primary cells (REPC), rumen tissues, and Madin-Darby bovine kidney epithelial cells (MDBK). We demonstrated that each chromatin state exhibited specific enrichment for sequence ontology, transcription, methylation, trait-associated variants, gene expression-associated variants, selection signatures, and evolutionarily conserved elements, implying distinct biological functions. After butyrate treatments, we observed that the weak enhancers and flanking active transcriptional start sites (TSS) were the most dynamic chromatin states, occurred concomitantly with significant alterations in gene expression and DNA methylation, which was significantly associated with heifer conception rate and stature economic traits. CONCLUSION Our results demonstrate the crucial role of functional genome annotation for understanding genome regulation, complex trait variation, and adaptive evolution in livestock. Using butyrate to induce the dynamics of the epigenomic landscape, we were able to establish the correlation among nutritional elements, chromatin states, gene activities, and phenotypic outcomes.
Collapse
Affiliation(s)
- Lingzhao Fang
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD 20742 USA
| | - Shuli Liu
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193 China
| | - Mei Liu
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
- College of Animal Science and Technology, Shaanxi Key Laboratory of Agricultural Molecular Biology, Northwest A&F University, Yangling, 712100 Shaanxi China
| | - Xiaolong Kang
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
- College of Agriculture, Ningxia University, Yinchuan, 750021 China
| | - Shudai Lin
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
- Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture, College of Animal Science of South China Agricultural University, Guangzhou, 510642 China
| | - Bingjie Li
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
| | - Erin E. Connor
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
| | - Ransom L. Baldwin
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
| | - Albert Tenesa
- The Roslin Institute, University of Edinburgh, Edinburgh, EH4 2XU UK
| | - Li Ma
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD 20742 USA
| | - George E. Liu
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
| | - Cong-jun Li
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705 USA
| |
Collapse
|
33
|
Fang L, Zhou Y, Liu S, Jiang J, Bickhart DM, Null DJ, Li B, Schroeder SG, Rosen BD, Cole JB, Van Tassell CP, Ma L, Liu GE. Integrating Signals from Sperm Methylome Analysis and Genome-Wide Association Study for a Better Understanding of Male Fertility in Cattle. EPIGENOMES 2019; 3:epigenomes3020010. [PMID: 34968233 PMCID: PMC8594688 DOI: 10.3390/epigenomes3020010] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2019] [Revised: 05/03/2019] [Accepted: 05/11/2019] [Indexed: 01/18/2023] Open
Abstract
Decreased male fertility is a big concern in both human society and the livestock industry. Sperm DNA methylation is commonly believed to be associated with male fertility. However, due to the lack of accurate male fertility records (i.e., limited mating times), few studies have investigated the comprehensive impacts of sperm DNA methylation on male fertility in mammals. In this study, we generated 10 sperm DNA methylomes and performed a preliminary correlation analysis between signals from sperm DNA methylation and signals from large-scale (n = 27,214) genome-wide association studies (GWAS) of 35 complex traits (including 12 male fertility-related traits). We detected genomic regions, which experienced DNA methylation alterations in sperm and were associated with aging and extreme fertility phenotypes (e.g., sire-conception rate or SCR). In dynamic hypomethylated regions (HMRs) and partially methylated domains (PMDs), we found genes (e.g., HOX gene clusters and microRNAs) that were involved in the embryonic development. We demonstrated that genomic regions, which gained rather than lost methylations during aging, and in animals with low SCR were significantly and selectively enriched for GWAS signals of male fertility traits. Our study discovered 16 genes as the potential candidate markers for male fertility, including SAMD5 and PDE5A. Collectively, this initial effort supported a hypothesis that sperm DNA methylation may contribute to male fertility in cattle and revealed the usefulness of functional annotations in enhancing biological interpretation and genomic prediction for complex traits and diseases.
Collapse
Affiliation(s)
- Lingzhao Fang
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705, USA
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD 20742, USA
| | - Yang Zhou
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction, Education Ministry of China, Huazhong Agricultural University, Wuhan 430070, China
| | - Shuli Liu
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705, USA
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture & National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Jicai Jiang
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD 20742, USA
| | - Derek M. Bickhart
- Dairy Forage Research Center, Agricultural Research Service, USDA, Madison, WI 53718, USA
| | - Daniel J. Null
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705, USA
| | - Bingjie Li
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705, USA
| | - Steven G. Schroeder
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705, USA
| | - Benjamin D. Rosen
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705, USA
| | - John B. Cole
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705, USA
| | - Curtis P. Van Tassell
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705, USA
| | - Li Ma
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD 20742, USA
- Correspondence: (L.M.); (G.E.L.); Tel.: +1-301-405-1389 (L.M.); +1-301-504-9843 (G.E.L.); Fax: +1-301-504-8414 (G.E.L.)
| | - George E. Liu
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD 20705, USA
- Correspondence: (L.M.); (G.E.L.); Tel.: +1-301-405-1389 (L.M.); +1-301-504-9843 (G.E.L.); Fax: +1-301-504-8414 (G.E.L.)
| |
Collapse
|
34
|
Fang L, Zhou Y, Liu S, Jiang J, Bickhart DM, Null DJ, Li B, Schroeder SG, Rosen BD, Cole JB, Van Tassell CP, Ma L, Liu GE. Comparative analyses of sperm DNA methylomes among human, mouse and cattle provide insights into epigenomic evolution and complex traits. Epigenetics 2019; 14:260-276. [PMID: 30810461 PMCID: PMC6557555 DOI: 10.1080/15592294.2019.1582217] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
Sperm DNA methylation is crucial for fertility and viability of offspring but epigenome evolution in mammals is largely understudied. By comparing sperm DNA methylomes and large-scale genome-wide association study (GWAS) signals between human and cattle, we aimed to examine the DNA methylome evolution and its associations with complex phenotypes in mammals. Our analysis revealed that genes with conserved non-methylated promoters (e.g., ANKS1A and WNT7A) among human and cattle were involved in common system and embryo development, and enriched for GWAS signals of body conformation traits in both species, while genes with conserved hypermethylated promoters (e.g., TCAP and CD80) were engaged in immune responses and highlighted by immune-related traits. On the other hand, genes with human-specific hypomethylated promoters (e.g., FOXP2 and HYDIN) were engaged in neuron system development and enriched for GWAS signals of brain-related traits, while genes with cattle-specific hypomethylated promoters (e.g., LDHB and DGAT2) mainly participated in lipid storage and metabolism. We validated our findings using sperm-retained nucleosome, preimplantation transcriptome, and adult tissue transcriptome data, as well as sequence evolutionary features, including motif binding sites, mutation rates, recombination rates and evolution signatures. In conclusion, our results demonstrate important roles of epigenome evolution in shaping the genetic architecture underlying complex phenotypes, hence enhance signal prioritization in GWAS and provide valuable information for human neurological disorders and livestock genetic improvement.
Collapse
Affiliation(s)
- Lingzhao Fang
- a Animal Genomics and Improvement Laboratory, BARC , Agricultural Research Service, USDA , Beltsville , MD , USA.,b Department of Animal and Avian Sciences , University of Maryland , College Park , MD , USA
| | - Yang Zhou
- c Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction, Education Ministry of China , Huazhong Agricultural University , Wuhan , Hubei , China
| | - Shuli Liu
- a Animal Genomics and Improvement Laboratory, BARC , Agricultural Research Service, USDA , Beltsville , MD , USA.,d Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture & National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology , China Agricultural University , Beijing , China
| | - Jicai Jiang
- b Department of Animal and Avian Sciences , University of Maryland , College Park , MD , USA
| | - Derek M Bickhart
- e Dairy Forage Research Center , Agricultural Research Service, USDA , Madison , WI , USA
| | - Daniel J Null
- a Animal Genomics and Improvement Laboratory, BARC , Agricultural Research Service, USDA , Beltsville , MD , USA
| | - Bingjie Li
- a Animal Genomics and Improvement Laboratory, BARC , Agricultural Research Service, USDA , Beltsville , MD , USA
| | - Steven G Schroeder
- a Animal Genomics and Improvement Laboratory, BARC , Agricultural Research Service, USDA , Beltsville , MD , USA
| | - Benjamin D Rosen
- a Animal Genomics and Improvement Laboratory, BARC , Agricultural Research Service, USDA , Beltsville , MD , USA
| | - John B Cole
- a Animal Genomics and Improvement Laboratory, BARC , Agricultural Research Service, USDA , Beltsville , MD , USA
| | - Curtis P Van Tassell
- a Animal Genomics and Improvement Laboratory, BARC , Agricultural Research Service, USDA , Beltsville , MD , USA
| | - Li Ma
- b Department of Animal and Avian Sciences , University of Maryland , College Park , MD , USA
| | - George E Liu
- a Animal Genomics and Improvement Laboratory, BARC , Agricultural Research Service, USDA , Beltsville , MD , USA
| |
Collapse
|
35
|
Fang L, Jiang J, Li B, Zhou Y, Freebern E, Vanraden PM, Cole JB, Liu GE, Ma L. Genetic and epigenetic architecture of paternal origin contribute to gestation length in cattle. Commun Biol 2019; 2:100. [PMID: 30886909 PMCID: PMC6418173 DOI: 10.1038/s42003-019-0341-6] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2018] [Accepted: 02/06/2019] [Indexed: 12/19/2022] Open
Abstract
The length of gestation can affect offspring health and performance. Both maternal and fetal effects contribute to gestation length; however, paternal contributions to gestation length remain elusive. Using genome-wide association study (GWAS) in 27,214 Holstein bulls with millions of gestation records, here we identify nine paternal genomic loci associated with cattle gestation length. We demonstrate that these GWAS signals are enriched in pathways relevant to embryonic development, and in differentially methylated regions between sperm samples with long and short gestation length. We reveal that gestation length shares genetic and epigenetic architecture in sperm with calving ability, body depth, and conception rate. While several candidate genes are detected in our fine-mapping analysis, we provide evidence indicating ZNF613 as a promising candidate for cattle gestation length. Collectively, our findings support that the paternal genome and epigenome can impact gestation length potentially through regulation of the embryonic development. Lingzhao Fang et al. studied the paternal genetic variants that affect gestational length in cattle. They found that paternal genes from pathways involved in embryonic development were associated with gestation length, and that these were often found in differentially methylated regions of the genome.
Collapse
Affiliation(s)
- Lingzhao Fang
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD, 20742, USA.,Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD, 20705, USA
| | - Jicai Jiang
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD, 20742, USA
| | - Bingjie Li
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD, 20705, USA
| | - Yang Zhou
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction, Education Ministry of China, Huazhong Agricultural University, 430070, Wuhan, Hubei, China
| | - Ellen Freebern
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD, 20742, USA
| | - Paul M Vanraden
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD, 20705, USA
| | - John B Cole
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD, 20705, USA
| | - George E Liu
- Animal Genomics and Improvement Laboratory, BARC, Agricultural Research Service, USDA, Beltsville, MD, 20705, USA
| | - Li Ma
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD, 20742, USA.
| |
Collapse
|
36
|
Abstract
Mastitis is a prevalent and costly disease on dairy farms. Improved management and hygiene can reduce the risk of infection by contagious or environmental pathogens, and genetic selection can confer permanent improvement in mastitis resistance. National veterinary recording systems in the Nordic countries have allowed direct selection for sire families with low incidence of clinical mastitis for 3 decades, whereas other countries have practiced indirect selection for lower somatic cell count. Recently, pooling of producer-recorded data from on-farm herd management software programs has enabled selection for reduced incidence of clinical mastitis in the United States and other leading dairy countries.
Collapse
Affiliation(s)
- Kent A Weigel
- Department of Dairy Science, University of Wisconsin-Madison, 1675 Observatory Drive, Madison, WI 53706-1205, USA.
| | - George E Shook
- Department of Dairy Science, University of Wisconsin-Madison, 1675 Observatory Drive, Madison, WI 53706-1205, USA
| |
Collapse
|
37
|
Genomic-polygenic and polygenic predictions for milk yield, fat yield, and age at first calving in Thai multibreed dairy population using genic and functional sets of genotypes. Livest Sci 2019. [DOI: 10.1016/j.livsci.2018.11.008] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
|
38
|
Pegolo S, Bergamaschi M, Gasperi F, Biasioli F, Cecchinato A, Bittante G. Integrated PTR-ToF-MS, GWAS and biological pathway analyses reveal the contribution of cow's genome to cheese volatilome. Sci Rep 2018; 8:17002. [PMID: 30451907 PMCID: PMC6242841 DOI: 10.1038/s41598-018-35323-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2018] [Accepted: 11/01/2018] [Indexed: 11/21/2022] Open
Abstract
Volatile organic compounds (VOCs) are small molecules that contribute to the distinctive flavour of cheese which is an important attribute for consumer acceptability. To investigate whether cow's genetic background might contribute to cheese volatilome, we carried out genome-wide association studies (GWAS) and pathway-based analyses for 173 spectrometric peaks tentatively associated with several VOCs obtained from proton-transfer-reaction mass spectrometry (PTR-ToF-MS) analyses of 1,075 model cheeses produced using raw whole-milk from Brown Swiss cows. Overall, we detected 186 SNPs associated with 120 traits, several of which mapped close to genes involved in protein (e.g. CSN3, GNRHR and FAM169A), fat (e.g. AGPAT3, SCD5, and GPAM) and carbohydrate (e.g. B3GNT2, B4GALT1, and PHKB) metabolism. Gene set enrichment analysis showed that pathways connected with proteolysis/amino acid metabolism (purine and nitrogen metabolism) as well as fat metabolism (long-term potentiation) and mammary gland function (tight junction) were overrepresented. Our results provide the first evidence of a putative link between cow's genes and cheese flavour and offer new insights into the role of potential candidate loci and the biological functions contributing to the cheese volatilome.
Collapse
Affiliation(s)
- Sara Pegolo
- Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padua, Viale dell'Università 16, 35020, Legnaro, Padua, Italy.
| | - Matteo Bergamaschi
- Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padua, Viale dell'Università 16, 35020, Legnaro, Padua, Italy
| | - Flavia Gasperi
- Department of Food Quality and Nutrition, Research and Innovation Centre, Fondazione Edmund Mach (FEM), Via E. Mach 1, 38010, San Michele all'Adige, TN, Italy
| | - Franco Biasioli
- Department of Food Quality and Nutrition, Research and Innovation Centre, Fondazione Edmund Mach (FEM), Via E. Mach 1, 38010, San Michele all'Adige, TN, Italy
| | - Alessio Cecchinato
- Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padua, Viale dell'Università 16, 35020, Legnaro, Padua, Italy
| | - Giovanni Bittante
- Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padua, Viale dell'Università 16, 35020, Legnaro, Padua, Italy
| |
Collapse
|
39
|
Okada D, Endo S, Matsuda H, Ogawa S, Taniguchi Y, Katsuta T, Watanabe T, Iwaisaki H. An intersection network based on combining SNP coassociation and RNA coexpression networks for feed utilization traits in Japanese Black cattle. J Anim Sci 2018; 96:2553-2566. [PMID: 29762780 DOI: 10.1093/jas/sky170] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2018] [Accepted: 05/11/2018] [Indexed: 11/12/2022] Open
Abstract
Genome-wide association studies (GWAS) of quantitative traits have detected numerous genetic associations, but they encounter difficulties in pinpointing prominent candidate genes and inferring gene networks. The present study used a systems genetics approach integrating GWAS results with external RNA-expression data to detect candidate gene networks in feed utilization and growth traits of Japanese Black cattle, which are matters of concern. A SNP coassociation network was derived from significant correlations between SNPs with effects estimated by GWAS across 7 phenotypic traits. The resulting network genes contained significant numbers of annotations related to the traits. Using bovine transcriptome data from a public database, an RNA coexpression network was inferred based on the similarity of expression patterns across different tissues. An intersection network was then generated by superimposing the SNP and RNA networks and extracting shared interactions. This intersection network contained 4 tissue-specific modules: nervous system, reproductive system, muscular system, and glands. To characterize the structure (topographical properties) of the 3 networks, their scale-free properties were evaluated, which revealed that the intersection network was the most scale-free. In the subnetwork containing the most connected transcription factors (URI1, ROCK2, and ETV6), most genes were widely expressed across tissues, and genes previously shown to be involved in the traits were found. Results indicated that the current approach might be used to construct a gene network that better reflects biological information, providing encouragement for the genetic dissection of economically important quantitative traits.
Collapse
Affiliation(s)
- Daigo Okada
- Faculty of Agriculture, Kyoto University, Kyoto, Japan
| | - Satoko Endo
- Graduate School of Agriculture, Kyoto University, Kyoto, Japan
| | | | | | - Yukio Taniguchi
- Graduate School of Agriculture, Kyoto University, Kyoto, Japan
| | | | - Toshio Watanabe
- National Livestock Breeding Center, Nishigo, Fukushima, Japan.,Shirakawa Institute of Animal Genetics, Japan Livestock Technology Association, Nishigo, Fukushima, Japan
| | | |
Collapse
|
40
|
MicroRNA-guided prioritization of genome-wide association signals reveals the importance of microRNA-target gene networks for complex traits in cattle. Sci Rep 2018; 8:9345. [PMID: 29921979 PMCID: PMC6008395 DOI: 10.1038/s41598-018-27729-y] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2018] [Accepted: 06/06/2018] [Indexed: 02/07/2023] Open
Abstract
MicroRNAs (miRNA) are key modulators of gene expression and so act as putative fine-tuners of complex phenotypes. Here, we hypothesized that causal variants of complex traits are enriched in miRNAs and miRNA-target networks. First, we conducted a genome-wide association study (GWAS) for seven functional and milk production traits using imputed sequence variants (13~15 million) and >10,000 animals from three dairy cattle breeds, i.e., Holstein (HOL), Nordic red cattle (RDC) and Jersey (JER). Second, we analyzed for enrichments of association signals in miRNAs and their miRNA-target networks. Our results demonstrated that genomic regions harboring miRNA genes were significantly (P < 0.05) enriched with GWAS signals for milk production traits and mastitis, and that enrichments within miRNA-target gene networks were significantly higher than in random gene-sets for the majority of traits. Furthermore, most between-trait and across-breed correlations of enrichments with miRNA-target networks were significantly greater than with random gene-sets, suggesting pleiotropic effects of miRNAs. Intriguingly, genes that were differentially expressed in response to mammary gland infections were significantly enriched in the miRNA-target networks associated with mastitis. All these findings were consistent across three breeds. Collectively, our observations demonstrate the importance of miRNAs and their targets for the expression of complex traits.
Collapse
|
41
|
Rohde PD, Østergaard S, Kristensen TN, Sørensen P, Loeschcke V, Mackay TFC, Sarup P. Functional Validation of Candidate Genes Detected by Genomic Feature Models. G3 (BETHESDA, MD.) 2018; 8:1659-1668. [PMID: 29519937 PMCID: PMC5940157 DOI: 10.1534/g3.118.200082] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/29/2018] [Accepted: 03/07/2018] [Indexed: 12/11/2022]
Abstract
Understanding the genetic underpinnings of complex traits requires knowledge of the genetic variants that contribute to phenotypic variability. Reliable statistical approaches are needed to obtain such knowledge. In genome-wide association studies, variants are tested for association with trait variability to pinpoint loci that contribute to the quantitative trait. Because stringent genome-wide significance thresholds are applied to control the false positive rate, many true causal variants can remain undetected. To ameliorate this problem, many alternative approaches have been developed, such as genomic feature models (GFM). The GFM approach tests for association of set of genomic markers, and predicts genomic values from genomic data utilizing prior biological knowledge. We investigated to what degree the findings from GFM have biological relevance. We used the Drosophila Genetic Reference Panel to investigate locomotor activity, and applied genomic feature prediction models to identify gene ontology (GO) categories predictive of this phenotype. Next, we applied the covariance association test to partition the genomic variance of the predictive GO terms to the genes within these terms. We then functionally assessed whether the identified candidate genes affected locomotor activity by reducing gene expression using RNA interference. In five of the seven candidate genes tested, reduced gene expression altered the phenotype. The ranking of genes within the predictive GO term was highly correlated with the magnitude of the phenotypic consequence of gene knockdown. This study provides evidence for five new candidate genes for locomotor activity, and provides support for the reliability of the GFM approach.
Collapse
Affiliation(s)
- Palle Duun Rohde
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830 Tjele, Denmark
- The Lundbeck Foundation Initiative for Integrative Psychiatric Research, iPSYCH, 8000 Aarhus, Denmark
- Center for Integrative Sequencing, Aarhus University, 8000 Aarhus, Denmark
| | - Solveig Østergaard
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830 Tjele, Denmark
| | - Torsten Nygaard Kristensen
- Section for Genetics, Ecology and Evolution, Department of Bioscience, Aarhus University, 8000 Aarhus, Denmark
- Section for Biology and Environmental Science, Department of Chemistry and Bioscience, Aalborg University, 9220 Aalborg, Denmark
| | - Peter Sørensen
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830 Tjele, Denmark
| | - Volker Loeschcke
- Section for Genetics, Ecology and Evolution, Department of Bioscience, Aarhus University, 8000 Aarhus, Denmark
| | - Trudy F C Mackay
- Department of Biological Sciences, North Carolina State University, Raleigh, North Carolina 27695
- Program in Genetics, North Carolina State University, Raleigh, North Carolina 27695
- W. M. Keck Center for Behavioral Biology, North Carolina State University, Raleigh, North Carolina 27695
| | - Pernille Sarup
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830 Tjele, Denmark
| |
Collapse
|
42
|
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
|
43
|
Tholey A, Taylor NL, Heazlewood JL, Bendixen E. We Are Not Alone: The iMOP Initiative and Its Roles in a Biology- and Disease-Driven Human Proteome Project. J Proteome Res 2017; 16:4273-4280. [DOI: 10.1021/acs.jproteome.7b00408] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Andreas Tholey
- Systematic Proteome Research & Bioanalytics, Institute for Experimental Medicine, Christian-Albrechts-Universität zu Kiel, 24105 Kiel, Germany
| | - Nicolas L. Taylor
- Australian
Research Council Centre of Excellence in Plant Energy Biology, School
of Molecular Sciences and Institute of Agriculture, The University of Western Australia, Crawley, Western Australia 6009, Australia
| | - Joshua L. Heazlewood
- School
of BioSciences, The University of Melbourne, Parkville, Victoria 3010, Australia
| | - Emøke Bendixen
- Department
of Molecular Biology and Genetics, Faculty of Science and Technology, Aarhus University, 8000 Aarhus, Denmark
| |
Collapse
|
44
|
Rohde PD, Gaertner B, Ward K, Sørensen P, Mackay TFC. Genomic Analysis of Genotype-by-Social Environment Interaction for Drosophila melanogaster Aggressive Behavior. Genetics 2017; 206:1969-1984. [PMID: 28550016 PMCID: PMC5560801 DOI: 10.1534/genetics.117.200642] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2017] [Accepted: 05/22/2017] [Indexed: 02/06/2023] Open
Abstract
Human psychiatric disorders such as schizophrenia, bipolar disorder, and attention-deficit/hyperactivity disorder often include adverse behaviors including increased aggressiveness. Individuals with psychiatric disorders often exhibit social withdrawal, which can further increase the probability of conducting a violent act. Here, we used the inbred, sequenced lines of the Drosophila Genetic Reference Panel (DGRP) to investigate the genetic basis of variation in male aggressive behavior for flies reared in a socialized and socially isolated environment. We identified genetic variation for aggressive behavior, as well as significant genotype-by-social environmental interaction (GSEI); i.e., variation among DGRP genotypes in the degree to which social isolation affected aggression. We performed genome-wide association (GWA) analyses to identify genetic variants associated with aggression within each environment. We used genomic prediction to partition genetic variants into gene ontology (GO) terms and constituent genes, and identified GO terms and genes with high prediction accuracies in both social environments and for GSEI. The top predictive GO terms significantly increased the proportion of variance explained, compared to prediction models based on all segregating variants. We performed genomic prediction across environments, and identified genes in common between the social environments that turned out to be enriched for genome-wide associated variants. A large proportion of the associated genes have previously been associated with aggressive behavior in Drosophila and mice. Further, many of these genes have human orthologs that have been associated with neurological disorders, indicating partially shared genetic mechanisms underlying aggression in animal models and human psychiatric disorders.
Collapse
Affiliation(s)
- Palle Duun Rohde
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830 Tjele, Denmark
- iPSYCH, The Lundbeck Foundation Initiative for Integrative Psychiatric Research, 8000 Aarhus, Denmark
- ISEQ, Center for Integrative Sequencing, Aarhus University, 8000 Aarhus, Denmark
| | - Bryn Gaertner
- Department of Biological Sciences, North Carolina State University, Raleigh, North Carolina 27695
- Program in Genetics, North Carolina State University, Raleigh, North Carolina 27695
- W.M. Keck Center for Behavioral Biology, North Carolina State University, Raleigh, North Carolina 27695
| | - Kirsty Ward
- Department of Biological Sciences, North Carolina State University, Raleigh, North Carolina 27695
- Program in Genetics, North Carolina State University, Raleigh, North Carolina 27695
- W.M. Keck Center for Behavioral Biology, North Carolina State University, Raleigh, North Carolina 27695
| | - Peter Sørensen
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Aarhus University, 8830 Tjele, Denmark
| | - Trudy F C Mackay
- Department of Biological Sciences, North Carolina State University, Raleigh, North Carolina 27695
- Program in Genetics, North Carolina State University, Raleigh, North Carolina 27695
- W.M. Keck Center for Behavioral Biology, North Carolina State University, Raleigh, North Carolina 27695
| |
Collapse
|