Detecting Differentially Expressed Genes Research Articles

Theoretical framework for the difference of two negative binomial distributions and its application in comparative analysis of sequencing data.

High-throughput sequencing (HTS) technologies have been instrumental in investigating biological questions at the bulk and single-cell levels. Comparative analysis of two HTS data sets often relies on testing the statistical significance for the difference of two negative binomial distributions (DOTNB). Although negative binomial distributions are well studied, the theoretical results for DOTNB remain largely unexplored. Here, we derive basic analytical results for DOTNB and examine its asymptotic properties. As a state-of-the-art application of DOTNB, we introduce DEGage, a computational method for detecting differentially expressed genes (DEGs) in scRNA-seq data. DEGage calculates the mean of the sample-wise differences of gene expression levels as the test statistic and determines significant differential expression by computing the P-value with DOTNB. Extensive validation using simulated and real scRNA-seq data sets demonstrates that DEGage outperforms five popular DEG analysis tools: DEGseq2, DEsingle, edgeR, Monocle3, and scDD. DEGage is robust against high dropout levels and exhibits superior sensitivity when applied to balanced and imbalanced data sets, even with small sample sizes. We utilize DEGage to analyze prostate cancer scRNA-seq data sets and identify marker genes for 17 cell types. Furthermore, we apply DEGage to scRNA-seq data sets of mouse neurons with and without fear memory and reveal eight potential memory-related genes overlooked in previous analyses. The theoretical results and supporting software for DOTNB can be widely applied to comparative analyses of dispersed count data in HTS and broad research questions.

Detecting differentially expressed genes from RNA-seq data using fuzzy clustering.

A two-group comparison test is generally performed on RNA sequencing data to detect differentially expressed genes (DEGs). However, the accuracy of this method is low due to the small sample size. To address this, we propose a method using fuzzy clustering that artificially generates data with expression patterns similar to those of DEGs to identify genes that are highly likely to be classified into the same cluster as the initial cluster data. The proposed method is advantageous in that it does not perform any test. Furthermore, a certain level of accuracy can be maintained even when the sample size is biased, and we show that such a situation may improve the accuracy of the proposed method. We compared the proposed method with the conventional method using simulations. In the simulations, we changed the sample size and difference between the expression levels of group 1 and group 2 in the DEGs to obtain the desired accuracy of the proposed method. The results show that the proposed method is superior in all cases under the conditions simulated. We also show that the effect of the difference between group 1 and group 2 on the accuracy is more prominent when the sample size is biased.

Meta-analysis of RNA-seq studies with an adaptive weight and truncation p-value combination test

For RNA-seq studies, an important task is to identify genes that are differentially expressed between groups. Nevertheless, the identification of differentially expressed (DE) genes from a single study can be challenging as the small number of samples and the high sensitivity of data perturbations. Combining p-values from multiple RNA-seq studies can increase the statistical power and the accuracy in detecting DE genes. In this paper, we propose a weight and truncation p-value combination test for meta-analyzing RNA-seq studies. We show that with proper weight and truncation parameters, our new test has a higher statistical power over the existing tests, when the genes are weakly expressed in most studies and the sample size is unbalanced. We then present simulations and apply to two real data analysis. Finally, we provide some practical guidance for our new test and some existing p-value combination tests.

Differentially expressed heterogeneous overdispersion genes testing for count data.

The mRNA-seq data analysis is a powerful technology for inferring information from biological systems of interest. Specifically, the sequenced RNA fragments are aligned with genomic reference sequences, and we count the number of sequence fragments corresponding to each gene for each condition. A gene is identified as differentially expressed (DE) if the difference in its count numbers between conditions is statistically significant. Several statistical analysis methods have been developed to detect DE genes based on RNA-seq data. However, the existing methods could suffer decreasing power to identify DE genes arising from overdispersion and limited sample size, where overdispersion refers to the empirical phenomenon that the variance of read counts is larger than the mean of read counts. We propose a new differential expression analysis procedure: heterogeneous overdispersion genes testing (DEHOGT) based on heterogeneous overdispersion modeling and a post-hoc inference procedure. DEHOGT integrates sample information from all conditions and provides a more flexible and adaptive overdispersion modeling for the RNA-seq read count. DEHOGT adopts a gene-wise estimation scheme to enhance the detection power of differentially expressed genes when the number of replicates is limited as long as the number of conditions is large. DEHOGT is tested on the synthetic RNA-seq read count data and outperforms two popular existing methods, DESeq2 and EdgeR, in detecting DE genes. We apply the proposed method to a test dataset using RNAseq data from microglial cells. DEHOGT tends to detect more differently expressed genes potentially related to microglial cells under different stress hormones treatments.

Transcriptome analysis reveals steroid hormones biosynthesis pathway involved in abdominal fat deposition in broilers

Excessive abdominal fat deposition reduces the feed efficiency and increase the cost of production in broilers. Therefore, it is an important task for poultry breeders to breed broilers with low abdominal fat. Abdominal fat deposition is a highly complex biological process, and its molecular basis remains elusive. In this study, we performed transcriptome analysis to compare gene expression profiles at different stages of abdominal fat deposition to identify the key genes and pathways involved in abdominal fat accumulation. We found that abdominal fat weight (AFW) increased gradually from day 35 (D35) to 91 (D91), and then decreased at day 119 (D119). Accordingly, after detecting differentially expressed genes (DEGs) by comparing gene expression profiles at D35 vs. D63 and D35 vs. D91, and identifying gene modules associated with fat deposition by weighted gene co-expression network analysis (WGCNA), we performed intersection analysis of the detected DEGs and WGCNA gene modules and identified 394 and 435 intersecting genes, respectively. The results of the Gene Ontology (GO) functional enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses showed that the steroid hormone biosynthesis and insulin signaling pathways were co-enriched in all intersecting genes, steroid hormones have been shown that regulated insulin signaling pathway, indicating the importance of the steroid hormone biosynthesis pathway in the development of broiler abdominal fat. We then identified 6 hub genes (ACTB, SOX9, RHOBTB2, PDLIM3, NEDD9, and DOCK4) related to abdominal fat deposition. Further analysis also revealed that there were direct interactions between 6 hub genes. SOX9 has been shown to bind to proteins required for steroid hormone receptor binding, and RHOBTB2 indirectly regulates the steroid hormones biosynthesis through cyclin factor, and ultimately affect fat deposition. Our results suggest that the genes RHOBTB2 and SOX9 play an important role in fat deposition in broilers, by regulating steroid hormone synthesis. These findings provide new targets and directions for further studies on the mechanisms of fat deposition in chicken.

Microarray Expression Profile of Myricetin-Treated THP-1 Macrophages Exhibits Alterations in Atherosclerosis-Related Regulator Molecules and LXR/RXR Pathway.

Atherosclerosis is a chronic inflammation characterized by macrophage infiltration, lipid deposition, and arterial wall thickening. Prevention of atherosclerosis by nutraceuticals is gaining attention. Myricetin, a dietary flavonol, is claimed to possess anti-atherosclerosis properties. We studied myricetin's effect on the atherosclerosis-associated molecular mechanism. Cytotoxicity and proliferation testing to check the viability of myricetin-treated THP-1 macrophages and monocyte migration study in the presence and absence of myricetin was performed. The whole transcriptome analysis was conducted using the Affymetrix microarray platform. The Partek genomics suite for detecting differentially expressed genes (DEGs) and ingenuity pathway analysis was used to identify canonical pathways. Cytotoxicity assays exhibited no significant toxicity in THP-1 macrophages treated with different myricetin concentrations (10-200 μM). Genome-wide expression profiling revealed 58 DEGs (53 upregulated and 5 downregulated) in myricetin-treated THP-1 macrophages. Pathway analysis revealed inhibition of LXR/RXR activation and angiogenesis inhibition by thrombospondin-1 and activated phagocytosis in myricetin-treated THP-1 macrophages. The cytotoxicity assay shows myricetin as a safe phytochemical. In vitro and in silico pathway studies on THP-1 macrophages showed that they can inhibit THP-1 monocyte migration and alter the cholesterol efflux mediated via LXR/RXR signaling. Therefore, myricetin could help in the prevention of cell infiltration in atherosclerotic plaque with reduced risk of stroke or brain damage.

Revealing a novel contributing landscape of ferroptosis-related genes in Parkinson’s disease

Transcriptomics studies have yielded great insights into disease processes by detecting differentially expressed genes (DEGs). In this study, due to the high heritability of Parkinson’s disease (PD), we performed bioinformatics analyses on nine transcriptomic datasets regarding substantia nigra from Gene Expression Omnibus database, including seven microarray datasets and two next-generation sequencing datasets. As a result, between age-matched PD patients and normal control, we identified 630 DEGs, of which 22 hub DEGs involved in PD or ferroptosis were found to be associated with each other at the transcriptional level and protein-protein interaction network, suggesting their high correlations among these hub genes. Moreover, 16 DEGs were singled out due to their comparable AUC (>0.6) in random forest classifiers, including seven PD-related genes (MAP4K4, LRP10, UCHL1, PAM, RIT2, SNCA, GCH1) and nine ferroptosis-related genes (GCH1, DDIT4, RGS4, MAPK9, CAV1, RELA, DUSP1, ATP6V1G2, ATF4 and ISCU). Furthermore, to probe the potential of those hub genes in predicting the PD progression and survival, we constructed a Cox model featured by an eight-gene signature, including four PD-related genes (SNCA, UCHL1, LRP10, and GCH1) and four ferroptosis-related genes (DDIT4, RGS4, RELA, and CAV1), and validated it successful in an independent dataset, indicating that it would be an effective tool for clinical research to predict PD progression. In conclusion, ferroptosis-related DEGs identified in this study were closely correlated with the known PD-related genes, revealing the involvement of ferroptosis in the development of PD. This study presented the potential of several ferroptosis-related genes as novel clinical biomarkers for PD.

Characterization and Duodenal Transcriptome Analysis of Chinese Beef Cattle With Divergent Feed Efficiency Using RNA-Seq.

Residual feed intake (RFI) is an important measure of feed efficiency for agricultural animals. Factors associated with cattle RFI include physiology, dietary factors, and the environment. However, a precise genetic mechanism underlying cattle RFI variations in duodenal tissue is currently unavailable. The present study aimed to identify the key genes and functional pathways contributing to variance in cattle RFI phenotypes using RNA sequencing (RNA-seq). Six bulls with extremely high or low RFIs were selected for detecting differentially expressed genes (DEGs) by RNA-seq, followed by conducting GO, KEGG enrichment, protein-protein interaction (PPI), and co-expression network (WGCNA, n = 10) analysis. A total of 380 differentially expressed genes was obtained from high and low RFI groups, including genes related to energy metabolism (ALDOA, HADHB, INPPL1), mitochondrial function (NDUFS1, RFN4, CUL1), and feed intake behavior (CCK). Two key sub-networks and 26 key genes were detected using GO analysis of DEGs and PPI analysis, such as TPM1 and TPM2, which are involved in mitochondrial pathways and protein synthesis. Through WGCNA, a gene network was built, and genes were sorted into 27 modules, among which the blue (r = 0.72, p = 0.03) and salmon modules (r = −0.87, p = 0.002) were most closely related with RFI. DEGs and genes from the main sub-networks and closely related modules were largely involved in metabolism; oxidative phosphorylation; glucagon, ribosome, and N-glycan biosynthesis, and the MAPK and PI3K-Akt signaling pathways. Through WGCNA, five key genes, including FN1 and TPM2, associated with the biological regulation of oxidative processes and skeletal muscle development were identified. Taken together, our data suggest that the duodenum has specific biological functions in regulating feed intake. Our findings provide broad-scale perspectives for identifying potential pathways and key genes involved in the regulation of feed efficiency in beef cattle.

A scaling-free minimum enclosing ball method to detect differentially expressed genes for RNA-seq data

BackgroundIdentifying differentially expressed genes between the same or different species is an urgent demand for biological and medical research. For RNA-seq data, systematic technical effects and different sequencing depths are usually encountered when conducting experiments. Normalization is regarded as an essential step in the discovery of biologically important changes in expression. The present methods usually involve normalization of the data with a scaling factor, followed by detection of significant genes. However, more than one scaling factor may exist because of the complexity of real data. Consequently, methods that normalize data by a single scaling factor may deliver suboptimal performance or may not even work.The development of modern machine learning techniques has provided a new perspective regarding discrimination between differentially expressed (DE) and non-DE genes. However, in reality, the non-DE genes comprise only a small set and may contain housekeeping genes (in same species) or conserved orthologous genes (in different species). Therefore, the process of detecting DE genes can be formulated as a one-class classification problem, where only non-DE genes are observed, while DE genes are completely absent from the training data.ResultsIn this study, we transform the problem to an outlier detection problem by treating DE genes as outliers, and we propose a scaling-free minimum enclosing ball (SFMEB) method to construct a smallest possible ball to contain the known non-DE genes in a feature space. The genes outside the minimum enclosing ball can then be naturally considered to be DE genes. Compared with the existing methods, the proposed SFMEB method does not require data normalization, which is particularly attractive when the RNA-seq data include more than one scaling factor. Furthermore, the SFMEB method could be easily extended to different species without normalization.ConclusionsSimulation studies demonstrate that the SFMEB method works well in a wide range of settings, especially when the data are heterogeneous or biological replicates. Analysis of the real data also supports the conclusion that the SFMEB method outperforms other existing competitors. The R package of the proposed method is available at https://bioconductor.org/packages/MEB.

RNAflow: An Effective and Simple RNA-Seq Differential Gene Expression Pipeline Using Nextflow.

RNA-Seq enables the identification and quantification of RNA molecules, often with the aim of detecting differentially expressed genes (DEGs). Although RNA-Seq evolved into a standard technique, there is no universal gold standard for these data’s computational analysis. On top of that, previous studies proved the irreproducibility of RNA-Seq studies. Here, we present a portable, scalable, and parallelizable Nextflow RNA-Seq pipeline to detect DEGs, which assures a high level of reproducibility. The pipeline automatically takes care of common pitfalls, such as ribosomal RNA removal and low abundance gene filtering. Apart from various visualizations for the DEG results, we incorporated downstream pathway analysis for common species as Homo sapiens and Mus musculus. We evaluated the DEG detection functionality while using qRT-PCR data serving as a reference and observed a very high correlation of the logarithmized gene expression fold changes.

Integrative Analysis of Axolotl Gene Expression Data from Regenerative and Wound Healing Limb Tissues

Axolotl (Ambystoma mexicanum) is a urodele amphibian endowed with remarkable regenerative capacities manifested in scarless wound healing and restoration of amputated limbs, which makes it a powerful experimental model for regenerative biology and medicine. Previous studies have utilized microarrays and RNA-Seq technologies for detecting differentially expressed (DE) genes in different phases of the axolotl limb regeneration. However, sufficient consistency may be lacking due to statistical limitations arising from intra-laboratory analyses. This study aims to bridge such gaps by performing an integrative analysis of publicly available microarray and RNA-Seq data from axolotl limb samples having comparable study designs using the “merging” method. A total of 351 genes were found DE in regenerative samples compared to the control in data of both technologies, showing an adjusted p-value < 0.01 and log fold change magnitudes >1. Downstream analyses illustrated consistent correlations of the directionality of DE genes within and between data of both technologies, as well as concordance with the literature on regeneration related biological processes. qRT-PCR analysis validated the observed expression level differences of five of the top DE genes. Future studies may benefit from the utilized concept and approach for enhanced statistical power and robust discovery of biomarkers of regeneration.

Detection of differentially expressed genes in discrete single-cell RNA sequencing data using a hurdle model with correlated random effects.

Single-cell RNA sequencing (scRNA-seq) technologies are revolutionary tools allowing researchers to examine gene expression at the level of a single cell. Traditionally, transcriptomic data have been analyzed from bulk samples, masking the heterogeneity now seen across individual cells. Even within the same cellular population, genes can be highly expressed in some cells but not expressed (or lowly expressed) in others. Therefore, the computational approaches used to analyze bulk RNA sequencing data are not appropriate for the analysis of scRNA-seq data. Here, we present a novel statistical model for high dimensional and zero-inflated scRNA-seq count data to identify differentially expressed (DE) genes across cell types. Correlated random effects are employed based on an initial clustering of cells to capture the cell-to-cell variability within treatment groups. Moreover, this model is flexible and can be easily adapted to an independent random effect structure if needed. We apply our proposed methodology to both simulated and real data and compare results to other popular methods designed for detecting DE genes. Due to the hurdle model's ability to detect differences in the proportion of cells expressed and the average expression level (among the expressed cells), our methods naturally identify some genes as DE that other methods do not, and we demonstrate with real data that these uniquely detected genes are associated with similar biological processes and functions.

Comparative analysis of differential gene expression analysis tools for single-cell RNA sequencing data

BackgroundThe analysis of single-cell RNA sequencing (scRNAseq) data plays an important role in understanding the intrinsic and extrinsic cellular processes in biological and biomedical research. One significant effort in this area is the detection of differentially expressed (DE) genes. scRNAseq data, however, are highly heterogeneous and have a large number of zero counts, which introduces challenges in detecting DE genes. Addressing these challenges requires employing new approaches beyond the conventional ones, which are based on a nonzero difference in average expression. Several methods have been developed for differential gene expression analysis of scRNAseq data. To provide guidance on choosing an appropriate tool or developing a new one, it is necessary to evaluate and compare the performance of differential gene expression analysis methods for scRNAseq data.ResultsIn this study, we conducted a comprehensive evaluation of the performance of eleven differential gene expression analysis software tools, which are designed for scRNAseq data or can be applied to them. We used simulated and real data to evaluate the accuracy and precision of detection. Using simulated data, we investigated the effect of sample size on the detection accuracy of the tools. Using real data, we examined the agreement among the tools in identifying DE genes, the run time of the tools, and the biological relevance of the detected DE genes.ConclusionsIn general, agreement among the tools in calling DE genes is not high. There is a trade-off between true-positive rates and the precision of calling DE genes. Methods with higher true positive rates tend to show low precision due to their introducing false positives, whereas methods with high precision show low true positive rates due to identifying few DE genes. We observed that current methods designed for scRNAseq data do not tend to show better performance compared to methods designed for bulk RNAseq data. Data multimodality and abundance of zero read counts are the main characteristics of scRNAseq data, which play important roles in the performance of differential gene expression analysis methods and need to be considered in terms of the development of new methods.

Classical and Bayesian random-effects meta-analysis models with sample quality weights in gene expression studies

BackgroundRandom-effects (RE) models are commonly applied to account for heterogeneity in effect sizes in gene expression meta-analysis. The degree of heterogeneity may differ due to inconsistencies in sample quality. High heterogeneity can arise in meta-analyses containing poor quality samples. We applied sample-quality weights to adjust the study heterogeneity in the DerSimonian and Laird (DSL) and two-step DSL (DSLR2) RE models and the Bayesian random-effects (BRE) models with unweighted and weighted data, Gibbs and Metropolis-Hasting (MH) sampling algorithms, weighted common effect, and weighted between-study variance. We evaluated the performance of the models through simulations and illustrated application of the methods using Alzheimer’s gene expression datasets.ResultsSample quality adjusting within study variance (wP6) models provided an appropriate reduction of differentially expressed (DE) genes compared to other weighted functions in classical RE models. The BRE model with a uniform(0,1) prior was appropriate for detecting DE genes as compared to the models with other prior distributions. The precision of DE gene detection in the heterogeneous data was increased with the DSLR2wP6 weighted model compared to the DSLwP6 weighted model. Among the BRE weighted models, the wP6weighted- and unweighted-data models and both Gibbs- and MH-based models performed similarly. The wP6 weighted common-effect model performed similarly to the unweighted model in the homogeneous data, but performed worse in the heterogeneous data. The wP6weighted data were appropriate for detecting DE genes with high precision, while the wP6weighted between-study variance models were appropriate for detecting DE genes with high overall accuracy. Without the weight, when the number of genes in microarray increased, the DSLR2 performed stably, while the overall accuracy of the BRE model was reduced. When applying the weighted models in the Alzheimer’s gene expression data, the number of DE genes decreased in all metadata sets with the DSLR2wP6weighted and the wP6weighted between study variance models. Four hundred and forty-six DE genes identified by the wP6weighted between study variance model could be potentially down-regulated genes that may contribute to good classification of Alzheimer’s samples.ConclusionsThe application of sample quality weights can increase precision and accuracy of the classical RE and BRE models; however, the performance of the models varied depending on data features, levels of sample quality, and adjustment of parameter estimates.

DiPhiSeq: robust comparison of expression levels on RNA-Seq data with large sample sizes.

In the analysis of RNA-Seq data, detecting differentially expressed (DE) genes has been a hot research area in recent years and many methods have been proposed. DE genes show different average expression levels in different sample groups, and thus can be important biological markers. While generally very successful, these methods need to be further tailored and improved for cancerous data, which often features quite diverse expression in the samples from the cancer group, and this diversity is much larger than that in the control group. We propose a statistical method that can detect not only genes that show different average expressions, but also genes that show different diversities of expressions in different groups. These 'differentially dispersed' genes can be important clinical markers. Our method uses a redescending penalty on the quasi-likelihood function, and thus has superior robustness against outliers and other noise. Simulations and real data analysis demonstrate that DiPhiSeq outperforms existing methods in the presence of outliers, and identifies unique sets of genes. DiPhiSeq is publicly available as an R package on CRAN: https://cran.r-project.org/package=DiPhiSeq. Supplementary data are available at Bioinformatics online.

SigEMD: A powerful method for differential gene expression analysis in single-cell RNA sequencing data.

Differential gene expression analysis is one of the significant efforts in single cell RNA sequencing (scRNAseq) analysis to discover the specific changes in expression levels of individual cell types. Since scRNAseq exhibits multimodality, large amounts of zero counts, and sparsity, it is different from the traditional bulk RNA sequencing (RNAseq) data. The new challenges of scRNAseq data promote the development of new methods for identifying differentially expressed (DE) genes. In this study, we proposed a new method, SigEMD, that combines a data imputation approach, a logistic regression model and a nonparametric method based on the Earth Mover's Distance, to precisely and efficiently identify DE genes in scRNAseq data. The regression model and data imputation are used to reduce the impact of large amounts of zero counts, and the nonparametric method is used to improve the sensitivity of detecting DE genes from multimodal scRNAseq data. By additionally employing gene interaction network information to adjust the final states of DE genes, we further reduce the false positives of calling DE genes. We used simulated datasets and real datasets to evaluate the detection accuracy of the proposed method and to compare its performance with those of other differential expression analysis methods. Results indicate that the proposed method has an overall powerful performance in terms of precision in detection, sensitivity, and specificity.

A novel joint analysis framework improves identification of differentially expressed genes in cross disease transcriptomic analysis

MotivationDetecting differentially expressed (DE) genes between disease and normal control group is one of the most common analyses in genome-wide transcriptomic data. Since most studies don’t have a lot of samples, researchers have used meta-analysis to group different datasets for the same disease. Even then, in many cases the statistical power is still not enough. Taking into account the fact that many diseases share the same disease genes, it is desirable to design a statistical framework that can identify diseases’ common and specific DE genes simultaneously to improve the identification power.ResultsWe developed a novel empirical Bayes based mixture model to identify DE genes in specific study by leveraging the shared information across multiple different disease expression data sets. The effectiveness of joint analysis was demonstrated through comprehensive simulation studies and two real data applications. The simulation results showed that our method consistently outperformed single data set analysis and two other meta-analysis methods in identification power. In real data analysis, overall our method demonstrated better identification power in detecting DE genes and prioritized more disease related genes and disease related pathways than single data set analysis. Over 150% more disease related genes are identified by our method in application to Huntington’s disease. We expect that our method would provide researchers a new way of utilizing available data sets from different diseases when sample size of the focused disease is limited.

SRC and TP53 play critical role in low-grade dysplasia colorectal mucosa transformation into cancer.

Determination of crucial genes that are involved in onset and progress of dysplasia of colorectal mucosa is the aim of this study. Management of dysplasia as one of the risk factors of colon cancer is very challenging. Molecular studies could be helpful in this matter. Here, the transcriptome profile of low-grade dysplasia in colon tissue in comparison with normal one is studied by protein-protein interaction (PPI) network analysis. For detecting differentially expressed genes (DEGs) of dysplasia lesion, the data was downloaded from the gene chip GSE31106, platform GPL1261, GSM770092-94 as normal colorectal mucosa group and GSM770098-100 as low-grade dysplasia colorectal mucosa from the Gene Expression Omnibus database (GEO). The expression profile is evaluated by GEO2R and a network of DEGs is constructed and analyzed by Cytoscape algorithms. The findings indicate that a PPI network analysis of 113 DEGs is consist of 8 nodes that 6 of them are common with inflammation state. Only SRC and TP53 were recognized as the specific makers for dysplasia. In this respect, a subnetwork of these two genes introduce a panel of 8 nodes consist of HRAS, MYC, PIK3CA, PIK3CB, PIK3CD, PIK3CG, SRC, and TP53. It can be concluded that SRC and TP53 may play prominent role in dysplasia pathogenicity after running validation tests.

XBSeq2: a fast and accurate quantification of differential expression and differential polyadenylation

BackgroundRNA sequencing (RNA-seq) is a high throughput technology that profiles gene expression in a genome-wide manner. RNA-seq has been mainly used for testing differential expression (DE) of transcripts between two conditions and has recently been used for testing differential alternative polyadenylation (APA). In the past, many algorithms have been developed for detecting differentially expressed genes (DEGs) from RNA-seq experiments, including the one we developed, XBSeq, which paid special attention to the context-specific background noise that is ignored in conventional gene expression quantification and DE analysis of RNA-seq data.ResultsWe present several major updates in XBSeq2, including alternative statistical testing and parameter estimation method for detecting DEGs, capacity to directly process alignment files and methods for testing differential APA usage. We evaluated the performance of XBSeq2 against several other methods by using simulated datasets in terms of area under the receiver operating characteristic (ROC) curve (AUC), number of false discoveries and statistical power. We also benchmarked different methods concerning execution time and computational memory consumed. Finally, we demonstrated the functionality of XBSeq2 by using a set of in-house generated clear cell renal carcinoma (ccRCC) samples.ConclusionsWe present several major updates to XBSeq. By using simulated datasets, we demonstrated that, overall, XBSeq2 performs equally well as XBSeq in terms of several statistical metrics and both perform better than DESeq2 and edgeR. In addition, XBSeq2 is faster in speed and consumes much less computational memory compared to XBSeq, allowing users to evaluate differential expression and APA events in parallel. XBSeq2 is available from Bioconductor: http://bioconductor.org/packages/XBSeq/

Powerful differential expression analysis incorporating network topology for next-generation sequencing data.

RNA-seq has become the technology of choice for interrogating the transcriptome. However, most methods for RNA-seq differential expression (DE) analysis do not utilize prior knowledge of biological networks to detect DE genes. With the increased availability and quality of biological network databases, methods that can utilize this prior knowledge are needed and will offer biologists with a viable, more powerful alternative when analyzing RNA-seq data. We propose a three-state Markov Random Field (MRF) method that utilizes known biological pathways and interaction to improve sensitivity and specificity and therefore reducing false discovery rates (FDRs) when detecting differentially expressed genes from RNA-seq data. The method requires normalized count data (e.g. in Fragments or Reads Per Kilobase of transcript per Million mapped reads (FPKM/RPKM) format) as its input and it is implemented in an R package pathDESeq available from Github. Simulation studies demonstrate that our method outperforms the two-state MRF model for various sample sizes. Furthermore, for a comparable FDR, it has better sensitivity than DESeq, EBSeq, edgeR and NOISeq. The proposed method also picks more top Gene Ontology terms and KEGG pathways terms when applied to real dataset from colorectal cancer and hepatocellular carcinoma studies, respectively. Overall, these findings clearly highlight the power of our method relative to the existing methods that do not utilize prior knowledge of biological network. As an R package at https://github.com/MalathiSIDona/pathDESeq. install_github("MalathiSIDona/pathDESeq",build_vignettes = TRUE). After installation, type vignette("pathDESeq") to access the vignette. a.salim@latrobe.edu.au. Supplementary data are available at Bioinformatics online.