Search

Mutation density patterns reveal unique biological properties of specific genomic regions and shed light on the mechanisms of carcinogenesis. Although previous studies reported insightful mutation density patterns associated with certain genomic regions such as transcription start sites and DNA replication origins, a tool that can systematically investigate mutational spatial patterns is still lacking. Thus, we developed MutDens, a bioinformatic tool for comprehensive analysis of mutation density patterns around genomic features, namely, genomic positions, in humans and model species. By scanning the bidirectional vicinity regions of given positions, MutDens systematically characterizes the mutation density for single-base substitution mutational classes after adjusting for total mutation burden and local nucleotide proportion. Analysis results using MutDens not only verified the previously reported transcriptional strand bias around transcription start sites and replicative strand bias around DNA replication origins, but also identified novel mutation density patterns around other genomics features, such as enhancers and retrotransposon insertion polymorphism sites. To our knowledge, MutDens is the first tool that systematically calculates, examines, and compares mutation density patterns, thus providing a valuable avenue for investigating the mutational landscapes associated with important genomic features.

Background: The metabolome is the best representation of cancer phenotypes. Gene expression can be considered a confounding covariate affecting metabolite levels. Data integration across metabolomics and genomics to establish the biological relevance of cancer metabolism is challenging. This study aimed to eliminate the confounding effect of metabolic gene expression to reflect actual metabolite levels in microsatellite instability (MSI) cancers. Methods: In this study, we propose a new strategy using covariate-adjusted tensor classification in high dimensions (CATCH) models to integrate metabolite and metabolic gene expression data to classify MSI and microsatellite stability (MSS) cancers. We used datasets from the Cancer Cell Line Encyclopedia (CCLE) phase II project and treated metabolomic data as tensor predictors and data on gene expression of metabolic enzymes as confounding covariates. Results: The CATCH model performed well, with high accuracy (0.82), sensitivity (0.66), specificity (0.88), precision (0.65), and F1 score (0.65). Seven metabolite features adjusted for metabolic gene expression, namely, 3-phosphoglycerate, 6-phosphogluconate, cholesterol ester, lysophosphatidylethanolamine (LPE), phosphatidylcholine, reduced glutathione, and sarcosine, were found in MSI cancers. Only one metabolite, Hippurate, was present in MSS cancers. The gene expression of phosphofructokinase 1 (PFKP), which is involved in the glycolytic pathway, was related to 3-phosphoglycerate. ALDH4A1and GPT2were associated with sarcosine. LPE was associated with the expression of CHPT1, which is involved in lipid metabolism. The glycolysis, nucleotide, glutamate, and lipid metabolic pathways were enriched in MSI cancers. Conclusions: We propose an effective CATCH model for predicting MSI cancer status. By controlling the confounding effect of metabolic gene expression, we identified cancer metabolic biomarkers and therapeutic targets. In addition, we provided the possible biology and genetics of MSI cancer metabolism.

Next generation sequencing technology has matured, and with its current affordability, will replace the SNP chip as the genotyping tool of choice. Even with the current affordability of NGS, large scale studies will require careful study design to reduce cost. In this study, we designed an experiment to assess the accuracy of allele frequency estimated from pooled sequencing data. We compared the allele frequency estimated from sequencing data with the allele frequency estimated from individual SNP chip data and observed high correlations between them. However, by calculating error rate, we found that many SNPs had their allele frequency estimated from sequencing data significantly different from allele frequency estimated from SNP chip data. In conclusion, we found correlation is not an ideal measurement for comparing allele frequencies. And for the purpose of estimating allele frequency, we do not recommend using pooling with NGS as a cheaper alternative to genotype each sample individually.

Sample size determination is an important issue in the experimental design of biomedical research. Because of the complexity of RNA-seq experiments, however, the field currently lacks a sample size method widely applicable to differential expression studies utilising RNA-seq technology. In this report, we propose several methods for sample size calculation for single-gene differential expression analysis of RNA-seq data under Poisson distribution. These methods are then extended to multiple genes, with consideration for addressing the multiple testing problem by controlling false discovery rate. Moreover, most of the proposed methods allow for closed-form sample size formulas with specification of the desired minimum fold change and minimum average read count, and thus are not computationally intensive. Simulation studies to evaluate the performance of the proposed sample size formulas are presented; the results indicate that our methods work well, with achievement of desired power. Finally, our sample size calculation methods are applied to three real RNA-seq data sets.