IsoDOT Detects Differential RNA-isoform Expression/Usage with respect to a Categorical or Continuous Covariate with High Sensitivity and Specificity.
Ontology highlight
ABSTRACT:
SUBMITTER: Sun W
PROVIDER: S-EPMC4662594 | biostudies-literature | 2015
REPOSITORIES: biostudies-literature
Publications
IsoDOT Detects Differential RNA-isoform Expression/Usage with respect to a Categorical or Continuous Covariate with High Sensitivity and Specificity.
Sun Wei W Liu Yufeng Y Crowley James J JJ Chen Ting-Hued TH Zhou Hua H Chu Haitao H Huang Shunping S Kuan Pei-Fen PF Li Yuan Y Miller Darla R DR Shaw Ginger D GD Wu Yichao Y Zhabotynsky Vasyl V McMillan Leonard L Zou Fei F Sullivan Patrick F PF de Villena Fernando Pardo-Manuel FP
Journal of the American Statistical Association 20151107 511
We have developed a statistical method named IsoDOT to assess differential isoform expression (DIE) and differential isoform usage (DIU) using RNA-seq data. Here isoform usage refers to relative isoform expression given the total expression of the corresponding gene. IsoDOT performs two tasks that cannot be accomplished by existing methods: to test DIE/DIU with respect to a continuous covariate, and to test DIE/DIU for one case versus one control. The latter task is not an uncommon situation in ...[more]
Similar Datasets
Project description:Multiple imputation (MI) has become popular for analyses with missing data in medical research. The standard implementation of MI is based on the assumption of data being missing at random (MAR). However, for missing data generated by missing not at random mechanisms, MI performed assuming MAR might not be satisfactory. For an incomplete variable in a given data set, its corresponding population marginal distribution might also be available in an external data source. We show how this information can be readily utilised in the imputation model to calibrate inference to the population by incorporating an appropriately calculated offset termed the "calibrated-? adjustment." We describe the derivation of this offset from the population distribution of the incomplete variable and show how, in applications, it can be used to closely (and often exactly) match the post-imputation distribution to the population level. Through analytic and simulation studies, we show that our proposed calibrated-? adjustment MI method can give the same inference as standard MI when data are MAR, and can produce more accurate inference under two general missing not at random missingness mechanisms. The method is used to impute missing ethnicity data in a type 2 diabetes prevalence case study using UK primary care electronic health records, where it results in scientifically relevant changes in inference for non-White ethnic groups compared with standard MI. Calibrated-? adjustment MI represents a pragmatic approach for utilising available population-level information in a sensitivity analysis to explore potential departures from the MAR assumption.
| S-EPMC6492126 | biostudies-literature
Project description:BACKGROUND: Zebrafish embryos are transcriptionally silent until activation of the zygotic genome during the 10th cell cycle. Onset of transcription is followed by cellular and morphological changes involving cell speciation and gastrulation. Previous genome-wide surveys of transcriptional changes only assessed gene expression levels; however, recent studies have shown the necessity to map isoform-specific transcriptional changes. Here, we perform isoform discovery and quantification on transcriptome sequences from before and after zebrafish zygotic genome activation (ZGA). RESULTS: We identify novel isoforms and isoform switches during ZGA for genes related to cell adhesion, pluripotency and DNA methylation. Isoform switching events include alternative splicing and changes in transcriptional start sites and in 3' untranslated regions. New isoforms are identified even for well-characterized genes such as pou5f1, sall4 and dnmt1. Genes involved in cell-cell interactions such as f11r and magi1 display isoform switches with alterations of coding sequences. We also detect over 1000 transcripts that acquire a longer 3' terminal exon when transcribed by the zygote compared to their maternal transcript counterparts. ChIP-sequencing data mapped onto skipped exon events reveal a correlation between histone H3K36 trimethylation peaks and skipped exons, suggesting epigenetic marks being part of alternative splicing regulation. CONCLUSIONS: The novel isoforms and isoform switches reported here include regulators of transcriptional, cellular and morphological changes taking place around ZGA. Our data display an array of isoform-related functional changes and represent a valuable resource complementary to existing early embryo transcriptomes.
| S-EPMC3747860 | biostudies-literature
Project description:BACKGROUND:Most genes in mammals generate several transcript isoforms that differ in stability and translational efficiency through alternative splicing. Such alternative splicing can be tissue- and developmental stage-specific, and such specificity is sometimes associated with disease. Thus, detecting differential isoform usage for a gene between tissues or cell lines/types (differences in the fraction of total expression of a gene represented by the expression of each of its isoforms) is potentially important for cell and developmental biology. RESULTS:We present a new method IUTA that is designed to test each gene in the genome for differential isoform usage between two groups of samples. IUTA also estimates isoform usage for each gene in each sample as well as averaged across samples within each group. IUTA is the first method to formulate the testing problem as testing for equal means of two probability distributions under the Aitchison geometry, which is widely recognized as the most appropriate geometry for compositional data (vectors that contain the relative amount of each component comprising the whole). Evaluation using simulated data showed that IUTA was able to provide test results for many more genes than was Cuffdiff2 (version 2.2.0, released in Mar. 2014), and IUTA performed better than Cuffdiff2 for the limited number of genes that Cuffdiff2 did analyze. When applied to actual mouse RNA-Seq datasets from six tissues, IUTA identified 2,073 significant genes with clear patterns of differential isoform usage between a pair of tissues. IUTA is implemented as an R package and is available at http://www.niehs.nih.gov/research/resources/software/biostatistics/iuta/index.cfm. CONCLUSIONS:Both simulation and real-data results suggest that IUTA accurately detects differential isoform usage. We believe that our analysis of RNA-seq data from six mouse tissues represents the first comprehensive characterization of isoform usage in these tissues. IUTA will be a valuable resource for those who study the roles of alternative transcripts in cell development and disease.
| S-EPMC4195885 | biostudies-literature
Project description:In population pharmacokinetic analyses, missing categorical data are often encountered. We evaluated several methods of performing covariate analyses with partially missing categorical covariate data. Missing data methods consisted of discarding data (DROP), additional effect parameter for the group with missing data (EXTRA), and mixture methods in which the mixing probability was fixed to the observed fraction of categories (MIX(obs)), based on the likelihood of the concentration data (MIX(conc)), or combined likelihood of observed covariate data and concentration data (MIX(joint)). Simulations were implemented to study bias and imprecision of the methods in datasets with equal-sized and unbalanced category ratios for a binary covariate as well as datasets with non-random missingness (MNAR). Additionally, the performance and feasibility of implementation was assessed in two real datasets. At either low (10%) or high (50%) levels of missingness, all methods performed similarly well. Performance was similar for situations with unbalanced datasets (3:1 covariate distribution) and balanced datasets. In the MNAR scenario, the MIX methods showed a higher bias in the estimation of CL and covariate effect than EXTRA. All methods could be applied to real datasets, except DROP. All methods perform similarly at the studied levels of missingness, but the DROP and EXTRA methods provided less bias than the mixture methods in the case of MNAR. However, EXTRA was associated with inflated type I error rates of covariate selection, while DROP handled data inefficiently.
| S-EPMC3385822 | biostudies-other
Project description:Neuromyelitis Optica (NMO) is an autoimmune demyelinating disease, characterized by the presence of autoantibody (NMO-IgG) to Aquaporin-4 (AQP4). NMO-IgG identification supports NMO diagnosis and several diagnostic tests have been developed, but their sensitivity is too variable, and some assay show low sensitivity. This impairs correct diagnosis of NMO. By cell based assay (CBA) we here evaluate the efficacy of different strategies to express AQP4 in mammalian cells in terms of: a) AQP4 translation initiation signals; b) AQP4 isoforms (M1 and M23) and fluorescent tag position; c) NMO serum concentration and AQP4 degradation. Our results demonstrate that when using AQP4-M1, the nucleotide in position -3 of the AUG greatly affects the AQP4-M1/M23 protein ratio, NMO-IgG binding, and consequently test sensitivity. Test sensitivity was highest with M23 expressing cells (97.5%) and only 27.5% with AQP4-M1. The fluorescent tag added to the N-terminus of AQP4-M23 considerably affected the NMO-IgG binding, and test sensitivity, due to disruption of AQP4 suprastructures. Furthermore, sera used at high concentration resulted in AQP4 degradation which affected test sensitivity. To further evaluate the reliability of the M23 based CBA test, samples of one NMO patient collected during about 2 years clinical follow-up were tested. The results of serum titer correlated with disease activity and treatment response. In conclusion, we provide a molecular explanation for the contrasting CBA test data reported and suggest the use of M23 with a C-terminus fluorescent tag as the proper test for NMO diagnosis.
| S-EPMC3829826 | biostudies-literature
Project description:Data generated from a system of interest typically consists of measurements on many covariate features and possibly multiple response features across all subjects in a designated ensemble. Such data is naturally represented by one response-matrix against one covariate-matrix. A matrix lattice is an advantageous platform for simultaneously accommodating heterogeneous data types: continuous, discrete and categorical, and exploring hidden dependency among/between features and subjects. After each feature being individually renormalized with respect to its own histogram, the categorical version of mutual conditional entropy is evaluated for all pairs of response and covariate features according to the combinatorial information theory. Then, by applying Data Could Geometry (DCG) algorithmic computations on such a mutual conditional entropy matrix, multiple synergistic feature-groups are partitioned. Distinct synergistic feature-groups embrace distinct structures of dependency. The explicit details of dependency among members of synergistic features are seen through mutliscale compositions of blocks computed by a computing paradigm called Data Mechanics. We then propose a categorical pattern matching approach to establish a directed associative linkage: from the patterned response dependency to serial structured covariate dependency. The graphic display of such a directed associative linkage is termed an information flow and the degrees of association are evaluated via tree-to-tree mutual conditional entropy. This new universal way of discovering system knowledge is illustrated through five data sets. In each case, the emergent visible heterogeneity is an organization of discovered knowledge.
| S-EPMC6006982 | biostudies-literature
Project description:RNA-seq has been a boon to the quantitative analysis of transcriptomes. A notable application is the detection of changes in transcript usage between experimental conditions. For example, discovery of pathological alternative splicing may allow the development of new treatments or better management of patients. From an analysis perspective, there are several ways to approach RNA-seq data to unravel differential transcript usage, such as annotation-based exon-level counting, differential analysis of the percentage spliced in, or quantitative analysis of assembled transcripts. The goal of this research is to compare and contrast current state-of-the-art methods, and to suggest improvements to commonly used work flows.We assess the performance of representative work flows using synthetic data and explore the effect of using non-standard counting bin definitions as input to DEXSeq, a state-of-the-art inference engine. Although the canonical counting provided the best results overall, several non-canonical approaches were as good or better in specific aspects and most counting approaches outperformed the evaluated event- and assembly-based methods. We show that an incomplete annotation catalog can have a detrimental effect on the ability to detect differential transcript usage in transcriptomes with few isoforms per gene and that isoform-level prefiltering can considerably improve false discovery rate control.Count-based methods generally perform well in the detection of differential transcript usage. Controlling the false discovery rate at the imposed threshold is difficult, particularly in complex organisms, but can be improved by prefiltering the annotation catalog.
| S-EPMC4729156 | biostudies-literature
Project description:High-throughput sequencing technologies, including RNA-seq, have made it possible to move beyond gene expression analysis to study transcriptional events including alternative splicing and gene fusions. Furthermore, recent studies in cancer have suggested the importance of identifying transcriptionally altered loci as biomarkers for improved prognosis and therapy. While many statistical methods have been proposed for identifying novel transcriptional events with RNA-seq, nearly all rely on contrasting known classes of samples, such as tumor and normal. Few tools exist for the unsupervised discovery of such events without class labels. In this paper, we present SigFuge for identifying genomic loci exhibiting differential transcription patterns across many RNA-seq samples. SigFuge combines clustering with hypothesis testing to identify genes exhibiting alternative splicing, or differences in isoform expression. We apply SigFuge to RNA-seq cohorts of 177 lung and 279 head and neck squamous cell carcinoma samples from the Cancer Genome Atlas, and identify several cases of differential isoform usage including CDKN2A, a tumor suppressor gene known to be inactivated in a majority of lung squamous cell tumors. By not restricting attention to known sample stratifications, SigFuge offers a novel approach to unsupervised screening of genetic loci across RNA-seq cohorts. SigFuge is available as an R package through Bioconductor.
| S-EPMC4132703 | biostudies-literature
Project description:Zebrafish embryos are transcriptional silent until activation of the zygotic genome during the 10th cell cycle. Onset of transcription is followed by cellular and morphological changes involving cell speciation and gastrulation. Previous genome-wide surveys of transcriptional changes only assessed gene expression levels; however, recent studies have shown the necessity to map isoform-specific transcriptional changes. Here we perform isoform discovery and quantification on transcriptome sequences from before and after zebrafish zygotic genome activation (ZGA). We identify novel isoforms and isoform switches during ZGA for genes related to cell adhesion, pluripotency and DNA methylation. Isoform switching events include alternative splicing and changes in transcriptional start sites and in 3’ untranslated regions. New isoforms are identified even for well-characterized genes such as pou5f1, sall4 and dnmt1. Genes involved in cell-cell interactions such as f11r and magi1 display isoform switches with alterations of coding sequences. We also detect over 1000 transcripts that acquire a longer 3’ terminal exon when transcribed zygotically relative to the maternal transcript counterparts. ChIP-seq data mapped onto skipped exon events reveals a correlation between histone H3K36trimethylation peaks and the skipped exons, suggesting epigenetic marks being part of alternative splicing regulation. The novel isoforms and isoform switches reported here include regulators of the transcriptional, cellular and morphological changes taking place around ZGA. Our data display an array of isoform-related functional changes and represent a valuable resource complementary to existing early embryo transcriptomes. Examination H3K36me3 in zebrafish whole embryos at the Post-MBT stage
2012-11-01 | E-GEOD-40629 | biostudies-arrayexpress
Project description:Zebrafish embryos are transcriptional silent until activation of the zygotic genome during the 10th cell cycle. Onset of transcription is followed by cellular and morphological changes involving cell speciation and gastrulation. Previous genome-wide surveys of transcriptional changes only assessed gene expression levels; however, recent studies have shown the necessity to map isoform-specific transcriptional changes. Here we perform isoform discovery and quantification on transcriptome sequences from before and after zebrafish zygotic genome activation (ZGA). We identify novel isoforms and isoform switches during ZGA for genes related to cell adhesion, pluripotency and DNA methylation. Isoform switching events include alternative splicing and changes in transcriptional start sites and in 3’ untranslated regions. New isoforms are identified even for well-characterized genes such as pou5f1, sall4 and dnmt1. Genes involved in cell-cell interactions such as f11r and magi1 display isoform switches with alterations of coding sequences. We also detect over 1000 transcripts that acquire a longer 3’ terminal exon when transcribed zygotically relative to the maternal transcript counterparts. ChIP-seq data mapped onto skipped exon events reveals a correlation between histone H3K36trimethylation peaks and the skipped exons, suggesting epigenetic marks being part of alternative splicing regulation. The novel isoforms and isoform switches reported here include regulators of the transcriptional, cellular and morphological changes taking place around ZGA. Our data display an array of isoform-related functional changes and represent a valuable resource complementary to existing early embryo transcriptomes.
2012-11-01 | GSE40629 | GEO