1
|
Muterko A, Kiseleva A, Salina E. A Transcriptome Response of Bread Wheat ( Triticum aestivum L.) to a 5B Chromosome Substitution from Wild Emmer. PLANTS (BASEL, SWITZERLAND) 2024; 13:1514. [PMID: 38891322 PMCID: PMC11174853 DOI: 10.3390/plants13111514] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/26/2024] [Revised: 05/16/2024] [Accepted: 05/29/2024] [Indexed: 06/21/2024]
Abstract
Over the years, alien chromosome substitution has attracted the attention of geneticists and breeders as a rich source of remarkable genetic diversity for improvement in narrowly adapted wheat cultivars. One of the problems encountered along this way is the coadaptation and realization of the genome of common wheat against the background of the introduced genes. Here, using RNA-Seq, we assessed a transcriptome response of hexaploid wheat Triticum aestivum L. (cultivar Chinese Spring) to a 5B chromosome substitution with its homolog from wild emmer (tetraploid wheat T. dicoccoides Koern) and discuss how complete the physiological compensation for this alien chromatin introgression is. The main signature of the transcriptome in the substituted line was a sharp significant drop of activity before the beginning of the photoperiod with a gradual increase up to overexpression in the middle of the night. The differential expression altered almost all biological processes and pathways tested. Because in most cases, the differential expression or its fold change were modest, and this was only a small proportion of the expressed transcriptome, the physiological compensation of the 5B chromosome substitution in common wheat seemed overall satisfactory, albeit not completely. No over- or under-representation of differential gene expression was found in specific chromosomes, implying that local structural changes in the genome can trigger a global transcriptome response.
Collapse
Affiliation(s)
- Alexandr Muterko
- Institute of Cytology and Genetics SB RAS, 10 Akad. Lavrentyeva Avenue, Novosibirsk 630090, Russia
| | | | | |
Collapse
|
2
|
Stansfield JC, Cresswell KG, Dozmorov MG. multiHiCcompare: joint normalization and comparative analysis of complex Hi-C experiments. Bioinformatics 2020; 35:2916-2923. [PMID: 30668639 DOI: 10.1093/bioinformatics/btz048] [Citation(s) in RCA: 41] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2018] [Revised: 12/14/2018] [Accepted: 01/17/2019] [Indexed: 12/17/2022] Open
Abstract
MOTIVATION With the development of chromatin conformation capture technology and its high-throughput derivative Hi-C sequencing, studies of the three-dimensional interactome of the genome that involve multiple Hi-C datasets are becoming available. To account for the technology-driven biases unique to each dataset, there is a distinct need for methods to jointly normalize multiple Hi-C datasets. Previous attempts at removing biases from Hi-C data have made use of techniques which normalize individual Hi-C datasets, or, at best, jointly normalize two datasets. RESULTS Here, we present multiHiCcompare, a cyclic loess regression-based joint normalization technique for removing biases across multiple Hi-C datasets. In contrast to other normalization techniques, it properly handles the Hi-C-specific decay of chromatin interaction frequencies with the increasing distance between interacting regions. multiHiCcompare uses the general linear model framework for comparative analysis of multiple Hi-C datasets, adapted for the Hi-C-specific decay of chromatin interaction frequencies. multiHiCcompare outperforms other methods when detecting a priori known chromatin interaction differences from jointly normalized datasets. Applied to the analysis of auxin-treated versus untreated experiments, and CTCF depletion experiments, multiHiCcompare was able to recover the expected epigenetic and gene expression signatures of loss of chromatin interactions and reveal novel insights. AVAILABILITY AND IMPLEMENTATION multiHiCcompare is freely available on GitHub and as a Bioconductor R package https://bioconductor.org/packages/multiHiCcompare. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- John C Stansfield
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA, USA
| | - Kellen G Cresswell
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA, USA
| | - Mikhail G Dozmorov
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA, USA
| |
Collapse
|
3
|
Xu T, Demmer RT, Li G. Zero-inflated Poisson factor model with application to microbiome read counts. Biometrics 2020; 77:91-101. [PMID: 32277466 DOI: 10.1111/biom.13272] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2019] [Revised: 03/17/2020] [Accepted: 03/30/2020] [Indexed: 11/29/2022]
Abstract
Dimension reduction of high-dimensional microbiome data facilitates subsequent analysis such as regression and clustering. Most existing reduction methods cannot fully accommodate the special features of the data such as count-valued and excessive zero reads. We propose a zero-inflated Poisson factor analysis model in this paper. The model assumes that microbiome read counts follow zero-inflated Poisson distributions with library size as offset and Poisson rates negatively related to the inflated zero occurrences. The latent parameters of the model form a low-rank matrix consisting of interpretable loadings and low-dimensional scores that can be used for further analyses. We develop an efficient and robust expectation-maximization algorithm for parameter estimation. We demonstrate the efficacy of the proposed method using comprehensive simulation studies. The application to the Oral Infections, Glucose Intolerance, and Insulin Resistance Study provides valuable insights into the relation between subgingival microbiome and periodontal disease.
Collapse
Affiliation(s)
- Tianchen Xu
- Department of Biostatistics, Mailman School of Public Health, Columbia University, New York
| | - Ryan T Demmer
- Division of Epidemiology, School of Public Health, University of Minnesota, Minneapolis, Minnesota
| | - Gen Li
- Department of Biostatistics, Mailman School of Public Health, Columbia University, New York
| |
Collapse
|
4
|
Baik B, Yoon S, Nam D. Benchmarking RNA-seq differential expression analysis methods using spike-in and simulation data. PLoS One 2020; 15:e0232271. [PMID: 32353015 PMCID: PMC7192453 DOI: 10.1371/journal.pone.0232271] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2019] [Accepted: 04/10/2020] [Indexed: 12/27/2022] Open
Abstract
Benchmarking RNA-seq differential expression analysis methods using spike-in and simulated RNA-seq data has often yielded inconsistent results. The spike-in data, which were generated from the same bulk RNA sample, only represent technical variability, making the test results less reliable. We compared the performance of 12 differential expression analysis methods for RNA-seq data, including recent variants in widely used software packages, using both RNA spike-in and simulation data for negative binomial (NB) model. Performance of edgeR, DESeq2, and ROTS was particularly different between the two benchmark tests. Then, each method was tested under most extensive simulation conditions especially demonstrating the large impacts of proportion, dispersion, and balance of differentially expressed (DE) genes. DESeq2, a robust version of edgeR (edgeR.rb), voom with TMM normalization (voom.tmm) and sample weights (voom.sw) showed an overall good performance regardless of presence of outliers and proportion of DE genes. The performance of RNA-seq DE gene analysis methods substantially depended on the benchmark used. Based on the simulation results, suitable methods were suggested under various test conditions.
Collapse
Affiliation(s)
- Bukyung Baik
- Department of Biological Sciences, Ulsan National Institute of Science and Technology, Ulsan, Republic of Korea
| | - Sora Yoon
- Department of Biological Sciences, Ulsan National Institute of Science and Technology, Ulsan, Republic of Korea
| | - Dougu Nam
- Department of Biological Sciences, Ulsan National Institute of Science and Technology, Ulsan, Republic of Korea
- Department of Mathematical Sciences, Ulsan National Institute of Science and Technology, Ulsan, Republic of Korea
- * E-mail:
| |
Collapse
|
5
|
Färkkilä A, Gulhan DC, Casado J, Jacobson CA, Nguyen H, Kochupurakkal B, Maliga Z, Yapp C, Chen YA, Schapiro D, Zhou Y, Graham JR, Dezube BJ, Munster P, Santagata S, Garcia E, Rodig S, Lako A, Chowdhury D, Shapiro GI, Matulonis UA, Park PJ, Hautaniemi S, Sorger PK, Swisher EM, D'Andrea AD, Konstantinopoulos PA. Immunogenomic profiling determines responses to combined PARP and PD-1 inhibition in ovarian cancer. Nat Commun 2020; 11:1459. [PMID: 32193378 PMCID: PMC7081234 DOI: 10.1038/s41467-020-15315-8] [Citation(s) in RCA: 168] [Impact Index Per Article: 42.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2019] [Accepted: 02/26/2020] [Indexed: 11/09/2022] Open
Abstract
Combined PARP and immune checkpoint inhibition has yielded encouraging results in ovarian cancer, but predictive biomarkers are lacking. We performed immunogenomic profiling and highly multiplexed single-cell imaging on tumor samples from patients enrolled in a Phase I/II trial of niraparib and pembrolizumab in ovarian cancer (NCT02657889). We identify two determinants of response; mutational signature 3 reflecting defective homologous recombination DNA repair, and positive immune score as a surrogate of interferon-primed exhausted CD8 + T-cells in the tumor microenvironment. Presence of one or both features associates with an improved outcome while concurrent absence yields no responses. Single-cell spatial analysis reveals prominent interactions of exhausted CD8 + T-cells and PD-L1 + macrophages and PD-L1 + tumor cells as mechanistic determinants of response. Furthermore, spatial analysis of two extreme responders shows differential clustering of exhausted CD8 + T-cells with PD-L1 + macrophages in the first, and exhausted CD8 + T-cells with cancer cells harboring genomic PD-L1 and PD-L2 amplification in the second.
Collapse
Affiliation(s)
- Anniina Färkkilä
- Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, MA, 02215, USA.,Harvard Medical School, 25 Shattuck Street, Boston, MA, 02115, USA.,Research Program in Systems Oncology, University of Helsinki, Haartmaninkatu 8, 00014, Helsinki, Finland.,Laboratory of Systems Pharmacology, Harvard Medical School, Boston, 200 Longwood Avenue, MA, 02115, USA
| | - Doga C Gulhan
- Harvard Medical School, 25 Shattuck Street, Boston, MA, 02115, USA
| | - Julia Casado
- Research Program in Systems Oncology, University of Helsinki, Haartmaninkatu 8, 00014, Helsinki, Finland
| | - Connor A Jacobson
- Laboratory of Systems Pharmacology, Harvard Medical School, Boston, 200 Longwood Avenue, MA, 02115, USA
| | - Huy Nguyen
- Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, MA, 02215, USA
| | - Bose Kochupurakkal
- Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, MA, 02215, USA
| | - Zoltan Maliga
- Laboratory of Systems Pharmacology, Harvard Medical School, Boston, 200 Longwood Avenue, MA, 02115, USA
| | - Clarence Yapp
- Laboratory of Systems Pharmacology, Harvard Medical School, Boston, 200 Longwood Avenue, MA, 02115, USA
| | - Yu-An Chen
- Laboratory of Systems Pharmacology, Harvard Medical School, Boston, 200 Longwood Avenue, MA, 02115, USA
| | - Denis Schapiro
- Laboratory of Systems Pharmacology, Harvard Medical School, Boston, 200 Longwood Avenue, MA, 02115, USA
| | - Yinghui Zhou
- TESARO: A GSK company, 1000 Winter Street, Waltham, MA, 02451, USA
| | - Julie R Graham
- TESARO: A GSK company, 1000 Winter Street, Waltham, MA, 02451, USA
| | - Bruce J Dezube
- TESARO: A GSK company, 1000 Winter Street, Waltham, MA, 02451, USA
| | - Pamela Munster
- Helen Diller Family Comprehensive Cancer Center, 1450 3rd Street, San Francisco, CA, 94158, USA
| | - Sandro Santagata
- Brigham and Women's Hospital, Laboratory for Systems Pharmacology, 75 Francis Street, Boston, MA, 02115, USA
| | - Elizabeth Garcia
- Department of Pathology, Brigham and Women's Hospital, Harvard Medical School, 75 Francis Street, Boston, MA, 02115, USA
| | - Scott Rodig
- Department of Pathology, Brigham and Women's Hospital, Harvard Medical School, 75 Francis Street, Boston, MA, 02115, USA
| | - Ana Lako
- Department of Pathology, Brigham and Women's Hospital, Harvard Medical School, 75 Francis Street, Boston, MA, 02115, USA
| | - Dipanjan Chowdhury
- Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, MA, 02215, USA
| | - Geoffrey I Shapiro
- Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, MA, 02215, USA
| | - Ursula A Matulonis
- Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, MA, 02215, USA
| | - Peter J Park
- Department of Biomedical Informatics, Harvard Medical School, 25 Shattuck Street, Boston, MA, 02115, USA
| | - Sampsa Hautaniemi
- Research Program in Systems Oncology, University of Helsinki, Haartmaninkatu 8, 00014, Helsinki, Finland
| | - Peter K Sorger
- Laboratory of Systems Pharmacology, Harvard Medical School, Boston, 200 Longwood Avenue, MA, 02115, USA
| | | | - Alan D D'Andrea
- Dana-Farber Cancer Institute, 450 Brookline Avenue, Boston, MA, 02215, USA.
| | | |
Collapse
|
6
|
Probabilistic cell typing enables fine mapping of closely related cell types in situ. Nat Methods 2019; 17:101-106. [PMID: 31740815 PMCID: PMC6949128 DOI: 10.1038/s41592-019-0631-4] [Citation(s) in RCA: 122] [Impact Index Per Article: 24.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2018] [Accepted: 10/08/2019] [Indexed: 12/24/2022]
Abstract
Understanding the function of a tissue requires knowing the spatial organization of its constituent cell types. In the cerebral cortex, single-cell RNA sequencing (scRNA-seq) has revealed the genome-wide expression patterns that define its many, closely related neuronal types, but cannot reveal their spatial arrangement. Here we introduce probabilistic cell typing by in situ sequencing (pciSeq), an approach that leverages prior scRNA-seq classification to identify cell types using multiplexed in situ RNA detection. We applied this method by mapping the inhibitory neurons of hippocampal area CA1, for which ground truth is available from extensive prior work identifying their laminar organization. Our method identified these closely related classes in a spatial arrangement matching ground truth, and further identified multiple classes of isocortical pyramidal cell in a pattern matching their known organization. This method will allow identifying the spatial organization of fine cell types across the brain and other tissues.
Collapse
|
7
|
Ullah I, Paul S, Hong Z, Wang YG. Significance tests for analyzing gene expression data with small sample sizes. Bioinformatics 2019; 35:3996-4003. [PMID: 30874796 DOI: 10.1093/bioinformatics/btz189] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2018] [Revised: 02/20/2019] [Accepted: 03/13/2019] [Indexed: 11/12/2022] Open
Abstract
MOTIVATION Under two biologically different conditions, we are often interested in identifying differentially expressed genes. It is usually the case that the assumption of equal variances on the two groups is violated for many genes where a large number of them are required to be filtered or ranked. In these cases, exact tests are unavailable and the Welch's approximate test is most reliable one. The Welch's test involves two layers of approximations: approximating the distribution of the statistic by a t-distribution, which in turn depends on approximate degrees of freedom. This study attempts to improve upon Welch's approximate test by avoiding one layer of approximation. RESULTS We introduce a new distribution that generalizes the t-distribution and propose a Monte Carlo based test that uses only one layer of approximation for statistical inferences. Experimental results based on extensive simulation studies show that the Monte Carol based tests enhance the statistical power and performs better than Welch's t-approximation, especially when the equal variance assumption is not met and the sample size of the sample with a larger variance is smaller. We analyzed two gene-expression datasets, namely the childhood acute lymphoblastic leukemia gene-expression dataset with 22 283 genes and Golden Spike dataset produced by a controlled experiment with 13 966 genes. The new test identified additional genes of interest in both datasets. Some of these genes have been proven to play important roles in medical literature. AVAILABILITY AND IMPLEMENTATION R scripts and the R package mcBFtest is available in CRAN and to reproduce all reported results are available at the GitHub repository, https://github.com/iullah1980/MCTcodes. SUPPLEMENTARY INFORMATION Supplementary data is available at Bioinformatics online.
Collapse
Affiliation(s)
- Insha Ullah
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, QLD, Australia
| | - Sudhir Paul
- Department of Mathematics and Statistics, University of Windsor, Windsor, ON, Canada
| | - Zhenjie Hong
- College of Mathematics and Physics, Wenzhou University, Wenzhou, Zhejiang, China
| | - You-Gan Wang
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, QLD, Australia
| |
Collapse
|
8
|
Lo C, Marculescu R. MetaNN: accurate classification of host phenotypes from metagenomic data using neural networks. BMC Bioinformatics 2019; 20:314. [PMID: 31216991 PMCID: PMC6584521 DOI: 10.1186/s12859-019-2833-2] [Citation(s) in RCA: 40] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023] Open
Abstract
BACKGROUND Microbiome profiles in the human body and environment niches have become publicly available due to recent advances in high-throughput sequencing technologies. Indeed, recent studies have already identified different microbiome profiles in healthy and sick individuals for a variety of diseases; this suggests that the microbiome profile can be used as a diagnostic tool in identifying the disease states of an individual. However, the high-dimensional nature of metagenomic data poses a significant challenge to existing machine learning models. Consequently, to enable personalized treatments, an efficient framework that can accurately and robustly differentiate between healthy and sick microbiome profiles is needed. RESULTS In this paper, we propose MetaNN (i.e., classification of host phenotypes from Metagenomic data using Neural Networks), a neural network framework which utilizes a new data augmentation technique to mitigate the effects of data over-fitting. CONCLUSIONS We show that MetaNN outperforms existing state-of-the-art models in terms of classification accuracy for both synthetic and real metagenomic data. These results pave the way towards developing personalized treatments for microbiome related diseases.
Collapse
Affiliation(s)
- Chieh Lo
- Department of Electrical and Computer Engineering, Carnegie Mellon University, 5000 Forbes Ave., Pittsburgh, PA, USA
| | - Radu Marculescu
- Department of Electrical and Computer Engineering, Carnegie Mellon University, 5000 Forbes Ave., Pittsburgh, PA, USA.
| |
Collapse
|
9
|
Dietary Fat, but Not Protein or Carbohydrate, Regulates Energy Intake and Causes Adiposity in Mice. Cell Metab 2018; 28:415-431.e4. [PMID: 30017356 DOI: 10.1016/j.cmet.2018.06.010] [Citation(s) in RCA: 175] [Impact Index Per Article: 29.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/14/2017] [Revised: 02/25/2018] [Accepted: 06/14/2018] [Indexed: 12/30/2022]
Abstract
The impacts of different macronutrients on body weight regulation remain unresolved, with different studies suggesting increased dietary fat, increased carbohydrates (particularly sugars), or reduced protein may all stimulate overconsumption and drive obesity. We exposed C57BL/6 mice to 29 different diets varying from 8.3% to 80% fat, 10% to 80% carbohydrate, 5% to 30% protein, and 5% to 30% sucrose. Only increased dietary fat content was associated with elevated energy intake and adiposity. This response was associated with increased gene expression in the 5-HT receptors, and the dopamine and opioid signaling pathways in the hypothalamus. We replicated the core findings in four other mouse strains (DBA/2, BALB/c, FVB, and C3H). Mice regulate their food consumption primarily to meet an energy rather than a protein target, but this system can be over-ridden by hedonic factors linked to fat, but not sucrose, consumption.
Collapse
|
10
|
JUM is a computational method for comprehensive annotation-free analysis of alternative pre-mRNA splicing patterns. Proc Natl Acad Sci U S A 2018; 115:E8181-E8190. [PMID: 30104386 DOI: 10.1073/pnas.1806018115] [Citation(s) in RCA: 44] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Alternative pre-mRNA splicing (AS) greatly diversifies metazoan transcriptomes and proteomes and is crucial for gene regulation. Current computational analysis methods of AS from Illumina RNA-sequencing data rely on preannotated libraries of known spliced transcripts, which hinders AS analysis with poorly annotated genomes and can further mask unknown AS patterns. To address this critical bioinformatics problem, we developed a method called the junction usage model (JUM) that uses a bottom-up approach to identify, analyze, and quantitate global AS profiles without any prior transcriptome annotations. JUM accurately reports global AS changes in terms of the five conventional AS patterns and an additional "composite" category composed of inseparable combinations of conventional patterns. JUM stringently classifies the difficult and disease-relevant pattern of intron retention (IR), reducing the false positive rate of IR detection commonly seen in other annotation-based methods to near-negligible rates. When analyzing AS in RNA samples derived from Drosophila heads, human tumors, and human cell lines bearing cancer-associated splicing factor mutations, JUM consistently identified approximately twice the number of novel AS events missed by other methods. Computational simulations showed JUM exhibits a 1.2 to 4.8 times higher true positive rate at a fixed cutoff of 5% false discovery rate. In summary, JUM provides a framework and improved method that removes the necessity for transcriptome annotations and enables the detection, analysis, and quantification of AS patterns in complex metazoan transcriptomes with superior accuracy.
Collapse
|
11
|
Wang D, Pentzold S, Kunert M, Groth M, Brandt W, Pasteels JM, Boland W, Burse A. A subset of chemosensory genes differs between two populations of a specialized leaf beetle after host plant shift. Ecol Evol 2018; 8:8055-8075. [PMID: 30250684 PMCID: PMC6145003 DOI: 10.1002/ece3.4246] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2018] [Revised: 04/13/2018] [Accepted: 04/17/2018] [Indexed: 01/26/2023] Open
Abstract
Due to its fundamental role in shaping host selection behavior, we have analyzed the chemosensory repertoire of Chrysomela lapponica. This specialized leaf beetle evolved distinct populations which shifted from the ancestral host plant, willow (Salix sp., Salicaceae), to birch (Betula rotundifolia, Betulaceae). We identified 114 chemosensory candidate genes in adult C. lapponica: 41 olfactory receptors (ORs), eight gustatory receptors, 17 ionotropic receptors, four sensory neuron membrane proteins, 32 odorant binding proteins (OBPs), and 12 chemosensory proteins (CSP) by RNA-seq. Differential expression analyses in the antennae revealed significant upregulation of one minus-C OBP (Clap OBP27) and one CSP (Clap CSP12) in the willow feeders. In contrast, one OR (Clap OR17), four minus-C OBPs (Clap OBP02, 07, 13, 20), and one plus-C OBP (Clap OBP32) were significantly upregulated in birch feeders. The differential expression pattern in the legs was more complex. To narrow down putative ligands acting as cues for host discrimination, the relative abundance and diversity of volatiles of the two host plant species were analyzed. In addition to salicylaldehyde (willow-specific), both plant species differed mainly in their emission rate of terpenoids such as (E,E)-α-farnesene (high in willow) or 4,8-dimethylnona-1,3,7-triene (high in birch). Qualitatively, the volatiles were similar between willow and birch leaves constituting an "olfactory bridge" for the beetles. Subsequent structural modeling of the three most differentially expressed OBPs and docking studies using 22 host volatiles indicated that ligands bind with varying affinity. We suggest that the evolution of particularly minus-C OBPs and ORs in C. lapponica facilitated its host plant shift via chemosensation of the phytochemicals from birch as novel host plant.
Collapse
Affiliation(s)
- Ding Wang
- Department of Bioorganic ChemistryMax Planck Institute for Chemical EcologyJenaGermany
| | - Stefan Pentzold
- Department of Bioorganic ChemistryMax Planck Institute for Chemical EcologyJenaGermany
| | - Maritta Kunert
- Department of Bioorganic ChemistryMax Planck Institute for Chemical EcologyJenaGermany
| | - Marco Groth
- Leibniz Institute on Aging – Fritz Lipmann InstituteJenaGermany
| | | | | | - Wilhelm Boland
- Department of Bioorganic ChemistryMax Planck Institute for Chemical EcologyJenaGermany
| | - Antje Burse
- Department of Bioorganic ChemistryMax Planck Institute for Chemical EcologyJenaGermany
| |
Collapse
|
12
|
Harris KD, Hochgerner H, Skene NG, Magno L, Katona L, Bengtsson Gonzales C, Somogyi P, Kessaris N, Linnarsson S, Hjerling-Leffler J. Classes and continua of hippocampal CA1 inhibitory neurons revealed by single-cell transcriptomics. PLoS Biol 2018; 16:e2006387. [PMID: 29912866 PMCID: PMC6029811 DOI: 10.1371/journal.pbio.2006387] [Citation(s) in RCA: 168] [Impact Index Per Article: 28.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2018] [Revised: 07/03/2018] [Accepted: 05/22/2018] [Indexed: 01/19/2023] Open
Abstract
Understanding any brain circuit will require a categorization of its constituent neurons. In hippocampal area CA1, at least 23 classes of GABAergic neuron have been proposed to date. However, this list may be incomplete; additionally, it is unclear whether discrete classes are sufficient to describe the diversity of cortical inhibitory neurons or whether continuous modes of variability are also required. We studied the transcriptomes of 3,663 CA1 inhibitory cells, revealing 10 major GABAergic groups that divided into 49 fine-scale clusters. All previously described and several novel cell classes were identified, with three previously described classes unexpectedly found to be identical. A division into discrete classes, however, was not sufficient to describe the diversity of these cells, as continuous variation also occurred between and within classes. Latent factor analysis revealed that a single continuous variable could predict the expression levels of several genes, which correlated similarly with it across multiple cell types. Analysis of the genes correlating with this variable suggested it reflects a range from metabolically highly active faster-spiking cells that proximally target pyramidal cells to slower-spiking cells targeting distal dendrites or interneurons. These results elucidate the complexity of inhibitory neurons in one of the simplest cortical structures and show that characterizing these cells requires continuous modes of variation as well as discrete cell classes.
Collapse
Affiliation(s)
- Kenneth D. Harris
- University College London Institute of Neurology, London, United Kingdom
- University College London Department of Neuroscience, Physiology and Pharmacology, London, United Kingdom
| | - Hannah Hochgerner
- Division of Molecular Neurobiology, Department of Medical Biochemistry and Biophysics, Karolinska Institutet, Stockholm, Sweden
| | - Nathan G. Skene
- University College London Institute of Neurology, London, United Kingdom
- Division of Molecular Neurobiology, Department of Medical Biochemistry and Biophysics, Karolinska Institutet, Stockholm, Sweden
| | - Lorenza Magno
- University College London Wolfson Institute for Biomedical Research, London, United Kingdom
| | - Linda Katona
- Department of Pharmacology, University of Oxford, Oxford, United Kingdom
| | - Carolina Bengtsson Gonzales
- Division of Molecular Neurobiology, Department of Medical Biochemistry and Biophysics, Karolinska Institutet, Stockholm, Sweden
| | - Peter Somogyi
- Department of Pharmacology, University of Oxford, Oxford, United Kingdom
| | - Nicoletta Kessaris
- University College London Wolfson Institute for Biomedical Research, London, United Kingdom
| | - Sten Linnarsson
- Division of Molecular Neurobiology, Department of Medical Biochemistry and Biophysics, Karolinska Institutet, Stockholm, Sweden
| | - Jens Hjerling-Leffler
- Division of Molecular Neurobiology, Department of Medical Biochemistry and Biophysics, Karolinska Institutet, Stockholm, Sweden
| |
Collapse
|
13
|
Xia Y, Sun J, Chen DG. Introductory Overview of Statistical Analysis of Microbiome Data. STATISTICAL ANALYSIS OF MICROBIOME DATA WITH R 2018. [DOI: 10.1007/978-981-13-1534-3_3] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
|
14
|
Ran D, Daye ZJ. Gene expression variability and the analysis of large-scale RNA-seq studies with the MDSeq. Nucleic Acids Res 2017; 45:e127. [PMID: 28535263 PMCID: PMC5737414 DOI: 10.1093/nar/gkx456] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2016] [Revised: 04/10/2017] [Accepted: 05/19/2017] [Indexed: 12/21/2022] Open
Abstract
Rapidly decreasing cost of next-generation sequencing has led to the recent availability of large-scale RNA-seq data, that empowers the analysis of gene expression variability, in addition to gene expression means. In this paper, we present the MDSeq, based on the coefficient of dispersion, to provide robust and computationally efficient analysis of both gene expression means and variability on RNA-seq counts. The MDSeq utilizes a novel reparametrization of the negative binomial to provide flexible generalized linear models (GLMs) on both the mean and dispersion. We address challenges of analyzing large-scale RNA-seq data via several new developments to provide a comprehensive toolset that models technical excess zeros, identifies outliers efficiently, and evaluates differential expressions at biologically interesting levels. We evaluated performances of the MDSeq using simulated data when the ground truths are known. Results suggest that the MDSeq often outperforms current methods for the analysis of gene expression mean and variability. Moreover, the MDSeq is applied in two real RNA-seq studies, in which we identified functionally relevant genes and gene pathways. Specifically, the analysis of gene expression variability with the MDSeq on the GTEx human brain tissue data has identified pathways associated with common neurodegenerative disorders when gene expression means were conserved.
Collapse
Affiliation(s)
- Di Ran
- Mel and Enid Zuckerman College of Public Health, The University of Arizona, Tucson, AZ 85724, USA
| | | |
Collapse
|
15
|
Lopez-Maestre H, Brinza L, Marchet C, Kielbassa J, Bastien S, Boutigny M, Monnin D, Filali AE, Carareto CM, Vieira C, Picard F, Kremer N, Vavre F, Sagot MF, Lacroix V. SNP calling from RNA-seq data without a reference genome: identification, quantification, differential analysis and impact on the protein sequence. Nucleic Acids Res 2016; 44:e148. [PMID: 27458203 PMCID: PMC5100560 DOI: 10.1093/nar/gkw655] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2015] [Accepted: 07/11/2016] [Indexed: 11/14/2022] Open
Abstract
SNPs (Single Nucleotide Polymorphisms) are genetic markers whose precise identification is a prerequisite for association studies. Methods to identify them are currently well developed for model species, but rely on the availability of a (good) reference genome, and therefore cannot be applied to non-model species. They are also mostly tailored for whole genome (re-)sequencing experiments, whereas in many cases, transcriptome sequencing can be used as a cheaper alternative which already enables to identify SNPs located in transcribed regions. In this paper, we propose a method that identifies, quantifies and annotates SNPs without any reference genome, using RNA-seq data only. Individuals can be pooled prior to sequencing, if not enough material is available from one individual. Using pooled human RNA-seq data, we clarify the precision and recall of our method and discuss them with respect to other methods which use a reference genome or an assembled transcriptome. We then validate experimentally the predictions of our method using RNA-seq data from two non-model species. The method can be used for any species to annotate SNPs and predict their impact on the protein sequence. We further enable to test for the association of the identified SNPs with a phenotype of interest.
Collapse
Affiliation(s)
- Hélène Lopez-Maestre
- Université de Lyon, F-69000, Lyon; Université Lyon 1; CNRS, UMR5558, Laboratoire de Biométrie et Biologie Evolutive, F-69622 Villeurbanne, France.,EPI ERABLE - Inria Grenoble, Rhône-Alpes
| | - Lilia Brinza
- PT Génomique et Transcriptomique, BIOASTER, Lyon, France
| | - Camille Marchet
- Université de Rennes, F-35000 Rennes; équipe GenScale, IRISA, Rennes
| | - Janice Kielbassa
- Synergie-Lyon-Cancer, Universite Lyon 1, Centre Leon Berard, Lyon, France
| | - Sylvère Bastien
- Université de Lyon, F-69000, Lyon; Université Lyon 1; CNRS, UMR5558, Laboratoire de Biométrie et Biologie Evolutive, F-69622 Villeurbanne, France.,EPI ERABLE - Inria Grenoble, Rhône-Alpes
| | - Mathilde Boutigny
- Université de Lyon, F-69000, Lyon; Université Lyon 1; CNRS, UMR5558, Laboratoire de Biométrie et Biologie Evolutive, F-69622 Villeurbanne, France.,EPI ERABLE - Inria Grenoble, Rhône-Alpes
| | - David Monnin
- Université de Lyon, F-69000, Lyon; Université Lyon 1; CNRS, UMR5558, Laboratoire de Biométrie et Biologie Evolutive, F-69622 Villeurbanne, France
| | - Adil El Filali
- Université de Lyon, F-69000, Lyon; Université Lyon 1; CNRS, UMR5558, Laboratoire de Biométrie et Biologie Evolutive, F-69622 Villeurbanne, France
| | - Claudia Marcia Carareto
- Department of Biology, UNESP - São Paulo State University, São José do Rio Preto, São Paulo, Brazil
| | - Cristina Vieira
- Université de Lyon, F-69000, Lyon; Université Lyon 1; CNRS, UMR5558, Laboratoire de Biométrie et Biologie Evolutive, F-69622 Villeurbanne, France.,EPI ERABLE - Inria Grenoble, Rhône-Alpes
| | - Franck Picard
- Université de Lyon, F-69000, Lyon; Université Lyon 1; CNRS, UMR5558, Laboratoire de Biométrie et Biologie Evolutive, F-69622 Villeurbanne, France
| | - Natacha Kremer
- Université de Lyon, F-69000, Lyon; Université Lyon 1; CNRS, UMR5558, Laboratoire de Biométrie et Biologie Evolutive, F-69622 Villeurbanne, France
| | - Fabrice Vavre
- Université de Lyon, F-69000, Lyon; Université Lyon 1; CNRS, UMR5558, Laboratoire de Biométrie et Biologie Evolutive, F-69622 Villeurbanne, France.,EPI ERABLE - Inria Grenoble, Rhône-Alpes
| | - Marie-France Sagot
- Université de Lyon, F-69000, Lyon; Université Lyon 1; CNRS, UMR5558, Laboratoire de Biométrie et Biologie Evolutive, F-69622 Villeurbanne, France.,EPI ERABLE - Inria Grenoble, Rhône-Alpes
| | - Vincent Lacroix
- Université de Lyon, F-69000, Lyon; Université Lyon 1; CNRS, UMR5558, Laboratoire de Biométrie et Biologie Evolutive, F-69622 Villeurbanne, France .,EPI ERABLE - Inria Grenoble, Rhône-Alpes
| |
Collapse
|
16
|
Konishi T. Parametric analysis of RNA-seq expression data. Genes Cells 2016; 21:639-47. [PMID: 27198878 DOI: 10.1111/gtc.12372] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/25/2015] [Accepted: 04/04/2016] [Indexed: 11/30/2022]
Abstract
Various methods had been introduced for normalization and comparison of RNA-seq count data. However, they lacked objectivity because they based on ad hoc assumptions that were never verified their appropriateness. Here, we introduced a method that assumes parsimony models on data distribution; the assumptions were verified according to exploratory data analysis. As was expected, count data were lognormally distributed. The level of noise in recent data appeared to be much higher than those of microarrays. Still, the appropriate distribution model would improve certainty and accuracy of normalization, by finding out the reliable range of data. Primary cause of noise was not the principle of the methodology; that is, each read is a trial that which transcript is read. Rather, the cause would be overlooking of transcripts, and the overlooking occurred more often among lower range of data. To find out genes likely to be overlooked, number of replications would be more important than read depth, which will not prevent overlooking. Both signal and noise in the reliable range of data were distributed normally, showing the suitability to use generalized linear model to evaluate differences in expression levels. In the framework, normalized data can be compared and combined freely beyond studies.
Collapse
Affiliation(s)
- Tomokazu Konishi
- Faculty of Bioresource Sciences, Akita Prefectural University, Akita, 010-0195, Japan
| |
Collapse
|
17
|
Farahmand S, Goliaei S, Ansari-Pour N, Razaghi-Moghadam Z. GTA: a game theoretic approach to identifying cancer subnetwork markers. MOLECULAR BIOSYSTEMS 2016; 12:818-25. [DOI: 10.1039/c5mb00684h] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
The identification of genetic markers (e.g. genes, pathways and subnetworks) for cancer has been one of the most challenging research areas in recent years.
Collapse
Affiliation(s)
- S. Farahmand
- Research Laboratory of Computational Biology
- Faculty of New Sciences and Technology
- University of Tehran
- Tehran
- Iran
| | - S. Goliaei
- Research Laboratory of Computational Biology
- Faculty of New Sciences and Technology
- University of Tehran
- Tehran
- Iran
| | - N. Ansari-Pour
- Faculty of New Sciences and Technology
- University of Tehran
- Tehran
- Iran
- School of Biological Sciences
| | - Z. Razaghi-Moghadam
- Faculty of New Sciences and Technology
- University of Tehran
- Tehran
- Iran
- School of Biological Sciences
| |
Collapse
|
18
|
Zhao G, Li Q, Wang IM, Liu X, Fang X, Zhang XD. An effective analytic method for detecting tissue-specific genes in RNA-seq experiments. Pharmacogenomics 2015; 16:1769-79. [PMID: 26554622 DOI: 10.2217/pgs.15.118] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
AIM To develop an analytic method for identifying tissue-specific (TS) genes from RNA-seq data. MATERIALS & METHODS Based on a negative binomial distribution, we develop a statistical method containing consecutive procedures incorporating data variability from replicates in each tissue. RESULTS Simulations show that our approach can effectively identify at least 94% of the truly TS genes if the sample size is 3 and at least 84% of the TS genes detected by our method are truly TS genes. We illustrated the utility of our method in an in-house RNA-seq project and produced sensible results. CONCLUSION Our approach not only directly works on discrete data but also naturally incorporates data variability. It works effectively for detecting TS genes.
Collapse
Affiliation(s)
- Guoqing Zhao
- School of Mathematical Sciences, Peking University, Beijing 100871, China.,BARDS,MSD R&D (China), Beijing 100015, China
| | - Qiao Li
- BARDS,MSD R&D (China), Beijing 100015, China
| | - I-Ming Wang
- Discovery Pharmacogenomics, Merck Research Laboratories, West Point, PA 19486, USA
| | - Xiaoqiao Liu
- Scientific Informatics, MSD R&D (China), Beijing 100015, China
| | - Xiangzhong Fang
- School of Mathematical Sciences, Peking University, Beijing 100871, China
| | | |
Collapse
|
19
|
Nguyen Y, Nettleton D, Liu H, Tuggle CK. Detecting Differentially Expressed Genes with RNA-seq Data Using Backward Selection to Account for the Effects of Relevant Covariates. JOURNAL OF AGRICULTURAL BIOLOGICAL AND ENVIRONMENTAL STATISTICS 2015; 20:577-597. [PMID: 26660449 PMCID: PMC4666287 DOI: 10.1007/s13253-015-0226-1] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/14/2015] [Accepted: 09/17/2015] [Indexed: 11/25/2022]
Abstract
A common challenge in analysis of transcriptomic data is to identify differentially expressed genes, i.e., genes whose mean transcript abundance levels differ across the levels of a factor of scientific interest. Transcript abundance levels can be measured simultaneously for thousands of genes in multiple biological samples using RNA sequencing (RNA-seq) technology. Part of the variation in RNA-seq measures of transcript abundance may be associated with variation in continuous and/or categorical covariates measured for each experimental unit or RNA sample. Ignoring relevant covariates or modeling the effects of irrelevant covariates can be detrimental to identifying differentially expressed genes. We propose a backward selection strategy for selecting a set of covariates whose effects are accounted for when searching for differentially expressed genes. We illustrate our approach through the analysis of an RNA-seq study intended to identify genes differentially expressed between two lines of pigs divergently selected for residual feed intake. We use simulation to show the advantages of our backward selection procedure over alternative strategies that either ignore or adjust for all measured covariates.
Collapse
Affiliation(s)
- Yet Nguyen
- Department of Statistics, Iowa State University, Ames, IA 50010 USA ; Institute of Mathematics, VAST, Hanoi, Vietnam
| | - Dan Nettleton
- Department of Statistics, Iowa State University, Ames, IA 50010 USA
| | - Haibo Liu
- Department of Animal Science, Iowa State University, Ames, IA 50010 USA
| | | |
Collapse
|
20
|
Yu F, Chen MH, Kuo L, Talbott H, Davis JS. Confident difference criterion: a new Bayesian differentially expressed gene selection algorithm with applications. BMC Bioinformatics 2015; 16:245. [PMID: 26250443 PMCID: PMC4527130 DOI: 10.1186/s12859-015-0664-3] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2014] [Accepted: 07/07/2015] [Indexed: 11/16/2022] Open
Abstract
Background Recently, the Bayesian method becomes more popular for analyzing high dimensional gene expression data as it allows us to borrow information across different genes and provides powerful estimators for evaluating gene expression levels. It is crucial to develop a simple but efficient gene selection algorithm for detecting differentially expressed (DE) genes based on the Bayesian estimators. Results In this paper, by extending the two-criterion idea of Chen et al. (Chen M-H, Ibrahim JG, Chi Y-Y. A new class of mixture models for differential gene expression in DNA microarray data. J Stat Plan Inference. 2008;138:387–404), we propose two new gene selection algorithms for general Bayesian models and name these new methods as the confident difference criterion methods. One is based on the standardized differences between two mean expression values among genes; the other adds the differences between two variances to it. The proposed confident difference criterion methods first evaluate the posterior probability of a gene having different gene expressions between competitive samples and then declare a gene to be DE if the posterior probability is large. The theoretical connection between the proposed first method based on the means and the Bayes factor approach proposed by Yu et al. (Yu F, Chen M-H, Kuo L. Detecting differentially expressed genes using alibrated Bayes factors. Statistica Sinica. 2008;18:783–802) is established under the normal-normal-model with equal variances between two samples. The empirical performance of the proposed methods is examined and compared to those of several existing methods via several simulations. The results from these simulation studies show that the proposed confident difference criterion methods outperform the existing methods when comparing gene expressions across different conditions for both microarray studies and sequence-based high-throughput studies. A real dataset is used to further demonstrate the proposed methodology. In the real data application, the confident difference criterion methods successfully identified more clinically important DE genes than the other methods. Conclusion The confident difference criterion method proposed in this paper provides a new efficient approach for both microarray studies and sequence-based high-throughput studies to identify differentially expressed genes. Electronic supplementary material The online version of this article (doi:10.1186/s12859-015-0664-3) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Fang Yu
- Department of Biostatistics, University of Nebraska Medical Center, Omaha, 68198-4350, NE, USA.
| | - Ming-Hui Chen
- Department of Statistics, University of Connecticut, Storrs, 06269-4120, CT, USA.
| | - Lynn Kuo
- Department of Statistics, University of Connecticut, Storrs, 06269-4120, CT, USA.
| | - Heather Talbott
- Department of Biochemistry and Molecular Biology and Department of Obstetrics and Gynecology, University of Nebraska Medical Center, Omaha, 68198-5870, NE, USA.
| | - John S Davis
- VA Nebraska-Western Iowa Health Care System and Department of Obstetrics and Gynecology, University of Nebraska Medical Center, Omaha, 68198-3255, NE, USA.
| |
Collapse
|
21
|
Farkas MH, Au ED, Sousa ME, Pierce EA. RNA-Seq: Improving Our Understanding of Retinal Biology and Disease. Cold Spring Harb Perspect Med 2015; 5:a017152. [PMID: 25722474 PMCID: PMC4561396 DOI: 10.1101/cshperspect.a017152] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Over the past several years, rapid technological advances have allowed for a dramatic increase in our knowledge and understanding of the transcriptional landscape, because of the ability to study gene expression in greater depth and with more detail than previously possible. To this end, RNA-Seq has quickly become one of the most widely used methods for studying transcriptomes of tissues and individual cells. Unlike previously favored analysis methods, RNA-Seq is extremely high-throughput, and is not dependent on an annotated transcriptome, laying the foundation for novel genetic discovery. Additionally, RNA-Seq derived transcriptomes provide a basis for widening the scope of research to identify potential targets in the treatment of retinal disease.
Collapse
Affiliation(s)
- Michael H Farkas
- Ocular Genomics Institute, Department of Ophthalmology, Massachusetts Eye and Ear Infirmary, Harvard Medical School, Boston, Massachusetts 02114
| | - Elizabeth D Au
- Ocular Genomics Institute, Department of Ophthalmology, Massachusetts Eye and Ear Infirmary, Harvard Medical School, Boston, Massachusetts 02114
| | - Maria E Sousa
- Ocular Genomics Institute, Department of Ophthalmology, Massachusetts Eye and Ear Infirmary, Harvard Medical School, Boston, Massachusetts 02114
| | - Eric A Pierce
- Ocular Genomics Institute, Department of Ophthalmology, Massachusetts Eye and Ear Infirmary, Harvard Medical School, Boston, Massachusetts 02114
| |
Collapse
|
22
|
Vavoulis DV, Francescatto M, Heutink P, Gough J. DGEclust: differential expression analysis of clustered count data. Genome Biol 2015; 16:39. [PMID: 25853652 PMCID: PMC4365804 DOI: 10.1186/s13059-015-0604-6] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2014] [Accepted: 02/03/2015] [Indexed: 11/10/2022] Open
Abstract
We present a statistical methodology, DGEclust, for differential expression analysis of digital expression data. Our method treats differential expression as a form of clustering, thus unifying these two concepts. Furthermore, it simultaneously addresses the problem of how many clusters are supported by the data and uncertainty in parameter estimation. DGEclust successfully identifies differentially expressed genes under a number of different scenarios, maintaining a low error rate and an excellent control of its false discovery rate with reasonable computational requirements. It is formulated to perform particularly well on low-replicated data and be applicable to multi-group data. DGEclust is available at http://dvav.github.io/dgeclust/.
Collapse
Affiliation(s)
| | - Margherita Francescatto
- />Genome Biology of Neurodegenerative Diseases, Deutsches Zentrum für Neurodegenerative Erkrankungen, Tübingen, Germany
| | - Peter Heutink
- />Genome Biology of Neurodegenerative Diseases, Deutsches Zentrum für Neurodegenerative Erkrankungen, Tübingen, Germany
| | - Julian Gough
- />Department of Computer Science, University of Bristol, Bristol, UK
| |
Collapse
|
23
|
Zhang H, Xu J, Jiang N, Hu X, Luo Z. PLNseq: a multivariate Poisson lognormal distribution for high-throughput matched RNA-sequencing read count data. Stat Med 2015; 34:1577-89. [PMID: 25641202 DOI: 10.1002/sim.6449] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2014] [Revised: 12/17/2014] [Accepted: 01/17/2015] [Indexed: 01/28/2023]
Abstract
High-throughput RNA-sequencing (RNA-seq) technology provides an attractive platform for gene expression analysis. In many experimental settings, RNA-seq read counts are measured from matched samples or taken from the same subject under multiple treatment conditions. The induced correlation therefore should be evaluated and taken into account in deriving tests of differential expression. We proposed a novel method 'PLNseq', which uses a multivariate Poisson lognormal distribution to model matched read count data. The correlation is directly modeled through Gaussian random effects, and inferences are made by likelihood methods. A three-stage numerical algorithm is developed to estimate unknown parameters and conduct differential expression analysis. Results using simulated data demonstrate that our method performs reasonably well in terms of parameter estimation, DE analysis power, and robustness. PLNseq also has better control of FDRs than the benchmarks edgeR and DESeq2 in the situations where the correlation is different across the genes but can still be accurately estimated. Furthermore, direct evaluation of correlation through PLNseq enables us to develop a new and more powerful test for DE analysis. Application to a lung cancer study is provided to illustrate the practical utilities of our method. An R package implementing the method is also publicly available.
Collapse
Affiliation(s)
- Hong Zhang
- Department of Biostatistics and Computational Biology, School of Life Sciences, Fudan University, China
| | | | | | | | | |
Collapse
|
24
|
Luo H, Li J, Chia BKH, Robson P, Nagarajan N. The importance of study design for detecting differentially abundant features in high-throughput experiments. Genome Biol 2014; 15:527. [PMID: 25517037 PMCID: PMC4253014 DOI: 10.1186/s13059-014-0527-7] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2014] [Accepted: 11/03/2014] [Indexed: 12/20/2022] Open
Abstract
High-throughput assays, such as RNA-seq, to detect differential abundance are widely used. Variable performance across statistical tests, normalizations, and conditions leads to resource wastage and reduced sensitivity. EDDA represents a first, general design tool for RNA-seq, Nanostring, and metagenomic analysis, that rationally selects tests, predicts performance, and plans experiments to minimize resource wastage. Case studies highlight EDDA’s ability to model single-cell RNA-seq, suggesting ways to reduce sequencing costs up to five-fold and improving metagenomic biomarker detection through improved test selection. EDDA’s novel mode-based normalization for detecting differential abundance improves robustness by 10% to 20% and precision by up to 140%.
Collapse
|
25
|
Williams AG, Thomas S, Wyman SK, Holloway AK. RNA-seq Data: Challenges in and Recommendations for Experimental Design and Analysis. ACTA ACUST UNITED AC 2014; 83:11.13.1-20. [PMID: 25271838 DOI: 10.1002/0471142905.hg1113s83] [Citation(s) in RCA: 59] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
RNA-seq is widely used to determine differential expression of genes or transcripts as well as identify novel transcripts, identify allele-specific expression, and precisely measure translation of transcripts. Thoughtful experimental design and choice of analysis tools are critical to ensure high-quality data and interpretable results. Important considerations for experimental design include number of replicates, whether to collect paired-end or single-end reads, sequence length, and sequencing depth. Common analysis steps in all RNA-seq experiments include quality control, read alignment, assigning reads to genes or transcripts, and estimating gene or transcript abundance. Our aims are two-fold: to make recommendations for common components of experimental design and assess tool capabilities for each of these steps. We also test tools designed to detect differential expression, since this is the most widespread application of RNA-seq. We hope that these analyses will help guide those who are new to RNA-seq and will generate discussion about remaining needs for tool improvement and development.
Collapse
|
26
|
Gu X, Zou Y, Huang W, Shen L, Arendsee Z, Su Z. Phylogenomic distance method for analyzing transcriptome evolution based on RNA-seq data. Genome Biol Evol 2014; 5:1746-53. [PMID: 23940099 PMCID: PMC3787673 DOI: 10.1093/gbe/evt121] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
Thanks to the microarray technology, our understanding of transcriptome evolution at the genome level has been considerably advanced in the past decade. Yet, further investigation was challenged by several technical limitations of this technology. Recent innovation of next-generation sequencing, particularly the invention of RNA-seq technology, has shed insightful lights on resolving this problem. Though a number of statistical and computational methods have been developed to analyze RNA-seq data, the analytical framework specifically designed for evolutionary genomics remains an open question. In this article we develop a new method for estimating the genome expression distance from the RNA-seq data, which has explicit interpretations under the model of gene expression evolution. Moreover, this distance measure takes the data overdispersion, gene length variation, and sequencing depth variation into account so that it can be applied to multiple genomes from different species. Using mammalian RNA-seq data as example, we demonstrated that this expression distance is useful in phylogenomic analysis.
Collapse
Affiliation(s)
- Xun Gu
- State Key Laboratory of Genetic Engineering and MOE Key Laboratory of Contemporary Anthropology, School of Life Sciences, Fudan University, Shanghai, China
| | | | | | | | | | | |
Collapse
|
27
|
McMurdie PJ, Holmes S. Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol 2014; 10:e1003531. [PMID: 24699258 PMCID: PMC3974642 DOI: 10.1371/journal.pcbi.1003531] [Citation(s) in RCA: 1689] [Impact Index Per Article: 168.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2013] [Accepted: 02/03/2014] [Indexed: 02/07/2023] Open
Abstract
Current practice in the normalization of microbiome count data is inefficient in the statistical sense. For apparently historical reasons, the common approach is either to use simple proportions (which does not address heteroscedasticity) or to use rarefying of counts, even though both of these approaches are inappropriate for detection of differentially abundant species. Well-established statistical theory is available that simultaneously accounts for library size differences and biological variability using an appropriate mixture model. Moreover, specific implementations for DNA sequencing read count data (based on a Negative Binomial model for instance) are already available in RNA-Seq focused R packages such as edgeR and DESeq. Here we summarize the supporting statistical theory and use simulations and empirical data to demonstrate substantial improvements provided by a relevant mixture model framework over simple proportions or rarefying. We show how both proportions and rarefied counts result in a high rate of false positives in tests for species that are differentially abundant across sample classes. Regarding microbiome sample-wise clustering, we also show that the rarefying procedure often discards samples that can be accurately clustered by alternative methods. We further compare different Negative Binomial methods with a recently-described zero-inflated Gaussian mixture, implemented in a package called metagenomeSeq. We find that metagenomeSeq performs well when there is an adequate number of biological replicates, but it nevertheless tends toward a higher false positive rate. Based on these results and well-established statistical theory, we advocate that investigators avoid rarefying altogether. We have provided microbiome-specific extensions to these tools in the R package, phyloseq. The term microbiome refers to the ecosystem of microbes that live in a defined environment. The decreasing cost and increasing speed of DNA sequencing technology has recently provided scientists with affordable and timely access to the genes and genomes of microbiomes that inhabit our planet and even our own bodies. In these investigations many microbiome samples are sequenced at the same time on the same DNA sequencing machine, but often result in total numbers of sequences per sample that are vastly different. The common procedure for addressing this difference in sequencing effort across samples – different library sizes – is to either (1) base analyses on the proportional abundance of each species in a library, or (2) rarefy, throw away sequences from the larger libraries so that all have the same, smallest size. We show that both of these normalization methods can work when comparing obviously-different whole microbiomes, but that neither method works well when comparing the relative proportions of each bacterial species across microbiome samples. We show that alternative methods based on a statistical mixture model perform much better and can be easily adapted from a separate biological sub-discipline, called RNA-Seq analysis.
Collapse
Affiliation(s)
- Paul J. McMurdie
- Statistics Department, Stanford University, Stanford, California, United States of America
| | - Susan Holmes
- Statistics Department, Stanford University, Stanford, California, United States of America
- * E-mail:
| |
Collapse
|
28
|
Oh S, Song S, Dasgupta N, Grabowski G. The analytical landscape of static and temporal dynamics in transcriptome data. Front Genet 2014; 5:35. [PMID: 24600473 PMCID: PMC3929947 DOI: 10.3389/fgene.2014.00035] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2013] [Accepted: 01/30/2014] [Indexed: 12/16/2022] Open
Abstract
Interpreting gene expression profiles often involves statistical analysis of large numbers of differentially expressed genes, isoforms, and alternative splicing events at either static or dynamic spectrums. Reduced sequencing costs have made feasible dense time-series analysis of gene expression using RNA-seq; however, statistical methods in the context of temporal RNA-seq data are poorly developed. Here we will review current methods for identifying temporal changes in gene expression using RNA-seq, which are limited to static pairwise comparisons of time points and which fail to account for temporal dependencies in gene expression patterns. We also review recently developed very few number of temporal dynamic RNA-seq specific methods. Application and development of RNA-specific temporal dynamic methods have been continuously under the development, yet, it is still in infancy. We fully cover microarray specific temporal methods and transcriptome studies in initial digital technology (e.g., SAGE) between traditional microarray and new RNA-seq.
Collapse
Affiliation(s)
- Sunghee Oh
- Division of Human Genetics, Department of Pediatrics, Cincinnati Children's Hospital Medical Center Cincinnati, OH, USA
| | - Seongho Song
- Department of Mathematical Sciences, McMicken College of Arts and Sciences, University of Cincinnati Cincinnati, OH, USA
| | - Nupur Dasgupta
- Division of Human Genetics, Department of Pediatrics, Cincinnati Children's Hospital Medical Center Cincinnati, OH, USA
| | - Gregory Grabowski
- Division of Human Genetics, Department of Pediatrics, Cincinnati Children's Hospital Medical Center Cincinnati, OH, USA
| |
Collapse
|
29
|
Saha KK. Inference concerning a common dispersion of several treatment groups in the analysis of over/underdispersed count data. Biom J 2014; 56:441-60. [DOI: 10.1002/bimj.201300105] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2013] [Revised: 09/24/2013] [Accepted: 10/25/2013] [Indexed: 11/10/2022]
Affiliation(s)
- Krishna K. Saha
- Department of Mathematical Sciences; Central Connecticut State University; New Britain CT 06050 USA
| |
Collapse
|
30
|
Chung LM, Ferguson JP, Zheng W, Qian F, Bruno V, Montgomery RR, Zhao H. Differential expression analysis for paired RNA-Seq data. BMC Bioinformatics 2013; 14:110. [PMID: 23530607 PMCID: PMC3663822 DOI: 10.1186/1471-2105-14-110] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2012] [Accepted: 03/01/2013] [Indexed: 11/18/2022] Open
Abstract
Background RNA-Seq technology measures the transcript abundance by generating sequence reads and counting their frequencies across different biological conditions. To identify differentially expressed genes between two conditions, it is important to consider the experimental design as well as the distributional property of the data. In many RNA-Seq studies, the expression data are obtained as multiple pairs, e.g., pre- vs. post-treatment samples from the same individual. We seek to incorporate paired structure into analysis. Results We present a Bayesian hierarchical mixture model for RNA-Seq data to separately account for the variability within and between individuals from a paired data structure. The method assumes a Poisson distribution for the data mixed with a gamma distribution to account variability between pairs. The effect of differential expression is modeled by two-component mixture model. The performance of this approach is examined by simulated and real data. Conclusions In this setting, our proposed model provides higher sensitivity than existing methods to detect differential expression. Application to real RNA-Seq data demonstrates the usefulness of this method for detecting expression alteration for genes with low average expression levels or shorter transcript length.
Collapse
Affiliation(s)
- Lisa M Chung
- Department of Biostatistics, Yale School of Public Health, New Haven, CT, USA.
| | | | | | | | | | | | | |
Collapse
|
31
|
Abstract
RNA-seq is a powerful tool for the study of alternative splicing and other forms of alternative isoform expression. Understanding the regulation of these processes requires sensitive and specific detection of differential isoform abundance in comparisons between conditions, cell types, or tissues. We present DEXSeq, a statistical method to test for differential exon usage in RNA-seq data. DEXSeq uses generalized linear models and offers reliable control of false discoveries by taking biological variation into account. DEXSeq detects with high sensitivity genes, and in many cases exons, that are subject to differential exon usage. We demonstrate the versatility of DEXSeq by applying it to several data sets. The method facilitates the study of regulation and function of alternative exon usage on a genome-wide scale. An implementation of DEXSeq is available as an R/Bioconductor package.
Collapse
Affiliation(s)
- Simon Anders
- European Molecular Biology Laboratory, 69111 Heidelberg, Germany.
| | | | | |
Collapse
|
32
|
McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res 2012. [PMID: 22287627 DOI: 10.1093/nar/gks042%jnucleicacidsresearch] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/18/2023] Open
Abstract
A flexible statistical framework is developed for the analysis of read counts from RNA-Seq gene expression studies. It provides the ability to analyse complex experiments involving multiple treatment conditions and blocking variables while still taking full account of biological variation. Biological variation between RNA samples is estimated separately from the technical variation associated with sequencing technologies. Novel empirical Bayes methods allow each gene to have its own specific variability, even when there are relatively few biological replicates from which to estimate such variability. The pipeline is implemented in the edgeR package of the Bioconductor project. A case study analysis of carcinoma data demonstrates the ability of generalized linear model methods (GLMs) to detect differential expression in a paired design, and even to detect tumour-specific expression changes. The case study demonstrates the need to allow for gene-specific variability, rather than assuming a common dispersion across genes or a fixed relationship between abundance and variability. Genewise dispersions de-prioritize genes with inconsistent results and allow the main analysis to focus on changes that are consistent between biological replicates. Parallel computational approaches are developed to make non-linear model fitting faster and more reliable, making the application of GLMs to genomic data more convenient and practical. Simulations demonstrate the ability of adjusted profile likelihood estimators to return accurate estimators of biological variability in complex situations. When variation is gene-specific, empirical Bayes estimators provide an advantageous compromise between the extremes of assuming common dispersion or separate genewise dispersion. The methods developed here can also be applied to count data arising from DNA-Seq applications, including ChIP-Seq for epigenetic marks and DNA methylation analyses.
Collapse
Affiliation(s)
- Davis J McCarthy
- Bioinformatics Division, The Walter and Eliza Hall Institute of Medical Research, 1G Royal Parade, Parkville, Victoria 3052, Australia
| | | | | |
Collapse
|
33
|
Kvam VM, Liu P, Si Y. A comparison of statistical methods for detecting differentially expressed genes from RNA-seq data. AMERICAN JOURNAL OF BOTANY 2012; 99:248-56. [PMID: 22268221 DOI: 10.3732/ajb.1100340] [Citation(s) in RCA: 156] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/18/2023]
Abstract
RNA-Seq technologies are quickly revolutionizing genomic studies, and statistical methods for RNA-seq data are under continuous development. Timely review and comparison of the most recently proposed statistical methods will provide a useful guide for choosing among them for data analysis. Particular interest surrounds the ability to detect differential expression (DE) in genes. Here we compare four recently proposed statistical methods, edgeR, DESeq, baySeq, and a method with a two-stage Poisson model (TSPM), through a variety of simulations that were based on different distribution models or real data. We compared the ability of these methods to detect DE genes in terms of the significance ranking of genes and false discovery rate control. All methods compared are implemented in freely available software. We also discuss the availability and functions of the currently available versions of these software.
Collapse
Affiliation(s)
- Vanessa M Kvam
- Department of Statistics, Iowa State University, Snedecor Hall, Ames, Iowa 50011-1210, USA
| | | | | |
Collapse
|
34
|
McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res 2012; 40:4288-97. [PMID: 22287627 PMCID: PMC3378882 DOI: 10.1093/nar/gks042] [Citation(s) in RCA: 3366] [Impact Index Per Article: 280.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
A flexible statistical framework is developed for the analysis of read counts from RNA-Seq gene expression studies. It provides the ability to analyse complex experiments involving multiple treatment conditions and blocking variables while still taking full account of biological variation. Biological variation between RNA samples is estimated separately from the technical variation associated with sequencing technologies. Novel empirical Bayes methods allow each gene to have its own specific variability, even when there are relatively few biological replicates from which to estimate such variability. The pipeline is implemented in the edgeR package of the Bioconductor project. A case study analysis of carcinoma data demonstrates the ability of generalized linear model methods (GLMs) to detect differential expression in a paired design, and even to detect tumour-specific expression changes. The case study demonstrates the need to allow for gene-specific variability, rather than assuming a common dispersion across genes or a fixed relationship between abundance and variability. Genewise dispersions de-prioritize genes with inconsistent results and allow the main analysis to focus on changes that are consistent between biological replicates. Parallel computational approaches are developed to make non-linear model fitting faster and more reliable, making the application of GLMs to genomic data more convenient and practical. Simulations demonstrate the ability of adjusted profile likelihood estimators to return accurate estimators of biological variability in complex situations. When variation is gene-specific, empirical Bayes estimators provide an advantageous compromise between the extremes of assuming common dispersion or separate genewise dispersion. The methods developed here can also be applied to count data arising from DNA-Seq applications, including ChIP-Seq for epigenetic marks and DNA methylation analyses.
Collapse
Affiliation(s)
- Davis J McCarthy
- Bioinformatics Division, The Walter and Eliza Hall Institute of Medical Research, 1G Royal Parade, Parkville, Victoria 3052, Australia
| | | | | |
Collapse
|
35
|
Li J, Tibshirani R. Finding consistent patterns: a nonparametric approach for identifying differential expression in RNA-Seq data. Stat Methods Med Res 2011; 22:519-36. [PMID: 22127579 DOI: 10.1177/0962280211428386] [Citation(s) in RCA: 284] [Impact Index Per Article: 21.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
We discuss the identification of features that are associated with an outcome in RNA-Sequencing (RNA-Seq) and other sequencing-based comparative genomic experiments. RNA-Seq data takes the form of counts, so models based on the normal distribution are generally unsuitable. The problem is especially challenging because different sequencing experiments may generate quite different total numbers of reads, or 'sequencing depths'. Existing methods for this problem are based on Poisson or negative binomial models: they are useful but can be heavily influenced by 'outliers' in the data. We introduce a simple, non-parametric method with resampling to account for the different sequencing depths. The new method is more robust than parametric methods. It can be applied to data with quantitative, survival, two-class or multiple-class outcomes. We compare our proposed method to Poisson and negative binomial-based methods in simulated and real data sets, and find that our method discovers more consistent patterns than competing methods.
Collapse
Affiliation(s)
- Jun Li
- 1Department of Statistics, Stanford University, Stanford, CA 94305, USA
| | | |
Collapse
|
36
|
Li J, Witten DM, Johnstone IM, Tibshirani R. Normalization, testing, and false discovery rate estimation for RNA-sequencing data. Biostatistics 2011; 13:523-38. [PMID: 22003245 DOI: 10.1093/biostatistics/kxr031] [Citation(s) in RCA: 214] [Impact Index Per Article: 16.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
We discuss the identification of genes that are associated with an outcome in RNA sequencing and other sequence-based comparative genomic experiments. RNA-sequencing data take the form of counts, so models based on the Gaussian distribution are unsuitable. Moreover, normalization is challenging because different sequencing experiments may generate quite different total numbers of reads. To overcome these difficulties, we use a log-linear model with a new approach to normalization. We derive a novel procedure to estimate the false discovery rate (FDR). Our method can be applied to data with quantitative, two-class, or multiple-class outcomes, and the computation is fast even for large data sets. We study the accuracy of our approaches for significance calculation and FDR estimation, and we demonstrate that our method has potential advantages over existing methods that are based on a Poisson or negative binomial model. In summary, this work provides a pipeline for the significance analysis of sequencing data.
Collapse
Affiliation(s)
- Jun Li
- Department of Statistics, Stanford University, Stanford, CA 94305, USA.
| | | | | | | |
Collapse
|
37
|
Jacob BJ, Baxter NN, Moineddin R, Sutradhar R, Del Giudice L, Urbach DR. Social disparities in the use of colonoscopy by primary care physicians in Ontario. BMC Gastroenterol 2011; 11:102. [PMID: 21955593 PMCID: PMC3206464 DOI: 10.1186/1471-230x-11-102] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/28/2010] [Accepted: 09/28/2011] [Indexed: 11/18/2022] Open
Abstract
Background It is unclear if all persons in Ontario have equal access to colonoscopy. This research was designed to describe long-term trends in the use of colonoscopy by primary care physicians (PCPs) in Ontario, and to determine whether PCP characteristics influence the use of colonoscopy. Methods We conducted a population-based retrospective study of PCPs in Ontario between the years 1996-2005. Using administrative data we identified a screen-eligible group of patients aged 50-74 years in Ontario. These patients were linked to the PCP who provided the most continuous care to them during each year. We determined the use of any colonoscopy among these patients. We calculated the rate of colonoscopy for each PCP as the number of patients undergoing colonoscopies per 100 screen eligible patients. Negative binomial regression was used to identify factors associated with the rate of colonoscopy, using generalized estimating equations to account for clustering of patients within PCPs. Results Between 7,955 and 8,419 PCPs in Ontario per year (median age 43 years) had at least 10 eligible patients in their practices. The use of colonoscopy by PCPs increased sharply in Ontario during the study period, from a median rate of 1.51 [inter quartile range (IQR) 0.57-2.62] per 100 screen eligible patients in 1996 to 4.71 (IQR 2.70-7.53) in 2005. There was substantial variation between PCPs in their use of colonoscopy. PCPs who were Canadian medical graduates and with more years of experience were more likely to use colonoscopy after adjusting for their patient characteristics. PCPs were more likely to use colonoscopy if their patient populations were predominantly women, older, had more illnesses, and if their patients resided in less marginalized neighborhoods (lower unemployment, fewer immigrants, higher income, higher education, and higher English/French fluency). Conclusions There is substantial variation in the use of colonoscopy by PCPs, and this variation has increased as the overall use of colonoscopy increased over time. PCPs whose patients were more marginalized were less likely to use colonoscopy, suggesting that there are inequities in access.
Collapse
Affiliation(s)
- Binu J Jacob
- Clinical Decision Making & Health Care, Toronto General Hospital, 200 Elizabeth Street, Toronto, Ontario M5G2C4, Canada.
| | | | | | | | | | | |
Collapse
|
38
|
Wu ZJ, Meyer CA, Choudhury S, Shipitsin M, Maruyama R, Bessarabova M, Nikolskaya T, Sukumar S, Schwartzman A, Liu JS, Polyak K, Liu XS. Gene expression profiling of human breast tissue samples using SAGE-Seq. Genome Res 2010; 20:1730-9. [PMID: 21045080 PMCID: PMC2989999 DOI: 10.1101/gr.108217.110] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2010] [Accepted: 09/24/2010] [Indexed: 12/29/2022]
Abstract
We present a powerful application of ultra high-throughput sequencing, SAGE-Seq, for the accurate quantification of normal and neoplastic mammary epithelial cell transcriptomes. We develop data analysis pipelines that allow the mapping of sense and antisense strands of mitochondrial and RefSeq genes, the normalization between libraries, and the identification of differentially expressed genes. We find that the diversity of cancer transcriptomes is significantly higher than that of normal cells. Our analysis indicates that transcript discovery plateaus at 10 million reads/sample, and suggests a minimum desired sequencing depth around five million reads. Comparison of SAGE-Seq and traditional SAGE on normal and cancerous breast tissues reveals higher sensitivity of SAGE-Seq to detect less-abundant genes, including those encoding for known breast cancer-related transcription factors and G protein-coupled receptors (GPCRs). SAGE-Seq is able to identify genes and pathways abnormally activated in breast cancer that traditional SAGE failed to call. SAGE-Seq is a powerful method for the identification of biomarkers and therapeutic targets in human disease.
Collapse
Affiliation(s)
- Zhenhua Jeremy Wu
- Department of Biostatistics and Computational Biology Dana-Farber Cancer Institute, Boston, Massachusetts 02115, USA
- Harvard School of Public Health, Boston, Massachusetts 02115, USA
| | - Clifford A. Meyer
- Department of Biostatistics and Computational Biology Dana-Farber Cancer Institute, Boston, Massachusetts 02115, USA
- Harvard School of Public Health, Boston, Massachusetts 02115, USA
| | - Sibgat Choudhury
- Department of Medical Oncology, Dana-Farber Cancer Institute, Boston, Massachusetts 02115, USA
- Department of Medicine, Brigham and Women's Hospital, Boston, Massachusetts 02115, USA
- Department of Medicine, Harvard Medical School, Boston, Massachusetts 02115, USA
| | - Michail Shipitsin
- Department of Medical Oncology, Dana-Farber Cancer Institute, Boston, Massachusetts 02115, USA
- Department of Medicine, Brigham and Women's Hospital, Boston, Massachusetts 02115, USA
- Department of Medicine, Harvard Medical School, Boston, Massachusetts 02115, USA
| | - Reo Maruyama
- Department of Medical Oncology, Dana-Farber Cancer Institute, Boston, Massachusetts 02115, USA
- Department of Medicine, Brigham and Women's Hospital, Boston, Massachusetts 02115, USA
- Department of Medicine, Harvard Medical School, Boston, Massachusetts 02115, USA
| | - Marina Bessarabova
- Vavilov Institute for General Genetics, Russian Academy of Sciences, Moscow 119331, Russia
| | - Tatiana Nikolskaya
- Vavilov Institute for General Genetics, Russian Academy of Sciences, Moscow 119331, Russia
| | | | - Armin Schwartzman
- Department of Biostatistics and Computational Biology Dana-Farber Cancer Institute, Boston, Massachusetts 02115, USA
- Harvard School of Public Health, Boston, Massachusetts 02115, USA
| | - Jun S. Liu
- Department of Statistics, Harvard University, Science Center 715, Cambridge, Massachusetts 02138, USA
| | - Kornelia Polyak
- Department of Medical Oncology, Dana-Farber Cancer Institute, Boston, Massachusetts 02115, USA
- Department of Medicine, Brigham and Women's Hospital, Boston, Massachusetts 02115, USA
- Department of Medicine, Harvard Medical School, Boston, Massachusetts 02115, USA
| | - X. Shirley Liu
- Department of Biostatistics and Computational Biology Dana-Farber Cancer Institute, Boston, Massachusetts 02115, USA
- Harvard School of Public Health, Boston, Massachusetts 02115, USA
| |
Collapse
|
39
|
Hardcastle TJ, Kelly KA. baySeq: empirical Bayesian methods for identifying differential expression in sequence count data. BMC Bioinformatics 2010; 11:422. [PMID: 20698981 PMCID: PMC2928208 DOI: 10.1186/1471-2105-11-422] [Citation(s) in RCA: 499] [Impact Index Per Article: 35.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2010] [Accepted: 08/10/2010] [Indexed: 11/24/2022] Open
Abstract
Background High throughput sequencing has become an important technology for studying expression levels in many types of genomic, and particularly transcriptomic, data. One key way of analysing such data is to look for elements of the data which display particular patterns of differential expression in order to take these forward for further analysis and validation. Results We propose a framework for defining patterns of differential expression and develop a novel algorithm, baySeq, which uses an empirical Bayes approach to detect these patterns of differential expression within a set of sequencing samples. The method assumes a negative binomial distribution for the data and derives an empirically determined prior distribution from the entire dataset. We examine the performance of the method on real and simulated data. Conclusions Our method performs at least as well, and often better, than existing methods for analyses of pairwise differential expression in both real and simulated data. When we compare methods for the analysis of data from experimental designs involving multiple sample groups, our method again shows substantial gains in performance. We believe that this approach thus represents an important step forward for the analysis of count data from sequencing experiments.
Collapse
Affiliation(s)
- Thomas J Hardcastle
- Department of Plant Sciences, University of Cambridge, Downing Street, Cambridge, UK.
| | | |
Collapse
|
40
|
Yu F, Chen MH, Kuo L, Huang P, Yang W. Bayesian hierarchical modeling and selection of differentially expressed genes for the EST data. Biometrics 2010; 67:142-50. [PMID: 20560937 DOI: 10.1111/j.1541-0420.2010.01447.x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Expressed sequence tag (EST) sequencing is a one-pass sequencing reading of cloned cDNAs derived from a certain tissue. The frequency of unique tags among different unbiased cDNA libraries is used to infer the relative expression level of each tag. In this article, we propose a hierarchical multinomial model with a nonlinear Dirichlet prior for the EST data with multiple libraries and multiple types of tissues. A novel hierarchical prior is developed and the properties of the proposed prior are examined. An efficient Markov chain Monte Carlo algorithm is developed for carrying out the posterior computation. We also propose a new selection criterion for detecting which genes are differentially expressed between two tissue types. Our new method with the new gene selection criterion is demonstrated via several simulations to have low false negative and false positive rates. A real EST data set is used to motivate and illustrate the proposed method.
Collapse
Affiliation(s)
- Fang Yu
- Department of Biostatistics, University of Nebraska Medical Center, Omaha, Nebraska 68198-4350,
| | | | | | | | | |
Collapse
|
41
|
Abstract
Next-generation sequencing technologies are quickly becoming the preferred approach for characterizing and quantifying entire genomes. Even though data produced from these technologies are proving to be the most informative of any thus far, very little attention has been paid to fundamental design aspects of data collection and analysis, namely sampling, randomization, replication, and blocking. We discuss these concepts in an RNA sequencing framework. Using simulations we demonstrate the benefits of collecting replicated RNA sequencing data according to well known statistical designs that partition the sources of biological and technical variation. Examples of these designs and their corresponding models are presented with the goal of testing differential expression.
Collapse
|
42
|
Silva IT, Vêncio RZN, Oliveira TYK, Molfetta GA, Silva WA. ProbFAST: Probabilistic functional analysis system tool. BMC Bioinformatics 2010; 11:161. [PMID: 20353576 PMCID: PMC2868004 DOI: 10.1186/1471-2105-11-161] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2009] [Accepted: 03/30/2010] [Indexed: 01/01/2023] Open
Abstract
Background The post-genomic era has brought new challenges regarding the understanding of the organization and function of the human genome. Many of these challenges are centered on the meaning of differential gene regulation under distinct biological conditions and can be performed by analyzing the Multiple Differential Expression (MDE) of genes associated with normal and abnormal biological processes. Currently MDE analyses are limited to usual methods of differential expression initially designed for paired analysis. Results We proposed a web platform named ProbFAST for MDE analysis which uses Bayesian inference to identify key genes that are intuitively prioritized by means of probabilities. A simulated study revealed that our method gives a better performance when compared to other approaches and when applied to public expression data, we demonstrated its flexibility to obtain relevant genes biologically associated with normal and abnormal biological processes. Conclusions ProbFAST is a free accessible web-based application that enables MDE analysis on a global scale. It offers an efficient methodological approach for MDE analysis of a set of genes that are turned on and off related to functional information during the evolution of a tumor or tissue differentiation. ProbFAST server can be accessed at http://gdm.fmrp.usp.br/probfast.
Collapse
Affiliation(s)
- Israel T Silva
- Department of Genetics, Faculty of Medicine, University of São Paulo, Ribeirão Preto, Brazil.
| | | | | | | | | |
Collapse
|
43
|
Bullard JH, Purdom E, Hansen KD, Dudoit S. Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments. BMC Bioinformatics 2010; 11:94. [PMID: 20167110 PMCID: PMC2838869 DOI: 10.1186/1471-2105-11-94] [Citation(s) in RCA: 1124] [Impact Index Per Article: 80.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2009] [Accepted: 02/18/2010] [Indexed: 11/21/2022] Open
Abstract
Background High-throughput sequencing technologies, such as the Illumina Genome Analyzer, are powerful new tools for investigating a wide range of biological and medical questions. Statistical and computational methods are key for drawing meaningful and accurate conclusions from the massive and complex datasets generated by the sequencers. We provide a detailed evaluation of statistical methods for normalization and differential expression (DE) analysis of Illumina transcriptome sequencing (mRNA-Seq) data. Results We compare statistical methods for detecting genes that are significantly DE between two types of biological samples and find that there are substantial differences in how the test statistics handle low-count genes. We evaluate how DE results are affected by features of the sequencing platform, such as, varying gene lengths, base-calling calibration method (with and without phi X control lane), and flow-cell/library preparation effects. We investigate the impact of the read count normalization method on DE results and show that the standard approach of scaling by total lane counts (e.g., RPKM) can bias estimates of DE. We propose more general quantile-based normalization procedures and demonstrate an improvement in DE detection. Conclusions Our results have significant practical and methodological implications for the design and analysis of mRNA-Seq experiments. They highlight the importance of appropriate statistical methods for normalization and DE inference, to account for features of the sequencing platform that could impact the accuracy of results. They also reveal the need for further research in the development of statistical and computational methods for mRNA-Seq.
Collapse
Affiliation(s)
- James H Bullard
- Division of Biostatistics, University of California, Berkeley, Berkeley, CA, USA.
| | | | | | | |
Collapse
|
44
|
Zaretzki RL, Gilchrist MA, Briggs WM, Armagan A. Bias correction and Bayesian analysis of aggregate counts in SAGE libraries. BMC Bioinformatics 2010; 11:72. [PMID: 20128916 PMCID: PMC2829012 DOI: 10.1186/1471-2105-11-72] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2009] [Accepted: 02/03/2010] [Indexed: 12/02/2022] Open
Abstract
Background Tag-based techniques, such as SAGE, are commonly used to sample the mRNA pool of an organism's transcriptome. Incomplete digestion during the tag formation process may allow for multiple tags to be generated from a given mRNA transcript. The probability of forming a tag varies with its relative location. As a result, the observed tag counts represent a biased sample of the actual transcript pool. In SAGE this bias can be avoided by ignoring all but the 3' most tag but will discard a large fraction of the observed data. Taking this bias into account should allow more of the available data to be used leading to increased statistical power. Results Three new hierarchical models, which directly embed a model for the variation in tag formation probability, are proposed and their associated Bayesian inference algorithms are developed. These models may be applied to libraries at both the tag and aggregate level. Simulation experiments and analysis of real data are used to contrast the accuracy of the various methods. The consequences of tag formation bias are discussed in the context of testing differential expression. A description is given as to how these algorithms can be applied in that context. Conclusions Several Bayesian inference algorithms that account for tag formation effects are compared with the DPB algorithm providing clear evidence of superior performance. The accuracy of inferences when using a particular non-informative prior is found to depend on the expression level of a given gene. The multivariate nature of the approach easily allows both univariate and joint tests of differential expression. Calculations demonstrate the potential for false positive and negative findings due to variation in tag formation probabilities across samples when testing for differential expression.
Collapse
Affiliation(s)
- Russell L Zaretzki
- Department of Statistics, Operations, and Management Science, The University of Tennessee, 331 Stokely Management Center, Knoxville, TN 37996, USA.
| | | | | | | |
Collapse
|
45
|
Tang D, Zhu Q, Yang F. A Poisson-based adaptive affinity propagation clustering for SAGE data. Comput Biol Chem 2010; 34:63-70. [DOI: 10.1016/j.compbiolchem.2009.11.001] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2009] [Revised: 11/06/2009] [Accepted: 11/13/2009] [Indexed: 11/28/2022]
|
46
|
Quon G, Morris Q. ISOLATE: a computational strategy for identifying the primary origin of cancers using high-throughput sequencing. Bioinformatics 2009; 25:2882-9. [PMID: 19542156 PMCID: PMC2781747 DOI: 10.1093/bioinformatics/btp378] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
Motivation: One of the most deadly cancer diagnoses is the carcinoma of unknown primary origin. Without the knowledge of the site of origin, treatment regimens are limited in their specificity and result in high mortality rates. Though supervised classification methods have been developed to predict the site of origin based on gene expression data, they require large numbers of previously classified tumors for training, in part because they do not account for sample heterogeneity, which limits their application to well-studied cancers. Results: We present ISOLATE, a new statistical method that simultaneously predicts the primary site of origin of cancers and addresses sample heterogeneity, while taking advantage of new high-throughput sequencing technology that promises to bring higher accuracy and reproducibility to gene expression profiling experiments. ISOLATE makes predictions de novo, without having seen any training expression profiles of cancers with identified origin. Compared with previous methods, ISOLATE is able to predict the primary site of origin, de-convolve and remove the effect of sample heterogeneity and identify differentially expressed genes with higher accuracy, across both synthetic and clinical datasets. Methods such as ISOLATE are invaluable tools for clinicians faced with carcinomas of unknown primary origin. Availability: ISOLATE is available for download at: http://morrislab.med.utoronto.ca/software Contact:gerald.quon@utoronto.ca; quaid.morris@utoronto.ca Supplementary information:Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Gerald Quon
- Department of Computer Science, University of Toronto, Toronto, Canada.
| | | |
Collapse
|
47
|
Abstract
BACKGROUND Metagenomics is the study of the genomic content of an environmental sample of microbes. Advances in the through-put and cost-efficiency of sequencing technology is fueling a rapid increase in the number and size of metagenomic datasets being generated. Bioinformatics is faced with the problem of how to handle and analyze these datasets in an efficient and useful way. One goal of these metagenomic studies is to get a basic understanding of the microbial world both surrounding us and within us. One major challenge is how to compare multiple datasets. Furthermore, there is a need for bioinformatics tools that can process many large datasets and are easy to use. RESULTS This article describes two new and helpful techniques for comparing multiple metagenomic datasets. The first is a visualization technique for multiple datasets and the second is a new statistical method for highlighting the differences in a pairwise comparison. We have developed implementations of both methods that are suitable for very large datasets and provide these in Version 3 of our standalone metagenome analysis tool MEGAN. CONCLUSION These new methods are suitable for the visual comparison of many large metagenomes and the statistical comparison of two metagenomes at a time. Nevertheless, more work needs to be done to support the comparative analysis of multiple metagenome datasets. AVAILABILITY Version 3 of MEGAN, which implements all ideas presented in this article, can be obtained from our web site at: www-ab.informatik.uni-tuebingen.de/software/megan. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Suparna Mitra
- Center for Bioinformatics ZBIT, Tübingen University, Sand 14, 72076 Tübingen, Germany.
| | | | | |
Collapse
|
48
|
White JR, Nagarajan N, Pop M. Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput Biol 2009; 5:e1000352. [PMID: 19360128 PMCID: PMC2661018 DOI: 10.1371/journal.pcbi.1000352] [Citation(s) in RCA: 1116] [Impact Index Per Article: 74.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2008] [Accepted: 03/09/2009] [Indexed: 12/18/2022] Open
Abstract
Numerous studies are currently underway to characterize the microbial communities inhabiting our world. These studies aim to dramatically expand our understanding of the microbial biosphere and, more importantly, hope to reveal the secrets of the complex symbiotic relationship between us and our commensal bacterial microflora. An important prerequisite for such discoveries are computational tools that are able to rapidly and accurately compare large datasets generated from complex bacterial communities to identify features that distinguish them.We present a statistical method for comparing clinical metagenomic samples from two treatment populations on the basis of count data (e.g. as obtained through sequencing) to detect differentially abundant features. Our method, Metastats, employs the false discovery rate to improve specificity in high-complexity environments, and separately handles sparsely-sampled features using Fisher's exact test. Under a variety of simulations, we show that Metastats performs well compared to previously used methods, and significantly outperforms other methods for features with sparse counts. We demonstrate the utility of our method on several datasets including a 16S rRNA survey of obese and lean human gut microbiomes, COG functional profiles of infant and mature gut microbiomes, and bacterial and viral metabolic subsystem data inferred from random sequencing of 85 metagenomes. The application of our method to the obesity dataset reveals differences between obese and lean subjects not reported in the original study. For the COG and subsystem datasets, we provide the first statistically rigorous assessment of the differences between these populations. The methods described in this paper are the first to address clinical metagenomic datasets comprising samples from multiple subjects. Our methods are robust across datasets of varied complexity and sampling level. While designed for metagenomic applications, our software can also be applied to digital gene expression studies (e.g. SAGE). A web server implementation of our methods and freely available source code can be found at http://metastats.cbcb.umd.edu/.
Collapse
Affiliation(s)
- James Robert White
- Applied Mathematics and Scientific Computation Program, Center for Bioinformatics and Computational Biology, University of Maryland, College Park, Maryland, United States of America
| | - Niranjan Nagarajan
- Center for Bioinformatics and Computational Biology, University of Maryland, College Park, Maryland, United States of America
| | - Mihai Pop
- Department of Computer Science, Center for Bioinformatics and Computational Biology, University of Maryland, College Park, Maryland, United States of America
- * E-mail:
| |
Collapse
|
49
|
Wang H, Zheng H, Azuaje F. Clustering-based approaches to SAGE data mining. BioData Min 2008; 1:5. [PMID: 18822151 PMCID: PMC2553774 DOI: 10.1186/1756-0381-1-5] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2008] [Accepted: 07/17/2008] [Indexed: 11/12/2022] Open
Abstract
Serial analysis of gene expression (SAGE) is one of the most powerful tools for global gene expression profiling. It has led to several biological discoveries and biomedical applications, such as the prediction of new gene functions and the identification of biomarkers in human cancer research. Clustering techniques have become fundamental approaches in these applications. This paper reviews relevant clustering techniques specifically designed for this type of data. It places an emphasis on current limitations and opportunities in this area for supporting biologically-meaningful data mining and visualisation.
Collapse
Affiliation(s)
- Haiying Wang
- School of Computing and Mathematics, University of Ulster, Newtownabbey, BT37 0QB, Co, Antrim, Northern Ireland, UK.
| | | | | |
Collapse
|
50
|
Schaaf GJ, van Ruissen F, van Kampen A, Kool M, Ruijter JM. Statistical comparison of two or more SAGE libraries: one tag at a time. Methods Mol Biol 2008; 387:151-68. [PMID: 18287630 DOI: 10.1007/978-1-59745-454-4_12] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Several statistical tests have been introduced for the comparison of serial analysis of gene expression (SAGE) libraries to quantitatively analyze the differential expression of genes. As each SAGE library is only one measurement, the necessary information on biological variation or experimental precision is lacking. Therefore, each test includes its own approach to derive such a variance measure from the data set or a theoretical distribution. Because the confidence in tag counts depends on the library size, a test between two or more libraries should be based on original tag counts. When groups of libraries are compared, the test should determine that the proportion of a specific tag in all libraries is the same (null hypothesis), but also offer the possibility to detect specific differences between individual libraries and groups of libraries. The Z-test and the G-test encompass these characteristics and are described for the comparison of two libraries and (two or more) groups of libraries, respectively.
Collapse
Affiliation(s)
- Gerben J Schaaf
- Department of Human Genetics, Academic Medical Center, Amsterdam, The Netherlands
| | | | | | | | | |
Collapse
|