Genome-wide identification of directed gene networks using large-scale population genomics data.
Ontology highlight
ABSTRACT:
SUBMITTER: Luijk R
PROVIDER: S-EPMC6079029 | biostudies-literature | 2018 Aug
REPOSITORIES: biostudies-literature
Publications
Genome-wide identification of directed gene networks using large-scale population genomics data.
Luijk René R Dekkers Koen F KF van Iterson Maarten M Arindrarto Wibowo W Claringbould Annique A Hop Paul P Boomsma Dorret I DI van Duijn Cornelia M CM van Greevenbroek Marleen M J MMJ Veldink Jan H JH Wijmenga Cisca C Franke Lude L 't Hoen Peter A C PAC Jansen Rick R van Meurs Joyce J Mei Hailiang H Slagboom P Eline PE Heijmans Bastiaan T BT van Zwet Erik W EW
Nature communications 20180806 1
Identification of causal drivers behind regulatory gene networks is crucial in understanding gene function. Here, we develop a method for the large-scale inference of gene-gene interactions in observational population genomics data that are both directed (using local genetic instruments as causal anchors, akin to Mendelian Randomization) and specific (by controlling for linkage disequilibrium and pleiotropy). Analysis of genotype and whole-blood RNA-sequencing data from 3072 individuals identifi ...[more]
Similar Datasets
Project description:Principal component analysis (PCA) is routinely used to analyze genome-wide single-nucleotide polymorphism (SNP) data, for detecting population structure and potential outliers. However, the size of SNP datasets has increased immensely in recent years and PCA of large datasets has become a time consuming task. We have developed flashpca, a highly efficient PCA implementation based on randomized algorithms, which delivers identical accuracy in extracting the top principal components compared with existing tools, in substantially less time. We demonstrate the utility of flashpca on both HapMap3 and on a large Immunochip dataset. For the latter, flashpca performed PCA of 15,000 individuals up to 125 times faster than existing tools, with identical results, and PCA of 150,000 individuals using flashpca completed in 4 hours. The increasing size of SNP datasets will make tools such as flashpca essential as traditional approaches will not adequately scale. This approach will also help to scale other applications that leverage PCA or eigen-decomposition to substantially larger datasets.
| S-EPMC3981753 | biostudies-literature
Project description:Genome-wide association studies (GWAS) have successfully identified about 70 genomic loci associated with breast cancer. Owing to the complexity of linkage disequilibrium and environmental exposures in different populations, it is essential to perform regional GWAS for better risk prediction. This study aimed to investigate the genetic architecture and to assess common genetic risk model of breast cancer with 6,669 breast cancer patients and 21,930 female controls in the Japanese population. This GWAS identified 11 genomic loci that surpass genome-wide significance threshold of P < 5.0 × 10-8 with nine previously reported loci and two novel loci that include rs9862599 on 3q13.11 (ALCAM) and rs75286142 on 21q22.12 (CLIC6-RUNX1). Validation study was carried out with 981 breast cancer cases and 1,394 controls from the Aichi Cancer Center. Pathway analyses of GWAS signals identified association of dopamine receptor medicated signaling and protein amino acid deacetylation with breast cancer. Weighted genetic risk score showed that individuals who were categorized in the highest risk group are approximately 3.7 times more likely to develop breast cancer compared to individuals in the lowest risk group. This well-powered GWAS is a representative study to identify SNPs that are associated with breast cancer in the Japanese population.
| S-EPMC6874604 | biostudies-literature
Project description:In the era of CRISPR gene editing and genetic screening, there is an increasing demand for quick and reliable nucleic acid extraction pipelines for rapid genotyping of large and diverse sample sets. Despite continuous improvements of current workflows, the handling-time and material costs per sample remain major limiting factors. Here we present a robust method for low-cost DIY-pipet tips addressing these needs; i.e. using a cellulose filter disc inserted into a regular pipet tip. These filter-in-tips allow for a rapid, stand-alone four-step genotyping workflow by simply binding the DNA contained in the primary lysate to the cellulose filter, washing it in water and eluting it directly into the buffer for the downstream application (e.g. PCR). This drastically cuts down processing time to maximum 30 seconds per sample, with the potential for parallelizing and automation. We show the ease and sensitivity of our procedure by genotyping genetically modified medaka (Oryzias latipes) and zebrafish (Danio rerio) embryos (targeted by CRISPR/Cas9 knock-out and knock-in) in a 96-well plate format. The robust isolation and detection of multiple alleles of various abundancies in a mosaic genetic background allows phenotype-genotype correlation already in the injected generation, demonstrating the reliability and sensitivity of the protocol. Our method is applicable across kingdoms to samples ranging from cells to tissues i. e. plant seedlings, adult flies, mouse cell culture and tissue as well as adult fish fin-clips.
| S-EPMC6400387 | biostudies-literature
Project description:Single-locus analysis is often used to analyze genome-wide association (GWA) data, but such analysis is subject to severe multiple comparisons adjustment. Multivariate logistic regression is proposed to fit a multi-locus model for case-control data. However, when the sample size is much smaller than the number of single-nucleotide polymorphisms (SNPs) or when correlation among SNPs is high, traditional multivariate logistic regression breaks down. To accommodate the scale of data from a GWA while controlling for collinearity and overfitting in a high dimensional predictor space, we propose a variable selection procedure using Bayesian logistic regression. We explored a connection between Bayesian regression with certain priors and L1 and L2 penalized logistic regression. After analyzing large number of SNPs simultaneously in a Bayesian regression, we selected important SNPs for further consideration. With much fewer SNPs of interest, problems of multiple comparisons and collinearity are less severe. We conducted simulation studies to examine probability of correctly selecting disease contributing SNPs and applied developed methods to analyze Genetic Analysis Workshop 16 North American Rheumatoid Arthritis Consortium data.
| S-EPMC2795912 | biostudies-literature
Project description:More and more large-scale imaging genetic studies are being widely conducted to collect a rich set of imaging, genetic, and clinical data to detect putative genes for complexly inherited neuropsychiatric and neurodegenerative disorders. Several major big-data challenges arise from testing genome-wide (NC>12 million known variants) associations with signals at millions of locations (NV~10(6)) in the brain from thousands of subjects (n~10(3)). The aim of this paper is to develop a Fast Voxelwise Genome Wide Association analysiS (FVGWAS) framework to efficiently carry out whole-genome analyses of whole-brain data. FVGWAS consists of three components including a heteroscedastic linear model, a global sure independence screening (GSIS) procedure, and a detection procedure based on wild bootstrap methods. Specifically, for standard linear association, the computational complexity is O (nNVNC) for voxelwise genome wide association analysis (VGWAS) method compared with O ((NC+NV)n(2)) for FVGWAS. Simulation studies show that FVGWAS is an efficient method of searching sparse signals in an extremely large search space, while controlling for the family-wise error rate. Finally, we have successfully applied FVGWAS to a large-scale imaging genetic data analysis of ADNI data with 708 subjects, 193,275voxels in RAVENS maps, and 501,584 SNPs, and the total processing time was 203,645s for a single CPU. Our FVGWAS may be a valuable statistical toolbox for large-scale imaging genetic analysis as the field is rapidly advancing with ultra-high-resolution imaging and whole-genome sequencing.
| S-EPMC4554832 | biostudies-literature
Project description:MotivationMuch of a cell's regulatory response to changing environments occurs at the transcriptional level. Particularly in higher organisms, transcription factors (TFs), microRNAs and epigenetic modifications can combine to form a complex regulatory network. Part of this system can be modeled as a collection of regulatory modules: co-regulated genes, the conditions under which they are co-regulated and sequence-level regulatory motifs.ResultsWe present the Combinatorial Algorithm for Expression and Sequence-based Cluster Extraction (COALESCE) system for regulatory module prediction. The algorithm is efficient enough to discover expression biclusters and putative regulatory motifs in metazoan genomes (>20,000 genes) and very large microarray compendia (>10,000 conditions). Using Bayesian data integration, it can also include diverse supporting data types such as evolutionary conservation or nucleosome placement. We validate its performance using a functional evaluation of co-clustered genes, known yeast and Escherichea coli TF targets, synthetic data and various metazoan data compendia. In all cases, COALESCE performs as well or better than current biclustering and motif prediction tools, with high accuracy in functional and TF/target assignments and zero false positives on synthetic data. COALESCE provides an efficient and flexible platform within which large, diverse data collections can be integrated to predict metazoan regulatory networks.AvailabilitySource code (C++) is available at http://function.princeton.edu/sleipnir, and supporting data and a web interface are provided at http://function.princeton.edu/coalesce.Contactogt@cs.princeton.edu; hcoller@princeton.edu.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC2788929 | biostudies-literature
Project description:Technological developments in large-scale biological experiments, coupled with bioinformatics tools, have opened the doors to computational approaches for the global analysis of whole genomes. This has provided the opportunity to look at genes within their context in the cell. The integration of vast amounts of data generated by these technologies provides a strategy for identifying potential drug targets within microbial pathogens, the causative agents of infectious diseases. As proteins are druggable targets, functional interaction networks between proteins are used to identify proteins essential to the survival, growth, and virulence of these microbial pathogens. Here we have integrated functional genomics data to generate functional interaction networks between Mycobacterium tuberculosis proteins and carried out computational analyses to dissect the functional interaction network produced for identifying drug targets using network topological properties. This study has provided the opportunity to expand the range of potential drug targets and to move towards optimal target-based strategies.
| S-EPMC3235424 | biostudies-literature
Project description:MotivationGenome-wide datasets produced for association studies have dramatically increased in size over the past few years, with modern datasets commonly including millions of variants measured in dozens of thousands of individuals. This increase in data size is a major challenge severely slowing down genomic analyses, leading to some software becoming obsolete and researchers having limited access to diverse analysis tools.ResultsHere we present two R packages, bigstatsr and bigsnpr, allowing for the analysis of large scale genomic data to be performed within R. To address large data size, the packages use memory-mapping for accessing data matrices stored on disk instead of in RAM. To perform data pre-processing and data analysis, the packages integrate most of the tools that are commonly used, either through transparent system calls to existing software, or through updated or improved implementation of existing methods. In particular, the packages implement fast and accurate computations of principal component analysis and association studies, functions to remove single nucleotide polymorphisms in linkage disequilibrium and algorithms to learn polygenic risk scores on millions of single nucleotide polymorphisms. We illustrate applications of the two R packages by analyzing a case-control genomic dataset for celiac disease, performing an association study and computing polygenic risk scores. Finally, we demonstrate the scalability of the R packages by analyzing a simulated genome-wide dataset including 500?000 individuals and 1 million markers on a single desktop computer.Availability and implementationhttps://privefl.github.io/bigstatsr/ and https://privefl.github.io/bigsnpr/.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC6084588 | biostudies-literature
Project description:With the rapid advancement of high-throughput sequencing technology, it is now possible to identify individual gene families from genomes on a large scale in order to study their functions. WRKY transcription factors are a key class of regulators that regulate plant growth and abiotic stresses. Here, a total of 74 WRKY genes were identified from Dendrobium officinale Kimura et Migo genome. Based on the genome-wide analysis, an in-depth analysis of gene structure and conserved motif was performed. The phylogenetic analysis indicated that DoWRKYs could be classified into three main groups: I, II, and III, with group II divided into five subgroups: II-a, II-b, II-c, II-d, and II-e. The sequence alignment indicated that these WRKY transcriptional factors contained a highly conserved WRKYGQK heptapeptide. The localization analysis of chromosomes showed that WRKY genes were irregularly distributed across several chromosomes of D. officinale. These genes comprised diverse patterns in both number and species, and there were certain distinguishing motifs among subfamilies. Moreover, the phylogenetic tree and chromosomal location results indicated that DoWRKYs may have undergone a widespread genome duplication event. Based on an evaluation of expression profiles, we proposed that DoWRKY5, 54, 57, 21, etc. may be involved in the transcriptional regulation of the JA signaling pathway. These results provide a scientific reference for the study of DoWRKY family genes.
| S-EPMC9868916 | biostudies-literature
Project description:It is well documented that drug responses are related to Absorption, Distribution, Metabolism, and Excretion (ADME) characteristics of individual patients. Several studies have identified genetic variability in pharmacogenes, that are either directly responsible for or are associated with ADME, giving rise to individualized treatments. Our objective was to provide a comprehensive overview of pharmacogenetic variation in the Saudi population. We mined next generation sequencing (NGS) data from 11,889 unrelated Saudi nationals, to determine the presence and frequencies of known functional SNP variants in 8 clinically relevant pharmacogenes (CYP2C9, CYP2C19, CYP3A5, CYP4F2, VKORC1, DPYD, TPMT and NUDT15), recommended by the Clinical Pharmacogenetics Implementation Consortium (CPIC), and collectively identified 82 such star alleles. Functionally significant pharmacogenetic variants were prevalent especially in CYP genes (excluding CYP3A5), with 10-44.4% of variants predicted to be inactive or to have decreased activity. In CYP3A5, inactive alleles (87.5%) were the most common. Only 1.8%, 0.7% and 0.7% of NUDT15, TPMT and DPYD variants respectively, were predicted to affect gene activity. In contrast, VKORC1 was found functionally, to be highly polymorphic with 53.7% of Saudi individuals harboring variants predicted to result in decreased activity and 31.3% having variants leading to increased metabolic activity. Furthermore, among the 8 pharmacogenes studied, we detected six rare variants with an aggregated frequency of 1.1%, that among several other ethnicities, were uniquely found in Saudi population. Similarly, within our cohort, the 8 pharmacogenes yielded forty-six novel variants predicted to be deleterious. Based upon our findings, 99.2% of individuals from the Saudi population carry at least one actionable pharmacogenetic variant.
| S-EPMC8797234 | biostudies-literature