High-throughput Expression Data Research Articles

Co-regulatory network analysis of the main secondary metabolite (SM) biosynthesis in Crocus sativus L.

Saffron (Crocus sativus L.) is being embraced as the most important medicinal plant and the commercial source of saffron spice. Despite the beneficial economic and medicinal properties of saffron, the regulatory mechanism of the correlation of TFs and genes related to the biosynthesis of the apocarotenoids pathway is less obvious. Realizing these regulatory hierarchies of gene expression networks related to secondary metabolites production events is the main challenge owing to the complex and extensive interactions between the genetic behaviors. Recently, high throughput expression data have been highly feasible for constructing co-regulation networks to reveal the regulated processes and identifying novel candidate hub genes in response to complex processes of the biosynthesis of secondary metabolites. Herein, we performed Weighted Gene Co-expression Network Analysis (WGCNA), a systems biology method, to identify 11 regulated modules and hub TFs related to secondary metabolites. Three specialized modules were found in the apocarotenoids pathway. Several hub TFs were identified in notable modules, including MADS, C2H2, ERF, bZIP, HD-ZIP, and zinc finger protein MYB and HB, which were potentially associated with apocarotenoid biosynthesis. Furthermore, the expression levels of six hub TFs and six co-regulated genes of apocarotenoids were validated with RT-qPCR. The results confirmed that hub TFs specially MADS, C2H2, and ERF had a high correlation (P < 0.05) and a positive effect on genes under their control in apocarotenoid biosynthesis (CCD2, GLT2, and ADH) among different C. sativus ecotypes in which the metabolite contents were assayed. Promoter analysis of the co-expressed genes of the modules involved in apocarotenoids biosynthesis pathway suggested that not only are the genes co-expressed, but also share common regulatory motifs specially related to hub TFs of each module and that they may describe their common regulation. The result can be used to engineer valuable secondary metabolites of C. sativus by manipulating the hub regulatory TFs.

Keyword-based Scoring for The Prioritization of Differentially Expressed Genes Using Gene Ontology Annotations and Network Propagation

Computational gene prioritization provides a basis for the utilization of high-throughput expression data; however, methods for prioritization considering relevance to both biological processes and selected keywords are lacking. In this study, a method was developed to rank differentially expressed genes (DEGs) by utilizing the Gene Ontology (GO) database for keyword searches and propagating the results through a protein–protein interaction network. A scoring system that effectively biases the scores of genes relevant to given keywords for avoidance and preference was established. This scoring method was combined with scoring based on expression characteristic groups (ECGs) with network propagation to obtain a final combined score (cScore). The performance of the new approach was evaluated using a rat middle cerebral artery occlusion (MCAO) dataset, revealing that the method more effectively filtered out DEGs compared with conventional methods based on both significance and fold change values, excluding 76% of genes in average, while retaining genes of interest. Further improvements, including addressing the inability of downward accumulation to balance terms with conflicting hits and the limited number of databases utilized for scoring, can further improve the performance of the method. Overall, the newly developed method can improve the interpretation of DEG data.

Analysis and Prediction of Significant Genes in Gastric Cancer Based on TCGA Database.

Use TCGA data to screen the significantly expressed genes in gastric cancer,then analyze the biological functions of characteristic genes.To identify potential targets for the diagnosis, prognosis and therapeutic monitoring of gastric cancer. Download the high-throughput gastric cancer patient tissue mRNA expression data from the TCGA database, and used the edgeR software package to perform gene differences using the Rstudio software according to the |logFC|>1 and p<0.05 standards analysis, then performed enrichment analysis and constructed protein interaction network of the differently expressed genes to find the significant genes.Finally,analyzed the relationship between significant genes and gastric cancer patient prognosis. The RNA-seq expression data of 375 cases of gastric cancer and 32 adjacent tissues were obtained from the TCGA database. The differential analysis yielded 4320 differential genes, of which 2718 were highly expressed and 1602 were low expressed. Through enrichment analysis, it is found that the differential genes are mainly related to the neural active ligand-receptor interaction, calcium signaling pathway, Protein digestion and absorption, chemical carcinogenesis, and cytochrome P450 metabolism of foreign compounds. Further survival analysis screened out the five significant genes of ALB, AFP, IGFBP1, APOH, and TF are related to the prognosis of gastric cancer. The five significant genes screened by TCGA data may be related to the prognosis of gastric cancer and are expected to become prognostic biomarkers.

Hub genes and pathways in gastric cancer: A comparison between studies that used normal tissues adjacent to the tumour and studies that used healthy tissues as calibrator.

Several bioinformatics studies have been performed on high-throughput expression data to determine the cellular pathways and hub genes affected by Gastric cancer (GC). However, these studies differ in using a healthy tissue or normal tissue adjacent to the tumour (NAT) as calibrator tissues. This study was designed to find how using healthy or NAT tissues as calibrator tissues could affect pathway enrichment data and hub genes in GC. Two gene expression datasets with NAT tissues (GSE79973 and GSE118916) and one dataset with healthy tissues (GSE54129) were downloaded and processed by the limma package to screen the differentially expressed genes (DEGs). Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) enrichment analysis were performed by the Enrichr online tool. Protein-protein interaction network construction, module analysis, and hub genes selection were performed by Cytoscape software, Molecular Complex Detection plugin, and cytoHubba plugin, respectively. The gene expression profiling interactive analysis web server was used to analyse RNA sequencing expression data from The Cancer Genome Atlas Program. The Kaplan-Meier plotter was used to perform survival analysis. Our results showed that some KEGG and GO pathways were shared between studies with NAT and the study with healthy tissues. However, some terms, especially inflammation-related terms, were missed when NAT tissues were used as calibrator tissues. Also, only FN1 and COL1A1 are common hub genes between DEGs of the studies with NAT and healthy tissues. Since hub genes are usually extracted and suggested as candidate targets for GC diagnosis, prognosis, or treatment, selecting healthy or NAT tissues may affect the hub genes selection.

In silico identification of potential miRNAs -mRNA inflammatory networks implicated in the pathogenesis of COVID-19

COVID-19 has been found to affect the expression profile of several mRNAs and miRNAs, leading to dysregulation of a number of signaling pathways, particularly those related to inflammatory responses. In the current study, a systematic biology procedure was used for the analysis of high-throughput expression data from blood specimens of COVID-19 and healthy individuals. Differentially expressed miRNAs in blood specimens of COVID-19 vs. healthy specimens were then identified to construct and analyze miRNA–mRNA networks and predict key miRNAs and genes in inflammatory pathways. Our results showed that 171 miRNAs were expressed as outliers in box plot and located in the critical areas according to our statistical analysis. Among them, 8 miRNAs, namely miR-1275, miR-4429, miR-4489, miR-6721-5p, miR-5010-5p, miR-7110-5p, miR-6804-5p and miR-6881–3p were found to affect expression of key genes in NF-KB, JAK/STAT and MAPK signaling pathways implicated in COVID-19 pathogenesis. In addition, our results predicted that 25 genes involved in above-mentioned inflammatory pathways were targeted not only by these 8 miRNAs but also by other obtained miRNAs (163 miRNAs). The results of the current in silico study represent candidate targets for further studies in COVID-19.

Akut ve Kronik Miyeloid Lösemili Hastalarda Prognostik miRNA İmzasının Biyoinformatik Analiz ile Karşılaştırılması

Aim: In this study, differentially expressed miRNA profiles were determined using high-throughput expression data from samples of AML and CML patients to identify miRNAs involved in the therapeutic response.&#x0D; Material and Methods: miRNA microarray datasets GSE142699 and GSE90773 were downloaded via the GEO database and analysis was performed with the online analysis tool GEO2R. Data no. GSE142699 was made with 24 control and 24 newly diagnosed AML patients, data no. GSE90773 was made with 8 control and 10 newly diagnosed CML patients. After the analysis, they were grouped according to fold change (FC) values and p

Clustering Ensemble Technology Based on Granular Computing to Extract Cervical Cancer Predictors

Background. Cervical cancer is the most common gynecological malignancy, and its incidence has tended to be younger in recent years. Through the analysis of high-throughput expression data, the identification of key genes in cancer and healthy individuals as predictors of cervical cancer is of great significance for the early detection and early treatment of cervical cancer. Method. Granular computing is a concept and computing paradigm to deal with problems through information granulation, and the process of granulation can be realized by means of clustering. Based on this, this paper proposes an AB method to obtain representative elements in a multiattribute data system. First, the evaluation index FHEI of the clustering structure is introduced, and Algorithm 1 is designed to obtain the optimal clustering structure of each attribute of the data system and use it as the base cluster. Secondly, based on the clustering ensemble technology of granular computing, Algorithm 2 is designed with the help of the concept of information entropy. The algorithm takes the base cluster as the input to obtain the optimal ensemble clustering structure. Finally, using the nearest center principle, the representative elements of each class in the optimal ensemble clustering structure are obtained. Results. In this paper, the differentially expressed genes (DEGs) are screened out by using the gene expression data of cervical cancer, and the scores of the four interaction relationships among the DEGs are used as a multiattribute data system and input into the AB method. The five representative elements obtained are RTTN, SAMD10, ZNF207, WAC, and METTL14, which are the predictors of cervical cancer. The classification accuracy of these predictors is as high as 98.82%. This paper also conducts a comparative study between the AB method and other classical methods on six independent gene expression datasets. The results show that the number of predictors obtained by the AB method is small but has a high classification accuracy in the classification of patient samples.

DeepReGraph co-clusters temporal gene expression and cis-regulatory elements through heterogeneous graph representation learning

This work presents DeepReGraph, a novel method for co-clustering genes and cis-regulatory elements (CREs) into candidate regulatory networks. Gene expression data, as well as data from three CRE activity markers from a publicly available dataset of mouse fetal heart tissue, were used for DeepReGraph concept proofing. In this study we used open chromatin accessibility from ATAC-seq experiments, as well as H3K27ac and H3K27me3 histone marks as CREs activity markers. However, this method can be executed with other sets of markers. We modelled all data sources as a heterogeneous graph and adapted a state-of-the-art representation learning algorithm to produce a low-dimensional and easy-to-cluster embedding of genes and CREs. Deep graph auto-encoders and an adaptive-sparsity generative model are the algorithmic core of DeepReGraph. The main contribution of our work is the design of proper combination rules for the heterogeneous gene expression and CRE activity data and the computational encoding of well-known gene expression regulatory mechanisms into a suitable objective function for graph embedding. We showed that the co-clusters of genes and CREs in the final embedding shed light on developmental regulatory mechanisms in mouse fetal-heart tissue. Such clustering could not be achieved by using only gene expression data. Function enrichment analysis proves that the genes in the co-clusters are involved in distinct biological processes. The enriched transcription factor binding sites in CREs prioritize the candidate transcript factors which drive the temporal changes in gene expression. Consequently, we conclude that DeepReGraph could foster hypothesis-driven tissue development research from high-throughput expression and epigenomic data. Full source code and data are available on the DeepReGraph GitHub project.

Boolean factor graph model for biological systems: the yeast cell-cycle network

BackgroundThe desire to understand genomic functions and the behavior of complex gene regulatory networks has recently been a major research focus in systems biology. As a result, a plethora of computational and modeling tools have been proposed to identify and infer interactions among biological entities. Here, we consider the general question of the effect of perturbation on the global dynamical network behavior as well as error propagation in biological networks to incite research pertaining to intervention strategies.ResultsThis paper introduces a computational framework that combines the formulation of Boolean networks and factor graphs to explore the global dynamical features of biological systems. A message-passing algorithm is proposed for this formalism to evolve network states as messages in the graph. In addition, the mathematical formulation allows us to describe the dynamics and behavior of error propagation in gene regulatory networks by conducting a density evolution (DE) analysis. The model is applied to assess the network state progression and the impact of gene deletion in the budding yeast cell cycle. Simulation results show that our model predictions match published experimental data. Also, our findings reveal that the sample yeast cell-cycle network is not only robust but also consistent with real high-throughput expression data. Finally, our DE analysis serves as a tool to find the optimal values of network parameters for resilience against perturbations, especially in the inference of genetic graphs.ConclusionOur computational framework provides a useful graphical model and analytical tools to study biological networks. It can be a powerful tool to predict the consequences of gene deletions before conducting wet bench experiments because it proves to be a quick route to predicting biologically relevant dynamic properties without tunable kinetic parameters.

Identification of RCC Subtype-Specific microRNAs-Meta-Analysis of High-Throughput RCC Tumor microRNA Expression Data.

Simple SummaryIn the majority of renal cancer cases, the disease course is non-symptomatic which frequently leads to late diagnosis of disease. Currently, there are no molecular tools dedicated to the detection and monitoring of renal cancer. Our study aimed to investigate changes in microRNA (miRNA) expression in tissue samples of renal cancer patients. We performed meta-analysis using results of 14 high-throughput studies (both, NGS and microarrays) and as a result, selected a group of miRNAs deregulated in renal cancer and its subtypes. Later, the expression changes of selected miRNA were validated in an independent sample set. We confirmed that the investigation of miRNA expression might be potentially applicable in the detection and monitoring of renal cancer and its subtypes.Renal cell carcinoma (RCC) is one of the most common cancers worldwide with a nearly non-symptomatic course until the advanced stages of the disease. RCC can be distinguished into three subtypes: papillary (pRCC), chromophobe (chRCC) and clear cell renal cell carcinoma (ccRCC) representing up to 75% of all RCC cases. Detection and RCC monitoring tools are limited to standard imaging techniques, in combination with non-RCC specific morphological and biochemical read-outs. RCC subtype identification relays mainly on results of pathological examination of tumor slides. Molecular, clinically applicable and ideally non-invasive tools aiding RCC management are still non-existent, although molecular characterization of RCC is relatively advanced. Hence, many research efforts concentrate on the identification of molecular markers that will assist with RCC sub-classification and monitoring. Due to stability and tissue-specificity miRNAs are promising candidates for such biomarkers. Here, we performed a meta-analysis study, utilized seven NGS and seven microarray RCC studies in order to identify subtype-specific expression of miRNAs. We concentrated on potentially oncocytoma-specific miRNAs (miRNA-424-5p, miRNA-146b-5p, miRNA-183-5p, miRNA-218-5p), pRCC-specific (miRNA-127-3p, miRNA-139-5p) and ccRCC-specific miRNAs (miRNA-200c-3p, miRNA-362-5p, miRNA-363-3p and miRNA-204-5p, 21-5p, miRNA-224-5p, miRNA-155-5p, miRNA-210-3p) and validated their expression in an independent sample set. Additionally, we found ccRCC-specific miRNAs to be differentially expressed in ccRCC tumor according to Fuhrman grades and identified alterations in their isoform composition in tumor tissue. Our results revealed that changes in the expression of selected miRNA might be potentially utilized as a tool aiding ccRCC subclass discrimination and we propose a miRNA panel aiding RCC subtype distinction.

Machine learning integrated ensemble of feature selection methods followed by survival analysis for predicting breast cancer subtype specific miRNA biomarkers

Breast cancer is the second leading cancer type among females. In this regard, it is found that microRNAs play an important role by regulating the gene expressions at the post-transcriptional phase. However, identification of the most influencing miRNAs in breast cancer subtypes is a challenging task, while the recent advancement in Next Generation Sequencing techniques allows analyzing high throughput expression data of miRNAs. Thus, we have conducted this research with the help of NGS data of breast cancer in order to identify the most significant miRNA biomarkers. The selected miRNA biomarkers are highly associated with the multiple breast cancer subtypes. For this purpose, a two-phase technique, called Machine Learning Integrated Ensemble of Feature Selection Methods, followed by survival analysis, is proposed. In the first phase, we have selected the best among seven machine learning techniques based on classification accuracy using the entire set of features (in this case miRNAs). Subsequently, eight different feature selection methods are used separately in order to rank the features and validate each set of top features using the selected machine learning technique by considering a multi-class classification task of the breast cancer subtypes. In the second phase, based on the classification accuracy values, the top features from each feature selection method are considered to make an ensemble to provide further categorization of the miRNAs as 8*, 7* up to 1*. The 8* miRNAs provide the highest average classification accuracy of 86% after 10-fold cross-validation. Thereafter, 27 miRNAs are identified from the list that is confined within 8* to 4* miRNAs based on their importance in survival for breast cancer subtypes using Cox regression based survival analysis. Moreover, expression analysis, regulatory network analysis, protein-protein interaction analysis, KEGG pathway and gene ontology enrichment analysis are performed in order to validate biological significance of the proposed solution. Additionally, we have prepared a miRNA-protein-drug interaction network to identify possible drug for the selected miRNAs. Thus, our findings may be considered during a clinical trial for the treatment of breast cancer patients.

Mining Protein Expression Databases Using Network Meta-Analysis.

Public databases featuring original, raw data from "Omics" experiments enable researchers to perform meta-analyses by combining either the raw data or the summarized results of several independent studies. In proteomics, high-throughput protein expression data is measured by diverse techniques such as mass spectrometry, 2-D gel electrophoresis or protein arrays yielding data of different scales. Therefore, direct data merging can be problematic, and combining the summarized data of the individual studies can be advantageous. A special form of meta-analysis is network meta-analysis, where studies with different settings of experimental groups can be combined. However, all studies must be linked by one experimental group that has to appear in each study. Usually that is the control group. Then, a study network is formed and indirect statistical inferences can also be made between study groups that appear not in each of the studies.In this chapter, we describe the working principle of and available software for network meta-analysis. The applicability to high-throughput protein expression data is demonstrated in an example from breast cancer research. We also describe the special challenges when applying this method.

A Metabolic Gene Signature to Predict Overall Survival in Head and Neck Squamous Cell Carcinoma.

Background Head and neck squamous cell carcinoma (HNSCC) is a common malignancy that emanates from the lips, mouth, paranasal sinuses, oropharynx, larynx, nasopharynx, and from other pharyngeal cancers. The availability of high-throughput expression data has made it possible to use global gene expression data to analyze the relationship between metabolic-related gene expression and clinical outcomes in HNSCC patients. Method In this study, we used RNA sequencing (RNA-seq) data from the cancer genome atlas (TCGA), with validation in the GEO dataset to profile the metabolic microenvironment and define potential biomarkers for metabolic therapy. Results We extracted data for 529 patients and 327 metabolic genes (198 upregulated and 129 downregulated genes) in the TCGA database. Carbonic anhydrase 9 (CA9) and CA6 had the largest logFCs in the upregulated and downregulated genes, respectively. Our Cox regression model data showed 51 prognostic-related genes with lysocardiolipin acyltransferase 1 (LCLAT1) and choline dehydrogenase (CHDH) being associated with the highest risk (HR = 1.144, 95% CI = 1.044 ~ 1.251) and the lowest risk (HR = 0.580, 95% CI = 0.400 ~ 0.839) in HNSCC, respectively. We next used the ROC curve to evaluate whether the differentially expressed metabolic-related genes could serve as early predictors of HNSCC. The findings showed an AUC of 0.745 and 0.618 in the TCGA and GEO analysis, respectively. Besides, the ability for the genes to predict clinicopathological HNSCC status was analyzed and the data showed that the AUC for age, gender, grade, stage, T, M, and N was 0.520, 0.495, 0.568, 0.606, 0.577, 0.476, and 0.673, respectively, in the TCGA dataset. On the other hand, the AUC for age, gender, stage, T, M, N, smoking, and HPV16-pos was 0.599, 0.531, 0.622, 0.606, 0.616, 0.550, 0.614, 0.519, and 0.397, respectively, in the GEO dataset. Conclusion Taken together, our study unearths a novel metabolic gene signature for the prediction of HNSCC prognosis based on the TCGA dataset. Our signature might point out the metabolic microenvironment disorders and provides potential treatment targets and prognostic biomarkers.

GCNG: graph convolutional networks for inferring gene interaction from spatial transcriptomics data

Most methods for inferring gene-gene interactions from expression data focus on intracellular interactions. The availability of high-throughput spatial expression data opens the door to methods that can infer such interactions both within and between cells. To achieve this, we developed Graph Convolutional Neural networks for Genes (GCNG). GCNG encodes the spatial information as a graph and combines it with expression data using supervised training. GCNG improves upon prior methods used to analyze spatial transcriptomics data and can propose novel pairs of extracellular interacting genes. The output of GCNG can also be used for downstream analysis including functional gene assignment.Supporting website with software and data: https://github.com/xiaoyeye/GCNG.

Bioinformatic Analysis Identifies Potential Key Genes in the Pathogenesis of Melanoma

Melanoma is the deadliest skin tumor and is prone to distant metastases. The incidence of melanoma has increased rapidly in the past few decades, and current trends indicate that this growth is continuing. This study was aimed to explore the molecular mechanisms of melanoma pathogenesis and discover underlying pathways and genes associated with melanoma. We used high-throughput expression data to study differential expression profiles of related genes in melanoma. The differentially expressed genes (DEGs) of melanoma in GSE15605, GSE46517, GSE7553, and the Cancer Genome Atlas (TCGA) datasets were analyzed. Differentially expressed genes (DEGs) were identified by paired t-test. Then the DEGs were performed cluster and principal component analyses and protein–protein interaction (PPI) network construction. After that, we analyzed the differential genes through bioinformatics and got hub genes. Finally, the expression of hub genes was confirmed in the TCGA databases and collected patient tissue samples. Total 144 up-regulated DEGs and 16 down-regulated DEGs were identified. A total of 17 gene ontology analysis (GO) terms and 11 pathways were closely related to melanoma. Pathway of pathways in cancer was enriched in 8 DEGs, such as junction plakoglobin (JUP) and epidermal growth factor receptor (EGFR). In the PPI networks, 9 hub genes were obtained, such as loricrin (LOR), filaggrin (FLG), keratin 5 (KRT5), corneodesmosin (CDSN), desmoglein 1 (DSG1), desmoglein 3 (DSG3), keratin 1 (KRT1), involucrin (IVL), and EGFR. The pathway of pathways in cancer and its enriched DEGs may play important roles in the process of melanoma. The hub genes of DEGs may become promising melanoma candidate genes. Five key genes FLG, DSG1, DSG3, IVL, and EGFR were identified in the TCGA database and melanoma tissues. The results suggested that FLG, DSG1, DSG3, IVL, and EGFR might play important roles and potentially be valuable in the prognosis and treatment of melanoma. These hub genes might well have clinical significance as diagnostic markers.

Non-Coding RNAs in the Brain-Heart Axis: The Case of Parkinson's Disease.

Parkinson’s disease (PD) is a complex and heterogeneous disorder involving multiple genetic and environmental influences. Although a wide range of PD risk factors and clinical markers for the symptomatic motor stage of the disease have been identified, there are still no reliable biomarkers available for the early pre-motor phase of PD and for predicting disease progression. High-throughput RNA-based biomarker profiling and modeling may provide a means to exploit the joint information content from a multitude of markers to derive diagnostic and prognostic signatures. In the field of PD biomarker research, currently, no clinically validated RNA-based biomarker models are available, but previous studies reported several significantly disease-associated changes in RNA abundances and activities in multiple human tissues and body fluids. Here, we review the current knowledge of the regulation and function of non-coding RNAs in PD, focusing on microRNAs, long non-coding RNAs, and circular RNAs. Since there is growing evidence for functional interactions between the heart and the brain, we discuss the benefits of studying the role of non-coding RNAs in organ interactions when deciphering the complex regulatory networks involved in PD progression. We finally review important concepts of harmonization and curation of high throughput datasets, and we discuss the potential of systems biomedicine to derive and evaluate RNA biomarker signatures from high-throughput expression data.

ICDA-CGR: Identification of circRNA-disease associations based on Chaos Game Representation.

Found in recent research, tumor cell invasion, proliferation, or other biological processes are controlled by circular RNA. Understanding the association between circRNAs and diseases is an important way to explore the pathogenesis of complex diseases and promote disease-targeted therapy. Most methods, such as k-mer and PSSM, based on the analysis of high-throughput expression data have the tendency to think functionally similar nucleic acid lack direct linear homology regardless of positional information and only quantify nonlinear sequence relationships. However, in many complex diseases, the sequence nonlinear relationship between the pathogenic nucleic acid and ordinary nucleic acid is not much different. Therefore, the analysis of positional information expression can help to predict the complex associations between circRNA and disease. To fill up this gap, we propose a new method, named iCDA-CGR, to predict the circRNA-disease associations. In particular, we introduce circRNA sequence information and quantifies the sequence nonlinear relationship of circRNA by Chaos Game Representation (CGR) technology based on the biological sequence position information for the first time in the circRNA-disease prediction model. In the cross-validation experiment, our method achieved 0.8533 AUC, which was significantly higher than other existing methods. In the validation of independent data sets including circ2Disease, circRNADisease and CRDD, the prediction accuracy of iCDA-CGR reached 95.18%, 90.64% and 95.89%. Moreover, in the case studies, 19 of the top 30 circRNA-disease associations predicted by iCDA-CGR on circRDisease dataset were confirmed by newly published literature. These results demonstrated that iCDA-CGR has outstanding robustness and stability, and can provide highly credible candidates for biological experiments.

Phylogenomic Analysis of R2R3 MYB Transcription Factors in Sorghum and their Role in Conditioning Biofuel Syndrome.

Background Large scale cultivation of sorghum for food, feed, and biofuel requires concerted efforts for engineering multipurpose cultivars with optimised agronomic traits. Due to their vital role in regulating the biosynthesis of phenylpropanoid-derived compounds, biomass composition, biotic, and abiotic stress response, R2R3-MYB family transcription factors are ideal targets for improving environmental resilience and economic value of sorghum.Methods We used diverse computational biology tools to survey the sorghum genome to identify R2R3-MYB transcription factors followed by their structural and phylogenomic analysis. We used in-house generated as well as publicly available high throughput expression data to analyse the R2R3 expression patterns in various sorghum tissue types.Results We have identified a total of 134 R2R3-MYB genes from sorghum and developed a framework to predict gene functions. Collating information from the physical location, duplication, structural analysis, orthologous sequences, phylogeny, and expression patterns revealed the role of duplications in clade-wise expansion of the R2R3-MYB family as well as intra-clade functional diversification. Using publicly available and in-house generated RNA sequencing data, we provide MYB candidates for conditioning biofuel syndrome by engineering phenylpropanoid biosynthesis and sugar signalling pathways in sorghum.Conclusion The results presented here are pivotal to prioritize MYB genes for functional validation and optimize agronomic traits in sorghum.

Integrating high-throughput microRNA and mRNA expression data to identify risk mRNA signature for pancreatic cancer prognosis.

Pancreatic cancer is a malignancy of the digestive system characterized by poor prognosis. A number of prognostic messenger RNA (mRNA) signatures have been identified by using the high-throughput expression profiles. MicroRNAs (miRNA) play a critical role in regulating multiple cellular functions. However, no such integrated analysis of miRNAs and mRNAs for studying the prognostic mechanisms of pancreatic cancer has been reported. In this study, we first identified prognostic mRNAs and miRNAs based on The Cancer Genome Atlas datasets, and then performed an enrichment analysis to explore the underlying biological mechanisms involved in pancreatic cancer prognosis at the mRNA level. Furthermore, we performed an integrated analysis of mRNAs and miRNAs to identify prognostic subpathways, which were closely associated with pancreatic cancer genes and tumor hallmarks and involved in hypoxia, oxidative phosphyorylation and xenobiotic metabolisms. Meanwhile, we performed a random walk algorithm based on global network, prognostic mRNAs and miRNAs, and identified top risk mRNAs as the prognostic signature. Finally, an independent testing set was used to confirm the predictive power of the top mRNA signature, and most of these genes involved were known oncogenes. In conclusion, we performed a series of integrated analyses by comprehensively exploring pancreatic cancer prognosis and systematically optimized the prognostic signature for clinical use.

Normalization Methods for the Analysis of Unbalanced Transcriptome Data: A Review.

Dozens of normalization methods for correcting experimental variation and bias in high-throughput expression data have been developed during the last two decades. Up to 23 methods among them consider the skewness of expression data between sample states, which are even more than the conventional methods, such as loess and quantile. From the perspective of reference selection, we classified the normalization methods for skewed expression data into three categories, data-driven reference, foreign reference, and entire gene set. We separately introduced and summarized these normalization methods designed for gene expression data with global shift between compared conditions, including both microarray and RNA-seq, based on the reference selection strategies. To our best knowledge, this is the most comprehensive review of available preprocessing algorithms for the unbalanced transcriptome data. The anatomy and summarization of these methods shed light on the understanding and appropriate application of preprocessing methods.