Dataset Information

On the substructure controls in rare variant analysis: Principal components or variance components?

ABSTRACT: Recent studies showed that population substructure (PS) can have more complex impact on rare variant tests and that similarity-based collapsing tests (e.g., SKAT) may suffer more severely by PS than burden-based tests. In this work, we evaluate the performance of SKAT coupling with principal components (PC) or variance components (VC) based PS correction methods. We consider confounding effects caused by PS including stratified populations, admixed populations, and spatially distributed nongenetic risk; we investigate which types of variants (e.g., common, less frequent, rare, or all variants) should be used to effectively control for confounding effects. We found that (i) PC-based methods can account for confounding effects in most scenarios except for admixture, although the number of sufficient PCs depends on the PS complexity and the type of variants used. (ii) PCs based on all variants (i.e., common + less frequent + rare) tend to require equal or fewer sufficient PCs and often achieve higher power than PCs based on other variant types. (iii) VC-based methods can effectively adjust for confounding in all scenarios (even for admixture), though the type of variants should be used to construct VC may vary. (iv) VC based on all variants works consistently in all scenarios, though its power may be sometimes lower than VC based on other variant types. Given that the best-performed method and which variants to use depend on the underlying unknown confounding mechanisms, a robust strategy is to perform SKAT analyses using VC-based methods based on all variants.

SUBMITTER: Luo Y

PROVIDER: S-EPMC5851819 | biostudies-literature | 2018 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

On the substructure controls in rare variant analysis: Principal components or variance components?

Luo Yiwen Y Maity Arnab A Wu Michael C MC Smith Chris C Duan Qing Q Li Yun Y Tzeng Jung-Ying JY

Genetic epidemiology 20171226 3

Recent studies showed that population substructure (PS) can have more complex impact on rare variant tests and that similarity-based collapsing tests (e.g., SKAT) may suffer more severely by PS than burden-based tests. In this work, we evaluate the performance of SKAT coupling with principal components (PC) or variance components (VC) based PS correction methods. We consider confounding effects caused by PS including stratified populations, admixed populations, and spatially distributed nongenet ...[more]

PMID: 29280188

Similar Datasets

Project description:BackgroundAccurate inference of genetic ancestry is of fundamental interest to many biomedical, forensic, and anthropological research areas. Genetic ancestry memberships may relate to genetic disease risks. In a genome association study, failing to account for differences in genetic ancestry between cases and controls may also lead to false-positive results. Although a number of strategies for inferring and taking into account the confounding effects of genetic ancestry are available, applying them to large studies (tens thousands samples) is challenging. The goal of this study is to develop an approach for inferring genetic ancestry of samples with unknown ancestry among closely related populations and to provide accurate estimates of ancestry for application to large-scale studies.MethodsIn this study we developed a novel distance-based approach, Ancestry Inference using Principal component analysis and Spatial analysis (AIPS) that incorporates an Inverse Distance Weighted (IDW) interpolation method from spatial analysis to assign individuals to population memberships.ResultsWe demonstrate the benefits of AIPS in analyzing population substructure, specifically related to the four most commonly used tools EIGENSTRAT, STRUCTURE, fastSTRUCTURE, and ADMIXTURE using genotype data from various intra-European panels and European-Americans. While the aforementioned commonly used tools performed poorly in inferring ancestry from a large number of subpopulations, AIPS accurately distinguished variations between and within subpopulations.ConclusionsOur results show that AIPS can be applied to large-scale data sets to discriminate the modest variability among intra-continental populations as well as for characterizing inter-continental variation. The method we developed will protect against spurious associations when mapping the genetic basis of a disease. Our approach is more accurate and computationally efficient method for inferring genetic ancestry in the large-scale genetic studies.

Project description:To decipher gene regulatory networks, we used systematic suppression of 97 transcription factors (TFs) and 7 other genes with shRNA in mouse embryonic stem (ES) cells, followed by global gene expression profiling. Meta-analysis of these data together with the earlier data obtained by the induction of 50 TFs and the existing genome-wide data on TF binding sites identified the sets of regulated target genes for 23 TFs. Principal component analysis shows different roles of two groups of TFs that are active in ES cells: Pou5f1 and Sox2 support the expression of their target genes and prevent the upregulation of trophectoderm-related genes, whereas Esrrb, Sall4, Nanog, Gbx2, Grhl2, Mtf2, Aff1, Tcfap4, and Cdc5l support the expression of targets of Esrrb, including glycolysis genes, and prevent upregulation of targets of Trp53 and Polycomb TFs. If TFs from the second group are downregulated while Pou5f1 and Sox2 are still active, then the cell state changes towards epiblast lineages. Sequences for shRNA were designed to target 3 untranslated region of genes (Supplementary Table S8). Gene expression change was checked with real time qPCR (Supplementary Figure S7) and Westerm blot. ES cells (ES[MC1R(20)], passage 20) were cultured without feeders and were co-transfected with 1.6 g of shRNA expression vector and 0.4 g of pPyCAG-EGFP-IP carrying expression cassettes for puromycin resistant genes and EGFP using Effectene (Qiagen). Transfected cells were cultured in presence of 1.5 g/ml of puromycin and were harvested at 72 h after transfection. Mock cells were treated with transfection reagent without DNA and cultured in the absence of puromycin. Experiments were done in 3 replications with 2 of them used for gene expression profiling with microarrays. Total RNA was isolated by TRIzol (Invitrogen) after 2 days. Cy3-CTP labeled sample targets were prepared with total RNA by Low RNA Input Fluorescent Linear Amplification Kit (Agilent). Cy5-CTP labeled reference target was Stratagene Universal Mouse Reference RNA. Note that the processed data in the paper, which is also attached as GSE26520_Table_data.txt.gz, is slightly different from the values columns in each sample. The original processed data are normalized with a batch method, as the new batches of arrays added in the set. The value columns in this submission reflect full normalization as described in the data processing fields in each sample.

Dataset Information

On the substructure controls in rare variant analysis: Principal components or variance components?

Publications

On the substructure controls in rare variant analysis: Principal components or variance components?

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets