Search

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]

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]

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]

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]

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]

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]

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]

Colorectal cancer (CRC) is one of the leading causes of mortality and morbidity in the world. Being a heterogeneous disease, cancer therapy and prognosis represent a significant challenge to medical care. The molecular information improves the accuracy with which patients are classified and treated since similar pathologies may show different clinical outcomes and other responses to treatment. However, the high dimensionality of gene expression data makes the selection of novel genes a problematic task. We propose TCox, a novel penalization function for Cox models, which promotes the selection of genes that have distinct correlation patterns in normal vs. tumor tissues. We compare TCox to other regularized survival models, Elastic Net, HubCox, and OrphanCox. Gene expression and clinical data of CRC and normal (TCGA) patients are used for model evaluation. Each model is tested 100 times. Within a specific run, eighteen of the features selected by TCox are also selected by the other survival regression models tested, therefore undoubtedly being crucial players in the survival of colorectal cancer patients. Moreover, the TCox model exclusively selects genes able to categorize patients into significant risk groups. Our work demonstrates the ability of the proposed weighted regularizer TCox to disclose novel molecular drivers in CRC survival by accounting for correlation-based network information from both tumor and normal tissue. The results presented support the relevance of network information for biomarker identification in high-dimensional gene expression data and foster new directions for the development of network-based feature selection methods in precision oncology. [ABSTRACT FROM AUTHOR]

Due to the complexity of the pathological mechanisms of neurodegenerative diseases, traditional differentially-expressed gene selection methods cannot detect disease-associated genes accurately. Recent studies have shown that consensus-guided unsupervised feature selection (CGUFS) performs well in feature selection for identifying disease-associated genes. Since the random initialization of the feature selection matrix in CGUFS results in instability of the final disease-associated gene set, for the purposes of this study we proposed an ensemble method based on CGUFS—namely, ensemble consensus-guided unsupervised feature selection (ECGUFS) in order to further improve the accuracy of disease-associated genes and the stability of feature gene sets. We also proposed a bagging integration strategy to integrate the results of CGUFS. Lastly, we conducted experiments with Huntington’s disease RNA sequencing (RNA-Seq) data and obtained the final feature gene set, where we detected 287 disease-associated genes. Enrichment analysis on these genes has shown that postsynaptic density and the postsynaptic membrane, synapse, and cell junction are all affected during the disease’s progression. However, ECGUFS greatly improved the accuracy of disease-associated gene prediction and the stability of the disease-associated gene set. We conducted a classification of samples with labels based on the linear support vector machine with 10-fold cross-validation. The average accuracy is 0.9, which suggests the effectiveness of the feature gene set. [ABSTRACT FROM AUTHOR]