Effect of reduced genomic representation on using runs of homozygosity for inbreeding characterization

Abstract

ABSTRACTRuns of homozygosity (ROHs) are proxy for genomic Identical-by-Descent segments and are increasingly used to measure individual inbreeding. ROHs analyses are mostly carried out on SNPs-arrays and whole-genome-sequencing data. Softwares recurrently used for their detection usually assume that genomic positions which have not been genotyped are non-variant. This might be true for whole-genome-sequencing data, but not for reduced genomic representations and can lead to spurious ROHs detection. We simulated the outputs of whole-genome-sequencing, two SNP-arrays and RAD-sequencing for three populations with different sizes. We compare the results of ROHs calling with two softwares: PLINK and RZooRoH. We demonstrate that to obtain meaningful estimates of inbreeding coefficients, RZooRoH requires fraction of genome seven times smaller compared to PLINK. When the SNP density is above 20 SNPs/Mb for PLINK and 3 SNPs/Mb for RZooRoH, ranks of ROHs-based inbreeding coefficients are conserved among individuals. With reduced genomic representations, ROHs distributions are consistently biased towards an underestimation of the total numbers of small and an overestimation of the total numbers of large ROHs, except for RZooRoH and high-density SNPs-arrays. We conclude that both ROHs-based inbreeding coefficients and ROHs distributions exact quantification are highly dependent on the fraction of genome sequenced and should thus be treated with caution. However, relative inbreeding estimates, such as comparison between individuals or populations, are reliable with reduced genomic representations providing that the fraction of genome sequenced is large enough. Consequently, we advise researchers working with reduced genomic data to use SNPs-independent measures or model-based ROHs calling methods for inbreeding estimations.