Uncovering MicroRNA and Transcription Factor Mediated Regulatory Networks in Glioblastoma

Computational inference of co-regulatory modules from transcription factors, MicroRNAs, and their targets using CanMod2

Gene regulators such as Transcription Factors (TFs) and microRNAs (miRNAs) regulate genes at the transcriptional and post-transcriptional levels, respectively. There is a complex interplay of regulatory patterns of TFs and miRNAs. Some TFs and miRNAs regulate the activity of their target genes individually, some co-regulate the activity of the same set of genes, some TFs regulate miRNA activity, and some miRNAs regulate TFs. As dysregulation in gene regulation can lead to various diseases like cancer, it is a significant problem to find the interplay among TFs, miRNAs, and their target genes. Here, we propose a computational pipeline, CanMod2, which infers modules of TFs, miRNAs, and their co-regulatory targets that are involved in a common biological process. In this work, we have introduced several algorithmic enhancements to the earlier version of CanMod2. We applied CanMod2 to five cancer types and analyzed the inferred modules extensively. Our results show that the inferred modules were enriched in cancer-related biological processes and pathways. The hub regulators that occur in many modules were among cancer-related genes and miRNAs. The inferred regulator-target interactions were significantly enriched in ground truth interactions. CanMod2 source code and documentation are publicly available at https://github.com/bozdaglab/CanMod2 .

Mechanistic Modeling the Role of MicroRNAs and Transcription Factors in Disease Progression

In this chapter, we illustrate the utilization of network analysis and mechanistic modeling, two potent branches of systems biology, to simplify the representation of intricate biological processes such as cell signaling, gene regulation, and metabolic pathways. Specifically, we demonstrate the application of a well-established method to generate a microRNA-transcription factor-gene regulatory feed-forward loop network extracted from the GEO dataset GSE163877. Furthermore, we outline a method for constructing a deterministic model using the LSODA method based on the sub-network. This model furnishes insights into the roles of crucial differentially expressed microRNAs and transcription factors in gene expression associated with Alzheimer's disease progression. Our analysis of the model reveals elevated kinetics of synthesis for EGR1, miR-6891, miR-4786, and LTBP1. The model suggests the linear upregulation of miR-8080, miR-3921, HSPB6, and downregulation MX2 gene. The rest of the miRNA, TFs, and genes shows a momentary variation in expression and if the system is undisturbed, they attain equilibrium. Thus, we elucidate how mechanistic modeling, along with perturbation studies and network analysis of expression data, can yield diverse insights into the trajectory of disease progression.

Investigating gene expression datasets of hippocampus tissue to discover Alzheimer's disease-associated molecular markers

BACKGROUND

Alzheimer's disease (AD) is an advancing neurodegenerative disorder distinguished by the formation of amyloid plaques and neurofibrillary tangles in the human brain. Nevertheless, the lack of peripheral biomarkers that can detect the development of AD remains a significant limitation.

OBJECTIVE

The main aim of this work was to discover the molecular markers associated with AD.

METHODS

We conducted a comprehensive microarray analysis of gene expression data from hippocampus tissue in AD patients and control samples using three microarray datasets (GSE1297, GSE28146, and GSE29378) collected from Gene Expression Omnibus (GEO). The datasets were pre-processed and normalized, revealing 346 significant genes, 103 of which were upregulated and 243 downregulated. The PPI network of significant genes was constructed to detect the top 50 hub genes, which were then further analyzed using Gene Ontology (GO) terms, Kyoto Encyclopedia of Genes and Genomes pathway (KEGG), and GSEA, revealing 47 key genes involved in AD-related pathways. These key genes were then subjected to feed forward loop (FFL) motif analysis for the prediction of transcriptional factors (TFs) and microRNAs (miRNAs) mediated gene regulatory networks.

RESULTS

The interaction of AD-associated TFs HNF4A, SPI1, EGR1, STAT3, and MYC and miRNAs hsa-miR-155-5p and hsa-miR-16-5p in the transcriptional and post-transcriptional events of 3 upregulated and 10 downregulated genes: H2AFZ, MCM3, MYO1C, AXIN1, CCND1, ETS2, MYH9, RELA, RHEB, SOCS3, TBL1X, TBP, TXNIP, and YWHAZ, respectively, has been identified. The miRNA/TF-mediated three types of the FFL motifs, i.e., miRNA-FFL, TF-FFL, and composite-FFL, were constructed, and seven common genes among these FFL were identified: CCND1, MYH9, SOCS3, RHEB, MYO1C, TXNIP, AXIN1, and TXNIP.

CONCLUSIONS

These findings may provide insights into the development of potential molecular markers for therapeutic management of AD.

GraphCVAE: Uncovering cell heterogeneity and therapeutic target discovery through residual and contrastive learning

Advancements in Spatial Transcriptomics (ST) technologies in recent years have transformed the analysis of tissue structure and function within spatial contexts. However, accurately identifying spatial domains remains challenging due to data sparsity and noise. Traditional clustering methods often fail to capture spatial dependencies, while spatial clustering methods struggle with batch effects and data integration. We introduce GraphCVAE, a model designed to enhance spatial domain identification by integrating spatial and morphological information, correcting batch effects, and managing heterogeneous data. GraphCVAE employs a multi-layer Graph Convolutional Network (GCN) and a variational autoencoder to improve the representation and integration of spatial information. Through contrastive learning, the model captures subtle differences between cell types and states. Extensive testing on various ST datasets demonstrates GraphCVAE's robustness and biological contributions. In the dorsolateral prefrontal cortex (DLPFC) dataset, it accurately delineates cortical layer boundaries. In glioblastoma, GraphCVAE reveals critical therapeutic targets such as TF and NFIB. In colorectal cancer, it explores the role of the extracellular matrix in colorectal cancer. The model's performance metrics consistently surpass existing methods, validating its effectiveness. GraphCVAE's advanced visualization capabilities further highlight its precision in resolving spatial structures, making it a powerful tool for spatial transcriptomics analysis and offering new insights into disease studies.

Use of microRNAs as Diagnostic, Prognostic, and Therapeutic Tools for Glioblastoma

Glioblastoma (GB) is the most aggressive and common type of cancer within the central nervous system (CNS). Despite the vast knowledge of its physiopathology and histology, its etiology at the molecular level has not been completely understood. Thus, attaining a cure has not been possible yet and it remains one of the deadliest types of cancer. Usually, GB is diagnosed when some symptoms have already been presented by the patient. This diagnosis is commonly based on a physical exam and imaging studies, such as computed tomography (CT) and magnetic resonance imaging (MRI), together with or followed by a surgical biopsy. As these diagnostic procedures are very invasive and often result only in the confirmation of GB presence, it is necessary to develop less invasive diagnostic and prognostic tools that lead to earlier treatment to increase GB patients' quality of life. Therefore, blood-based biomarkers (BBBs) represent excellent candidates in this context. microRNAs (miRNAs) are small, non-coding RNAs that have been demonstrated to be very stable in almost all body fluids, including saliva, serum, plasma, urine, cerebrospinal fluid (CFS), semen, and breast milk. In addition, serum-circulating and exosome-contained miRNAs have been successfully used to better classify subtypes of cancer at the molecular level and make better choices regarding the best treatment for specific cases. Moreover, as miRNAs regulate multiple target genes and can also act as tumor suppressors and oncogenes, they are involved in the appearance, progression, and even chemoresistance of most tumors. Thus, in this review, we discuss how dysregulated miRNAs in GB can be used as early diagnosis and prognosis biomarkers as well as molecular markers to subclassify GB cases and provide more personalized treatments, which may have a better response against GB. In addition, we discuss the therapeutic potential of miRNAs, the current challenges to their clinical application, and future directions in the field.

Investigating the link between miR-34a-5p and TLR6 signaling in sepsis-induced ARDS

Acute lung injury (ALI)/acute respiratory distress syndrome (ARDS) are lung complications diagnosed by impaired gaseous exchanges leading to mortality. From the diverse etiologies, sepsis is a prominent contributor to ALI/ARDS. In the present study, we retrieved sepsis-induced ARDS mRNA expression profile and identified 883 differentially expressed genes (DEGs). Next, we established an ARDS-specific weighted gene co-expression network (WGCN) and picked the blue module as our hub module based on highly correlated network properties. Later we subjected all hub module DEGs to form an ARDS-specific 3-node feed-forward loop (FFL) whose highest-order subnetwork motif revealed one TF (STAT6), one miRNA (miR-34a-5p), and one mRNA (TLR6). Thereafter, we screened a natural product library and identified three lead molecules that showed promising binding affinity against TLR6. We then performed molecular dynamics simulations to evaluate the stability and binding free energy of the TLR6-lead molecule complexes. Our results suggest these lead molecules may be potential therapeutic candidates for treating sepsis-induced ALI/ARDS. In-silico studies on clinical datasets for sepsis-induced ARDS indicate a possible positive interaction between miR-34a and TLR6 and an antagonizing effect on STAT6 to promote inflammation. Also, the translational study on septic mice lungs by IHC staining reveals a hike in the expression of TLR6. We report here that miR-34a actively augments the effect of sepsis on lung epithelial cell apoptosis. This study suggests that miR-34a promotes TLR6 to heighten inflammation in sepsis-induced ALI/ARDS.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03700-1.

H19X-encoded microRNAs induced by IL-4 in adipocyte precursors regulate proliferation to facilitate differentiation

Adipose tissue, an organ critical for systemic energy homeostasis, is influenced by type 2 immunity in its development and function. The type 2 cytokine interleukin (IL)-4 induces the proliferation of bipotential adipocyte precursors (APs) in white fat tissue and primes these cells for differentiation into beige adipocytes, which are specialized for thermogenesis. However, the underlying mechanisms have not yet been comprehensively examined. Here, we identified six microRNA (miRNA) genes upregulated upon IL-4 stimulation in APs, miR-322, miR-503, miR-351, miR-542, miR-450a, and miR-450b; these are encoded in the H19X locus of the genome. Their expression is positively regulated by the transcription factor Klf4, whose expression also increases upon IL-4 stimulation. These miRNAs shared a large set of target genes, of which 381 genes were downregulated in mRNA expression upon IL-4 stimulation and enriched in Wnt signaling pathways. Two genes with downregulated expression, Ccnd1 and Fzd6, were repressed by H19X-encoded miRNAs. Additionally, the Wnt signaling activator LiCl downregulated the expression of this group of miRNAs in APs, indicating that Wnt signaling-related genes and these miRNAs form a double-negative feedback regulatory loop. This miRNA/Wnt feedback regulation modulated the elevated proliferation of APs induced by IL-4 stimulation and contributed to priming them for beige adipocyte differentiation. Moreover, the aberrant expression of these miRNAs attenuates the differentiation of APs into beige adipocytes. Collectively, our results suggest that H19X-encoded miRNAs facilitate the transition of APs from proliferation to differentiation in the IL-4-mediated regulation.

Deciphering molecular mechanisms of SARS-CoV-2 pathogenesis and drug repurposing through GRN motifs: a comprehensive systems biology study

The world has recently been plagued by a new coronavirus infection called SARS-CoV-2. This virus may lead to severe acute respiratory syndrome followed by multiple organ failure. SARS-CoV-2 has approximately 80-90% genetic similarity to SARS-CoV. Given the limited omics data available for host response to the viruses (more limited data for SARS-CoV-2), we attempted to unveil the crucial molecular mechanisms underlying the SARS-CoV-2 pathogenesis by comparing its regulatory network motifs with SARS-CoV. We also attempted to identify the non-shared crucial molecules and their functions to predict the specific mechanisms for each infection and the processes responsible for their different manifestations. Deciphering the crucial shared and non-shared mechanisms at the molecular level and signaling pathways underlying both diseases may help shed light on their pathogenesis and pave the way for other new drug repurposing against COVID-19. We constructed the GRNs for host response to SARS-CoV and SARS-CoV-2 pathogens (in vitro) and identified the significant 3-node regulatory motifs by analyzing them topologically and functionally. We attempted to identify the shared and non-shared regulatory elements and signaling pathways between their host responses. Interestingly, our findings indicated that NFKB1, JUN, STAT1, FOS, KLF4, and EGR1 were the critical shared TFs between motif-related subnetworks in both SARS and COVID-1, which are considered genes with specific functions in the immune response. Enrichment analysis revealed that the NOD-like receptor signaling, TNF signaling, and influenza A pathway were among the first significant pathways shared between SARS and COVID-19 up-regulated DEGs networks, and the term "metabolic pathways" (hsa01100) among the down-regulated DEGs networks. WEE1, PMAIP1, and TSC22D2 were identified as the top three hubs specific to SARS. However, MYPN, SPRY4, and APOL6 were the tops specific to COVID-19 in vitro. The term "Complement and coagulation cascades" pathway was identified as the first top non-shared pathway for COVID-19 and the MAPK signaling pathway for SARS. We used the identified crucial DEGs to construct a drug-gene interaction network to propose some drug candidates. Zinc chloride, Fostamatinib, Copper, Tirofiban, Tretinoin, and Levocarnitine were the six drugs with higher scores in our drug-gene network analysis.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03518-x.

Identification of potential microRNA diagnostic panels and uncovering regulatory mechanisms in breast cancer pathogenesis

Early diagnosis of breast cancer (BC), as the most common cancer among women, increases the survival rate and effectiveness of treatment. MicroRNAs (miRNAs) control various cell behaviors, and their dysregulation is widely involved in pathophysiological processes such as BC development and progress. In this study, we aimed to identify potential miRNA biomarkers for early diagnosis of BC. We also proposed a consensus-based strategy to analyze the miRNA expression data to gain a deeper insight into the regulatory roles of miRNAs in BC initiation. Two microarray datasets (GSE106817 and GSE113486) were analyzed to explore the differentially expressed miRNAs (DEMs) in serum of BC patients and healthy controls. Utilizing multiple bioinformatics tools, six serum-based miRNA biomarkers (miR-92a-3p, miR-23b-3p, miR-191-5p, miR-141-3p, miR-590-5p and miR-190a-5p) were identified for BC diagnosis. We applied our consensus and integration approach to construct a comprehensive BC-specific miRNA-TF co-regulatory network. Using different combination of these miRNA biomarkers, two novel diagnostic models, consisting of miR-92a-3p, miR-23b-3p, miR-191-5p (model 1) and miR-92a-3p, miR-23b-3p, miR-141-3p, and miR-590-5p (model 2), were obtained from bioinformatics analysis. Validation analysis was carried out for the considered models on two microarray datasets (GSE73002 and GSE41922). The model based on similar network topology features, comprising miR-92a-3p, miR-23b-3p and miR-191-5p was the most promising model in the diagnosis of BC patients from healthy controls with 0.89 sensitivity, 0.96 specificity and area under the curve (AUC) of 0.98. These findings elucidate the regulatory mechanisms underlying BC and represent novel biomarkers for early BC diagnosis.

Pan-cancer analyses of classical protein tyrosine phosphatases and phosphatase-targeted therapy in cancer

Protein tyrosine phosphatases function in dephosphorylating target proteins to regulate signaling pathways that control a broad spectrum of fundamental physiological and pathological processes. Detailed knowledge concerning the roles of classical PTPs in human cancer merits in-depth investigation. We comprehensively analyzed the regulatory mechanisms and clinical relevance of classical PTPs in more than 9000 tumor patients across 33 types of cancer. The independent datasets and functional experiments were employed to validate our findings. We exhibited the extensive dysregulation of classical PTPs and constructed the gene regulatory network in human cancer. Moreover, we characterized the correlation of classical PTPs with both drug-resistant and drug-sensitive responses to anti-cancer drugs. To evaluate the PTP activity in cancer prognosis, we generated a PTPscore based on the expression and hazard ratio of classical PTPs. Our study highlights the notable role of classical PTPs in cancer biology and provides novel intelligence to improve potential therapeutic strategies based on pTyr regulation.

Integrated transcriptomic and regulatory network analyses uncovers the role of let-7b-5p, SPIB, and HLA-DPB1 in sepsis

Sepsis has affected millions of populations of all age groups, locations, and sexes worldwide. Immune systems, either innate or adaptive are dysregulated due to the infection. Various biomarkers are present to date, still sepsis is a primary cause of mortality. Globally, post-operative body infections can cause sepsis and septic shock in ICU. Abnormal antigen presentation to T-cells leads to a dysregulated immune system. miRNAs are sparkly evolved as biomarkers due to their high sensitivity and efficiency. In this work, we analyzed high-throughput mRNA data collected from Gene Expression Omnibus (GEO) and linked it to significant miRNAs and TFs using a network-based approach. Protein–protein interaction (PPI) network was constructed using sepsis-specific differentially expressed genes (DEGs) followed by enrichment analyses and hub module detection. Sepsis-linked decrease transcription of the classical HLA gene such as HLA-DPB1 and its interplay with miR-let-7b-5p and transcription factor SPIB was observed. This study helped to provide innovative targets for sepsis.

Differential Gene Expression and Weighted Correlation Network Dynamics in High-Throughput Datasets of Prostate Cancer

Precision oncology is an absolute need today due to the emergence of treatment resistance and heterogeneity among cancerous profiles. Target-propelled cancer therapy is one of the treasures of precision oncology which has come together with substantial medical accomplishment. Prostate cancer is one of the most common cancers in males, with tremendous biological heterogeneity in molecular and clinical behavior. The spectrum of molecular abnormalities and varying clinical patterns in prostate cancer suggest substantial heterogeneity among different profiles. To identify novel therapeutic targets and precise biomarkers implicated with prostate cancer, we performed a state-of-the-art bioinformatics study, beginning with analyzing high-throughput genomic datasets from The Cancer Genome Atlas (TCGA). Weighted gene co-expression network analysis (WGCNA) suggests a set of five dysregulated hub genes (MAF, STAT6, SOX2, FOXO1, and WNT3A) that played crucial roles in biological pathways associated with prostate cancer progression. We found overexpressed STAT6 and SOX2 and proposed them as candidate biomarkers and potential targets in prostate cancer. Furthermore, the alteration frequencies in STAT6 and SOX2 and their impact on the patients' survival were explored through the cBioPortal platform. The Kaplan-Meier survival analysis suggested that the alterations in the candidate genes were linked to the decreased overall survival of the patients. Altogether, the results signify that STAT6 and SOX2 and their genomic alterations can be explored in therapeutic interventions of prostate cancer for precision oncology, utilizing early diagnosis and target-propelled therapy.

Algorithmic reconstruction of glioblastoma network complexity

Glioblastoma is a complex disease that is difficult to treat. Network and data science offer alternative approaches to classical bioinformatics pipelines to study gene expression patterns from single-cell datasets, helping to distinguish genes associated with the control of differentiation and aggression. To identify the key molecular regulators of the networks driving glioblastoma/GSC and predict their cell fate dynamics, we applied a host of data theoretic techniques to gene expression patterns from pediatric and adult glioblastoma, and adult glioma-derived stem cells (GSCs). We identified eight transcription factors (OLIG1/2, TAZ, GATA2, FOXG1, SOX6, SATB2, and YY1) and four signaling genes (ATL3, MTSS1, EMP1, and TPT1) as coordinators of cell state transitions and, thus, clinically targetable putative factors differentiating pediatric and adult glioblastomas from adult GSCs. Our study provides strong evidence of complex systems approaches for inferring complex dynamics from reverse-engineering gene networks, bolstering the search for new clinically relevant targets in glioblastoma.

•

Complex cell fate attractors capture glioblastoma differentiation dynamics

•

Graph theoretic approaches decode master regulators of GBM glioblastoma cell fate decisions

•

Network dynamics of pediatric glioblastoma resemble adult GSCs

•

Transcriptional networks may help reprogram glioblastoma behavioral patterns

Crucial Roles of microRNA-16-5p and microRNA-27b-3p in Ameloblast Differentiation Through Regulation of Genes Associated With Amelogenesis Imperfecta

Amelogenesis imperfecta is a congenital disorder within a heterogeneous group of conditions characterized by enamel hypoplasia. Patients suffer from early tooth loss, social embarrassment, eating difficulties, and pain due to an abnormally thin, soft, fragile, and discolored enamel with poor aesthetics and functionality. The etiology of amelogenesis imperfecta is complicated by genetic interactions. To identify mouse amelogenesis imperfecta-related genes (mAIGenes) and their respective phenotypes, we conducted a systematic literature review and database search and found and curated 70 mAIGenes across all of the databases. Our pathway enrichment analysis indicated that these genes were enriched in tooth development-associated pathways, forming four distinct groups. To explore how these genes are regulated and affect the phenotype, we predicted microRNA (miRNA)-gene interaction pairs using our bioinformatics pipeline. Our miRNA regulatory network analysis pinpointed that miR-16-5p, miR-27b-3p, and miR-23a/b-3p were hub miRNAs. The function of these hub miRNAs was evaluated through ameloblast differentiation assays with/without the candidate miRNA mimics using cultured mouse ameloblast cells. Our results revealed that overexpression of miR-16-5p and miR-27b-3p, but not miR-23a/b-3p, significantly inhibited ameloblast differentiation through regulation of mAIGenes. Thus, our study shows that miR-16-5p and miR-27b-3p are candidate pathogenic miRNAs for amelogenesis imperfecta.

MicroRNAs Regulate TASK-1 and Are Linked to Myocardial Dilatation in Atrial Fibrillation

Background

Atrial fibrillation (AF) is the most common sustained cardiac arrhythmia. However, underlying molecular mechanisms are insufficiently understood. Previous studies suggested that microRNA (miRNA) dependent gene regulation plays an important role in the initiation and maintenance of AF. The 2‐pore‐domain potassium channel TASK‐1 (tandem of P domains in a weak inward rectifying K⁺ channel–related acid sensitive K⁺ channel 1) is an atrial‐specific ion channel that is upregulated in AF. Inhibition of TASK‐1 current prolongs the atrial action potential duration to similar levels as in patients with sinus rhythm. Here, we hypothesize that miRNAs might be responsible for the regulation of KCNK3 that encodes for TASK‐1.

Methods and Results

We selected miRNAs potentially regulating KCNK3 and studied their expression in atrial tissue samples obtained from patients with sinus rhythm, paroxysmal AF, or permanent/chronic AF. MiRNAs differentially expressed in AF were further investigated for their ability to regulate KCNK3 mRNA and TASK‐1 protein expression in human induced pluripotent stem cells, transfected with miRNA mimics or inhibitors. Thereby, we observed that miR‐34a increases TASK‐1 expression and current and further decreases the resting membrane potential of Xenopus laevis oocytes, heterologously expressing hTASK‐1. Finally, we investigated associations between miRNA expression in atrial tissues and clinical parameters of our patient cohort. A cluster containing AF stage, left ventricular end‐diastolic diameter, left ventricular end‐systolic diameter, left atrial diameter, atrial COL1A2 (collagen alpha‐2(I) chain), and TASK‐1 protein level was associated with increased expression of miR‐25, miR‐21, miR‐34a, miR‐23a, miR‐124, miR‐1, and miR‐29b as well as decreased expression of miR‐9 and miR‐485.

Conclusions

These results suggest an important pathophysiological involvement of miRNAs in the regulation of atrial expression of the TASK‐1 potassium channel in patients with atrial cardiomyopathy.

Deciphering Key Genes and miRNAs Associated With Hepatocellular Carcinoma via Network-Based Approach

Hepatocellular carcinoma (HCC)is a common type of liver cancer and has a high mortality world-widely. The diagnosis, prognoses, and therapeutics are very poor due to the unclear molecular mechanism of progression of the disease. To unveil the molecular mechanism of progression of HCC, we extract a large sample of mRNA expression levels from the GEO database where a total of 167 samples were used for study, and out of them, 115 samples were from HCC tumor tissue. This study aims to investigate the module of differentially expressed genes (DEGs)which are co-expressed only in HCC sample data but not in normal tissue samples. Thereafter, we identified the highly significant module of significant co-expressed genes and formed a PPI network for these genes. There were only six genes (namely, MSH3, DMC1, ALPP, IL10, ZNF223, and HSD17B7)obtained after analysis of the PPI network. Out of six only MSH3, DMC1, HSD17B7, and IL10 were found enriched in GO Term & Pathway enrichment analysis and these candidate genes were mainly involved in cellular process, metabolic and catalytic activity, which promote the development & progression of HCC. Lastly, the composite 3-node FFL reveals the driver miRNAs and TFs associated with our key genes.

Identification of microRNAs and gene regulatory networks in cleft lip common in humans and mice

The etiology of cleft lip with/without cleft palate (CL/P), one of the most frequent craniofacial birth defects worldwide, is complicated by contributions of both genetic and environmental factors. Understanding the etiology of these conditions is essential for developing preventive strategies. This study thus aims to identify regulatory networks of microRNAs (miRNAs), transcriptional factors (TFs) and non-TF genes associated with cleft lip (CL) that are conserved in humans and mice. Notably, we found that miR-27b, miR-133b, miR-205, miR-376b and miR-376c were involved in the regulation of CL-associated gene expression in both humans and mice. Among the candidate miRNAs, the overexpression of miR-27b, miR-133b and miR-205, but not miR-376b and miR-376c, significantly inhibited cell proliferation through suppression of CL-associated genes (miR-27b suppressed PAX9 and RARA; miR-133b suppressed FGFR1, PAX7, and SUMO1; and miR-205 suppressed PAX9 and RARA) in cultured human and mouse lip mesenchymal cells. Taken together, our results suggest that elevated expression of miR-27b, miR-133b and miR-205 may play a crucial role in CL through the suppression of genes associated with CL.

Comprehensive Integrative Analysis Reveals the Association of KLF4 with Macrophage Infiltration and Polarization in Lung Cancer Microenvironment

Macrophage polarization and infiltration to the tumor microenvironment (TME) is a critical determining factor for tumor progression. Macrophages are polarized into two states—M1 (pro-inflammatory, anti-tumorigenic and stimulated by LPS or IFN-γ) and M2 (anti-inflammatory pro-tumorigenic and stimulated by IL-4) phenotypes. Specifically, M2 macrophages enhance tumor cell growth and survival. Recent evidences suggest the pivotal role of microRNAs in macrophage polarization during the development of Non-small cell lung cancer (NSCLC), thus proposing a new therapeutic option to target lung cancer. In silico analysis determined cogent upregulation of KLF4, downregulation of IL-1β and miR-34a-5p in NSCLC tissues, consequently worsening the overall survival of NSCLC patients. We observed a significant association of KLF4 with macrophage infiltration and polarization in NSCLC. We found that KLF4 is critically implicated in M2 polarization of macrophages, which, in turn, promotes tumorigenesis. KLF4 expression correlated with miR-34a-5p and IL-1β in a feed-forward loop (FFL), both of which are implicated in immune regulation. Mechanistic overexpression of miR-34a-5p in macrophages (IL-4 stimulated) inhibits KLF4, along with downregulation of ARG1, REL-1MB (M2 macrophage specific markers), and upregulation of IL-1β, IL-6, (M1 macrophage specific markers), demonstrating macrophage polarization switch from M2 to M1 phenotype. Moreover, co-culture of these macrophages with NSCLC cells reduces their proliferation, wound healing, clonogenic capacity and enhanced NO-mediated apoptosis. Further, transfection of miR-34a-5p in NSCLC cells, also degrades KLF4, but enhances the expression of KLF4 regulated genes—IL-1β, IL-6 (pro-inflammatory mediators), which is further enhanced upon co-culture with IL-4 stimulated macrophages. Additionally, we observed a significant increase in i-NOS/NO content upon co-culture, suggesting polarization reversion of macrophages from M2 to M1, and eventually leading to anti-tumor effects. Our findings thus show a significant role of KLF4 in tumorigenesis and TAM polarization of NSCLC. However, miR-34a-5p mediated targeting of these molecular networks will provide a better therapeutic intervention for NSCLC.

miRGTF-net: Integrative miRNA-gene-TF network analysis reveals key drivers of breast cancer recurrence

Analysis of regulatory networks is a powerful framework for identification and quantification of intracellular interactions. We introduce miRGTF-net, a novel tool for construction of miRNA-gene-TF networks. We consider multiple transcriptional and post-transcriptional interaction types, including regulation of gene and miRNA expression by transcription factors, gene silencing by miRNAs, and co-expression of host genes with their intronic miRNAs. The underlying algorithm uses information on experimentally validated interactions as well as integrative miRNA/mRNA expression profiles in a given set of samples. The latter ensures simultaneous tissue-specificity and biological validity of interactions. We applied miRGTF-net to paired miRNA/mRNA-sequencing data of breast cancer samples from The Cancer Genome Atlas (TCGA). Together with topological analysis of the constructed network we showed that considered players can form reliable prognostic gene signatures for ER-positive breast cancer. A number of signatures demonstrated remarkably high accuracy on transcriptomic data obtained by both microarrays and RNA sequencing from several independent patient cohorts. Furthermore, an essential part of prognostic genes were identified as direct targets of transcription factor E2F1. The putative interplay between estrogen receptor alpha and E2F1 was suggested as a potential recurrence factor in patients treated with tamoxifen. Source codes of miRGTF-net are available at GitHub (https://github.com/s-a-nersisyan/miRGTF-net).

Construction of a TF-miRNA-gene feed-forward loop network predicts biomarkers and potential drugs for myasthenia gravis

Myasthenia gravis (MG) is an autoimmune disease and the most common type of neuromuscular disease. Genes and miRNAs associated with MG have been widely studied; however, the molecular mechanisms of transcription factors (TFs) and the relationship among them remain unclear. A TF–miRNA–gene network (TMGN) of MG was constructed by extracting six regulatory pairs (TF–miRNA, miRNA–gene, TF–gene, miRNA–TF, gene–gene and miRNA–miRNA). Then, 3/4/5-node regulatory motifs were detected in the TMGN. Then, the motifs with the highest Z-score, occurring as 3/4/5-node composite feed-forward loops (FFLs), were selected as statistically significant motifs. By merging these motifs together, we constructed a 3/4/5-node composite FFL motif-specific subnetwork (CFMSN). Then, pathway and GO enrichment analyses were performed to further elucidate the mechanism of MG. In addition, the genes, TFs and miRNAs in the CFMSN were also utilized to identify potential drugs. Five related genes, 3 TFs and 13 miRNAs, were extracted from the CFMSN. As the most important TF in the CFMSN, MYC was inferred to play a critical role in MG. Pathway enrichment analysis showed that the genes and miRNAs in the CFMSN were mainly enriched in pathways related to cancer and infections. Furthermore, 21 drugs were identified through the CFMSN, of which estradiol, estramustine, raloxifene and tamoxifen have the potential to be novel drugs to treat MG. The present study provides MG-related TFs by constructing the CFMSN for further experimental studies and provides a novel perspective for new biomarkers and potential drugs for MG.

Synergism of Proneurogenic miRNAs Provides a More Effective Strategy to Target Glioma Stem Cells

Simple Summary

miRNAs function as critical regulators of gene expression and have been defined as contributors of cancer phenotypes by acting as oncogenes or tumor suppressors. Based on these findings, miRNA-based therapies have been explored in the treatment of many different malignancies. The use of single miRNAs has faced some challenges and showed limited success. miRNAs cooperate to regulate distinct biological processes and pathways and, therefore, combination of related miRNAs could amplify the repression of oncogenic factors and the effect on cancer relevant pathways. We established that the combination of tumor suppressor miRNAs miR-124, miR-128, and miR-137 is much more effective than single miRNAs in disrupting proliferation and survival of glioma stem cells and neuroblastoma lines and promoting differentiation and response to radiation. Subsequent genomic analyses showed that other combinations of tumor suppressor miRNAs could be equally effective, and its use could provide new routes to target in special cancer-initiating cell populations.

Abstract

Tumor suppressor microRNAs (miRNAs) have been explored as agents to target cancer stem cells. Most strategies use a single miRNA mimic and present many disadvantages, such as the amount of reagent required and the diluted effect on target genes. miRNAs work in a cooperative fashion to regulate distinct biological processes and pathways. Therefore, we propose that miRNA combinations could provide more efficient ways to target cancer stem cells. We have previously shown that miR-124, miR-128, and miR-137 function synergistically to regulate neurogenesis. We used a combination of these three miRNAs to treat glioma stem cells and showed that this treatment was much more effective than single miRNAs in disrupting cell proliferation and survival and promoting differentiation and response to radiation. Transcriptomic analyses indicated that transcription regulation, angiogenesis, metabolism, and neuronal differentiation are among the main biological processes affected by transfection of this miRNA combination. In conclusion, we demonstrated the value of using combinations of neurogenic miRNAs to disrupt cancer phenotypes and glioma stem cell growth. The synergistic effect of these three miRNA amplified the repression of oncogenic factors and the effect on cancer relevant pathways. Future therapeutic approaches would benefit from utilizing miRNA combinations, especially when targeting cancer-initiating cell populations.

A developmental stage-specific network approach for studying dynamic co-regulation of transcription factors and microRNAs during craniofacial development

Craniofacial development is regulated through dynamic and complex mechanisms that involve various signaling cascades and gene regulations. Disruption of such regulations can result in craniofacial birth defects. Here, we propose the first developmental stage-specific network approach by integrating two crucial regulators, transcription factors (TFs) and microRNAs (miRNAs), to study their co-regulation during craniofacial development. Specifically, we used TFs, miRNAs and non-TF genes to form feed-forward loops (FFLs) using genomic data covering mouse embryonic days E10.5 to E14.5. We identified key novel regulators (TFs Foxm1, Hif1a, Zbtb16, Myog, Myod1 and Tcf7, and miRNAs miR-340-5p and miR-129-5p) and target genes (Col1a1, Sgms2 and Slc8a3) expression of which changed in a developmental stage-dependent manner. We found that the Wnt-FoxO-Hippo pathway (from E10.5 to E11.5), tissue remodeling (from E12.5 to E13.5) and miR-129-5p-mediated Col1a1 regulation (from E10.5 to E14.5) might play crucial roles in craniofacial development. Enrichment analyses further suggested their functions. Our experiments validated the regulatory roles of miR-340-5p and Foxm1 in the Wnt-FoxO-Hippo subnetwork, as well as the role of miR-129-5p in the miR-129-5p-Col1a1 subnetwork. Thus, our study helps understand the comprehensive regulatory mechanisms for craniofacial development.

Integrated analysis of microRNA and transcription factors in the bone marrow of patients with acute monocytic leukemia

Acutemonocytic leukemia (AMoL) is a distinct subtype of acute myeloid leukemia (AML) with poor prognosis. However, the molecular mechanisms and key regulators involved in the global regulation of gene expression levels in AMoL are poorly understood. In order to elucidate the role of microRNAs (miRNAs/miRs) and transcription factors (TFs) in AMoL pathogenesis at the network level, miRNA and TF expression level profiles were systematically analyzed by miRNA sequencing and TF array, respectively; this identified 285 differentially expressed miRNAs and 139 differentially expressed TFs in AMoL samples compared with controls. By combining expression level profile data and bioinformatics tools available for predicting TF and miRNA targets, a comprehensive AMoL-specific miRNA-TF-mediated regulatory network was constructed. A total of 26 miRNAs and 23 TFs were identified as hub nodes in the network. Among these hubs, miR-29b-3p, MYC, TP53 and NFKB1 were determined to be potential AMoL regulators, and were subsequently extracted to construct sub-networks. A hypothetical pathway model was also proposed for miR-29b-3p to reveal the potential co-regulatory mechanisms of miR-29b-3p, MYC, TP53 and NFKB1 in AMoL. The present study provided an effective approach to discover critical regulators via a comprehensive regulatory network in AMoL, in addition to enhancing understanding of the pathogenesis of this disease at the molecular level.

Notch1 in Cancer Therapy: Possible Clinical Implications and Challenges

The Notch family consists of four highly conserved transmembrane receptors. The release of the active intracellular domain requires the enzymatic activity of γ-secretase. Notch is involved in embryonic development and in many physiologic processes of normal cells, in which it regulates growth, apoptosis, and differentiation. Notch1, a member of the Notch family, is implicated in many types of cancer, including breast cancer (especially triple-negative breast cancer), leukemias, brain tumors, and many others. Notch1 is tightly connected to many signaling pathways that are therapeutically involved in tumorigenesis. Together, they impact apoptosis, proliferation, chemosensitivity, immune response, and the population of cancer stem cells. Notch1 inhibition can be achieved through various and diverse methods, the most common of which are the γ-secretase inhibitors, which produce a pan-Notch inhibition, or the use of Notch1 short interference RNA or Notch1 monoclonal antibodies, which produce a more specific blockade. Downregulation of Notch1 can be used alone or in combination with chemotherapy, which can achieve a synergistic effect and a decrease in chemoresistance. Targeting Notch1 in cancers that harbor high expression levels of Notch1 offers an addition to therapeutic strategies recruited for managing cancer. Considering available evidence, Notch1 offers a legitimate target that might be incorporated in future strategies for combating cancer. In this review, the possible clinical applications of Notch1 inhibition and the obstacles that hinder its clinical application are discussed. SIGNIFICANCE STATEMENT: Notch1 plays an important role in different types of cancer. Numerous approaches of Notch1 inhibition possess potential benefits in the management of various clinical aspects of cancer. The application of different Notch1 inhibition modalities faces many challenges.

Integrative Analysis of Regulatory Module Reveals Associations of Microgravity with Dysfunctions of Multi-body Systems and Tumorigenesis

Previous studies have demonstrated that microgravity could lead to health risks. The investigation of the molecular mechanisms from the aspect of systems biology has not been performed yet. Here, we integratively analyzed transcriptional and post-transcriptional regulations based on gene and miRNA expression profiles in human peripheral blood lymphocytes cultured in modeled microgravity. Two hundred and thirty dysregulated TF-miRNA (transcription factor and microRNA) feed-forward loops (FFLs) were identified in microgravity. The immune, cardiovascular, endocrine, nervous and skeletal system subnetworks were constructed according to the functions of dysregulated FFLs. Taking the skeletal system as an example, most of genes and miRNAs in the subnetwork were involved in bone loss. In addition, several drugs have been predicted to have potential to reduce bone loss, such as traditional Chinese medicines Emodin and Ginsenoside Rh2. Furthermore, we investigated the relationships between microgravity and 20 cancer types, and found that most of cancers might be promoted by microgravity. For example, rectum adenocarcinoma (READ) might be induced by microgravity through reducing antigen presentation and suppressing IgA-antibody-secreting cells' migration. Collectively, TF-miRNA FFL might provide a novel mechanism to elucidate the changes induced by microgravity, serve as drug targets to relieve microgravity effects, and give new insights to explore the relationships between microgravity and cancers.

Integrative mRNA/miRNA expression analysis in healing human gingiva

BACKGROUND

MicroRNAs (miRNAs) are implicated in the epigenetic regulation of complex biological processes. Their possible role in human oral wound healing, a process that differs from cutaneous wound healing by being faster and typically scar-free, has been unexplored. This report presents the miRNA expression profile of experimental human oral wounds and an integrative analysis of mRNA/miRNA expression.

METHODS

Nine healthy volunteers provided standardized normal and 5-day healing palatal biopsies, used for next generation miRNA and mRNA sequencing analysis, correlation and network analysis, real-time PCR (qPCR) and immunohistochemistry.

RESULTS

On average, 169 significantly regulated precursor miRNAs were detected, including 21 novel miRNAs, selectively confirmed by PCR. Hsa-miR-223-3p and hsa-miR-124-3p were, respectively, the most up- and downregulated miRNAs in healing gingiva. Hsa-miR-124-3p had the most predicted mRNA target interactions, with angiogenesis-related genes the most enriched. Correlation analysis showed the highest correlation between hsa-miR-181a-3p and SERPINB1; hsa-miR-223-5p and SLC2A3; hsa-miR-1301 and MS4A7. In addition, SERPINB1 mRNA had the most associations with differentially regulated miRNAs. IL33 was the only cytokine significantly correlated with miRNAs (ρ >  0.95). qPCR and immunohistochemistry verified the significant upregulation of SERPINB1 and IL33 in healing gingiva.

CONCLUSIONS

This study is the first to report on the miRNome of healing human gingiva and to provide an integrative analysis of miRNA/mRNA expression during human oral wound healing; the results offer novel insights into the participating molecular mechanisms and raise the possibility of SERPINB1 and IL-33 as potential wound healing therapeutic targets.

Temozolomide-Induced RNA Interactome Uncovers Novel LncRNA Regulatory Loops in Glioblastoma

Simple Summary

Glioblastoma (GBM) is the most aggressive brain tumor and most resistant to therapy. The identification of novel predictive biomarkers or targets to counteract chemoresistance, requires a better understanding of the GBM primary response to therapy. The aim of our study was to assess the molecular response of GBM to the standard of care chemotherapy by temozolomide (TMZ). We established a comprehensive map of gene expression changes after treatment and discovered that GBM cells elicit a coordinated gene expression program after chemotherapy that differs between sensitive and resistant cells. We found that a novel class of genes expressed as long non-coding RNAs (lncRNAs) is involved in gene regulatory circuits in GBM and could represent novel markers of GBM patient prognosis. By shedding light on the involvement of the non-coding genome in GBM, our results may provide new mechanistic insight on lncRNAs and their importance in chemoresistance.

Abstract

Resistance to chemotherapy by temozolomide (TMZ) is a major cause of glioblastoma (GBM) recurrence. So far, attempts to characterize factors that contribute to TMZ sensitivity have largely focused on protein-coding genes, and failed to provide effective therapeutic targets. Long noncoding RNAs (lncRNAs) are essential regulators of epigenetic-driven cell diversification, yet, their contribution to the transcriptional response to drugs is less understood. Here, we performed RNA-seq and small RNA-seq to provide a comprehensive map of transcriptome regulation upon TMZ in patient-derived GBM stem-like cells displaying different drug sensitivity. In a search for regulatory mechanisms, we integrated thousands of molecular associations stored in public databases to generate a background “RNA interactome”. Our systems-level analysis uncovered a coordinated program of TMZ response reflected by regulatory circuits that involve transcription factors, mRNAs, miRNAs, and lncRNAs. We discovered 22 lncRNAs involved in regulatory loops and/or with functional relevance in drug response and prognostic value in gliomas. Thus, the investigation of TMZ-induced gene networks highlights novel RNA-based predictors of chemosensitivity in GBM. The computational modeling used to identify regulatory circuits underlying drug response and prioritizing gene candidates for functional validation is applicable to other datasets.

H19, a Long Non-coding RNA, Mediates Transcription Factors and Target Genes through Interference of MicroRNAs in Pan-Cancer

Long non-coding RNAs (lncRNAs) have recently been found to be important in gene regulation. lncRNA H19 has been reported to play an oncogenic role in many human cancers. Its specific regulatory role is still elusive. In this study, we developed a novel analytic approach by integrating the synergistic regulation among lncRNAs (e.g., H19), transcription factors (TFs), target genes, and microRNAs (miRNAs) and then applied it to the pan-cancer expression datasets from The Cancer Genome Atlas (TCGA). Using linear regression models, we identified 88 H19-TF-gene co-regulatory triplets, in which 93% of the TF-gene pairs were related to cancer, indicating that our approach was effective to identify disease-related lncRNA-TF-gene co-regulation mechanisms. lncRNAs can function as miRNA sponges. Our further experiments found that H19 might regulate SP1-TGFBR2 through let-7b and miR-200b, ETS1-TGFBR2 through miR-29a and miR-200b, and STAT3-KLF11 through miR-17 in breast cancer cell lines. Our work suggests that miRNA-mediated lncRNA-TF-gene co-regulation is complicated yet important in cancer.

Construction of a potential microRNA, transcription factor and mRNA regulatory network in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is one of the most common malignancies worldwide and the third leading cause of cancer-related death. MicroRNAs and transcription factors (TFs) cooperate to regulate the same target gene, thus affecting the progression of HCC.

Methods

Differentially expressed miRNAs and mRNAs were screened. Functional enrichment analysis of these HCC-related mRNAs was performed, and a protein-protein interaction network was constructed. TFs that regulate these miRNAs and hub genes were also screened.

Results

Ten differentially upregulated miRNAs and 5 differentially downregulated miRNAs were screened. Additionally, 183 downregulated mRNAs and 303 upregulated mRNAs that are potentially bound to these differentially expressed miRNAs were identified. The Kyoto Encyclopedia of Genes and Genomes (KEGG) results showed that the differentially expressed mRNAs were significantly enriched in pathways in cancer, the Wnt signaling pathway, and the Rap1 signaling pathway. Then, 220 TFs were identified for 5 candidate genes of the downregulated mRNAs, and 258 TFs were identified for 9 candidate genes of the upregulated mRNAs. Finally, the 9 upregulated hub genes were related to higher overall survival (OS) in the low-expression group, and 4/5 downregulated hub genes were related to higher OS in the high-expression group.

Conclusions

This study constructed a potential regulatory network between candidate molecules and that need to be further verified. These regulatory relationships are expected to clarify the new molecular mechanisms of the occurrence and development of HCC.

Unravelling host-pathogen interactions: ceRNA network in SARS-CoV-2 infection (COVID-19)

•

Bioinformatics analysis of mouse mRNA expression dataset for presumptive SARS-CoV-2 targets.

•

Induction of ISGs-Isg15, Oasl1, Usp18 and Ddx58 with no marked changes in the expression of IFNs.

•

No induction of ACE2 and TMPRSS2, raising implications for host factor limitations.

•

Identification of ceRNA network including miR-124-3p, Ddx58, lncRNA (Gm26917) and circRNAs (Ppp1r10, C330019G07RiK).

•

Virus regulates the expression of lnc and circRNAs, acting as sponges for miR-124-3p targeting Ddx58.

COVID-19 is a lurking calamitous disease caused by an unusual virus, SARS-CoV-2, causing massive deaths worldwide. Nonetheless, explicit therapeutic drugs or clinically approved vaccines are not available for COVID-19. Thus, a comprehensive research is crucially needed to decode the pathogenic tools, plausible drug targets, committed to the development of efficient therapy. Host-pathogen interactions via host cellular components is an emerging field of research in this respect. miRNAs have been established as vital players in host-virus interactions. Moreover, viruses have the capability to manoeuvre the host miRNA networks according to their own obligations. Besides protein coding mRNAs, noncoding RNAs might also be targeted in infected cells and viruses can exploit the host miRNA network via ceRNA effect. We have predicted a ceRNA network involving one miRNA (miR-124-3p), one mRNA (Ddx58), one lncRNA (Gm26917) and two circRNAs (Ppp1r10, C330019G07RiK) in SARS-CoV infected cells. We have identified 4 DEGs-Isg15, Ddx58, Oasl1, Usp18 by analyzing a mRNA GEO dataset. There is no notable induction of IFNs and IFN-induced ACE2, significant receptor responsible for S-protein binding mediated viral entry. Pathway enrichment and GO analysis conceded the enrichment of pathways associated with interferon signalling and antiviral-mechanism by IFN-stimulated genes. Further, we have identified 3 noncoding RNAs, playing as potential ceRNAs to the genes associated with immune mechanisms. This integrative analysis has identified noncoding RNAs and their plausible targets, which could effectively enhance the understanding of molecular mechanisms associated with viral infection. However, validation of these targets is further corroborated to determine their therapeutic efficacy.

Transcriptomic analysis delineates potential signature genes and miRNAs associated with the pathogenesis of asthma

Asthma is a multifarious disease affecting several million people around the world. It has a heterogeneous risk architecture inclusive of both genetic and environmental factors. This heterogeneity can be utilised to identify differentially expressed biomarkers of the disease, which may ultimately aid in the development of more localized and molecularly targeted therapies. In this respect, our study complies with meta-analysis of microarray datasets containing mRNA expression profiles of both asthmatic and control patients, to identify the critical Differentially Expressed Genes (DEGs) involved in the pathogenesis of asthma. We found a total of 30 DEGs out of which 13 were involved in the pathway and functional enrichment analysis. Moreover, 5 DEGs were identified as the hub genes by network centrality-based analysis. Most hub genes were involved in protease/antiprotease pathways. Also, 26 miRNAs and 20 TFs having an association with these hub genes were found to be intricated in a 3-node miRNA Feed-Forward Loop. Out of these, miR-34b and miR-449c were identified as the key miRNAs regulating the expression of SERPINB2 gene and SMAD4 transcription factor. Thus, our study is suggestive of certain miRNAs and unexplored pathways which may pave a way to unravel critical therapeutic targets in asthma.

Critical microRNAs and regulatory motifs in cleft palate identified by a conserved miRNA-TF-gene network approach in humans and mice

Cleft palate (CP) is the second most common congenital birth defect. The etiology of CP is complicated, with involvement of various genetic and environmental factors. To investigate the gene regulatory mechanisms, we designed a powerful regulatory analytical approach to identify the conserved regulatory networks in humans and mice, from which we identified critical microRNAs (miRNAs), target genes and regulatory motifs (miRNA-TF-gene) related to CP. Using our manually curated genes and miRNAs with evidence in CP in humans and mice, we constructed miRNA and transcription factor (TF) co-regulation networks for both humans and mice. A consensus regulatory loop (miR17/miR20a-FOXE1-PDGFRA) and eight miRNAs (miR-140, miR-17, miR-18a, miR-19a, miR-19b, miR-20a, miR-451a and miR-92a) were discovered in both humans and mice. The role of miR-140, which had the strongest association with CP, was investigated in both human and mouse palate cells. The overexpression of miR-140-5p, but not miR-140-3p, significantly inhibited cell proliferation. We further examined whether miR-140 overexpression could suppress the expression of its predicted target genes (BMP2, FGF9, PAX9 and PDGFRA). Our results indicated that miR-140-5p overexpression suppressed the expression of BMP2 and FGF9 in cultured human palate cells and Fgf9 and Pdgfra in cultured mouse palate cells. In summary, our conserved miRNA-TF-gene regulatory network approach is effective in detecting consensus miRNAs, motifs, and regulatory mechanisms in human and mouse CP.

Subtype-specific signaling pathways and genomic aberrations associated with prognosis of glioblastoma

Background

A high heterogeneity and activation of multiple oncogenic pathways have been implicated in failure of targeted therapies in glioblastoma (GBM).

Methods

Using The Cancer Genome Atlas data, we identified subtype-specific prognostic core genes by a combined approach of genome-wide Cox regression and Gene Set Enrichment Analysis. The results were validated with 8 combined public datasets containing 608 GBMs. We further examined prognostic chromosome aberrations and mutations.

Results

In classical and mesenchymal subtypes, 2 receptor tyrosine kinases (RTKs) (MET and IGF1R), and the genes in RTK downstream pathways such as phosphatidylinositol-3 kinase (PI3K)/Akt/mammalian target of rapamycin (mTOR), and nuclear factor-kappaB (NF-kB), were commonly detected as prognostic core genes. Classical subtype-specific prognostic core genes included those in cell cycle, DNA repair, and the Janus kinase/signal transducers and activators of transcription (JAK-STAT) pathway. Immune-related genes were enriched in the prognostic genes showing negative promoter cytosine-phosphate-guanine (CpG) methylation/expression correlations. Mesenchymal subtype-specific prognostic genes were those related to mesenchymal cell movement, PI3K/Akt, mitogen-activated protein kinase (MAPK)/extracellular signal-regulated kinase (ERK), Wnt/β-catenin, and Wnt/Ca2+ pathways. In copy number alterations and mutations, 6p loss and TP53 mutation were associated with poor and good survival, respectively, in the classical subtype. In the mesenchymal subtype, patients with PIK3R1 or PCLO mutations showed poor prognosis. In the glioma CpG island methylator phenotype (G-CIMP) subtype, patients harboring 10q loss, 12p gain, or 14q loss exhibited poor survival. Furthermore, 10q loss was significantly associated with the recently recognized G-CIMP subclass showing relatively low CpG methylation and poor prognosis.

Conclusion

These subtype-specific alterations have promising potentials as new prognostic biomarkers and therapeutic targets combined with surrogate markers of GBM subtypes. However, considering the small number of events, the results of copy number alterations and mutations require further validations.

TFmiR2: constructing and analyzing disease-, tissue- and process-specific transcription factor and microRNA co-regulatory networks

SUMMARY

TFmiR2 is a freely available web server for constructing and analyzing integrated transcription factor (TF) and microRNA (miRNA) co-regulatory networks for human and mouse. TFmiR2 generates tissue- and biological process-specific networks for the set of deregulated genes and miRNAs provided by the user. Furthermore, the service can now identify key driver genes and miRNAs in the constructed networks by utilizing the graph theoretical concept of a minimum connected dominating set. These putative key players as well as the newly implemented four-node TF-miRNA motifs yield novel insights that may assist in developing new therapeutic approaches.

AVAILABILITY AND IMPLEMENTATION

The TFmiR2 web server is available at http://service.bioinformatik.uni-saarland.de/tfmir2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Micro RNA Molecules as Modulators of Treatment Resistance, Immune Checkpoints Controllers and Sensitive Biomarkers in Glioblastoma Multiforme

Based on genome sequencing, it is estimated that over 90% of genes stored in human genetic material are transcribed, but only 3% of them contain the information needed for the production of body proteins. This group also includes micro RNAs representing about 1%–3% of the human genome. Recent studies confirmed the hypothesis that targeting molecules called Immune Checkpoint (IC) open new opportunities to take control over glioblastoma multiforme (GBM). Detection of markers that indicate the presence of the cancer occupies a very important place in modern oncology. This function can be performed by both the cancer cells themselves as well as their components and other substances detected in the patients’ bodies. Efforts have been made for many years to find a suitable marker useful in the diagnosis and monitoring of gliomas, including glioblastoma.

MicroRNA and transcription factor co-regulatory networks and subtype classification of seminoma and non-seminoma in testicular germ cell tumors

Recent studies have revealed that feed-forward loops (FFLs) as regulatory motifs have synergistic roles in cellular systems and their disruption may cause diseases including cancer. FFLs may include two regulators such as transcription factors (TFs) and microRNAs (miRNAs). In this study, we extensively investigated TF and miRNA regulation pairs, their FFLs, and TF-miRNA mediated regulatory networks in two major types of testicular germ cell tumors (TGCT): seminoma (SE) and non-seminoma (NSE). Specifically, we identified differentially expressed mRNA genes and miRNAs in 103 tumors using the transcriptomic data from The Cancer Genome Atlas. Next, we determined significantly correlated TF-gene/miRNA and miRNA-gene/TF pairs with regulation direction. Subsequently, we determined 288 and 664 dysregulated TF-miRNA-gene FFLs in SE and NSE, respectively. By constructing dysregulated FFL networks, we found that many hub nodes (12 out of 30 for SE and 8 out of 32 for NSE) in the top ranked FFLs could predict subtype-classification (Random Forest classifier, average accuracy ≥90%). These hub molecules were validated by an independent dataset. Our network analysis pinpointed several SE-specific dysregulated miRNAs (miR-200c-3p, miR-25-3p, and miR-302a-3p) and genes (EPHA2, JUN, KLF4, PLXDC2, RND3, SPI1, and TIMP3) and NSE-specific dysregulated miRNAs (miR-367-3p, miR-519d-3p, and miR-96-5p) and genes (NR2F1 and NR2F2). This study is the first systematic investigation of TF and miRNA regulation and their co-regulation in two major TGCT subtypes.

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

BACKGROUND

Feed-forward loops (FFLs), consisting of miRNAs, transcription factors (TFs) and their common target genes, have been validated to be important for the initialization and development of complex diseases, including cancer. Esophageal Carcinoma (ESCA) and Stomach Adenocarcinoma (STAD) are two types of malignant tumors in the digestive tract. Understanding common and distinct molecular mechanisms of ESCA and STAD is extremely crucial.

RESULTS

In this paper, we presented a computational framework to explore common and distinct FFLs, and molecular biomarkers for ESCA and STAD. We identified FFLs by combining regulation pairs and RNA-seq data. Then we constructed disease-specific co-expression networks based on the FFLs identified. We also used random walk with restart (RWR) on disease-specific co-expression networks to prioritize candidate molecules. We identified 148 and 242 FFLs for these two types of cancer, respectively. And we found that one TF, E2F3 was related to ESCA, two genes, DTNA and KCNMA1 were related to STAD, while one TF ESR1 and one gene KIT were associated with both of the two types of cancer.

CONCLUSIONS

This proposed computational framework predicted disease-related biomolecules effectively and discovered the correlation between two types of cancers, which helped develop the diagnostic and therapeutic strategies of Esophageal Carcinoma and Stomach Adenocarcinoma.

Systems biology-based investigation of cooperating microRNAs as monotherapy or adjuvant therapy in cancer

MicroRNAs (miRNAs) are short, noncoding RNAs that regulate gene expression by suppressing mRNA translation and reducing mRNA stability. A miRNA can potentially bind many mRNAs, thereby affecting the expression of oncogenes and tumor suppressor genes as well as the activity of whole pathways. The promise of miRNA therapeutics in cancer is to harness this evolutionarily conserved mechanism for the coordinated regulation of gene expression, and thus restoring a normal cell phenotype. However, the promiscuous binding of miRNAs can provoke unwanted off-target effects, which are usually caused by high-dose single-miRNA treatments. Thus, it is desirable to develop miRNA therapeutics with increased specificity and efficacy. To achieve that, we propose the concept of miRNA cooperativity in order to exert synergistic repression on target genes, thus lowering the required total amount of miRNAs. We first review miRNA therapies in clinical application. Next, we summarize the knowledge on the molecular mechanism and biological function of miRNA cooperativity and discuss its application in cancer therapies. We then propose and discuss a systems biology approach to investigate miRNA cooperativity for the clinical setting. Altogether, we point out the potential of miRNA cooperativity to reduce off-target effects and to complement conventional, targeted, or immune-based therapies for cancer.

Angiotensin II-induced micro RNA-21 culprit for non-small-cell lung adenocarcinoma

Lung cancer is among the most complicated cancers, with an estimated 1.6 million deaths each year for both men and women. However, the proportion of lung cancer patients in developing nations has increased from 31% to 49.9% in the last two decades. There are two main subtypes of lung cancer, small-cell lung carcinoma and non-small-cell lung carcinoma (NSCLC), accounting for 15% and 85% of all lung cancer, respectively. Adenocarcinoma is the most common type of lung cancer in smokers and nonsmokers in men and women regardless of their age. Chemicals in cigarette smoke and nicotine enter our bloodstream and can then affect the entire body and finally lead to the activation of several important, pro-survival signaling pathways. The biologically active peptide of RAAS on overstimulation enhance Ang II mediates cell proliferation, fibrosis and inflammatory effects via AT1 receptor. Very few studies highlight the diagnostic and therapeutic potential of miRNAs with the EGFR-regulated miRNA-21.

Therapeutic Targeting of the Notch Pathway in Glioblastoma Multiforme

BACKGROUND

Glioblastoma (GBM) is the most common and deadly form of brain tumor. After standard treatment of resection, radiotherapy, and chemotherapy, the 5-year survival is <5%. In recent years, research has uncovered several potential targets within the Notch signaling pathway, which may lead to improved patient outcomes.

METHODS

A literature search was performed for articles containing the terms "Glioblastoma" and "Receptors, Notch" between 2003 and July 2015. Of the 62 articles retrieved, 46 met our criteria and were included in our review. Nine articles were identified from other sources and were subsequently included, leaving 55 articles reviewed.

RESULTS

Of the 55 articles reviewed, 47 used established human GBM cell lines. Seventeen articles used human GBM surgical samples. Forty-five of 48 articles that assessed Notch activity showed increased expression in GBM cell lines. Targeting the Notch pathway was carried out through Notch knockdown and overexpression and targeting δ-like ligand, Jagged, γ-secretase, ADAM10, ADAM17, and Mastermindlike protein 1. Arsenic trioxide, microRNAs, and several other compounds were shown to have an effect on the Notch pathway in GBM. Notch activity in GBM was also shown to be associated with hypoxia and certain cancer-related molecular pathways such as PI3K/AKT/mTOR and ERK/MAPK. Most articles concluded that Notch activity amplifies malignant characteristics in GBM and targeting this pathway can bring about amelioration of these effects.

CONCLUSIONS

Recent literature suggests targeting the Notch pathway has great potential for future therapies for GBM.

MicroRNA-655-3p and microRNA-497-5p inhibit cell proliferation in cultured human lip cells through the regulation of genes related to human cleft lip

Background

The etiology of cleft lip with or without cleft palate (CL/P), a common congenital birth defect, is complex and involves the contribution of genetic and environmental factors. Although many candidate genes have been identified, the regulation and interaction of these genes in CL/P remain unclear. In addition, the contribution of microRNAs (miRNAs), non-coding RNAs that regulate the expression of multiple genes, to the etiology of CL/P is largely unknown.

Methods

To identify the signatures of causative biological pathways for human CL/P, we conducted a systematic literature review for human CL/P candidate genes and subsequent bioinformatics analyses. Functional enrichment analyses of the candidate CL/P genes were conducted using the pathway databases GO and KEGG. The miRNA-mediated post-transcriptional regulation of the CL/P candidate genes was analyzed with miRanda, PITA, and TargetScan, and miRTarbase. Genotype-phenotype association analysis was conducted using GWAS. The functional significance of the candidate miRNAs was evaluated experimentally in cell proliferation and target gene regulation assays in human lip fibroblasts.

Results

Through an extensive search of the main biomedical databases, we mined 177 genes with mutations or association/linkage reported in individuals with CL/P, and considered them as candidate genes for human CL/P. The genotype-phenotype association study revealed that mutations in 12 genes (ABCA4, ADAM3A, FOXE1, IRF6, MSX2, MTHFR, NTN1, PAX7, TP63, TPM1, VAX1, and WNT9B) were significantly associated with CL/P. In addition, our bioinformatics analysis predicted 16 microRNAs (miRNAs) to be post-transcriptional regulators of CL/P genes. To validate the bioinformatics results, the top six candidate miRNAs (miR-124-3p, miR-369-3p, miR-374a-5p, miR-374b-5p, miR-497-5p, and miR-655-3p) were evaluated by cell proliferation/survival assays and miRNA-gene regulation assays in cultured human lip fibroblasts. We found that miR-497-5p and miR-655-3p significantly suppressed cell proliferation in these cells. Furthermore, the expression of the predicted miRNA-target genes was significantly downregulated by either miR-497-5p or miR-655-3p mimic.

Conclusion

Expression of miR-497-5p and miR-655-3p suppresses cell proliferation through the regulation of human CL/P-candidate genes. This study provides insights into the role of miRNAs in the etiology of CL/P and suggests possible strategies for the diagnosis of CL/P.

Electronic supplementary material

The online version of this article (10.1186/s12920-019-0535-2) contains supplementary material, which is available to authorized users.

Role of Notch Signaling Pathway in Glioblastoma Pathogenesis

Notch signaling is an evolutionarily conserved pathway that regulates important biological processes, such as cell proliferation, apoptosis, migration, self-renewal, and differentiation. In mammals, Notch signaling is composed of four receptors (Notch1–4) and five ligands (Dll1-3–4, Jagged1–2) that mainly contribute to the development and maintenance of the central nervous system (CNS). Neural stem cells (NSCs) are the starting point for neurogenesis and other neurological functions, representing an essential aspect for the homeostasis of the CNS. Therefore, genetic and functional alterations to NSCs can lead to the development of brain tumors, including glioblastoma. Glioblastoma remains an incurable disease, and the reason for the failure of current therapies and tumor relapse is the presence of a small subpopulation of tumor cells known as glioma stem cells (GSCs), characterized by their stem cell-like properties and aggressive phenotype. Growing evidence reveals that Notch signaling is highly active in GSCs, where it suppresses differentiation and maintains stem-like properties, contributing to Glioblastoma tumorigenesis and conventional-treatment resistance. In this review, we try to give a comprehensive view of the contribution of Notch signaling to Glioblastoma and its possible implication as a target for new therapeutic approaches.

Network-based identification of critical regulators as putative drivers of human cleft lip

BACKGROUND

Cleft lip (CL) is one of the most common congenital birth defects with complex etiology. While genome-wide association studies (GWAS) have made significant advances in our understanding of mutations and their related genes with potential involvement in the etiology of CL, it remains unknown how these genes are functionally regulated and interact with each other in lip development. Currently, identifying the disease-causing genes in human CL is urgently needed. So far, the causative CL genes have been largely undiscovered, making it challenging to design experiments to validate the functional influence of the mutations identified from large genomic studies such as CL GWAS.

RESULTS

Transcription factors (TFs) and microRNAs (miRNAs) are two important regulators in cellular system. In this study, we aimed to investigate the genetic interactions among TFs, miRNAs and the CL genes curated from the previous studies. We constructed miRNA-TF co-regulatory networks, from which the critical regulators as putative drivers in CL were examined. Based on the constructed networks, we identified ten critical hub genes with prior evidence in CL. Furthermore, the analysis of partitioned regulatory modules highlighted a number of biological processes involved in the pathology of CL, including a novel pathway "Signaling pathway regulating pluripotency of stem cells". Our subnetwork analysis pinpointed two candidate miRNAs, hsa-mir-27b and hsa-mir-497, activating the Wnt pathway that was associated with CL. Our results were supported by an independent gene expression dataset in CL.

CONCLUSIONS

This study represents the first regulatory network analysis of CL genes. Our work presents a global view of the CL regulatory network and a novel approach on investigating critical miRNAs, TFs and genes via combinatory regulatory networks in craniofacial development. The top genes and miRNAs will be important candidates for future experimental validation of their functions in CL.

BiModule: biclique modularity strategy for identifying transcription factor and microRNA co-regulatory modules

Systematic identification of gene regulatory modules can provide invaluable knowledge towards understanding aberrant transcriptional/post-transcriptional collaborative regulatory (co-regulatory) effects in cancer. Transcription factor (TF) and microRNA (miRNA) are known as two classes of prominent regulators that play crucial roles in gene regulation. Existing studies on gene regulatory modules identification mainly focused on the miRNA-mediated regulatory network, and few considered these two regulators in a co-occurring network. In this current study, we developed a computational method called BiModule for systematically identifying TF-miRNA co-regulatory modules. BiModule operates in two main stages: it first constructs a cancerspecific regulator-mRNA network and then identifies modules based on maximal bicliques by employing biclique modularity strategy, which is a novel flexible method for bipartite graph mining. We applied our model to a cervical cancer dataset. The results showed that the TF-miRNA co-regulatory modules identified by BiModule exhibit denser connections and stronger expression correlations than another existing related method. Moreover, the BiModule-modules exhibit high biological functional enrichment. In addition, based on Kaplan-Meier survival analysis, we found a number of modules with significant prognostic associations.

MicroRNA-365 suppressed cell proliferation and migration via targeting PAX6 in glioblastoma

MicroRNAs (miRNAs) act an important role in the progression of tumor. In this study, we showed that the serum expression of miR-365 was downregulated in the glioblastoma compared with in the healthy controls. We also demonstrated that miR-365 expression was downregulated in glioblastoma tissues compared with the adjacent normal tissues. Overexpression of miR-365 suppressed the glioblastoma cell proliferation and migration. Moreover, ectopic expression of miR-365 promoted the expression of Ecadherin while inhibited the expression of N-cadherin and Vimentin in U87 cell. Furthermore, we identified PAX6 as a direct target gene of miR-365 in U87 cell. Overexpression of miR-365 suppressed glioblastoma cell proliferation and migration and epithelial-to-mesenchymal transition through inhibiting PAX6 expression. These results suggested that miR-365 played a tumor suppressor in glioblastoma.

CMTCN: a web tool for investigating cancer-specific microRNA and transcription factor co-regulatory networks

Transcription factors (TFs) and microRNAs (miRNAs) are well-characterized trans-acting essential players in gene expression regulation. Growing evidence indicates that TFs and miRNAs can work cooperatively, and their dysregulation has been associated with many diseases including cancer. A unified picture of regulatory interactions of these regulators and their joint target genes would shed light on cancer studies. Although online resources developed to support probing of TF-gene and miRNA-gene interactions are available, online applications for miRNA-TF co-regulatory analysis, especially with a focus on cancers, are lacking. In light of this, we developed a web tool, namely CMTCN (freely available at http://www.cbportal.org/CMTCN), which constructs miRNA-TF co-regulatory networks and conducts comprehensive analyses within the context of particular cancer types. With its user-friendly provision of topological and functional analyses, CMTCN promises to be a reliable and indispensable web tool for biomedical studies.

Systematic identification and analysis of dysregulated miRNA and transcription factor feed-forward loops in hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is the most common genetic cardiovascular disease. Although some genes and miRNAs related with HCM have been studied, the molecular regulatory mechanisms between miRNAs and transcription factors (TFs) in HCM have not been systematically elucidated. In this study, we proposed a novel method for identifying dysregulated miRNA‐TF feed‐forward loops (FFLs) by integrating sample matched miRNA and gene expression profiles and experimentally verified interactions of TF‐target gene and miRNA‐target gene. We identified 316 dysregulated miRNA‐TF FFLs in HCM, which were confirmed to be closely related with HCM from various perspectives. Subpathway enrichment analysis demonstrated that the method was outperformed by the existing method. Furthermore, we systematically analysed the global architecture and feature of gene regulation by miRNAs and TFs in HCM, and the FFL composed of hsa‐miR‐17‐5p, FASN and STAT3 was inferred to play critical roles in HCM. Additionally, we identified two panels of biomarkers defined by three TFs (CEBPB, HIF1A, and STAT3) and four miRNAs (hsa‐miR‐155‐5p, hsa‐miR‐17‐5p, hsa‐miR‐20a‐5p, and hsa‐miR‐181a‐5p) in a discovery cohort of 126 samples, which could differentiate HCM patients from healthy controls with better performance. Our work provides HCM‐related dysregulated miRNA‐TF FFLs for further experimental study, and provides candidate biomarkers for HCM diagnosis and treatment.

MicroRNA-transcription factor network analysis reveals miRNAs cooperatively suppress RORA in oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) represents over 90% of oral cancer incidence, while its mechanisms of tumorigenesis remain poorly characterized. In this study, we applied RNA-seq and microRNA-seq methodologies in four pairs of cancer and adjacent normal tissues to profile the contribution of miRNAs to tumorigenesis-altered functional pathways by constructing a comprehensive miRNA-mediated mRNA regulatory network. There were 213 differentially expressed (DE) miRNAs and 2172 DE mRNAs with the involvement of negative miRNA-mRNA interactions identified by at least two pairs of cancerous tissues. GO analysis revealed that the upregulated microRNAs significantly contributed to a global down-regulation of a number of transcription factors (TFs) in OSCC. Among the negative regulatory networks between the selected miRNAs (133) and TFs (167), circadian rhythm genes (RORA, RORB, RORC, and CLOCK) simultaneously regulated by multiple microRNAs were of particular interest. For instance, RORA transcript was predicted to be targeted by 25 co-upregulated miRNAs, of which, miR-503-5p, miR-450b-5p, miR-27a-3p, miR-181a-5p and miR-183-5p were further validated to directly target RORA, resulting in a stronger effect on RORA suppression together. In addition, we showed that the mRNA and protein expression levels of RORα were significantly decreased in most OSCC samples, associated with advanced clinical stage and poor prognosis. RORα significantly suppressed the proliferation of OSCC cells in vitro and in vivo. Attenuated RORα decreased p53 protein expression and suppressed p53 phosphorylation activity. Altogether, our results strongly suggest the importance of the role of miRNAs in regulating the activity of circadian rhythm-related TFs network during OSCC tumorigenesis, and provide further clues to understand the clinical link between circadian rhythm and cancer therapy.

Molecular mechanism of Notch signaling with special emphasis on microRNAs: Implications for glioma

Glioma is the most aggressive primary brain tumor and is notorious for resistance to chemoradiotherapy. Although its associated mechanisms are still not completely understood, Notch signaling, an evolutionarily conserved pathway, appears to be the key processes involved. Nevertheless, its mechanisms are sophisticated, due to a variety of targets and signal pathways, especially microRNA. MicroRNAs, which are small noncoding regulatory RNA molecules, have been proposed as one of the key mechanisms in glioma pathogenesis. Among the known glioma associated microRNA, microRNA-129, microRNA-34 family, and microRNA-326 have been shown to influence the progress of glioma through Notch signaling. Evidence also indicates that recurrence is due to development or persistence of the glioma stem-like cells and active angiogenesis, which are tightly regulated by a variety of factors, including Notch signaling. In this review, we summarize the recent progress regarding the functional roles of Notch signaling in glioma, including Notch ligand, microRNA, intracellular crosstalk, glioma stem-like cells and active angiogenesis and explore their clinical implications as diagnostic or prognostic biomarkers and molecular therapeutic targets for glioma.