Date

2022-10

Author

Taç, İdil
Işıldak, Ulaş
Çubukcu, Hande
Karadavut, Damla
Vural, Kıvılcım Başak
Altınışık, Ezgi
Erdal, Yılmaz Selim
Somel, Mehmet
Özer, Füsun
Yet, İdil
Kılınç, Gülşah Merve

Metadata

Show full item record

This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

Item Usage Stats

175
views

83
downloads

Cite This

Ancient DNA is a field of study with genomic material obtained from biological samples that are highly damaged and not stored under special conditions. However, with advances in both sequencing techniques and bioinformatics methods such as imputation, which improve the quality of the data obtained, we have the opportunity to make maximum use of the material we have. When working with ancient genomes, which are mostly low-coverage ancient genomes, we cannot follow the pipeline we apply to modern genome datasets. While modern genomes can reliably perform diploid variant calling due to their high coverage, ancient genomes with <1 coverage often undergo a pseudo-haploidization procedure. In this study, we used different approaches using pseudo-haploidization or read depths to calculate and compare the allele frequencies of neutral and disease related SNPs in Neolithic Anatolian individuals. We calculated the minor allele frequency for each SNP using (i) pseudo-haploidization with random base selection using the ANGSD tool followed by pseudo-haploidization with the Plink program and (ii) diploid variant calling with samtools -mpileup followed by an algorithm that uses read depths and finds the values that maximize the log-likelihood calculated for each population using the binomial probability distribution. With the second method we use, our aim is to maximize the use of low-coverage ancient genome sequence data by using all the information we obtain about the locus of interest from variant calling without pseudo-haploidization in ancient DNA. One-way ANOVA test was used to test whether there was a significant difference between the methods (p value <2e-16 for neutral SNPs, <2e-16 for type 2 diabetes-related SNPs), Tukey test was performed as a post-hoc test; p values: method 1 vs method 2 (0.00) for neutral SNPs, (3.64e-11) for phenotype-related SNPs.

URI

https://hibit2022.ims.metu.edu.tr
https://hdl.handle.net/11511/101887

Collections

Graduate School of Informatics, Conference / Seminar