Search

Alzheimer's disease (AD) is the most common neurodegenerative disease, and it is currently untreatable. RNA sequencing (RNA‐Seq) is commonly used in the literature to identify AD‐associated molecular mechanisms by analysing changes in gene expression. RNA‐Seq data can also be used to detect genomic variants, enabling the identification of the genes with a higher load of deleterious variants in patients compared with controls. Here, we analysed AD RNA‐Seq datasets to obtain differentially expressed genes and genes with a higher load of pathogenic variants in AD, and we combined them in a single list. We mapped these genes on a human protein–protein interaction network to discover subnetworks perturbed by AD. Our results show that utilizing gene pathogenicity information from RNA‐Seq data positively contributes to the disclosure of AD‐related mechanisms. Moreover, dividing the discovered subnetworks into highly connected modules reveals a clearer picture of altered molecular pathways that, otherwise, would not be captured. Repeating the whole pipeline with human metabolic network genes led to results confirming the positive contribution of gene pathogenicity information and enabled a more detailed identification of altered metabolic pathways in AD. [ABSTRACT FROM AUTHOR]

Simple Summary: Bipolar disorder is a psychiatric condition in which one prominent symptom is the regular occurrence of mood swings. Its root cause remains unclear; thus, this study was intended to provide new insights to better understand the genetics underlying such a disorder. Using samples obtained from three regions in the brain, it was found that some molecular and cellular mechanisms are involved in the disruption of neurobiological processes. The dysregulated expression of certain genes eventually leads to neurogenesis and neurotransmission impairment events in humans. Furthermore, these genes significant in the onset of bipolar disorder were identified and evaluated for the presence of variants, which may be targeted to engineer better curative treatment strategies for the disorder. Bipolar disorder (BPD) is a serious psychiatric condition that is characterized by the frequent shifting of mood patterns, ranging from manic to depressive episodes. Although there are already treatment strategies that aim at regulating the manifestations of this disorder, its etiology remains unclear and continues to be a question of interest within the scientific community. The development of RNA sequencing techniques has provided newer and better approaches to studying disorders at the transcriptomic level. Hence, using RNA-seq data, we employed intramodular connectivity analysis and network pharmacology assessment of disease-associated variants to elucidate the biological pathways underlying the complex nature of BPD. This study was intended to characterize the expression profiles obtained from three regions in the brain, which are the nucleus accumbens (nAcc), the anterior cingulate cortex (AnCg), and the dorsolateral prefrontal cortex (DLPFC), provide insights into the specific roles of these regions in the onset of the disorder, and present potential targets for drug design and development. The nAcc was found to be highly associated with genes responsible for the deregulated transcription of neurotransmitters, while the DLPFC was greatly correlated with genes involved in the impairment of components crucial in neurotransmission. The AnCg did show association with some of the expressions, but the relationship was not as strong as the other two regions. Furthermore, disease-associated variants or single nucleotide polymorphisms (SNPs) were identified among the significant genes in BPD, which suggests the genetic interrelatedness of such a disorder and other mental illnesses. DRD2, GFRA2, and DCBLD1 were the genes with disease-associated variants expressed in the nAcc; ST8SIA2 and ADAMTS16 were the genes with disease-associated variants expressed in the AnCg; and FOXO3, ITGA9, CUBN, PLCB4, and RORB were the genes with disease-associated variants expressed in the DLPFC. Aside from unraveling the molecular and cellular mechanisms behind the expression of BPD, this investigation was envisioned to propose a new research pipeline in studying the transcriptome of psychiatric disorders to support and improve existing studies. [ABSTRACT FROM AUTHOR]

Purpose: To perform an extensive exploratory analysis to build a deeper insight into clinically relevant molecular biomarkers in Papillary, Follicular, and Anaplastic thyroid carcinomas (PTC, FTC, ATC). Methods: Thirteen Thyroid Cancer (THCA) datasets incorporating PTC, FTC, and ATC were derived from the Gene Expression Omnibus. Genes differentially expressed (DEGs) between THCA and normal were identified and subjected to GO and KEGG analyses. Multiple topological properties were harnessed and protein-protein interaction (PPI) networks were constructed to identify the hub genes followed by survival analysis and validation. Results: There were 70, 87, and 377 DEGs, and 23, 27, and 53 hub genes for PTC, FTC, and ATC samples, respectively. Survival analysis detected 39 overall and 49 relapse-free survival-relevant hub genes. Six hub genes, BCL2, FN1, ITPR1, LYVE1, NTRK2, TBC1D4, were found common to more than one THCA type. The most significant hub genes found in the study were: BCL2, CD44, DCN, FN1, IRS1, ITPR1, MFAP4, MKI67, NTRK2, PCLO, TGFA. The most enriched and significant GO terms were Melanocyte differentiation for PTC, Extracellular region for FTC, and Extracellular exosome for ATC. Prostate cancer for PTC was the most significantly enriched KEGG pathway. The results were validated using TCGA data. Conclusions: The findings unravel potential biomarkers and therapeutic targets of thyroid carcinomas. [ABSTRACT FROM AUTHOR]

Heat shock proteins (HSPs) are molecular chaperones that play essential roles in plant development and in response to various environmental stresses. Understanding R. delavayi HSP genes is of great importance since R. delavayi is severely affected by heat stress. In the present study, a total of 76 RdHSP genes were identified in the R. delavayi genome, which were divided into five subfamilies based on molecular weight and domain composition. Analyses of the chromosome distribution, gene structure, and conserved motif of the RdHSP family genes were conducted using bioinformatics analysis methods. Gene duplication analysis showed that 15 and 8 RdHSP genes were obtained and retained from the WGD/segmental duplication and tandem duplication, respectively. Cis-element analysis revealed the importance of RdHSP genes in plant adaptations to the environment. Moreover, the expression patterns of RdHSP family genes were investigated in R. delavayi treated with high temperature based on our RNA-seq data, which were further verified by qRT-PCR. Further analysis revealed that nine candidate genes, including six RdHSP20 subfamily genes (RdHSP20.4, RdHSP20.8, RdHSP20.6, RdHSP20.3, RdHSP20.10, and RdHSP20.15) and three RdHSP70 subfamily genes (RdHSP70.15, RdHSP70.21, and RdHSP70.16), might be involved in enhancing the heat stress tolerance. The subcellular localization of two candidate RdHSP genes (RdHSP20.8 and RdHSP20.6) showed that two candidate RdHSPs were expressed and function in the chloroplast and nucleus, respectively. These results provide a basis for the functional characterization of HSP genes and investigations on the molecular mechanisms of heat stress response in R. delavayi. [ABSTRACT FROM AUTHOR]

Background: Current RNA-seq analysis software for RNA-seq data tends to use similar parameters across different species without considering species-specific differences. However, the suitability and accuracy of these tools may vary when analyzing data from different species, such as humans, animals, plants, fungi, and bacteria. For most laboratory researchers lacking a background in information science, determining how to construct an analysis workflow that meets their specific needs from the array of complex analytical tools available poses a significant challenge. Results: By utilizing RNA-seq data from plants, animals, and fungi, it was observed that different analytical tools demonstrate some variations in performance when applied to different species. A comprehensive experiment was conducted specifically for analyzing plant pathogenic fungal data, focusing on differential gene analysis as the ultimate goal. In this study, 288 pipelines using different tools were applied to analyze five fungal RNA-seq datasets, and the performance of their results was evaluated based on simulation. This led to the establishment of a relatively universal and superior fungal RNA-seq analysis pipeline that can serve as a reference, and certain standards for selecting analysis tools were derived for reference. Additionally, we compared various tools for alternative splicing analysis. The results based on simulated data indicated that rMATS remained the optimal choice, although consideration could be given to supplementing with tools such as SpliceWiz. Conclusion: The experimental results demonstrate that, in comparison to the default software parameter configurations, the analysis combination results after tuning can provide more accurate biological insights. It is beneficial to carefully select suitable analysis software based on the data, rather than indiscriminately choosing tools, in order to achieve high-quality analysis results more efficiently. [ABSTRACT FROM AUTHOR]

Introduction: Our study undertakes a detailed exploration of gene expression dynamics within human lung organ tissue equivalents (OTEs) in response to Influenza A virus (IAV), Human metapneumovirus (MPV), and Parainfluenza virus type 3 (PIV3) infections. Through the analysis of RNA-Seq data from 19,671 genes, we aim to identify differentially expressed genes under various infection conditions, elucidating the complexities of virus-host interactions. Methods: We employ Generalized Linear Models (GLMs) with Quasi-Likelihood (QL) F-tests (GLMQL) and introduce the novel Magnitude-Altitude Score (MAS) and Relaxed Magnitude-Altitude Score (RMAS) algorithms to navigate the intricate landscape of RNA-Seq data. This approach facilitates the precise identification of potential biomarkers, highlighting the host's reliance on innate immune mechanisms. Our comprehensive methodological framework includes RNA extraction, library preparation, sequencing, and Gene Ontology (GO) enrichment analysis to interpret the biological significance of our findings. Results: The differential expression analysis unveils significant changes in gene expression triggered by IAV, MPV, and PIV3 infections. The MAS and RMAS algorithms enable focused identification of biomarkers, revealing a consistent activation of interferon-stimulated genes (e.g., IFIT1, IFIT2, IFIT3, OAS1) across all viruses. Our GO analysis provides deep insights into the host's defense mechanisms and viral strategies exploiting host cellular functions. Notably, changes in cellular structures, such as cilium assembly and mitochondrial ribosome assembly, indicate a strategic shift in cellular priorities. The precision of our methodology is validated by a 92% mean accuracy in classifying respiratory virus infections using multinomial logistic regression, demonstrating the superior efficacy of our approach over traditional methods. Discussion: This study highlights the intricate interplay between viral infections and host gene expression, underscoring the need for targeted therapeutic interventions. The stability and reliability of the MAS/RMAS ranking method, even under stringent statistical corrections, and the critical importance of adequate sample size for biomarker reliability are significant findings. Our comprehensive analysis not only advances our understanding of the host's response to viral infections but also sets a new benchmark for the identification of biomarkers, paving the way for the development of effective diagnostic and therapeutic strategies. [ABSTRACT FROM AUTHOR]

Although the pathogenesis of solar lentigo (SL) involves chronic ultraviolet (UV) exposure, cellular senescence, and upregulated melanogenesis, underlying molecular-level mechanisms associated with SL remain unclear. The aim of this study was to investigate the gene regulatory mechanisms intimately linked to inflammation in SL. Skin samples from patients with SL with or without histological inflammatory features were obtained. RNA-seq data from the samples were analyzed via multiple analysis approaches, including exploration of core inflammatory gene alterations, identifying functional pathways at both transcription and protein levels, comparison of inflammatory module (gene clusters) activation levels, and analyzing correlations between modules. These analyses disclosed specific core genes implicated in oxidative stress, especially the upregulation of nuclear factor kappa B in the inflammatory SLs, while genes associated with protective mechanisms, such as SLC6A9, were highly expressed in the non-inflammatory SLs. For inflammatory modules, Extracellular Immunity and Mitochondrial Innate Immunity were exclusively upregulated in the inflammatory SL. Analysis of protein–protein interactions revealed the significance of CXCR3 upregulation in the pathogenesis of inflammatory SL. In conclusion, the upregulation of stress response-associated genes and inflammatory pathways in response to UV-induced oxidative stress implies their involvement in the pathogenesis of inflammatory SL. [ABSTRACT FROM AUTHOR]

Background: Advances in sequencing technology and cost reduction have enabled an emergence of various statistical methods used in RNA-sequencing data, including the differential co-expression network analysis (or differential network analysis). A key benefit of this method is that it takes into consideration the interactions between or among genes and do not require an established knowledge in biological pathways. As of now, none of existing softwares can incorporate covariates that should be adjusted if they are confounding factors while performing the differential network analysis. Results: We develop an R package PRANA which a user can easily include multiple covariates. The main R function in this package leverages a novel pseudo-value regression approach for a differential network analysis in RNA-sequencing data. This software is also enclosed with complementary R functions for extracting adjusted p-values and coefficient estimates of all or specific variable for each gene, as well as for identifying the names of genes that are differentially connected (DC, hereafter) between subjects under biologically different conditions from the output. Conclusion: Herewith, we demonstrate the application of this package in a real data on chronic obstructive pulmonary disease. PRANA is available through the CRAN repositories under the GPL-3 license: https://cran.r-project.org/web/packages/PRANA/index.html. [ABSTRACT FROM AUTHOR]

Background. Pathological conditions may result in certain genes having expression variance that differs markedly from that of the control. Finding such genes from gene expression data can provide invaluable candidates for therapeutic intervention. Under the dominant paradigm for modeling RNA-Seq gene counts using the negative binomial model, tests of differential variability are challenging to develop, owing to dependence of the variance on the mean. Methods. Here, we describe clrDV, a statistical method for detecting genes that show differential variability between two populations. We present the skew-normal distribution for modeling gene-wise null distribution of centered log-ratio transformation of compositional RNA-seq data. Results. Simulation results show that clrDV has false discovery rate and probability of Type II error that are on par with or superior to existing methodologies. In addition, its run time is faster than its closest competitors, and remains relatively constant for increasing sample size per group. Analysis of a large neurodegenerative disease RNASeq dataset using clrDV successfully recovers multiple gene candidates that have been reported to be associated with Alzheimer's disease. [ABSTRACT FROM AUTHOR]

Matrix variate longitudinal discrete data can arise in transcriptomics studies when the data are collected for N genes at r conditions over t time points, and thus, each observation Y n for n = 1 , ... , N can be written as an r × t matrix. When dealing with such data, the number of parameters in the model can be greatly reduced by considering the matrix variate structure. The components of the covariance matrix then also provide a meaningful interpretation. In this work, a mixture of matrix variate Poisson-log normal distributions is introduced for clustering longitudinal read counts from RNA-seq studies. To account for the longitudinal nature of the data, a modified Cholesky-decomposition is utilized for a component of the covariance structure. Furthermore, a parsimonious family of models is developed by imposing constraints on elements of these decompositions. The models are applied to both real and simulated data, and it is demonstrated that the proposed approach can recover the underlying cluster structure. [ABSTRACT FROM AUTHOR]

Background: The differential network (DN) analysis identifies changes in measures of association among genes under two or more experimental conditions. In this article, we introduce a pseudo-value regression approach for network analysis (PRANA). This is a novel method of differential network analysis that also adjusts for additional clinical covariates. We start from mutual information criteria, followed by pseudo-value calculations, which are then entered into a robust regression model. Results: This article assesses the model performances of PRANA in a multivariable setting, followed by a comparison to dnapath and DINGO in both univariable and multivariable settings through variety of simulations. Performance in terms of precision, recall, and F1 score of differentially connected (DC) genes is assessed. By and large, PRANA outperformed dnapath and DINGO, neither of which is equipped to adjust for available covariates such as patient-age. Lastly, we employ PRANA in a real data application from the Gene Expression Omnibus database to identify DC genes that are associated with chronic obstructive pulmonary disease to demonstrate its utility. Conclusion: To the best of our knowledge, this is the first attempt of utilizing a regression modeling for DN analysis by collective gene expression levels between two or more groups with the inclusion of additional clinical covariates. By and large, adjusting for available covariates improves accuracy of a DN analysis. [ABSTRACT FROM AUTHOR]

High-throughput RNA-seq enables comprehensive analysis of the transcriptome for various purposes. However, this technology generally generates massive amounts of sequencing reads with a shorter read length. Consequently, fast, accurate, and flexible tools are needed for assembling raw RNA-seq data into full-length transcripts and quantifying their expression levels. In this protocol, we report TransBorrow, a novel transcriptome assembly software specifically designed for short RNA-seq reads. TransBorrow is employed in conjunction with a splice-aware alignment tool (e.g. Hisat2 and Star) and some other transcriptome assembly tools (e.g. StringTie, Cufflinks, and Scallop). The protocol encompasses all necessary steps, starting from downloading and processing raw sequencing data to assembling the full-length transcripts and quantifying their expressed abundances. The execution time of the protocol may vary depending on the sizes of processed datasets and computational platforms. [ABSTRACT FROM AUTHOR]

Selection of differentially expressed genes (DEGs) is a vital process to discover the causes of diseases. It has been shown that modelling of genomics data by considering relation among genes increases the predictive performance of methods compared to univariate analysis. However, there exist serious differences among most studies analyzing the same dataset for the reasons arising from the methods. Therefore, there is a strong need for easily accessible, user-friendly, and interactive tool to perform gene selection for RNA-seq data via machine learning algorithms simultaneously not to miss DEGs. We develop an open-source and freely available web-based tool for gene selection via machine learning algorithms that can deal with high performance computation. This tool includes six machine learning algorithms having different aspects. Moreover, the tool involves classical pre-processing steps; filtering, normalization, transformation, and univariate analysis. It also offers well-arranged graphical approaches; network plot, heatmap, venn diagram, and box-and-whisker plot. Gene ontology analysis is provided for both mRNA and miRNA DEGs. The implementation is carried out on Alzheimer RNA-seq data to demonstrate the use of this web-based tool. Eleven genes are suggested by at least two out of six methods. One of these genes, hsa-miR-148a-3p, might be considered as a new biomarker for Alzheimer's disease diagnosis. Kidney Chromophobe dataset is also analyzed to demonstrate the validity of GeneSelectML web tool on a different dataset. GeneSelectML is distinguished in that it simultaneously uses different machine learning algorithms for gene selection and can perform pre-processing, graphical representation, and gene ontology analyses on the same tool. This tool is freely available at www.softmed.hacettepe.edu.tr/GeneSelectML. [ABSTRACT FROM AUTHOR]

Background. Differential gene expression analysis using RNA sequencing technology (RNA-Seq) has become the most popular technique in transcriptome research. Although many R packages have been developed to analyze differentially expressed genes (DEGs), several evaluations have shown that no single DEG analysis method outperforms all others. The validity of DEG identification could be increased by using multiple methods and producing the consensus results. However, DEG analysis methods are complex and most of them require prior knowledge of a programming language or command-line shell. Users who do not have this knowledge need to invest time and effort to acquire it. Methods. We developed a novel web application called "bestDEG" to automatically analyze DEGs with different tools and compare the results. A differential expression (DE) analysis pipeline was created combining the edgeR, DESeq2, NOISeq, and EBSeq packages; selected because they use different statistical methods to identify DEGs. bestDEG was evaluated on human datasets from the MicroArray Quality Control (MAQC) project. Results. The performance of the bestDEG web application with the human datasets showed excellent results, and the consensus method outperformed the otherDEanalysis methods in terms of precision (94.71%) and specificity (97.01%). bestDEG is a rapid and efficient tool to analyze DEGs. With bestDEG, users can select DE analysis methods and parameters in the user-friendly web interface. bestDEG also provides a Venn diagram and a table of results. Moreover, the consensus method of this tool can maximize the precision or minimize the false discovery rate (FDR), which reduces the cost of gene expression validation by minimizing wet-lab experiments. [ABSTRACT FROM AUTHOR]

Increasingly amounts of biological data promote the development of various penalized regression models. This review discusses the recent advances in both linear and logistic regression models with penalization terms. This review is mainly focused on various penalized regression models, some of the corresponding optimization algorithms, and their applications in biological data. The pros and cons of different models in terms of response prediction, sample classification, network construction and feature selection are also reviewed. The performances of different models in a real-world RNA-seq dataset for breast cancer are explored. Finally, some future directions are discussed. [ABSTRACT FROM AUTHOR]

Acute myeloid leukemia (AML) is an aggressive type of leukemia, characterized by the accumulation of highly proliferative blasts with a disrupted myeloid differentiation program. Current treatments are ineffective for most patients, partly due to the genetic heterogeneity of AML. This is driven by genetically distinct leukemia stem cells, resulting in relapse even after most of the tumor cells are destroyed. Thus, personalized treatment approaches addressing cellular heterogeneity are urgently required. Reconstruction of Transcriptional regulatory Networks (RTN) is a tool for inferring transcriptional activity in patients with various diseases. In this study, we applied this method to transcriptome profiles of AML patients to test if it provided additional information for the interpretation of transcriptome data. We showed that when RTN results were added to RNA-seq results, superior clusters were formed, which were more homogenous and allowed the better separation of patients with low and high survival rates. We concluded that the external knowledge used for RTN analysis improved the ability of unsupervised machine learning to find meaningful patterns in the data. [ABSTRACT FROM AUTHOR]

Background: Esophageal cancer is a tumor type with high invasiveness and low prognosis. As immunotherapy has been shown to improve the prognosis of esophageal cancer patients, we were interested in the establishment of an immune-associated gene prognostic index to effectively predict the prognosis of patients. Methods: To establish the immune-related gene prognostic index of esophageal cancer (EC), we screened 363 upregulated and 83 downregulated immune-related genes that were differentially expressed in EC compared to normal tissues. By multivariate Cox regression and weighted gene coexpression network analysis (WGCNA), we built a prognostic model based on eight immune-related genes (IRGs). We confirmed the prognostic model in both TCGA and GEO cohorts and found that the low-risk group had better overall survival than the high-risk group. Results: In this study, we identified 363 upregulated IRGs and 83 downregulated IRGs. Next, we found a prognostic model that was constructed with eight IRGs (OSM, CEACAM8, HSPA6, HSP90AB1, PCSK2, PLXNA1, TRIB2, and HMGB3) by multivariate Cox regression analysis and WGCNA. According to the Kaplan–Meier survival analysis results, the model we constructed can predict the prognosis of patients with esophageal cancer. This result can be verified by the Gene Expression Omnibus (GEO). Patients were divided into two groups with different outcomes. IRGPI-low patients had better overall survival than IRGPI-high patients. Conclusion: Our findings indicated the potential value of the IRGPI risk model for predicting the prognosis of EC patients [ABSTRACT FROM AUTHOR]

Background: Salvia miltiorrhiza is an important traditional Chinese medicinal (TCM) plant and a model plant in the genetic study of TCM. A series of omics related to Danshen have been published. Integrating, managing, storing, and sharing data has become an urgent problem to be solved in S. miltiorrhiza genetic studies. Objectives: The genome database is the link for the exchange, acquisition, and use of different omics data between data producers and users, maximizing value and utilization of data. Methods: The genome database included DSS3 genome and five RNA-Seq data. The back-end performs data search and retrieval through the LAMP (Linux, Apache, MySQL, PHP) framework. Results: Here, we present SmGDB (S. miltiorrhiza genome database; http://8.140.162.85/), which houses the latest version of genome sequence and annotation data for S. miltiorrhiza, combining three unpublished RNA-Seq data from our group and two released RNA-Seq data. We also identified a novel gene cluster including seven CYP71D genes involved in the tanshinone synthesis pathway based on genome sequences and expression data. Besides, SmGDB provides user-friendly web interfaces for querying and browsing gene annotation, structure, location, and expression profiles for concerned genes. Popular bioinformatics tools such as 'BLAST', 'Search', 'Heatmap', 'JBrowse', etc., were also provided in SmGDB. Conclusions: SmGDB will provide utility for characterizing the structure of the S. miltiorrhiza genome and better understanding gene functions and biological processes underlying complex secondary metabolism in Danshen. [ABSTRACT FROM AUTHOR]

Alternative splicing (AS) is an essential co-transcriptional regulatory mechanism in eukaryotes. The accumulation of plant RNA-Seq data provides an unprecedented opportunity to investigate the global landscape of plant AS events. However, most existing AS identification tools were originally designed for animals, and their performance in plants was not rigorously benchmarked. In this work, we developed a simple and easy-to-use bioinformatics tool named ASTool for detecting AS events from plant RNA-Seq data. As an exon-based method, ASTool can detect 4 major AS types, including intron retention (IR), exon skipping (ES), alternative 5′ splice sites (A5SS), and alternative 3′ splice sites (A3SS). Compared with existing tools, ASTool revealed a favorable performance when tested in simulated RNA-Seq data, with both recall and precision values exceeding 95% in most cases. Moreover, ASTool also showed a competitive computational speed and consistent detection results with existing tools when tested in simulated or real plant RNA-Seq data. Considering that IR is the most predominant AS type in plants, ASTool allowed the detection and visualization of novel IR events based on known splice sites. To fully present the functionality of ASTool, we also provided an application example of ASTool in processing real RNA-Seq data of Arabidopsis in response to heat stress. [ABSTRACT FROM AUTHOR]

Considering human brain disorders, Major Depressive Disorder (MDD) is seen as a lethal disease in which a person goes to the extent of suicidal behavior. Physical detection of MDD patients is less precise but machine learning can aid in improved classification of disease. The present research included three RNA-seq data classes to classify DEGs and then train key gene data using a random forest machine learning method. The three classes in the sample are 29 CON (sudden death healthy control), 21 MDD-S (a Major Depressive Disorder Suicide) being included in the second group, and 9 MDD (non-suicides MDD) which are included in the third group. With PCA analysis, 99 key genes were obtained. 47.1% data variability is given by these 99 genes. The model training of 99 genes indicated improved classification. The RF classification model has an accuracy of 61.11% over test data and 97.56% over train data. It was also noticed that the RF method offered greater accuracy than the KNN method. 99 genes were annotated using DAVID and ClueGo packages. Some of the important pathways and function observed in the study were glutamatergic synapse, GABA receptor activation, long-term synaptic depression, and morphine addiction. Out Of 99 genes, four genes, namely DLGAP1, GNG2, GRIA1, and GRIA4, were found to be predominantly involved in the glutamatergic synapse pathway. Another substantial link was observed in the GABA receptor activation involving the following two genes, GABBR2 and GNG2. Also, the genes found responsible for long-term synaptic depression were GRIA1, MAPT, and PTEN. There was another finding of morphine addiction which comprises three genes, namely GABBR2, GNG2, and PDE4D. For massive datasets, this approach will act as the gold standard. The cases of CON, MDD, and MDD-S are physically distinct. There was dysregulation in the expression level of 12 genes. The 12 genes act as a possible biomarker for Major Depressive Disorder and open up a new path for depressed subjects to explore further. [ABSTRACT FROM AUTHOR]

Background: Recent development of bioinformatics tools for Next Generation Sequencing data has facilitated complex analyses and prompted large scale experimental designs for comparative genomics. When combined with the advances in network inference tools, this can lead to powerful methodologies for mining genomics data, allowing development of pipelines that stretch from sequence reads mapping to network inference. However, integrating various methods and tools available over different platforms requires a programmatic framework to fully exploit their analytic capabilities. Integrating multiple genomic analysis tools faces challenges from standardization of input and output formats, normalization of results for performing comparative analyses, to developing intuitive and easy to control scripts and interfaces for the genomic analysis pipeline. Results: We describe here NetSeekR, a network analysis R package that includes the capacity to analyze time series of RNA-Seq data, to perform correlation and regulatory network inferences and to use network analysis methods to summarize the results of a comparative genomics study. The software pipeline includes alignment of reads, differential gene expression analysis, correlation network analysis, regulatory network analysis, gene ontology enrichment analysis and network visualization of differentially expressed genes. The implementation provides support for multiple RNA-Seq read mapping methods and allows comparative analysis of the results obtained by different bioinformatics methods. Conclusion: Our methodology increases the level of integration of genomics data analysis tools to network inference, facilitating hypothesis building, functional analysis and genomics discovery from large scale NGS data. When combined with network analysis and simulation tools, the pipeline allows for developing systems biology methods using large scale genomics data. [ABSTRACT FROM AUTHOR]

Expression quantitative trait loci (eQTL) analysis has been widely used in interpreting disease-associated loci through correlating genetic variant loci with the expression of specific genes. RNA-sequencing (RNA-Seq), which can quantify gene expression at the genome-wide level, is often used in eQTL identification. Since different normalization methods of gene expression have substantial impacts on RNA-seq downstream analysis, it is of great necessity to systematically compare the effects of these methods on eQTL identification. Here, by using RNA-seq and genotype data of four different cancers in The Cancer Genome Atlas (TCGA) database, we comprehensively evaluated the effect of eight commonly used normalization methods on eQTL identification. Our results showed that the application of different methods could cause 20–30% differences in the final results of eQTL identification. Among these methods, COUNT, Median of Ratio (MED) and Trimmed Mean of M-values (TMM) generated similar results for identifying eQTLs, while Fragments Per Kilobase Million (FPKM) or RANK produced more differential results compared with other methods. Based on the accuracy and receiver operating characteristic (ROC) curve, the TMM method was found to be the optimal method for normalizing gene expression data in eQTLs analysis. In addition, we also evaluated the performance of different pairwise combinations of these methods. As a result, compared with single normalization methods, the combination of methods can not only identify more cis-eQTLs, but also improve the performance of the ROC curve. Overall, this study provides a comprehensive comparison of normalization methods for identifying eQTLs from RNA-seq data, and proposes some practical recommendations for diverse scenarios. [ABSTRACT FROM AUTHOR]

Copyright of Mycosystema is the property of Mycosystema Editorial Board and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

The tumor microenvironment is composed of tumor cells, stroma cells, immune cells, blood vessels, and other associated non-cancerous cells. Gene expression measurements on tumor samples are an average over cells in the microenvironment. However, research questions often seek answers about tumor cells rather than the surrounding non-tumor tissue. Previous studies have suggested that the tumor purity (TP)—the proportion of tumor cells in a solid tumor sample—has a confounding effect on differential expression (DE) analysis of high vs. low survival groups. We investigate three ways incorporating the TP information in the two statistical methods used for analyzing gene expression data, namely, differential network (DN) analysis and DE analysis. Analysis 1 ignores the TP information completely, Analysis 2 uses a truncated sample by removing the low TP samples, and Analysis 3 uses TP as a covariate in the underlying statistical models. We use three gene expression data sets related to three different cancers from the Cancer Genome Atlas (TCGA) for our investigation. The networks from Analysis 2 have greater amount of differential connectivity in the two networks than that from Analysis 1 in all three cancer datasets. Similarly, Analysis 1 identified more differentially expressed genes than Analysis 2. Results of DN and DE analyses using Analysis 3 were mostly consistent with those of Analysis 1 across three cancers. However, Analysis 3 identified additional cancer-related genes in both DN and DE analyses. Our findings suggest that using TP as a covariate in a linear model is appropriate for DE analysis, but a more robust model is needed for DN analysis. However, because true DN or DE patterns are not known for the empirical datasets, simulated datasets can be used to study the statistical properties of these methods in future studies. [ABSTRACT FROM AUTHOR]

Background: Identifying 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. Results: In 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. Conclusions: Simulation 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. [ABSTRACT FROM AUTHOR]

The Pan-Cancer Atlas consists of original sequencing data from various sources, provides the opportunity to perform systematic studies on the commonalities and differences between diverse cancers. The analysis for the pan-cancer dataset could help researchers to identify the key factors that could trigger cancer. In this paper, we present a novel pan-cancer classification method, referred to MI_DenseNetCAM, to identify a set of genes that can differentiate all tumor types accurately. First, the Mutual Information (MI) was utilized to eliminate noise and redundancy from the pan-cancer datasets. Then, the gene data was further converted to 2D images. Next, the DenseNet model was adopted as a classifier and the Guided Grad-CAM algorithm was applied to identify the key genes. Extensive experimental results on the public RNA-seq data sets with 33 different tumor types show that our method outperforms the other state-of-the-art classification methods. Moreover, gene analysis further demonstrated that the genes selected by our method were related to the corresponding tumor types. [ABSTRACT FROM AUTHOR]

Introduction (SNPs) are single nucleotide base variations, caused by transitions (C/T or G/A) or transversions (C/G, C/A, or T/A, T/G), in the same position between individual genomic DNA sequences. Single nucleotide polymorphisms have been applied as important molecular markers in genetics and breeding studies. About 40% of the Single nucleotide polymorphisms in the genes cause a change in an amino acid. The rapid advance of next generation sequencing provides a high-throughput means of SNP discovery. Transcriptome study can fill the gap between genotype and phenotype and help understanding the mechanisms from sequence to function. RNA sequencing (RNA-Seq) is a next generating sequencing based technology for studying of whole transcriptome and gene expression. It simultaneously enables study of transcriptomics sequences and very accurate quantitative gene expression (digital expression). Hence, these data are very suitable for high-throughput study of expression level of all transcribed genes and their SNPs (Single Nucleotide Polymorphism. Recently, RNA-Seq has also been used as an efficient and cost-effective method to systematically identify SNPs in transcribed regions in different species. A transcriptomics-based sequencing approach offers a cheaper alternative to identify a large number of polymorphisms and possibly to discover causative variants. Materials and Methods In this study, RNA-Seq data were used to SNP discovery in American Holstein (Bos taurus) and Pakistanian Cholistani (Bos indicus) cows. RNA-Seq data of 21,078,477 and 20940063 paired end reads with 75 bp length resulted from pooling of whole blood samples of 40 Holstein cows at the University of Wisconsin, Dairy Cattle Center, USA, and 45 Cholistani cows at Gujait Peer Farm, Bahawalpur, Punjab, Pakistan, respectively, obtained from SRA database in NCBI for Holstein cows (http://www.ncbi.nlm.nih.gov/sra/SRX317197) and Cholistani cows http://www.ncbi.nlm.nih.gov/sra/SRS454433). MRNA sequencing was run on Illumina Genome Analyzer IIx (Illumina Inc., San Diego, CA). Data were converted from Sra format to Fastq format by fastq-dump command from Ubuntu linux version of Sratoolkit 2.5.4-1. Data quality control was checked by FastQC (v0.11.3) likewise trimmed for linked adaptors and bad quality reads by Trimmomatic 0.33 Adaptors were considered according to sequencing instrument as default (TruSeq2-PE.fa) and the minimum read length was set at 50 bp. Trimmed reads were aligned on UMD3.1 reference genome (release 81) based on annotation data by Tophat2, which applies Bowtie2 as the aligner. The transcriptome was assembled by TopHat2 software in two cow's population by aligning and mapping the RNA-Seq reads on bovine reference genome. The SNPs were discovered by Samtools software. Results and Discussion After data editing, the removed and low quality reads in both breeds were almost equal and relatively low. The length of whole transcriptome assembled, for example 52798651 bases in Holstein, indicates around 2% of the whole genome (around 2.6 Mbp) expressed as mRNA. In Cholistani cows, read mapping rate for forward and reverse reads were 81.3 and 79.9%, respectively, and multiple alignments rate was about 9.4%. Overall read mapping was 80.6% and concordant pair alignment was 70.1%. In Holstein cows, read mapping rate for forward and reverse reads were 66.3 and 55.4%, respectively, and multiple alignments rate was about 7.2%,. Overall read mapping was 60.8% and concordant pair alignment was 51.3%. Results show that 50183 and 137954 SNPs were discovered on the assembled transcriptome of Holstein and Cholistani cow's samples, respectively, and 15308 SNPs were common in both breeds. No direct relation was found between the number of discovered SNPs and the chromosome length. Also 12 SNP types were identified including 4 transition and 8 transversion. The most commonly discovered SNP were transition, which were 70.6% in Cholistani and 69.6% in Holstein cows. The ratio of transition to transversion SNP (Ts / Tv) was 2.4 and 2.3 in Cholistani and Holstein cows, respectively. The number of discovered SNPs in Cholistani cows were approximately three times higher than Holstein cows. Because, for the alignment of both species used a same reference genome with Herford origin. Conclusion the expression difference between two alleles in a single-nucleotide position causes phenotype diversity and probably explains the large part of variances between these two bovine subspecies, especially in diversity, susceptibility to disease and parasites, tolerating environmental stress such as biological and non-biological stresses in different environmental conditions. While, differential gene expression analysis or even allelic specific expression in gene level may not be able to explain phenotype diversity. [ABSTRACT FROM AUTHOR]

Next-generation sequencing has emerged as an essential technology for the quantitative analysis of gene expression. In medical research, RNA sequencing (RNA-seq) data are commonly used to identify which type of disease a patient has. Because of the discrete nature of RNA-seq data, the existing statistical methods that have been developed for microarray data cannot be directly applied to RNA-seq data. Existing statistical methods usually model RNA-seq data by a discrete distribution, such as the Poisson, the negative binomial, or the mixture distribution with a point mass at zero and a Poisson distribution to further allow for data with an excess of zeros. Consequently, analytic tools corresponding to the above three discrete distributions have been developed: Poisson linear discriminant analysis (PLDA), negative binomial linear discriminant analysis (NBLDA), and zero-inflated Poisson logistic discriminant analysis (ZIPLDA). However, it is unclear what the real distributions would be for these classifications when applied to a new and real dataset. Considering that count datasets are frequently characterized by excess zeros and overdispersion, this paper extends the existing distribution to a mixture distribution with a point mass at zero and a negative binomial distribution and proposes a zero-inflated negative binomial logistic discriminant analysis (ZINBLDA) for classification. More importantly, we compare the above four classification methods from the perspective of model parameters, as an understanding of parameters is necessary for selecting the optimal method for RNA-seq data. Furthermore, we determine that the above four methods could transform into each other in some cases. Using simulation studies, we compare and evaluate the performance of these classification methods in a wide range of settings, and we also present a decision tree model created to help us select the optimal classifier for a new RNA-seq dataset. The results of the two real datasets coincide with the theory and simulation analysis results. The methods used in this work are implemented in the open-scource R scripts, with a source code freely available at https://github.com/FocusPaka/ZINBLDA. [ABSTRACT FROM AUTHOR]

Variations of gene expression levels play an important role in tumors. There are numerous methods to identify differentially expressed genes in high-throughput sequencing. Several algorithms endeavor to identify distinctive genetic patterns susceptable to particular diseases. Although these processes have been proved successful, the probability that the number of non-differentially expressed genes measured by false discovery rate (FDR) has a large standard deviation, and the misidentification rate (type I error) grows rapidly when the number of genes to be detected become larger. In this study we developed a new method, Unit Gamma Measurement (UGM), accounting for multiple hypotheses test statistics distribution, which could reduce the dependency problem. Simulated expression profile data and breast cancer RNA-Seq data were utilized to testify the accuracy of UGM. The results show that the number of non-differentially expressed genes identified by the UGM is very close to the real-evidence data, and the UGM also has a smaller standard error, range, quartile range and RMS error. In addition, the UGM can be used to screen many breast cancer-associated genes, such as BRCA1, BRCA2, PTEN, BRIP1, etc., provides better accuracy, robustness and efficiency, the method of identification differentially expressed genes in high-throughput sequencing. [ABSTRACT FROM AUTHOR]

Key message: An interologous PPI network of tea leaf is designed by developing reference transcriptome assembly and using experimentally validated PPIs in plants. Key regulatory proteins are proposed and potential TFs are predicted. Worldwide, tea (Camellia sinensis) is the most consumed beverage primarily due to the taste, flavour, and aroma of its newly formed leaves; and has been used as an important ingredient in several traditional medicinal systems because of its antioxidant properties. For this medicinally and commercially important plant, design principles of gene-regulatory and protein–protein interaction (PPI) networks at sub-cellular level are largely un-characterized. In this work, we report a tea leaf interologous PPI network (TeaLIPIN) consisting of 11,208 nodes and 197,820 interactions. A reference transcriptome assembly was first developed from all the 44 samples of 6 publicly available leaf transcriptomes (1,567,288,290 raw reads). By inferring the high-confidence interactions among potential proteins coded by these transcripts using known experimental information about PPIs in 14 plants, an interologous PPI network was constructed and its modular architecture was explored. Comparing this network with 10,000 realizations of two types of corresponding random networks (Erdős–Rényi and Barabási–Albert models) and examining over three network centrality metrics, we predict 2750 bottleneck proteins (having p values < 0.01). 247 of these are deduced to have transcription factor domains by in-house developed HMM models of known plant TFs and these were also mapped to the draft tea genome for searching their probable loci of origin. Co-expression analysis of the TeaLIPIN proteins was also performed and top ranking modules are elaborated. We believe that the proposed novel methodology can easily be adopted to develop and explore the PPI interactomes in other plant species by making use of the available transcriptomic data. [ABSTRACT FROM AUTHOR]

Golden pompano (Trachinotus ovatus) is an important aquaculture species in China and broadly cultivated in Asia–Pacific region. Our previous study found that T. ovatus had a better tolerance to the high level of carbohydrate in diet (360 g/kg) compared to Micropterus salmoides. Present study investigated the tolerance mechanism of T. ovatus to high level of dietary carbohydrate (PH‐360 g/kg of corn starch) on transcriptome level using Illumina sequencing compared to optimum level of dietary carbohydrate (PM‐240 g/kg corn starch). A total of 208,952,492 reads were obtained and assembled into 62,377 unigenes with an average length of 822 bp and a range of 200–16,206 bp. Transcriptome annotation revealed that 27,636 unigenes with a cut‐off e‐value of 10−5 were functionally annotated. Specifically, 7,903 unigenes were clustered into 24 functional categories by searching against the COG database, 12,161 unigenes to 84 GO terms, 26,835 unigenes to Nr database and 23,681 unigenes to SWISS‐PROT database. In addition, 13,739 unigenes were grouped into 240 KEGG pathways. There are 3,461 differentially expressed genes in PH compared to PM, in which 3,059 unigenes were significantly upregulated and 402 unigenes were significantly downregulated. Twelve key genes related to carbohydrate metabolism were further validated by quantitative real‐time RT–PCR, and the results were consistent with those of transcriptome profile.Under high level of dietary carbohydrate (PH‐360 g/kg), the results indicate that T. ovatus can adapt to the high level of dietary carbohydrate by enhancing glycolysis with the increased heaptic HK enzyme activity (p > 0.05), reducing gluconeogenesis with the decreased trend of PEPCK enzyme activity and maintaining the homeostasis of carbohydrate metabolism‐related hormones and hormone receptor gene expressions. Histopathological examination has shown that liver cavity area was significantly increased (p > 0.05), but no visible nutritional diseases are found in fish. This study provides a valuable understanding of the molecular responses of T. ovatus fed by the high level of dietary carbohydrate (PH‐360 g/kg). [ABSTRACT FROM AUTHOR]

Gene‐expression data is being widely used for various clinical research. It represents expression levels of thousands of genes across the various experimental conditions simultaneously. Mining conditions specific hub genes from gene‐expression data is a challenging task. Conditions specific hub genes signify the functional behaviour of bicluster across the subset of conditions and can act as prognostic or diagnostic markers of the diseases. In this study, the authors have introduced a new approach for identifying conditions specific hub genes from the RNA‐Seq data using a biclustering algorithm. In the proposed approach, efficient 'runibic' biclustering algorithm, the concept of gene co‐expression network and concept of protein–protein interaction network have been used for getting better performance. The result shows that the proposed approach extracts biologically significant conditions specific hub genes which play an important role in various biological processes and pathways. These conditions specific hub genes can be used as prognostic or diagnostic biomarkers. Conditions specific hub genes will be helpful to reduce the analysis time and increase the accuracy of further research. Also, they summarised application of the proposed approach to the drug discovery process. [ABSTRACT FROM AUTHOR]

Background: As RNA-seq becomes the assay of choice for measuring gene expression levels, differential expression analysis has received extensive attentions of researchers. To date, for the evaluation of DE methods, most attention has been paid on validity. Yet another important aspect of DE methods, stability, is overlooked and has not been studied to the best of our knowledge. Results: In this study, we empirically show the need of assessing stability of DE methods and propose a stability metric, called Area Under the Correlation curve (AUCOR), that generates the perturbed datasets by a mixture distribution and combines the information of similarities between sets of selected features from these perturbed datasets and the original dataset. Conclusion: Empirical results support that AUCOR can effectively rank the DE methods in terms of stability for given RNA-seq datasets. In addition, we explore how biological or technical factors from experiments and data analysis affect the stability of DE methods. AUCOR is implemented in the open-source R package AUCOR, with source code freely available at https://github.com/linbingqing/stableDE. [ABSTRACT FROM AUTHOR]

Background: Positively correlated with carcass weight and animal growth, the ribeye area (REA) and the backfat thickness (BFT) are economic important carcass traits, which impact directly on producer's payment. The selection of these traits has not been satisfactory since they are expressed later in the animal's life and multigene regulated. So, next-generation technologies have been applied in this area to improve animal's selection and better understand the molecular mechanisms involved in the development of these traits. Correlation network analysis, performed by tools like WGCNA (Weighted Correlation Network Analysis), has been used to explore gene-gene interactions and gene-phenotype correlations. Thus, this study aimed to identify putative candidate genes and metabolic pathways that regulate REA and BFT by constructing a gene co-expression network using WGCNA and RNA sequencing data, to better understand genetic and molecular variations behind these complex traits in Nelore cattle. Results: The gene co-expression network analysis, using WGCNA, were built using RNA-sequencing data normalized by transcript per million (TPM) from 43 Nelore steers. Forty-six gene clusters were constructed, between them, three were positively correlated (p-value< 0.1) to the BFT (Green Yellow, Ivory, and Light Yellow modules) and, one cluster was negatively correlated (p-value< 0.1) with REA (Salmon module). The enrichment analysis performed by DAVID and WebGestalt (FDR 5%) identified eight Gene Ontology (GO) terms and three KEGG pathways in the Green Yellow module, mostly associated with immune response and inflammatory mechanisms. The enrichment of the Salmon module demonstrated 19 GO terms and 21 KEGG pathways, related to muscle energy metabolism, lipid metabolism, muscle degradation, and oxidative stress diseases. The Ivory and Light yellow modules have not shown significant results in the enrichment analysis. Conclusion: With this study, we verified that inflammation and immune response pathways modulate the BFT trait. Energy and lipid metabolism pathways, highlighting fatty acid metabolism, were the central pathways associated with REA. Some genes, as RSAD2, EIF2AK2, ACAT1, and ACSL1 were considered as putative candidate related to these traits. Altogether these results allow us to a better comprehension of the molecular mechanisms that lead to muscle and fat deposition in bovine. [ABSTRACT FROM AUTHOR]

Cancer belongs to a class of highly aggressive diseases and a leading cause of death in the world. With more than 100 types of cancers, breast, lung and prostate cancer remain to be the most common types. To identify essential network markers (NMs) and therapeutic targets in these cancers, the authors present a novel approach which uses gene expression data from microarray and RNA‐seq platforms and utilises the results from this data to evaluate protein–protein interaction (PPI) network. Differentially expressed genes (DEGs) are extracted from microarray data using three different statistical methods in R, to produce a consistent set of genes. Also, DEGs are extracted from RNA‐seq data for the same three cancer types. DEG sets found to be common in both platforms are obtained at three fold change (FC) cut‐off levels to accurately identify the level of change in expression of these genes in all three cancers. A cancer network is built using PPI data characterising gene sets at log‐FC (LFC)>1, LFC>1.5 and LFC>2, and interconnection between principal hub nodes of these networks is observed. Resulting network of hubs at three FC levels highlights prime NMs with high confidence in multiple cancers as validated by Gene Ontology functional enrichment and maximal complete subgraphs from CFinder. [ABSTRACT FROM AUTHOR]

In bioinformatics, exon skipping (ES) event prediction is an essential part of alternative splicing (AS) event analysis. Although many methods have been developed to predict ES events, a solution has yet to be found. In this study, given the limitations of machine learning algorithms with RNA-Seq data or genome sequences, a new feature, called RS (RNA-seq and sequence) features, was constructed. These features include RNA-Seq features derived from the RNA-Seq data and sequence features derived from genome sequences. We propose a novel Rotation Forest classifier to predict ES events with the RS features (RotaF-RSES). To validate the efficacy of RotaF-RSES, a dataset from two human tissues was used, and RotaF-RSES achieved an accuracy of 98.4%, a specificity of 99.2%, a sensitivity of 94.1%, and an area under the curve (AUC) of 98.6%. When compared to the other available methods, the results indicate that RotaF-RSES is efficient and can predict ES events with RS features. [ABSTRACT FROM AUTHOR]

Background: Predicting disease-associated genes is helpful for understanding the molecular mechanisms during the disease progression. Since the pathological mechanisms of neurodegenerative diseases are very complex, traditional statistic-based methods are not suitable for identifying key genes related to the disease development. Recent studies have shown that the computational models with deep structure can learn automatically the features of biological data, which is useful for exploring the characteristics of gene expression during the disease progression. Results: In this paper, we propose a deep learning approach based on the restricted Boltzmann machine to analyze the RNA-seq data of Huntington's disease, namely stacked restricted Boltzmann machine (SRBM). According to the SRBM, we also design a novel framework to screen the key genes during the Huntington's disease development. In this work, we assume that the effects of regulatory factors can be captured by the hierarchical structure and narrow hidden layers of the SRBM. First, we select disease-associated factors with different time period datasets according to the differentially activated neurons in hidden layers. Then, we select disease-associated genes according to the changes of the gene energy in SRBM at different time periods. Conclusions: The experimental results demonstrate that SRBM can detect the important information for differential analysis of time series gene expression datasets. The identification accuracy of the disease-associated genes is improved to some extent using the novel framework. Moreover, the prediction precision of disease-associated genes for top ranking genes using SRBM is effectively improved compared with that of the state of the art methods. [ABSTRACT FROM AUTHOR]

As gene expression measurement technology is shifting from microarrays to sequencing, the statistical tools available for their analysis must be adapted since RNA-seq data are measured as counts. It has been proposed to model RNA-seq counts as continuous variables using nonparametric regression to account for their inherent heteroscedasticity. In this vein, we propose tcgsaseq, a principled, model-free, and efficient method for detecting longitudinal changes in RNA-seq gene sets defined a priori. The method identifies those gene sets whose expression varies over time, based on an original variance component score test accounting for both covariates and heteroscedasticity without assuming any specific parametric distribution for the (transformed) counts. We demonstrate that despite the presence of a nonparametric component, our test statistic has a simple form and limiting distribution, and both may be computed quickly. A permutation version of the test is additionally proposed for very small sample sizes. Applied to both simulated data and two real datasets, tcgsaseq is shown to exhibit very good statistical properties, with an increase in stability and power when compared to state-of-the-art methods ROAST (rotation gene set testing), edgeR, and DESeq2, which can fail to control the type I error under certain realistic settings. We have made the method available for the community in the R package tcgsaseq. [ABSTRACT FROM AUTHOR]

Background: Insecticide resistance is a substantial problem in controlling agricultural and medical pests. Detecting target site mutations is crucial to manage insecticide resistance. Though PCR-based methods have been widely used in this field, they are time-consuming and inefficient, and typically have a high false positive rate. Acetylcholinesterases (Ace) is the neural target of the widely used organophosphate (OP) and carbamate insecticides. However, there is not any software available to detect insecticide resistance associated mutations in RNA-Seq data at present. Results: A computational pipeline ACE was developed to detect resistance mutations of ace in insect RNA-Seq data. Known ace resistance mutations were collected and used as a reference. We constructed a Web server for ACE, and the standalone software in both Linux and Windows versions is available for download. ACE was used to analyse 971 RNA-Seq data from 136 studies in 7 insect pests. The mutation frequency of each RNA-Seq dataset was calculated. The results indicated that the resistance frequency was 30%-44% in an eastern Ugandan Anopheles population, thus suggesting this resistance-conferring mutation has reached high frequency in these mosquitoes in Uganda. Analyses of RNA-Seq data from the diamondback moth Plutella xylostella indicated that the G227A mutation was positively related with resistance levels to organophosphate or carbamate insecticides. The wasp Nasonia vitripennis had a low frequency of resistant reads (<5%), but the agricultural pests Chilo suppressalis and Bemisia tabaci had a high resistance frequency. All ace reads in the 30 B. tabaci RNA-Seq data were resistant reads, suggesting that insecticide resistance has spread to very high frequency in B. tabaci. Conclusions: To the best of our knowledge, the ACE pipeline is the first tool to detect resistance mutations from RNASeq data, and it facilitates the full utilization of large-scale genetic data obtained by using next-generation sequencing. [ABSTRACT FROM AUTHOR]

We present a framework for the simultaneous estimation of state and parameters of partially observed Boolean dynamical systems (POBDS). Simultaneous state and parameter estimation is achieved through the combined use of the Boolean Kalman filter and Boolean Kalman smoother, which provide the minimum mean-square error state estimators for the POBDS model, and maximum-likelihood (ML) parameter estimation; in the presence of continuous parameters, ML estimation is performed using the expectation–maximization algorithm. The performance of the proposed ML adaptive filter is demonstrated by numerical experiments with a POBDS model of gene regulatory networks observed through noisy next-generation sequencing (RNA-seq) time series data using the well-known p53-MDM2 negative-feedback loop gene regulatory model. [ABSTRACT FROM PUBLISHER]

In the past 5 years, RNA-Seq has become a powerful tool in transcriptome analysis even though computational methods dedicated to the analysis of high-throughput sequencing data are yet to be standardized. It is, however, now commonly accepted that the choice of a normalization procedure is an important step in such a process, for example in differential gene expression analysis. The present article highlights the similarities between three normalization methods: TMM from edgeR R package, RLE from DESeq2 R package, and MRN. Both TMM and DESeq2 are widely used for differential gene expression analysis. This paper introduces properties that show when these three methods will give exactly the same results. These properties are proven mathematically and illustrated by performing in silico calculations on a given RNA-Seq data set. [ABSTRACT FROM AUTHOR]

Next-generation sequencing technologies now constitute a method of choice to measure gene expression. Data to analyze are read counts, commonly modeled using negative binomial distributions. A relevant issue associated with this probabilistic framework is the reliable estimation of the overdispersion parameter, reinforced by the limited number of replicates generally observable for each gene. Many strategies have been proposed to estimate this parameter, but when differential analysis is the purpose, they often result in procedures based on plug-in estimates, and we show here that this discrepancy between the estimation framework and the testing framework can lead to uncontrolled type-I errors. Instead, we propose a mixture model that allows each gene to share information with other genes that exhibit similar variability. Three consistent statistical tests are developed for differential expression analysis. We show through a wide simulation study that the proposed method improves the sensitivity of detecting differentially expressed genes with respect to the common procedures, since it reaches the nominal value for the type-I error, while keeping elevate discriminative power between differentially and not differentially expressed genes. The method is finally illustrated on prostate cancer RNA-Seq data. [ABSTRACT FROM AUTHOR]

The identification of essential proteins is not only important for understanding organism structure on the molecular level, but also beneficial to drug-target detection and genetic disease prevention. Traditional methods often employ various centrality indices of static protein-protein interaction (PPI) networks and/or gene expression profiles to predict essential proteins. However, the prediction accuracy of most methods still has room to be further improved. In this study, we propose a strategy to increase the prediction accuracy of essential protein identification in three ways. Firstly, RNA-Seq datasets are employed to construct integrated dynamic PPI networks. Using a RNA-Seq dataset is expected to give more accurate predictions than using microarray gene expression profiles. Secondly, a novel integrated dynamic PPI network is constructed by considering both the co-expression pattern and the co-expression level of the RNA-Seq data. Thirdly, a novel two-step strategy is proposed to identify essential proteins from two known centrality indices. Numerical experiments have shown that the proposed strategy can increase the prediction accuracy dramatically, which can be generalized to many existing methods and centrality indices. [ABSTRACT FROM AUTHOR]

Background: RNA-sequencing (RNA-Seq) has become a popular tool for transcriptome profiling in mammals. However, accurate estimation of allele-specific expression (ASE) based on alignments of reads to the reference genome is challenging, because it contains only one allele on a mosaic haploid genome. Even with the information of diploid genome sequences, precise alignment of reads to the correct allele is difficult because of the high-similarity between the corresponding allele sequences. Results: We propose a Bayesian approach to estimate ASE from RNA-Seq data with diploid genome sequences. In the statistical framework, the haploid choice is modeled as a hidden variable and estimated simultaneously with isoform expression levels by variational Bayesian inference. Through the simulation data analysis, we demonstrate the effectiveness of the proposed approach in terms of identifying ASE compared to the existing approach. We also show that our approach enables better quantification of isoform expression levels compared to the existing methods, TIGAR2, RSEM and Cufflinks. In the real data analysis of the human reference lymphoblastoid cell line GM12878, some autosomal genes were identified as ASE genes, and skewed paternal X-chromosome inactivation in GM12878 was identified. Conclusions: The proposed method, called ASE-TIGAR, enables accurate estimation of gene expression from RNASeq data in an allele-specific manner. Our results show the effectiveness of utilizing personal genomic information for accurate estimation of ASE. An implementation of our method is available at http://nagasakilab.csml.org/ase-tigar. [ABSTRACT FROM AUTHOR]

Background: Recently, rapid improvements in technology and decrease in sequencing costs have made RNA-Seq a widely used technique to quantify gene expression levels. Various normalization approaches have been proposed, owing to the importance of normalization in the analysis of RNA-Seq data. A comparison of recently proposed normalization methods is required to generate suitable guidelines for the selection of the most appropriate approach for future experiments. Results: In this paper, we compared eight non-abundance (RC, UQ, Med, TMM, DESeq, Q, RPKM, and ERPKM) and two abundance estimation normalization methods (RSEM and Sailfish). The experiments were based on real Illumina high-throughput RNA-Seq of 35- and 76-nucleotide sequences produced in the MAQC project and simulation reads. Reads were mapped with human genome obtained from UCSC Genome Browser Database. For precise evaluation, we investigated Spearman correlation between the normalization results from RNA-Seq and MAQC qRT-PCR values for 996 genes. Based on this work, we showed that out of the eight non-abundance estimation normalization methods, RC, UQ, Med, TMM, DESeq, and Q gave similar normalization results for all data sets. For RNA-Seq of a 35-nucleotide sequence, RPKM showed the highest correlation results, but for RNA-Seq of a 76-nucleotide sequence, least correlation was observed than the other methods. ERPKM did not improve results than RPKM. Between two abundance estimation normalization methods, for RNA-Seq of a 35-nucleotide sequence, higher correlation was obtained with Sailfish than that with RSEM, which was better than without using abundance estimation methods. However, for RNA-Seq of a 76-nucleotide sequence, the results achieved by RSEM were similar to without applying abundance estimation methods, and were much better than with Sailfish. Furthermore, we found that adding a poly-A tail increased alignment numbers, but did not improve normalization results. Conclusion: Spearman correlation analysis revealed that RC, UQ, Med, TMM, DESeq, and Q did not noticeably improve gene expression normalization, regardless of read length. Other normalization methods were more efficient when alignment accuracy was low; Sailfish with RPKM gave the best normalization results. When alignment accuracy was high, RC was sufficient for gene expression calculation. And we suggest ignoring poly-A tail during differential gene expression analysis. [ABSTRACT FROM AUTHOR]

Markers associated with diagnosis, presentation and potential therapeutic targets have received widespread attention in ovarian cancer research in the past few years. However, the majority of these markers have been investigated individually, and the changes in expression and the association between them are rarely documented. Next-generation sequencing, also termed RNA-seq when the sequencing targets are cDNAs, can provide a whole blueprint of the transcriptome of a specific tissue. In the present study, RNA-seq data of human ovarian cancer samples were used to verify the expression of known markers and to identify the association between them. A total of 563 markers associated with ovarian cancer were retrieved from the database of the National Center of Biotechnology Information, and used as the target markers. The transcriptome of the ovarian tissue of four different tumors, containing tumor presentation and recurrence stages, were sequenced using the Illumina GAII platform. Approximately 85.97% markers were expressed of the total 563 markers, and the majority of them were involved in pathways associated with cancer, signaling and infection. In total, 85 markers were found to be aberrantly expressed in tumor cells from patients with ovarian cancer who had recurrences, including 33 upregulated markers at the recurrence stage. Therefore, they may have roles ovarian tumor due to their aberrant expression. Differentially expressed markers and the associations between them can be assessed by examining the RNA-seq data. These findings may provide novel information for further studies on ovarian cancer. [ABSTRACT FROM AUTHOR]

Background: Allelic specific expression (ASE) increases our understanding of the genetic control of gene expression and its links to phenotypic variation. ASE testing is implemented through binomial or beta-binomial tests of sequence read counts of alternative alleles at a cSNP of interest in heterozygous individuals. This requires prior ascertainment of the cSNP genotypes for all individuals. To meet the needs, we propose hidden Markov methods to call SNPs from next generation RNA sequence data when ASE possibly exists. Results: We propose two hidden Markov models (HMMs), HMM-ASE and HMM-NASE that consider or do not consider ASE, respectively, in order to improve genotyping accuracy. Both HMMs have the advantages of calling the genotypes of several SNPs simultaneously and allow mapping error which, respectively, utilize the dependence among SNPs and correct the bias due to mapping error. In addition, HMM-ASE exploits ASE information to further improve genotype accuracy when the ASE is likely to be present. Simulation results indicate that the HMMs proposed demonstrate a very good prediction accuracy in terms of controlling both the false discovery rate (FDR) and the false negative rate (FNR). When ASE is present, the HMM-ASE had a lower FNR than HMM-NASE, while both can control the false discovery rate (FDR) at a similar level. By exploiting linkage disequilibrium (LD), a real data application demonstrate that the proposed methods have better sensitivity and similar FDR in calling heterozygous SNPs than the VarScan method. Sensitivity and FDR are similar to that of the BCFtools and Beagle methods. The resulting genotypes show good properties for the estimation of the genetic parameters and ASE ratios. Conclusions: We introduce HMMs, which are able to exploit LD and account for the ASE and mapping errors, to simultaneously call SNPs from the next generation RNA sequence data. The method introduced can reliably call for cSNP genotypes even in the presence of ASE and under low sequencing coverage. As a byproduct, the proposed method is able to provide predictions of ASE ratios for the heterozygous genotypes, which can then be used for ASE testing. [ABSTRACT FROM AUTHOR]

Background: Allelic specific expression (ASE) increases our understanding of the genetic control of gene expression and its links to phenotypic variation. ASE testing is implemented through binomial or beta-binomial tests of sequence read counts of alternative alleles at a cSNP of interest in heterozygous individuals. This requires prior ascertainment of the cSNP genotypes for all individuals. To meet the needs, we propose hidden Markov methods to call SNPs from next generation RNA sequence data when ASE possibly exists. Results: We propose two hidden Markov models (HMMs), HMM-ASE and HMM-NASE that consider or do not consider ASE, respectively, in order to improve genotyping accuracy. Both HMMs have the advantages of calling the genotypes of several SNPs simultaneously and allow mapping error which, respectively, utilize the dependence among SNPs and correct the bias due to mapping error. In addition, HMM-ASE exploits ASE information to further improve genotype accuracy when the ASE is likely to be present. Simulation results indicate that the HMMs proposed demonstrate a very good prediction accuracy in terms of controlling both the false discovery rate (FDR) and the false negative rate (FNR). When ASE is present, the HMM-ASE had a lower FNR than HMM-NASE, while both can control the false discovery rate (FDR) at a similar level. By exploiting linkage disequilibrium (LD), a real data application demonstrate that the proposed methods have better sensitivity and similar FDR in calling heterozygous SNPs than the VarScan method. Sensitivity and FDR are similar to that of the BCFtools and Beagle methods. The resulting genotypes show good properties for the estimation of the genetic parameters and ASE ratios. Conclusions: We introduce HMMs, which are able to exploit LD and account for the ASE and mapping errors, to simultaneously call SNPs from the next generation RNA sequence data. The method introduced can reliably call for cSNP genotypes even in the presence of ASE and under low sequencing coverage. As a byproduct, the proposed method is able to provide predictions of ASE ratios for the heterozygous genotypes, which can then be used for ASE testing. [ABSTRACT FROM AUTHOR]

Differentially-expressed genes (DEGs) correlated to oral contraceptives (OCs) were identified by comparing the transcriptomes of cervical cancer patients who have taken OCs and those who have not. Their biological functions and relevance to clinical manifestations were investigated further in order to gain an understanding of the pathogenesis of cervical cancer and provide potential therapeutic targets. Level 3 RNA-sequencing (seq) data for cervical squamous-cell carcinoma and endocervical adenocarcinoma and the clinical information were downloaded from The Cancer Genome Atlas. The present study analyzed the RNA-seq data and information on OC use of 35 patients [OC users (n=18) and those who have never used OCs (n=17)]. Student's t-test was used in order to identify DEGs and the false discovery rate (FDR) was estimated by a Beta-Uniform Mixture model, which was adopted in multiple testing corrections. A functional enrichment analysis was performed with the Database for Annotation, Visualization and Integrated Discovery tool and BioCarta. A total of 80 DEGs were identified in OC users while FDR=0.3 was set as the cut-off value. The metabolic process and human telomerase RNA gene transcription were significantly upregulated in DEGs. Furthermore, secreted LY6/PLAUR domain containing 1 was identified to be correlated to the pathological response, while the synapse defective 1 Rho GTPase homolog 2 was found to be significantly associated with the histological grade and overall survival time. In conclusion, present study shed light on the effect of OC use on the oncogenesis of the cervix and may indicate novel approaches for a targeted therapy of cervical cancer. [ABSTRACT FROM AUTHOR]