Dataset Information

Exploring predictive biomarkers from clinical genome-wide association studies via multidimensional hierarchical mixture models.

ABSTRACT: Although the detection of predictive biomarkers is of particular importance for the development of accurate molecular diagnostics, conventional statistical analyses based on gene-by-treatment interaction tests lack sufficient statistical power for this purpose, especially in large-scale clinical genome-wide studies that require an adjustment for multiplicity of a huge number of tests. Here we demonstrate an alternative efficient multi-subgroup screening method using multidimensional hierarchical mixture models developed to overcome this issue, with application to stroke and breast cancer randomized clinical trials with genomic data. We show that estimated effect size distributions of single nucleotide polymorphisms (SNPs) associated with outcomes, which could provide clues for exploring predictive biomarkers, optimizing individualized treatments, and understanding biological mechanisms of diseases. Furthermore, using this method we detected three new SNPs that are associated with blood homocysteine levels, which are strongly associated with the risk of stroke. We also detected six new SNPs that are associated with progression-free survival in breast cancer patients.

SUBMITTER: Otani T

PROVIDER: S-EPMC6303260 | biostudies-literature | 2019 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Exploring predictive biomarkers from clinical genome-wide association studies via multidimensional hierarchical mixture models.

Otani Takahiro T Noma Hisashi H Sugasawa Shonosuke S Kuchiba Aya A Goto Atsushi A Yamaji Taiki T Kochi Yuta Y Iwasaki Motoki M Matsui Shigeyuki S Tsunoda Tatsuhiko T

European journal of human genetics : EJHG 20180910 1

Although the detection of predictive biomarkers is of particular importance for the development of accurate molecular diagnostics, conventional statistical analyses based on gene-by-treatment interaction tests lack sufficient statistical power for this purpose, especially in large-scale clinical genome-wide studies that require an adjustment for multiplicity of a huge number of tests. Here we demonstrate an alternative efficient multi-subgroup screening method using multidimensional hierarchical ...[more]

PMID: 30202041

Similar Datasets

Project description:Motivation:Genome-wide association studies (GWASs) have identified many genetic loci associated with complex traits. A substantial fraction of these identified loci is associated with multiple traits-a phenomena known as pleiotropy. Identification of pleiotropic associations can help characterize the genetic relationship among complex traits and can facilitate our understanding of disease etiology. Effective pleiotropic association mapping requires the development of statistical methods that can jointly model multiple traits with genome-wide single nucleic polymorphisms (SNPs) together. Results:We develop a joint modeling method, which we refer to as the integrative MApping of Pleiotropic association (iMAP). iMAP models summary statistics from GWASs, uses a multivariate Gaussian distribution to account for phenotypic correlation, simultaneously infers genome-wide SNP association pattern using mixture modeling and has the potential to reveal causal relationship between traits. Importantly, iMAP integrates a large number of SNP functional annotations to substantially improve association mapping power, and, with a sparsity-inducing penalty, is capable of selecting informative annotations from a large, potentially non-informative set. To enable scalable inference of iMAP to association studies with hundreds of thousands of individuals and millions of SNPs, we develop an efficient expectation maximization algorithm based on an approximate penalized regression algorithm. With simulations and comparisons to existing methods, we illustrate the benefits of iMAP in terms of both high association mapping power and accurate estimation of genome-wide SNP association patterns. Finally, we apply iMAP to perform a joint analysis of 48 traits from 31 GWAS consortia together with 40 tissue-specific SNP annotations generated from the Roadmap Project. Availability and implementation:iMAP is freely available at http://www.xzlab.org/software.html. Supplementary information:Supplementary data are available at Bioinformatics online.

Project description:We propose a new computational method for exploring chromatin structural organization based on Markov State Modelling of Hi-C data represented as an interaction network between genomic loci. A Markov process describes the random walk of a traveling probe in the corresponding energy landscape, mimicking the motion of a biomolecule involved in chromatin function. By studying the metastability of the associated Markov State Model upon annealing, the hierarchical structure of individual chromosomes is observed, and corresponding set of structural partitions is identified at each level of hierarchy. Then, the notion of effective interaction between partitions is derived, delineating the overall topology and architecture of chromosomes. Mapping epigenetic data on the graphs of intra-chromosomal effective interactions helps in understanding how chromosome organization facilitates its function. A sketch of whole-genome interactions obtained from the analysis of 539 partitions from all 23 chromosomes, complemented by distributions of gene expression regulators and epigenetic factors, sheds light on the structure-function relationships in chromatin, delineating chromosomal territories, as well as structural partitions analogous to topologically associating domains and active / passive epigenomic compartments. In addition to the overall genome architecture shown by effective interactions, the affinity between partitions of different chromosomes was analyzed as an indicator of the degree of association between partitions in functionally relevant genomic interactions. The overall static picture of whole-genome interactions obtained with the method presented in this work provides a foundation for chromatin structural reconstruction, for the modelling of chromatin dynamics, and for exploring the regulation of genome function. The algorithms used in this study are implemented in a freely available Python package ChromaWalker (https://bitbucket.org/ZhenWahTan/chromawalker).

Project description:BackgroundRecently mixed linear models are used to address the issue of "missing" heritability in traditional Genome-wide association studies (GWAS). The models assume that all single-nucleotide polymorphisms (SNPs) are associated with the phenotypes of interest. However, it is more common that only a small proportion of SNPs have significant effects on the phenotypes, while most SNPs have no or very small effects. To incorporate this feature, we propose an efficient Hierarchical Bayesian Model (HBM) that extends the existing mixed models to enforce automatic selection of significant SNPs. The HBM models the SNP effects using a mixture distribution of a point mass at zero and a normal distribution, where the point mass corresponds to those non-associative SNPs.ResultsWe estimate the HBM using Gibbs sampling. The estimation performance of our method is first demonstrated through two simulation studies. We make the simulation setups realistic by using parameters fitted on the Framingham Heart Study (FHS) data. The simulation studies show that our method can accurately estimate the proportion of SNPs associated with the simulated phenotype and identify these SNPs, as well as adapt to certain model mis-specification than the standard mixed models. In addition, we analyze data from the FHS and the Health and Retirement Study (HRS) to study the association between Body Mass Index (BMI) and SNPs on Chromosome 16, and replicate the identified genetic associations. The analysis of the FHS data identifies 0.3% SNPs on Chromosome 16 that affect BMI, including rs9939609 and rs9939973 on the FTO gene. These two SNPs are in strong linkage disequilibrium with rs1558902 (Rsq =0.901 for rs9939609 and Rsq =0.905 for rs9939973), which has been reported to be linked with obesity in previous GWAS. We then replicate the findings using the HRS data: the analysis finds 0.4% of SNPs associated with BMI on Chromosome 16. Furthermore, around 25% of the genes that are identified to be associated with BMI are common between the two studies.ConclusionsThe results demonstrate that the HBM and the associated estimation algorithm offer a powerful tool for identifying significant genetic associations with phenotypes of interest, among a large number of SNPs that are common in modern genetics studies.

Dataset Information

Exploring predictive biomarkers from clinical genome-wide association studies via multidimensional hierarchical mixture models.

Publications

Exploring predictive biomarkers from clinical genome-wide association studies via multidimensional hierarchical mixture models.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure