1
|
Ogawa S, Zoda A, Kagawa R, Obinata R. Comparing Methods to Select Candidates for Re-Genotyping to Impute Higher-Density Genotype Data in a Japanese Black Cattle Population: A Case Study. Animals (Basel) 2023; 13:ani13040638. [PMID: 36830425 PMCID: PMC9951718 DOI: 10.3390/ani13040638] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2022] [Revised: 02/04/2023] [Accepted: 02/10/2023] [Indexed: 02/15/2023] Open
Abstract
As optimization methods to identify the best animals for dense genotyping to construct a reference population for genotype imputation, the MCA and MCG methods, which use the pedigree-based additive genetic relationship matrix (A matrix) and the genomic relationship matrix (G matrix), respectively, have been proposed. We assessed the performance of MCA and MCG methods using 575 Japanese Black cows. Pedigree data were provided to trace back up to five generations to construct the A matrix with changing the pedigree depth from 1 to 5 (five MCA methods). Genotype information on 36,426 single-nucleotide polymorphisms was used to calculate the G matrix based on VanRaden's methods 1 and 2 (two MCG methods). The MCG always selected one cow per iteration, while MCA sometimes selected multiple cows. The number of commonly selected cows between the MCA and MCG methods was generally lower than that between different MCA methods or between different MCG methods. For the studied population, MCG appeared to be more reasonable than MCA in selecting cows as a reference population for higher-density genotype imputation to perform genomic prediction and a genome-wide association study.
Collapse
Affiliation(s)
- Shinichiro Ogawa
- Division of Meat Animal and Poultry Research, Institute of Livestock and Grassland Science, National Agriculture and Food Research Organization (NARO), Tsukuba 305-0901, Japan
- Correspondence: ; Tel.: +81-29-838-8627
| | - Atsushi Zoda
- Research and Development Group, Zen-Noh Embryo Transfer Center, Kamishihoro 080-1407, Japan
| | - Rino Kagawa
- Research and Development Group, Zen-Noh Embryo Transfer Center, Kamishihoro 080-1407, Japan
| | - Rui Obinata
- Research and Development Group, Zen-Noh Embryo Transfer Center, Kamishihoro 080-1407, Japan
| |
Collapse
|
2
|
Assessment of linkage disequilibrium patterns between structural variants and single nucleotide polymorphisms in three commercial chicken populations. BMC Genomics 2022; 23:193. [PMID: 35264116 PMCID: PMC8908679 DOI: 10.1186/s12864-022-08418-7] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2021] [Accepted: 02/24/2022] [Indexed: 12/29/2022] Open
Abstract
BACKGROUND Structural variants (SV) are causative for some prominent phenotypic traits of livestock as different comb types in chickens or color patterns in pigs. Their effects on production traits are also increasingly studied. Nevertheless, accurately calling SV remains challenging. It is therefore of interest, whether close-by single nucleotide polymorphisms (SNPs) are in strong linkage disequilibrium (LD) with SVs and can serve as markers. Literature comes to different conclusions on whether SVs are in LD to SNPs on the same level as SNPs to other SNPs. The present study aimed to generate a precise SV callset from whole-genome short-read sequencing (WGS) data for three commercial chicken populations and to evaluate LD patterns between the called SVs and surrounding SNPs. It is thereby the first study that assessed LD between SVs and SNPs in chickens. RESULTS The final callset consisted of 12,294,329 bivariate SNPs, 4,301 deletions (DEL), 224 duplications (DUP), 218 inversions (INV) and 117 translocation breakpoints (BND). While average LD between DELs and SNPs was at the same level as between SNPs and SNPs, LD between other SVs and SNPs was strongly reduced (DUP: 40%, INV: 27%, BND: 19% of between-SNP LD). A main factor for the reduced LD was the presence of local minor allele frequency differences, which accounted for 50% of the difference between SNP - SNP and DUP - SNP LD. This was potentially accompanied by lower genotyping accuracies for DUP, INV and BND compared with SNPs and DELs. An evaluation of the presence of tag SNPs (SNP in highest LD to the variant of interest) further revealed DELs to be slightly less tagged by WGS SNPs than WGS SNPs by other SNPs. This difference, however, was no longer present when reducing the pool of potential tag SNPs to SNPs located on four different chicken genotyping arrays. CONCLUSIONS The results implied that genomic variance due to DELs in the chicken populations studied can be captured by different SNP marker sets as good as variance from WGS SNPs, whereas separate SV calling might be advisable for DUP, INV, and BND effects.
Collapse
|
3
|
Li YF, Li YH, Su SS, Reif JC, Qi ZM, Wang XB, Wang X, Tian Y, Li DL, Sun RJ, Liu ZX, Xu ZJ, Fu GH, Ji YL, Chen QS, Liu JQ, Qiu LJ. SoySNP618K array: A high-resolution single nucleotide polymorphism platform as a valuable genomic resource for soybean genetics and breeding. JOURNAL OF INTEGRATIVE PLANT BIOLOGY 2022; 64:632-648. [PMID: 34914170 DOI: 10.1111/jipb.13202] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/05/2021] [Accepted: 12/05/2021] [Indexed: 05/13/2023]
Abstract
Innovations in genomics have enabled the development of low-cost, high-resolution, single nucleotide polymorphism (SNP) genotyping arrays that accelerate breeding progress and support basic research in crop science. Here, we developed and validated the SoySNP618K array (618,888 SNPs) for the important crop soybean. The SNPs were selected from whole-genome resequencing data containing 2,214 diverse soybean accessions; 29.34% of the SNPs mapped to genic regions representing 86.85% of the 56,044 annotated high-confidence genes. Identity-by-state analyses of 318 soybeans revealed 17 redundant accessions, highlighting the potential of the SoySNP618K array in supporting gene bank management. The patterns of population stratification and genomic regions enriched through domestication were highly consistent with previous findings based on resequencing data, suggesting that the ascertainment bias in the SoySNP618K array was largely compensated for. Genome-wide association mapping in combination with reported quantitative trait loci enabled fine-mapping of genes known to influence flowering time, E2 and GmPRR3b, and of a new candidate gene, GmVIP5. Moreover, genomic prediction of flowering and maturity time in 502 recombinant inbred lines was highly accurate (>0.65). Thus, the SoySNP618K array is a valuable genomic tool that can be used to address many questions in applied breeding, germplasm management, and basic crop research.
Collapse
Affiliation(s)
- Yan-Fei Li
- The National Key Facility for Crop Gene Resources and Genetic Improvement (NFCRI)/Key Laboratory of Crop Gene Resource and Germplasm Enhancement (MOA)/Key Laboratory of Soybean Biology (Beijing) (MOA), Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Ying-Hui Li
- The National Key Facility for Crop Gene Resources and Genetic Improvement (NFCRI)/Key Laboratory of Crop Gene Resource and Germplasm Enhancement (MOA)/Key Laboratory of Soybean Biology (Beijing) (MOA), Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Shan-Shan Su
- Beijing Compass Biotechnology Co. Ltd, Beijing, 102206, China
| | - Jochen C Reif
- Department of Breeding Research, Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), Gatersleben, 06466, Germany
| | - Zhao-Ming Qi
- Key Laboratory of Soybean Biology in Chinese Ministry of Education (Key Laboratory of Soybean Biology and Breeding/Genetics of Chinese Agriculture Ministry), Northeast Agricultural University, Harbin, 150030, China
| | - Xiao-Bo Wang
- School of Agronomy, Anhui Agricultural University, Hefei, 230036, China
| | - Xing Wang
- Xuzhou Institute of Agricultural Sciences of Xu-huai Region of Jiangsu, Xuzhou, 221131, China
| | - Yu Tian
- The National Key Facility for Crop Gene Resources and Genetic Improvement (NFCRI)/Key Laboratory of Crop Gene Resource and Germplasm Enhancement (MOA)/Key Laboratory of Soybean Biology (Beijing) (MOA), Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - De-Lin Li
- The National Key Facility for Crop Gene Resources and Genetic Improvement (NFCRI)/Key Laboratory of Crop Gene Resource and Germplasm Enhancement (MOA)/Key Laboratory of Soybean Biology (Beijing) (MOA), Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
- Department of Plant Genetics and Breeding, China Agricultural University, Beijing, 100193, China
| | - Ru-Jian Sun
- The National Key Facility for Crop Gene Resources and Genetic Improvement (NFCRI)/Key Laboratory of Crop Gene Resource and Germplasm Enhancement (MOA)/Key Laboratory of Soybean Biology (Beijing) (MOA), Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
- Key Laboratory of Soybean Biology in Chinese Ministry of Education (Key Laboratory of Soybean Biology and Breeding/Genetics of Chinese Agriculture Ministry), Northeast Agricultural University, Harbin, 150030, China
- Hulun Buir Institution of Agricultural Sciences, Zhalantun, Inner Mongolia, 021000, China
| | - Zhang-Xiong Liu
- The National Key Facility for Crop Gene Resources and Genetic Improvement (NFCRI)/Key Laboratory of Crop Gene Resource and Germplasm Enhancement (MOA)/Key Laboratory of Soybean Biology (Beijing) (MOA), Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Ze-Jun Xu
- Xuzhou Institute of Agricultural Sciences of Xu-huai Region of Jiangsu, Xuzhou, 221131, China
| | - Guang-Hui Fu
- Suzhou Academy of Agricultural Sciences, Suzhou, 234000, China
| | - Ya-Liang Ji
- Beijing Compass Biotechnology Co. Ltd, Beijing, 102206, China
| | - Qing-Shan Chen
- Key Laboratory of Soybean Biology in Chinese Ministry of Education (Key Laboratory of Soybean Biology and Breeding/Genetics of Chinese Agriculture Ministry), Northeast Agricultural University, Harbin, 150030, China
| | - Ji-Qiang Liu
- Beijing Compass Biotechnology Co. Ltd, Beijing, 102206, China
| | - Li-Juan Qiu
- The National Key Facility for Crop Gene Resources and Genetic Improvement (NFCRI)/Key Laboratory of Crop Gene Resource and Germplasm Enhancement (MOA)/Key Laboratory of Soybean Biology (Beijing) (MOA), Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| |
Collapse
|
4
|
de Faria DA, do Prado Paim T, Dos Santos CA, Paiva SR, Nogueira MB, McManus C. Selection signatures for heat tolerance in Brazilian horse breeds. Mol Genet Genomics 2022; 297:449-462. [PMID: 35150300 DOI: 10.1007/s00438-022-01862-w] [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] [Received: 03/23/2021] [Accepted: 01/24/2022] [Indexed: 12/14/2022]
Abstract
Since domestication, horse breeds have adapted to their environments and differentiated from one another. This paper uses two methods to detect selection signatures in 23 horse breeds, eight of which are Brazilian (610 animals), both cold-blooded and warm-blooded, from temperate and tropical regions. These animals were genotyped using the GGP Equine BeadChip and we analysed the data by Principal Component Analysis (PCA). The samples were separated into groups based on their geographical area of origin and PCA results studied. The genomic regions under selection were detected by hapFLK and PCAdapt methodologies, identifying six regions under selection with at least one Brazilian horse breed. These regions contain genes associated with heat tolerance, skin colour, body size, energy production/metabolism, genes involved in protein degradation/turnover/DNA repair, genes reducing the impact of oxidative stress/cellular repair, and transcriptional regulation. This work confirmed LCORL and NCAPG gene regions in previous studies associated with body size on Equine Chromosome Autosome 3 (ECA3). On the same ECA3, a region implicating genes linked to coat colour was identified, also previously related to heat stress. Regions with genes coding heat shock proteins were found on ECA1 and 2, and many candidate genes for oxidation-reduction which are a natural response to heat stress. However, a larger sample size and whole-genome SNPs are needed to understand better and identify new candidate regions as well as their functional relation with heat tolerance.
Collapse
Affiliation(s)
- Danielle Assis de Faria
- Faculdade de Agronomia e Veterinária, Instituto Central de Ciências, Campus Darcy Ribeiro, Universidade de Brasília, Asa Norte, Brasília, DF, 70910-900, Brazil
| | - Tiago do Prado Paim
- Instituto Federal de Educação, Ciência e Tecnologia Goiano, Rodovia Sul Goiana, Km 01, Zona Rural, Rio Verde, GO, 75901-970, Brazil
| | - Camila Alves Dos Santos
- Instituto Federal de Educação, Ciência e Tecnologia Goiano, Rodovia Sul Goiana, Km 01, Zona Rural, Rio Verde, GO, 75901-970, Brazil
| | - Samuel Rezende Paiva
- Embrapa Recursos Genéticos e Biotecnologia, Final W5 Norte, Brasília, DF, 70770-917, Brasil
| | - Marcelo Bchara Nogueira
- Faculdade de Agronomia e Veterinária, Instituto Central de Ciências, Campus Darcy Ribeiro, Universidade de Brasília, Asa Norte, Brasília, DF, 70910-900, Brazil
| | - Concepta McManus
- Departamento de Ciências Fisiológicas, Instituto de Biologia, Campus Darcy Ribeiro, Universidade de Brasília, Asa Norte, Brasília, DF, 70910-900, Brazil.
| |
Collapse
|