A computational framework for improving genetic variants identification from 5,061 sheep sequencing data

Journal of Animal Science and Biotechnology

Research

Shangqian Xie¹  Gabrielle Becker¹  Brenda M. Murdoch¹  Karissa Isaacs² 

[1] Department of Animal, Veterinary & Food Sciences, University of Idaho, Moscow, ID, USA;Superior Farms, California, USA;

关键词: Computational framework;    Genetic variants;    Multiple samples;    Sheep;   

DOI  :  10.1186/s40104-023-00923-3

 received in 2023-03-25, accepted in 2023-08-01,  发布年份 2023

来源: Springer

【 摘 要 】

BackgroundPan-genomics is a recently emerging strategy that can be utilized to provide a more comprehensive characterization of genetic variation. Joint calling is routinely used to combine identified variants across multiple related samples. However, the improvement of variants identification using the mutual support information from multiple samples remains quite limited for population-scale genotyping.ResultsIn this study, we developed a computational framework for joint calling genetic variants from 5,061 sheep by incorporating the sequencing error and optimizing mutual support information from multiple samples' data. The variants were accurately identified from multiple samples by using four steps: (1) Probabilities of variants from two widely used algorithms, GATK and Freebayes, were calculated by Poisson model incorporating base sequencing error potential; (2) The variants with high mapping quality or consistently identified from at least two samples by GATK and Freebayes were used to construct the raw high-confidence identification (rHID) variants database; (3) The high confidence variants identified in single sample were ordered by probability value and controlled by false discovery rate (FDR) using rHID database; (4) To avoid the elimination of potentially true variants from rHID database, the variants that failed FDR were reexamined to rescued potential true variants and ensured high accurate identification variants. The results indicated that the percent of concordant SNPs and Indels from Freebayes and GATK after our new method were significantly improved 12%−32% compared with raw variants and advantageously found low frequency variants of individual sheep involved several traits including nipples number (GPC5), scrapie pathology (PAPSS2), seasonal reproduction and litter size (GRM1), coat color (RAB27A), and lentivirus susceptibility (TMEM154).ConclusionThe new method used the computational strategy to reduce the number of false positives, and simultaneously improve the identification of genetic variants. This strategy did not incur any extra cost by using any additional samples or sequencing data information and advantageously identified rare variants which can be important for practical applications of animal breeding.

【 授权许可】

CC BY   
© Chinese Association of Animal Science and Veterinary Medicine 2023

【 图 表 】

Fig. 3

Fig. 1

Fig. 6

12864_2015_2198_Article_IEq1.gif

Fig. 5

12951_2017_270_Article_IEq8.gif

Fig. 1

13100_2023_302_Article_IEq2.gif

Fig. 1

Fig. 3

12951_2015_111_Article_IEq2.gif

Fig. 1

Fig. 5

12951_2015_155_Article_IEq25.gif

Fig. 2

Fig. 3

12951_2015_111_Article_IEq1.gif