Identification of gene-gene interaction using principal components.
Ontology highlight
ABSTRACT:
SUBMITTER: Li J
PROVIDER: S-EPMC2795980 | biostudies-literature | 2009 Dec
REPOSITORIES: biostudies-literature
Publications
Identification of gene-gene interaction using principal components.
Li Jia J Tang Rui R Biernacka Joanna M JM de Andrade Mariza M
BMC proceedings 20091215
After more than 200 genome-wide association studies, there have been some successful identifications of a single novel locus. Thus, the identification of single-nucleotide polymorphisms (SNP) with interaction effects is of interest. Using the Genetic Analysis Workshop 16 data from the North American Rheumatoid Arthritis Consortium, we propose an approach to screen for SNP-SNP interaction using a two-stage method and an approach for detecting gene-gene interactions using principal components. We ...[more]
Similar Datasets
Project description:The literature provides strong evidence that stock price values can be predicted from past price data. Principal component analysis (PCA) identifies a small number of principle components that explain most of the variation in a data set. This method is often used for dimensionality reduction and analysis of the data. In this paper, we develop a general method for stock price prediction using time-varying covariance information. To address the time-varying nature of financial time series, we assign exponential weights to the price data so that recent data points are weighted more heavily. Our proposed method involves a dimension-reduction operation constructed based on principle components. Projecting the noisy observation onto a principle subspace results in a well-conditioned problem. We illustrate our results based on historical daily price data for 150 companies from different market-capitalization categories. We compare the performance of our method to two other methods: Gauss-Bayes, which is numerically demanding, and moving average, a simple method often used by technical traders and researchers. We investigate the results based on mean squared error and directional change statistic of prediction, as measures of performance, and volatility of prediction as a measure of risk.
| S-EPMC7083277 | biostudies-literature
Project description:Heterogeneity in the phenotypic mean and variance across populations is often observed for complex traits. One way to understand heterogeneous phenotypes lies in uncovering heterogeneity in genetic effects. Previous studies on genetic heterogeneity across populations were typically based on discrete groups in populations stratified by different countries or cohorts, which ignored the difference of population characteristics for the individuals within each group and resulted in loss of information. Here, we introduce a novel concept of genotype-by-population (G × P) interaction where population is defined by the first and second ancestry principal components (PCs), which are less likely to be confounded with country/cohort-specific factors. We applied a reaction norm model fitting each of 70 complex traits with significant SNP-heritability and the PCs as covariates to examine G × P interactions across diverse populations including white British and other white Europeans from the UK Biobank (N = 22,229). Our results demonstrated a significant population genetic heterogeneity for behavioral traits such as age at first sexual intercourse and academic qualification. Our approach may shed light on the latent genetic architecture of complex traits that underlies the modulation of genetic effects across different populations.
| S-EPMC7186421 | biostudies-literature
Project description:BackgroundThere are many methods for analyzing microarray data that group together genes having similar patterns of expression over all conditions tested. However, in many instances the biologically important goal is to identify relatively small sets of genes that share coherent expression across only some conditions, rather than all or most conditions as required in traditional clustering; e.g. genes that are highly up-regulated and/or down-regulated similarly across only a subset of conditions. Equally important is the need to learn which conditions are the decisive ones in forming such gene sets of interest, and how they relate to diverse conditional covariates, such as disease diagnosis or prognosis.ResultsWe present a method for automatically identifying such candidate sets of biologically relevant genes using a combination of principal components analysis and information theoretic metrics. To enable easy use of our methods, we have developed a data analysis package that facilitates visualization and subsequent data mining of the independent sources of significant variation present in gene microarray expression datasets (or in any other similarly structured high-dimensional dataset). We applied these tools to two public datasets, and highlight sets of genes most affected by specific subsets of conditions (e.g. tissues, treatments, samples, etc.). Statistically significant associations for highlighted gene sets were shown via global analysis for Gene Ontology term enrichment. Together with covariate associations, the tool provides a basis for building testable hypotheses about the biological or experimental causes of observed variation.ConclusionWe provide an unsupervised data mining technique for diverse microarray expression datasets that is distinct from major methods now in routine use. In test uses, this method, based on publicly available gene annotations, appears to identify numerous sets of biologically relevant genes. It has proven especially valuable in instances where there are many diverse conditions (10's to hundreds of different tissues or cell types), a situation in which many clustering and ordering algorithms become problematic. This approach also shows promise in other topic domains such as multi-spectral imaging datasets.
| S-EPMC1501050 | biostudies-literature
Project description:Studies of complex human diseases and traits associated with candidate genes are potentially vulnerable to bias (confounding) due to population stratification and inbreeding, especially in admixed population. In GWAS, the principal components (PCs) method provides a global ancestry value per subject, allowing corrections for population stratification. However, these coefficients are typically estimated assuming unrelated individuals, and if family structure is present and ignored, such substructures may induce artifactual PCs. Extensions of the PCs method have been proposed by Konishi and Rao [Biometrika 1992;79:631-641], taking into account only siblings' relatedness, and by Oualkacha et al. [Stat Appl Genet Mol Biol 2012, DOI: 10.2202/1544-6115.1711], taking into account large pedigrees and high-dimensional phenotype data. In this work, we extend these methods to estimate the global individual ancestry coefficients from PCs derived from different variance component matrix estimators using SNPs from two simulated data sets and two real data sets: the GENOA sibship data consisting of European and African-American subjects and the Baependi Heart Study consisting of 80 extended Brazilian families, both with genotyping data from the Affymetrix 6.0 chip. Our results show that the family structure plays an important role in the estimation of the global individual ancestry value for extended pedigrees but not for sibships.
| S-EPMC4583840 | biostudies-literature
Project description:Functional principal components (FPC) analysis is widely used to decompose and express functional observations. Curve estimates implicitly condition on basis functions and other quantities derived from FPC decompositions; however these objects are unknown in practice. In this article, we propose a method for obtaining correct curve estimates by accounting for uncertainty in FPC decompositions. Additionally, pointwise and simultaneous confidence intervals that account for both model- and decomposition-based variability are constructed. Standard mixed model representations of functional expansions are used to construct curve estimates and variances conditional on a specific decomposition. Iterated expectation and variance formulas combine model-based conditional estimates across the distribution of decompositions. A bootstrap procedure is implemented to understand the uncertainty in principal component decomposition quantities. Our method compares favorably to competing approaches in simulation studies that include both densely and sparsely observed functions. We apply our method to sparse observations of CD4 cell counts and to dense white-matter tract profiles. Code for the analyses and simulations is publicly available, and our method is implemented in the R package refund on CRAN.
| S-EPMC3962763 | biostudies-literature
Project description:We propose a novel method for estimating population-level and subject-specific effects of covariates on the variability of functional data. We extend the functional principal components analysis framework by modeling the variance of principal component scores as a function of covariates and subject-specific random effects. In a setting where principal components are largely invariant across subjects and covariate values, modeling the variance of these scores provides a flexible and interpretable way to explore factors that affect the variability of functional data. Our work is motivated by a novel dataset from an experiment assessing upper extremity motor control, and quantifies the reduction in motion variance associated with skill learning.
| S-EPMC6223649 | biostudies-literature
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.
2013-03-14 | E-GEOD-26520 | biostudies-arrayexpress
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.
2013-03-14 | GSE26520 | GEO
Project description:BACKGROUND/AIMS:Association studies using unrelated individuals have become the most popular design for mapping complex traits. One of the major challenges of association mapping is avoiding spurious association due to population stratification. Principal component analysis (PCA) on genome-wide marker genotypes is one of the most popular population stratification control methods. It implicitly assumes that the markers are in linkage equilibrium, a condition that is rarely satisfied and that we plan to relax. METHODS:We carefully examined the impact of linkage disequilibrium (LD) on PCA, and proposed a simple modification of the standard PCA to automatically adjust for the correlations among markers. RESULTS:We demonstrated that LD patterns in genome-wide association datasets can distort the techniques for stratification control, showing 'subpopulations' reflecting localized LD phenomena rather than plausible population structure. We showed that the proposed method effectively removes the artifactual effect of LD patterns, and successfully recovers underlying population structure that is not apparent from standard PCA. CONCLUSION:PCA is highly influenced by sets of SNPs with high LD, obscuring the true population substructure. Our shrinkage PCA applies to all available markers, regardless of the LD patterns. The proposed method is easier to implement than most existing LD adjusted PCA methods.
| S-EPMC2912642 | biostudies-literature
Project description:The analysis of microbiome data has several technical challenges. In particular, count matrices contain a large proportion of zeros, some of which are biological, whereas others are technical. Furthermore, the measurements suffer from unequal sequencing depth, overdispersion, and data redundancy. These nuisance factors introduce substantial noise. We propose an accurate and robust method, mbDenoise, for denoising microbiome data. Assuming a zero-inflated probabilistic PCA (ZIPPCA) model, mbDenoise uses variational approximation to learn the latent structure and recovers the true abundance levels using the posterior, borrowing information across samples and taxa. mbDenoise outperforms state-of-the-art methods to extract the signal for downstream analyses.
| S-EPMC9011970 | biostudies-literature