Weighted likelihood inference of genomic autozygosity patterns in dense genotype data
Preprint
Abstract
Background: Genomic regions of autozygosity (ROA) arise when an individual is homozygous for haplotypes inherited identical-by-descent from ancestors shared by both parents. Over the past decade, they have gained importance for understanding evolutionary history and the genetic basis of complex diseases and traits. However, methods to detect ROA in dense genotype data have not evolved in step with advances in genome technology that now enable us to rapidly create large high-resolution genotype datasets, limiting our ability to investigate their constituent ROA patterns. Results: We report a weighted likelihood approach for identifying ROA in dense genotype data that accounts for autocorrelation among genotyped positions and the possibilities of unobserved mutation and recombination events, and variability in the confidence of individual genotype calls in whole genome sequence (WGS) data. Forward-time genetic simulations under two demographic scenarios that reflect situations where inbreeding and its effect on fitness are of interest suggest this approach is better powered than existing state-of-the-art methods to detect ROA at marker densities consistent with WGS and popular microarray genotyping platforms used in human and non-human studies. Moreover, we present evidence that suggests this approach is able to distinguish ROA arising via consanguinity from ROA arising via endogamy. Using subsets of The 1000 Genomes Project Phase 3 data we show that, relative to WGS, intermediate and long ROA are captured robustly with popular microarray platforms, while detection of short ROA is more variable and improves with marker density. Worldwide ROA patterns inferred from WGS data are found to accord well with those previously reported on the basis of microarray genotype data. Finally, we highlight the potential of this approach to detect genomic regions enriched for autozygosity signals in one group relative to another based upon comparisons of per-individual autozygosity likelihoods instead of inferred ROA frequencies. Conclusions: This weighted likelihood ROA detection approach can assist population- and disease-geneticists working with a wide variety of data types and species to explore ROA patterns and to identify genomic regions with differential ROA signals among groups, thereby advancing our understanding of evolutionary history and the role of recessive variation in phenotypic variation and disease.Keywords
All Related Versions
- Published version: BMC Genomics, 18 (1), 928.
This publication has 229 references indexed in Scilit:
- Genomic Patterns of Homozygosity in Worldwide Human PopulationsThe American Journal of Human Genetics, 2012
- Consanguinity and the risk of congenital heart diseaseAmerican Journal of Medical Genetics Part A, 2012
- Haplotype variation and genotype imputation in African populationsGenetic Epidemiology, 2011
- Inference of Unexpected Genetic Relatedness among Individuals in HapMap Phase IIIThe American Journal of Human Genetics, 2010
- Population structure and genome-wide patterns of variation in Ireland and BritainEuropean Journal of Human Genetics, 2010
- High-Resolution Detection of Identity by Descent in Unrelated IndividualsThe American Journal of Human Genetics, 2010
- The Genome-wide Patterns of Variation Expose Significant Substructure in a Founder PopulationThe American Journal of Human Genetics, 2008
- Genome‐wide autozygosity mapping in human populationsGenetic Epidemiology, 2008
- Runs of Homozygosity in European PopulationsThe American Journal of Human Genetics, 2008
- Natural Selection on Genes that Underlie Human Disease SusceptibilityCurrent Biology, 2008