Dataset Information

VIPR: A probabilistic algorithm for analysis of microbial detection microarrays

ABSTRACT: The development of rapid and sensitive assays capable of detecting a wide range of infectious agents is critical for the effective diagnosis of diseases that have multiple etiologies. In recent years, many microarray-based diagnostics have been developed to identify viruses present in clinical specimens in a highly parallel fashion. Unfortunately, the rate of development of algorithms to interpret data generated from such platforms has not been commensurate. In particular, none of the existing interpretive algorithms is capable of utilizing empirical training data in a Bayesian framework. We have developed an interpretive algorithm, VIPR (Viral Identification using a PRobabilistic algorithm), to capitalize on our ability to generate positive control data for analysis of microbial diagnostic arrays. To illustrate this approach, we have focused on the analysis of viruses that cause hemorrhagic fever (HF). To assess the efficacy of VIPR, we hybridized 33 viruses to 100 microarrays and applied our algorithm to this dataset. A microarray composed of nearly 15,000 oligonucleotides was designed using a custom viral taxonomy-based strategy. The performance of VIPR was assessed by performing a leave-one-out cross validation. VIPR was able to identity the infecting virus with an accuracy of 94%. VIPR outperformed previously described algorithms for the set of HF viruses tested. Bayesian interpretative algorithms such as VIPR should be considered for diagnostic microarray applications.

ORGANISM(S): Peribunyaviridae Mammarenavirus machupoense Ngari virus Yellow fever virus Omsk hemorrhagic fever virus unidentified Mammarenavirus juninense dengue virus type 3 Mammarenavirus amapariense La Crosse virus dengue virus type 4 Rift Valley fever virus Orthonairovirus haemorrhagiae dengue virus type 2 Mammarenavirus choriomeningitidis Kyasanur Forest disease virus Mammarenavirus lassaense Toscana virus Mammarenavirus ippyense Orthohantavirus seoulense Sudan ebolavirus Reston ebolavirus Mammarenavirus guanaritoense Mammarenavirus mopeiaense Mammarenavirus tacaribeense California encephalitis virus Zaire ebolavirus Mammarenavirus brazilense Orthohantavirus puumalaense Mammarenavirus praomyidis Rocio virus Filoviridae Ebola virus - Gabon (1994-1997) Orthohantavirus hantanense Lake Victoria marburgvirus - Angola2005 Flaviviridae Arenaviridae dengue virus type 1

PROVIDER: GSE21407 | GEO | 2011/04/19

SECONDARY ACCESSION(S): PRJNA126409

REPOSITORIES: GEO

ACCESS DATA

Similar Datasets

Project description:Single-cell analysis of the transcriptome deepens our understanding of an individual cell's contribution to its microenvironment. Using single-cell analysis to study complex biological processes requires state-of-the-art computational tools. Assessing similarity is highly important for bioinformatics algorithms in order to determine correlations between biological information. Similarity can appear by chance, particularly for low expressed entities. This is especially relevant in single cell RNA-seq (scRNA-seq) because the read counts obtained are lower compared to bulk RNA-sequencing and therefore classic bioinformatics tools are insufficient to obtain reproducible results. Recently, a Bayesian correlation scheme, that assigns low correlation values to correlations coming from low expressed genes, has been proposed to assess similarity for bulk RNA-seq and miRNA. This Bayesian method uses a prior distribution before using empirical evidence. Our goal was to extend the properties of this Bayesian correlation scheme to scRNA-seq data. We assessed 3 ways to compute similarity. First, we computed the similarity of each pair of genes over all cells. Second, we identified specific cell populations and computed the correlation in those specific cells. Third, we computed the similarity of each pair of genes over all clusters, by including the total mRNA expression in those cells. To study the effect of the number of cells on the method, we did not rely on simulated data, we generated 4 scRNA-seq mouse liver cell libraries with a varying number of input cells. Results: We show that Bayesian correlations are more reproducible than Pearson correlations in all the scenarios studied. Compared to Pearson correlations, Bayesian correlations have a smaller dependence on the number of input cells. We demonstrate that the Bayesian correlation algorithm assigns high similarity values to genes with a biological relevance in a specific population. Significance: Our results demonstrate that Bayesian correlation is a robust similarity measure for scRNA-seq datasets. The Bayesian method allows researchers to study similarity between pairs of genes without discarding low expressed entities and to minimize biasing the results by fake correlations. Taken together, using our method of Bayesian correlation the reproducibility of scRNA-seq experiments is increased significantly.

Project description:Gene transcription mediated by RNA polymerase II (pol-II) is a key step in gene expression. The dynamics of pol-II moving along the transcribed region influences the rate and timing of gene expression. In this work we present a probabilistic model of transcription dynamics which is fitted to pol-II occupancy time course data measured using ChIP-Seq. The model can be used to estimate transcription speed and to infer the temporal pol-II activity profile at the gene promoter. Model parameters are determined using either maximum likelihood estimation or via Bayesian inference using Markov chain Monte Carlo sampling. The Bayesian approach provides confidence intervals for parameter estimates and allows the use of priors that capture domain knowledge, e.g. the expected range of transcription speeds, based on previous experiments. The model describes the movement of pol-II down the gene body and can be used to identify the time of induction for transcriptionally engaged genes. By clustering the inferred promoter activity time profiles, we are able to determine which genes respond quickly to stimuli and group genes that share activity profiles and may therefore be co-regulated. We apply our methodology to biological data obtained using ChIP-seq to measure pol-II occupancy genome-wide when MCF-7 human breast cancer cells are treated with estradiol (E2). The transcription speeds we obtain agree with those obtained previously for smaller numbers of genes with the advantage that our approach can be applied genome- wide. We validate the biological significance of the pol-II promoter activity clusters by investigating cluster-specific transcription factor binding patterns and determining canonical pathway enrichment. MCF-7 cells were mock treated or with 10nM 17b-E2 to seven time points (5', 10', 20', 40', 80', 160', 320') Genome-wide identification of RNA polymerase II (RNAPII) occupancy at different time points following E2 induction of MCF-7 cells

Dataset Information

VIPR: A probabilistic algorithm for analysis of microbial detection microarrays

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure