Dataset Information

Non-gaussian distributions affect identification of expression patterns, functional annotation, and prospective classification in human cancer genomes.

ABSTRACT:

Introduction

Gene expression data is often assumed to be normally-distributed, but this assumption has not been tested rigorously. We investigate the distribution of expression data in human cancer genomes and study the implications of deviations from the normal distribution for translational molecular oncology research.

Methods

We conducted a central moments analysis of five cancer genomes and performed empiric distribution fitting to examine the true distribution of expression data both on the complete-experiment and on the individual-gene levels. We used a variety of parametric and nonparametric methods to test the effects of deviations from normality on gene calling, functional annotation, and prospective molecular classification using a sixth cancer genome.

Results

Central moments analyses reveal statistically-significant deviations from normality in all of the analyzed cancer genomes. We observe as much as 37% variability in gene calling, 39% variability in functional annotation, and 30% variability in prospective, molecular tumor subclassification associated with this effect.

Conclusions

Cancer gene expression profiles are not normally-distributed, either on the complete-experiment or on the individual-gene level. Instead, they exhibit complex, heavy-tailed distributions characterized by statistically-significant skewness and kurtosis. The non-Gaussian distribution of this data affects identification of differentially-expressed genes, functional annotation, and prospective molecular classification. These effects may be reduced in some circumstances, although not completely eliminated, by using nonparametric analytics. This analysis highlights two unreliable assumptions of translational cancer gene expression analysis: that "small" departures from normality in the expression data distributions are analytically-insignificant and that "robust" gene-calling algorithms can fully compensate for these effects.

SUBMITTER: Marko NF

PROVIDER: S-EPMC3485292 | biostudies-literature | 2012

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Non-gaussian distributions affect identification of expression patterns, functional annotation, and prospective classification in human cancer genomes.

Marko Nicholas F NF Weil Robert J RJ

PloS one 20121031 10

<h4>Introduction</h4>Gene expression data is often assumed to be normally-distributed, but this assumption has not been tested rigorously. We investigate the distribution of expression data in human cancer genomes and study the implications of deviations from the normal distribution for translational molecular oncology research.<h4>Methods</h4>We conducted a central moments analysis of five cancer genomes and performed empiric distribution fitting to examine the true distribution of expression d ...[more]

PMID: 23118863

Similar Datasets

Project description:Rhodopseudomonas palustris, a nonsulphur purple photosynthetic bacteria, has been extensively investigated for its metabolic versatility including ability to produce hydrogen gas from sunlight and biomass. The availability of the finished genome sequences of six R. palustris strains (BisA53, BisB18, BisB5, CGA009, HaA2 and TIE-1) combined with online bioinformatics software for integrated analysis presents new opportunities to determine the genomic basis of metabolic versatility and ecological lifestyles of the bacteria species. The purpose of this investigation was to compare the functional annotations available for multiple R. palustris genomes to identify annotations that can be further investigated for strain-specific or uniquely shared phenotypic characteristics. A total of 2,355 protein family Pfam domain annotations were clustered based on presence or absence in the six genomes. The clustering process identified groups of functional annotations including those that could be verified as strain-specific or uniquely shared phenotypes. For example, genes encoding water/glycerol transport were present in the genome sequences of strains CGA009 and BisB5, but absent in strains BisA53, BisB18, HaA2 and TIE-1. Protein structural homology modeling predicted that the two orthologous 240 aa R. palustris aquaporins have water-specific transport function. Based on observations in other microbes, the presence of aquaporin in R. palustris strains may improve freeze tolerance in natural conditions of rapid freezing such as nitrogen fixation at low temperatures where access to liquid water is a limiting factor for nitrogenase activation. In the case of adaptive loss of aquaporin genes, strains may be better adapted to survive in conditions of high-sugar content such as fermentation of biomass for biohydrogen production. Finally, web-based resources were developed to allow for interactive, user-defined selection of the relationship between protein family annotations and the R. palustris genomes.

Project description:BackgroundThe genome sequence data of more than 65985 species are publicly available as of October 2021 within the National Center for Biotechnology Information (NCBI) database alone and additional genome sequences are available in other databases and also continue to accumulate at a rapid pace. However, an error-free functional annotation of these genome is essential for the research communities to fully utilize these data in an optimum and efficient manner.ResultsAn analysis of proteome sequence data of 689 fungal species (7.15 million protein sequences) was conducted to identify the presence of functional annotation errors. Proteins associated with calcium signaling events, including calcium dependent protein kinases (CDPKs), calmodulins (CaM), calmodulin-like (CML) proteins, WRKY transcription factors, selenoproteins, and proteins associated with the terpene biosynthesis pathway, were targeted in the analysis. Gene associated with CDPKs and selenoproteins are known to be absent in fungal genomes. Our analysis, however, revealed the presence of proteins that were functionally annotated as CDPK proteins. However, InterproScan analysis indicated that none of the protein sequences annotated as "calcium dependent protein kinase" were found to encode calcium binding EF-hands at the regulatory domain. Similarly, none of a protein sequences annotated as a "selenocysteine" were found to contain a Sec (U) amino acid. Proteins annotated as CaM and CMLs also had significant discrepancies. CaM proteins should contain four calcium binding EF-hands, however, a range of 2-4 calcium binding EF-hands were present in the fungal proteins that were annotated as CaM proteins. Similarly, CMLs should possess four calcium binding EF-hands, but some of the CML annotated fungal proteins possessed either three or four calcium binding EF-hands. WRKY transcription factors are characterized by the presence of a WRKY domain and are confined to the plant kingdom. Several fungal proteins, however, were annotated as WRKY transcription factors, even though they did not contain a WRKY domain.ConclusionThe presence of functional annotation errors in fungal genome and proteome databases is of considerable concern and needs to be addressed in a timely manner.

Project description:In order to generate an atlas of the functional elements driving genome expression in domestic animals, the Functional Annotation of Animal Genome (FAANG) strategy was to sample many tissues from a few animals of different species, sexes, ages, and production stages. This article presents the collection of tissue samples for four species produced by two pilot projects, at INRAE (National Research Institute for Agriculture, Food and Environment) and the University of California, Davis. There were three mammals (cattle, goat, and pig) and one bird (chicken). It describes the metadata characterizing these reference sets (1) for animals with origin and selection history, physiological status, and environmental conditions; (2) for samples with collection site and tissue/cell processing; (3) for quality control; and (4) for storage and further distribution. Three sets are identified: set 1 comprises tissues for which collection can be standardized and for which representative aliquots can be easily distributed (liver, spleen, lung, heart, fat depot, skin, muscle, and peripheral blood mononuclear cells); set 2 comprises tissues requiring special protocols because of their cellular heterogeneity (brain, digestive tract, secretory organs, gonads and gametes, reproductive tract, immune tissues, cartilage); set 3 comprises specific cell preparations (immune cells, tracheal epithelial cells). Dedicated sampling protocols were established and uploaded in https://data.faang.org/protocol/samples. Specificities between mammals and chicken are described when relevant. A total of 73 different tissues or tissue sections were collected, and 21 are common to the four species. Having a common set of tissues will facilitate the transfer of knowledge within and between species and will contribute to decrease animal experimentation. Combining data on the same samples will facilitate data integration. Quality control was performed on some tissues with RNA extraction and RNA quality control. More than 5,000 samples have been stored with unique identifiers, and more than 4,000 were uploaded onto the Biosamples database, provided that standard ontologies were available to describe the sample. Many tissues have already been used to implement FAANG assays, with published results. All samples are available without restriction for further assays. The requesting procedure is described. Members of FAANG are encouraged to apply a range of molecular assays to characterize the functional status of collected samples and share their results, in line with the FAIR (Findable, Accessible, Interoperable, and Reusable) data principles.

Dataset Information

Non-gaussian distributions affect identification of expression patterns, functional annotation, and prospective classification in human cancer genomes.

Introduction

Methods

Results

Conclusions

Publications

Non-gaussian distributions affect identification of expression patterns, functional annotation, and prospective classification in human cancer genomes.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets