The exploration of disease-specific gene regulatory networks in esophageal carcinoma and stomach adenocarcinoma

Learning in Transcriptional Network Models: Computational Discovery of Pathway-Level Memory and Effective Interventions

Trainability, in any substrate, refers to the ability to change future behavior based on past experiences. An understanding of such capacity within biological cells and tissues would enable a particularly powerful set of methods for prediction and control of their behavior through specific patterns of stimuli. This top-down mode of control (as an alternative to bottom-up modification of hardware) has been extensively exploited by computer science and the behavioral sciences; in biology however, it is usually reserved for organism-level behavior in animals with brains, such as training animals towards a desired response. Exciting work in the field of basal cognition has begun to reveal degrees and forms of unconventional memory in non-neural tissues and even in subcellular biochemical dynamics. Here, we characterize biological gene regulatory circuit models and protein pathways and find them capable of several different kinds of memory. We extend prior results on learning in binary transcriptional networks to continuous models and identify specific interventions (regimes of stimulation, as opposed to network rewiring) that abolish undesirable network behavior such as drug pharmacoresistance and drug sensitization. We also explore the stability of created memories by assessing their long-term behavior and find that most memories do not decay over long time periods. Additionally, we find that the memory properties are quite robust to noise; surprisingly, in many cases noise actually increases memory potential. We examine various network properties associated with these behaviors and find that no one network property is indicative of memory. Random networks do not show similar memory behavior as models of biological processes, indicating that generic network dynamics are not solely responsible for trainability. Rational control of dynamic pathway function using stimuli derived from computational models opens the door to empirical studies of proto-cognitive capacities in unconventional embodiments and suggests numerous possible applications in biomedicine, where behavior shaping of pathway responses stand as a potential alternative to gene therapy.

A Complex Heterogeneous Network Model of Disease Regulated by Noncoding RNAs: A Case Study of Unstable Angina Pectoris

MicroRNAs (miRNAs) are important types of noncoding RNAs, and there is a lack of holistic and systematic understanding of the functions they play in disease. We proposed a research strategy, including two parts network analysis and network modelling, to analyze, model, and predict the regulatory network of miRNAs from a network perspective, using unstable angina pectoris as an example. In the network analysis section, we proposed the WGCNA & SimCluster method using both correlation and similarity to find hub miRNAs, and validation on two datasets showed better results than the methods using correlation or similarity alone. In the network modelling section, we used six knowledge graph or graph neural network models for link prediction of three types of edges and multilabel classification of two types of nodes. Comparative experiments showed that the RotatE model was a good model for link prediction, while the RGCN model was the best model for multilabel classification. Potential target genes were predicted for hub miRNAs and validation of hub miRNA-target gene interactions, target genes as biomarkers and target gene functions were performed using a three-step validation approach. In conclusion, our study provides a new strategy to analyze and model miRNA regulatory networks.

Next Generation Sequencing for Potential Regulated Genes and Micro-RNAs of Early Growth Response-1 in the Esophageal Squamous Cell Carcinoma

Esophageal cancer has a poor prognosis due to its aggressiveness and low survival rate. In Ease Asia, esophageal squamous cell carcinoma (ESCC) outnumbers esophageal adenocarcinoma (EAC). The ESCC patients still have high mortality despite modern surgical resection and neoadjuvant treatment. Determining patient and outcome prognostic factors is critical in ESCC treatment. In esophageal cancer, early growth response-1 (Egr-1) is a tumor suppressor gene, but the mechanism and associated genes are unknown. The study utilizes RNA interference method, the platform of Next Generation Sequencing (NGS) and bioinformatics analysis to investigate the influences after the Egr-1 gene slicing on the ESCC cells. The heat maps of differentially expressed mRNA and microRNAs were analyzed using the algorithm, Burrows-Wheller Aligner. The study showed that the expression of 51 mRNA and 26 microRNAs have significant changes in ESCC cells after Egr-1 knockdown. The KEGG enrichment analysis linked Egr-1-regulated genes and microRNAs. Egr-1 interactions with these genes and microRNAs may be important in tumor progression. In conclusions, this study provided the transcriptome patterns and relating pathway analysis for Egr-1 knockdown in ESCC cells. The mRNA and microRNAs altered by Egr-1 gene silencing might provide key information in the treatment of ESCC.

Identification of hub genes related to CD4+ memory T cell infiltration with gene co-expression network predicts prognosis and immunotherapy effect in colon adenocarcinoma

Background: CD4+ memory T cells (CD4+ MTCs), as an important part of the microenvironment affecting tumorigenesis and progression, have rarely been systematically analyzed. Our purpose was to comprehensively analyze the effect of CD4+ MTC infiltration on the prognosis of colon adenocarcinoma (COAD). Methods: Based on RNA-Seq data, weighted gene co-expression network analysis (WGCNA) was used to screen the CD4+ MTC infiltration genes most associated with colon cancer and then identify hub genes and construct a prognostic model using the least absolute shrinkage and selection operator algorithm (LASSO). Finally, survival analysis, immune efficacy analysis, and drug sensitivity analysis were performed to evaluate the role of the prognostic model in COAD. Results: We identified 929 differentially expressed genes (DEGs) associated with CD4+ MTCs and constructed a prognosis model based on five hub genes (F2RL2, TGFB2, DTNA, S1PR5, and MPP2) to predict overall survival (OS) in COAD. Kaplan-Meier analysis showed poor prognosis in the high-risk group, and the analysis of the hub gene showed that overexpression of TGFB2, DTNA, S1PR5, or MPP2 was associated with poor prognosis. Clinical prediction nomograms combining CD4+ MTC-related DEGs and clinical features were constructed to accurately predict OS and had high clinical application value. Immune efficacy and drug sensitivity analysis provide new insights for individualized treatment. Conclusion: We constructed a prognostic risk model to predict OS in COAD and analyzed the effects of risk score on immunotherapy efficacy or drug sensitivity. These studies have important clinical significance for individualized targeted therapy and prognosis.

Reconstruction of dynamic regulatory networks reveals signaling-induced topology changes associated with germ layer specification

Elucidating regulatory relationships between transcription factors (TFs) and target genes is fundamental to understanding how cells control their identity and behavior. Unfortunately, existing computational gene regulatory network (GRN) reconstruction methods are imprecise, computationally burdensome, and fail to reveal dynamic regulatory topologies. Here, we present Epoch, a reconstruction tool that uses single-cell transcriptomics to accurately infer dynamic networks. We apply Epoch to identify the dynamic networks underpinning directed differentiation of mouse embryonic stem cells (ESCs) guided by multiple signaling pathways, and we demonstrate that modulating these pathways drives topological changes that bias cell fate potential. We also find that Peg3 rewires the pluripotency network to favor mesoderm specification. By integrating signaling pathways with GRNs, we trace how Wnt activation and PI3K suppression govern mesoderm and endoderm specification, respectively. Finally, we identify regulatory circuits of patterning and axis formation that distinguish in vitro and in vivo mesoderm specification.

Prediction of Disease Genes Based on Stage-Specific Gene Regulatory Networks in Breast Cancer

Breast cancer is one of the most common malignant tumors in women, which seriously endangers women’s health. Great advances have been made over the last decades, however, most studies predict driver genes of breast cancer using biological experiments and/or computational methods, regardless of stage information. In this study, we propose a computational framework to predict the disease genes of breast cancer based on stage-specific gene regulatory networks. Firstly, we screen out differentially expressed genes and hypomethylated/hypermethylated genes by comparing tumor samples with corresponding normal samples. Secondly, we construct three stage-specific gene regulatory networks by integrating RNA-seq profiles and TF-target pairs, and apply WGCNA to detect modules from these networks. Subsequently, we perform network topological analysis and gene set enrichment analysis. Finally, the key genes of specific modules for each stage are screened as candidate disease genes. We obtain seven stage-specific modules, and identify 20, 12, and 22 key genes for three stages, respectively. Furthermore, 55%, 83%, and 64% of the genes are associated with breast cancer, for example E2F2, E2F8, TPX2, BUB1, and CKAP2L. So it may be of great importance for further verification by cancer experts.

Oncogenic microRNA-301b regulates tumor repressor dystrobrevin alpha to facilitate cell growth, invasion and migration in esophageal cancer

BACKGROUND

Esophageal cancer (EC) ranks the eighth in morbidity and the sixth in mortality around the whole world, which is an aggressive malignancy. To authenticate potential therapeutic targets for EC is therefore imperative. Although miR-301b might display changed expression in esophageal adenocarcinoma by utilizing Taqman miRNA profiling analysis, much less is known about the impact of miR-301b in EC.

METHODS AND RESULTS

By analyzing the data of 187 cancer tissues and 13 normal samples from TCGA database, we discovered that miR-301b was highly expressed in EC tissues. Then, RT-qPCR determined that miR-301b was up-regulated in EC cell lines (ECA109, JAR, TE-1 and OE33). Besides, miR-301b expression level was higher in ESCC cell line-TE-1 cells and lower in ESCC cell line-ECA109 cells compared to other EC cell lines. Hence, ECA109 cell line was used to up-regulate miR-301b expression while TE-1 cell line was applied to down-regulate miR-301b expression in the subsequent experiments. Additionally, OE33, as an ECA cell line, was applied to upregulate miR-301b expression to reflect the influence of miR-301b overexpression on EC progression. More interestingly, miR-301b appeared to act as a promoting effect on the proliferation of EC cells, which was tested by CCK8. Dystrobrevin alpha (DTNA) was a targeting gene of miR-301b, which was predicted by the websites of miRanda, miRWalk and TargetScan. Additionally, DTNA was low expressed in EC tissues and was an independent predictor of EC. Meanwhile, the low expression of DTNA was related to worse overall survival in EC patients. The Pearson correlation coefficient analyzed that DTNA expression was negatively correlated with miR-301b. Furthermore, RT-qPCR and western blotting assays ulteriorly indicated that DTNA was negatively modulated by miR-301b. The facilitating impact of miR-301b re-expression on ECA109 and OE33 cell growth, invasion and migration was receded by DTNA over-expression, whilst the repressive effect of miR-301b ablation on TE-1 cell growth, invasion and migration was inversed by DTNA silencing. Overexpression of miR-301b accelerated EC cell growth, migration and invasion through targeting DTNA.

CONCLUSIONS

Above all, we concluded that miR-301b was concerned with the progression of EC via regulating DTNA, suggesting that miR-301b and its target gene, DTNA, might serve as predictive biomarkers for EC therapy.

Gene regulatory networks exhibit several kinds of memory: quantification of memory in biological and random transcriptional networks

Gene regulatory networks (GRNs) process important information in developmental biology and biomedicine. A key knowledge gap concerns how their responses change over time. Hypothesizing long-term changes of dynamics induced by transient prior events, we created a computational framework for defining and identifying diverse types of memory in candidate GRNs. We show that GRNs from a wide range of model systems are predicted to possess several types of memory, including Pavlovian conditioning. Associative memory offers an alternative strategy for the biomedical use of powerful drugs with undesirable side effects, and a novel approach to understanding the variability and time-dependent changes of drug action. We find evidence of natural selection favoring GRN memory. Vertebrate GRNs overall exhibit more memory than invertebrate GRNs, and memory is most prevalent in differentiated metazoan cell networks compared with undifferentiated cells. Timed stimuli are a powerful alternative for biomedical control of complex in vivo dynamics without genomic editing or transgenes.

•

Gene regulatory networks' dynamics are modified by transient stimuli

•

GRNs have several different types of memory, including associative conditioning

•

Evolution favored GRN memory, and differentiated cells have the most memory capacity

•

Training GRNs offers a novel biomedical strategy not dependent on genetic rewiring

Prognostic prediction of carcinoma by a differential-regulatory-network-embedded deep neural network

The accurate prognostic prediction is essential for precise diagnosis and treatment of carcinoma. In addition to clinical survival prediction method, many computational methods based on transcriptomic data have been proposed to build the prediction models and study the prognosis of cancer patients. We propose a differential-regulatory-network-embedded deep neural network (DRE-DNN) method by integrating differential regulatory analysis based on gene co-expression network and deep neural network (DNN) method. From three public hepatocellular carcinoma (HCC) datasets, we derive differential regulatory network and embed regulatory information into DNN. By employing 1869 differential regulatory genes and survival data, we apply DRE-DNN to build a prediction model. We compare our method with the one which has all gene features in normal DNN, and results show that our method has better generalization ability and accuracy. We modify the normal DNN and develop an efficient method to predict prognosis of HCC from gene expression data. Our method decreases the inconsistence caused by the overfitting problem when the training sample size is small. DRE-DNN is also extendable for prognostic prediction of other cancers.

A 14-Methylation-Driven Differentially Expressed RNA as a Signature for Overall Survival Prediction in Patients with Uterine Corpus Endometrial Carcinoma

DNA methylation has been implicated as an important mechanism for the development of uterine corpus endometrial carcinoma (UCEC), indicating that methylation-driven genes may be potential biomarkers for survival prediction. In this study, we aimed to identify a new prognostic methylation signature for UCEC based on differentially expressed genes (DEGs) and long noncoding RNAs (lncRNAs) (DELs). Sample-matched RNA-sequencing and methylation-array data were downloaded from The Cancer Genome Atlas database, by analysis of which a total of 269 DEGs and 4 DELs were identified to be methylation driven. Least absolute shrinkage and selection operator analysis screened that 14 methylation-driven genes were significantly associated with overall survival (OS) and thus were used as a signature to establish a prognostic risk model. Based on the median threshold, the patients were divided into the low-risk and the high-risk groups, which showed significantly different survival periods under the Kaplan-Meier curve. The area under receiver operating characteristic curve (AUC) was 0.934, 0.919, and 0.952 for the training, validation, and entire cohort, respectively. Stratification analysis showed that the established risk model may add prognostic values to conventional clinical factors (age, neoplasm histologic grade, and clinical stage). A nomogram was constructed based on the risk model and clinical parameters, with the AUC of 0.978 and c-index of 0.8079. Database for Annotation, Visualization, and Integrated Discovery (DAVID) function enrichment and Human Protein Atlas (HPA) protein expression validation showed 5 of these 14 genes may be especially important for UCEC (hypermethylated lowly expressed: CCBE1, FOXL2, PHLDB2, and DTNA; hypomethylated highly expressed: CCNE1). Comparison with breast cancer in the methylation level indicated ABCA12, CCNE1, and CLRN3 may be specific methylation-driven genes for UCEC. LncRNA HCG11 may function by coexpressing with DTNA. In conclusion, this 14-DNA methylation signature combined with clinical factors may a potentially effective biomarker in predicting OS for UCEC patients.