Dataset Information

Inference of multiple-wave admixtures by length distribution of ancestral tracks.

ABSTRACT: The ancestral tracks in admixed genomes are valuable for population history inference. While a few methods have been developed to infer admixture history based on ancestral tracks, these methods suffer the same flaw: only population admixture history under some specific models can be inferred. In addition, the inference of history might be biased or even unreliable if the specific model deviates from the real situation. To address this problem, we firstly proposed a general discrete admixture model to describe the admixture history with multiple ancestral populations and multiple-wave admixtures. We next deduced the length distribution of ancestral tracks under the general discrete admixture model. We further developed a new method, MultiWaver, to explore multiple-wave admixture histories. Our method could automatically determine an optimal admixture model based on the length distribution of ancestral tracks, and estimate the corresponding parameters under this optimal model. Specifically, we used a likelihood ratio test (LRT) to determine the number of admixture waves, and implemented an expectation-maximization (EM) algorithm to estimate parameters. We used simulation studies to validate the reliability and effectiveness of our method. Finally, good performance was observed when our method was applied to real data sets of African Americans and Mexicans, and new insights were gained into the admixture history of Uyghurs and Hazaras.

SUBMITTER: Ni X

PROVIDER: S-EPMC5997750 | biostudies-literature | 2018 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Inference of multiple-wave admixtures by length distribution of ancestral tracks.

Ni Xumin X Yuan Kai K Yang Xiong X Feng Qidi Q Guo Wei W Ma Zhiming Z Xu Shuhua S

Heredity 20180123 1

The ancestral tracks in admixed genomes are valuable for population history inference. While a few methods have been developed to infer admixture history based on ancestral tracks, these methods suffer the same flaw: only population admixture history under some specific models can be inferred. In addition, the inference of history might be biased or even unreliable if the specific model deviates from the real situation. To address this problem, we firstly proposed a general discrete admixture mo ...[more]

PMID: 29358727

Similar Datasets

Project description:The complex correlation structure of a collection of orthologous DNA sequences is uniquely captured by the "ancestral recombination graph" (ARG), a complete record of coalescence and recombination events in the history of the sample. However, existing methods for ARG inference are computationally intensive, highly approximate, or limited to small numbers of sequences, and, as a consequence, explicit ARG inference is rarely used in applied population genomics. Here, we introduce a new algorithm for ARG inference that is efficient enough to apply to dozens of complete mammalian genomes. The key idea of our approach is to sample an ARG of [Formula: see text] chromosomes conditional on an ARG of [Formula: see text] chromosomes, an operation we call "threading." Using techniques based on hidden Markov models, we can perform this threading operation exactly, up to the assumptions of the sequentially Markov coalescent and a discretization of time. An extension allows for threading of subtrees instead of individual sequences. Repeated application of these threading operations results in highly efficient Markov chain Monte Carlo samplers for ARGs. We have implemented these methods in a computer program called ARGweaver. Experiments with simulated data indicate that ARGweaver converges rapidly to the posterior distribution over ARGs and is effective in recovering various features of the ARG for dozens of sequences generated under realistic parameters for human populations. In applications of ARGweaver to 54 human genome sequences from Complete Genomics, we find clear signatures of natural selection, including regions of unusually ancient ancestry associated with balancing selection and reductions in allele age in sites under directional selection. The patterns we observe near protein-coding genes are consistent with a primary influence from background selection rather than hitchhiking, although we cannot rule out a contribution from recurrent selective sweeps.

Project description:CHD7 is a member of the chromodomain helicase DNA binding domain family of ATP-dependent chromatin remodeling enzymes. De novo mutation of the CHD7 gene is a major cause of CHARGE syndrome, a genetic disease characterized by a complex constellation of birth defects. To gain insight to the function of CHD7, we mapped the distribution of the CHD7 protein on chromatin using the approach of chromatin immunoprecipitation on tiled microarrays (ChIP-chip). These studies were performed in human colorectal carcinoma cells, human neuroblastoma cells, and mouse embryonic stem (ES) cells before and after differentiation into neural precursor cells. The results indicate that CHD7 localizes to discrete locations along chromatin that are specific to each cell type, and that the cell-specific binding of CHD7 correlates with a subset of histone H3 methylated at lysine 4 (H3K4me). The CHD7 sites change concomitantly with H3K4me patterns during ES cell differentiation, suggesting that H3K4me is part of the epigenetic signature that defines lineage-specific association of CHD7 with specific sites on chromatin. Furthermore, the CHD7 sites are predominantly located distal to transcription start sites, most often contained within DNase hypersensitive sites, frequently conserved, and near genes expressed at relatively high levels. These features are similar to those of gene enhancer elements, raising the possibility that CHD7 functions in enhancer mediated transcription, and that the congenital anomalies in CHARGE syndrome are due to alterations in transcription of tissue-specific genes normally regulated by CHD7 during development. ChIP-chip experiments were performed for CHD7 and H3K4 mono-,di-, and trimethylation modifications in 4 cells types: human DLD1 and SH-SY5Y; mouse ES and differentiated neural precursor cells derived from mouse ES cells. Microarrays used in these experiments tiled all or subset of ENCODE regions (in mouse, analogous ENCODE regions were assayed). At least two biological replicates were performed for each CHD7 ChIP experiment; H3K4 ChIP's were performed once in each cell type.

Dataset Information

Inference of multiple-wave admixtures by length distribution of ancestral tracks.

Publications

Inference of multiple-wave admixtures by length distribution of ancestral tracks.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure