Target Regulatory Network Research Articles

Comprehensive analysis of the REST transcription factor regulatory networks in IDH mutant and IDH wild-type glioma cell lines and tumors

The RE1-silencing transcription factor (REST) acts either as a repressor or activator of transcription depending on the genomic and cellular context. REST is a key player in brain cell differentiation by inducing chromatin modifications, including DNA methylation, in a proximity of its binding sites. Its dysfunction may contribute to oncogenesis. Mutations in IDH1/2 significantly change the epigenome contributing to blockade of cell differentiation and glioma development. We aimed at defining how REST modulates gene activation and repression in the context of the IDH mutation-related phenotype in gliomas. We studied the effects of REST knockdown, genome wide occurrence of REST binding sites, and DNA methylation of REST motifs in IDH wild type and IDH mutant gliomas. We found that REST target genes, REST binding patterns, and TF motif occurrence proximal to REST binding sites differed in IDH wild-type and mutant gliomas. Among differentially expressed REST targets were genes involved in glial cell differentiation and extracellular matrix organization, some of which were differentially methylated at promoters or gene bodies. REST knockdown differently impacted invasion of the parental or IDH1 mutant glioma cells. The canonical REST-repressed gene targets showed significant correlation with the GBM NPC-like cellular state. Interestingly, results of REST or KAISO silencing suggested the interplay between these TFs in regulation of REST-activated and repressed targets. The identified gene regulatory networks and putative REST cooperativity with other TFs, such as KAISO, show distinct REST target regulatory networks in IDH-WT and IDH-MUT gliomas, without concomitant DNA methylation changes. We conclude that REST could be an important therapeutic target in gliomas.

Integrative analyses of ferroptosis and immune related biomarkers and the osteosarcoma associated mechanisms

Osteosarcoma (OS) is the most common primary malignant bone tumor with high metastatic potential and relapse risk. To study the regulatory mechanism of the OS microenvironment, a complex regulatory network involving the ferroptosis- and immune response-related genes remains to be established. In the present study, we determined the effect of a comprehensive evaluation system established on the basis of ferroptosis- and immune-related genes on the immune status, related biomarkers, prognosis, and the potential regulatory networks underlying OS based on the TARGET and Gene Expression Omnibus databases that contain information on OS patients by bioinformatics analyses. We first characterized individual ferroptosis scores and immune scores through gene set variation analysis (GSVA) against TARGET-OS datasets. We then identified differentially expressed genes by score groups. Weighted gene co-expression network analysis was performed to identify the most relevant ferroptosis-related and immune-related gene modules, which facilitated the identification of 327 ferroptosis gene and 306 immune gene candidates. A 4-gene (WAS, CORT, WNT16, and GLB1L2) signature was constructed and valuation using the least absolute shrinkage and selection operator-Cox regression models to effectively predict OS prognosis. The prediction efficiency was further validated by GSE39055. We stratified patients based on the prognostic scoring systems. Eight hub genes (namely CD3D, CD8A, CD3E, IL2, CD2, MYH6, MYH7, and MYL2) were identified, and TF–miRNA target regulatory networks were constructed. Furthermore, Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, gene set enrichment analysis, and GSVA were used to determine the signature’s potential pathways and biological functions, which showed that the hub genes were enriched in ferroptosis-associated biological functions and immune-associated molecular mechanisms. Thereafter, we investigated the proportion and infiltration extent of 22 infiltrating immune cells by using CIBERSORT, which revealed significant subgroup differences in CD8 + T cells, M0 macrophages, and M2 macrophages. In conclusion, we determined a new ferroptosis-related and immune-related gene signature for predicting OS patients’ prognosis and further explored the ferroptosis and immunity interactions during OS development, which provides insights into the exploration of molecular mechanisms and targeted therapies in patients with OS.

In silico study of novel niclosamide derivatives, SARS-CoV-2 nonstructural proteins catalytic residue-targeting small molecules drug candidates

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2)-mediated coronavirus disease 2019 (COVID-19) infection remains a global pandemic and health emergency with overwhelming social and economic impacts throughout the world. Therapeutics for COVID-19 are limited to only remdesivir; therefore, there is a need for combined, multidisciplinary efforts to develop new therapeutic molecules and explore the effectiveness of existing drugs against SARS-CoV-2. In the present study, we reported eight (SCOV-L-02, SCOV-L-09, SCOV-L-10, SCOV-L-11, SCOV-L-15, SCOV-L-18, SCOV-L-22, and SCOV-L-23) novel structurally related small-molecule derivatives of niclosamide (SCOV-L series) for their targeting potential against angiotensin-converting enzyme-2 (ACE2), type II transmembrane serine protease (TMPRSS2), and SARS-COV-2 nonstructural proteins (NSPs) including NSP5 (3CLpro), NSP3 (PLpro), and RdRp. Our correlation analysis suggested that ACE2 and TMPRSS2 modulate host immune response via regulation of immune-infiltrating cells at the site of tissue/organs entries. In addition, we identified someTMPRSS2 and ACE2 microRNAs target regulatory networks in SARS-CoV-2 infection and thus open up a new window for microRNAs-based therapy for the treatment of SARS-CoV-2 infection. Our in vitro study revealed that with the exception of SCOV-L-11 and SCOV-L-23 which were non-active, the SCOV-L series exhibited strict antiproliferative activities and non-cytotoxic effects against ACE2- and TMPRSS2-expressing cells. Our molecular docking for the analysis of receptor-ligand interactions revealed that SCOV-L series demonstrated high ligand binding efficacies (at higher levels than clinical drugs) against the ACE2, TMPRSS2, and SARS-COV-2 NSPs. SCOV-L-18, SCOV-L-15, and SCOV-L-09 were particularly found to exhibit strong binding affinities with three key SARS-CoV-2’s proteins: 3CLpro, PLpro, and RdRp. These compounds bind to the several catalytic residues of the proteins, and satisfied the criteria of drug-like candidates, having good adsorption, distribution, metabolism, excretion, and toxicity (ADMET) pharmacokinetic profile. Altogether, the present study suggests the therapeutic potential of SCOV-L series for preventing and managing SARs-COV-2 infection and are currently under detailed investigation in our lab.

High Expression MicroRNA-206 Inhibits the Growth of Tumor Cells in Human Malignant Fibrous Histiocytoma

Background: Malignant fibrous histiocytoma (MFH) is a common type of soft tissue sarcoma and a serious threat to human health. MFH often relapses locally after the curettage is related to the residual cancer stem cells (CSCs). Currently, the dysregulation of microRNA (miRNA) has been found to be closely related to the recurrence of CSCs. However, whether dysregulations of miRNAs exist in MFH, CSCs remained unknown.Methods: In this study, miRNAs in MFH CSCs and MFH common cells were examined by gene probe. Then, target genes and their functions involved in the signal pathway were predicted by the relevant database. Finally, the miRNAs’ target regulatory network was constructed. Furthermore, the miRNAs and target genes were identified by quantitative polymerase chain reaction, whereas miRNA analogs and antagonists were transfected in tumor cells to investigate cell proliferation ability further.Results: Results showed that a total of 47 miRNAs were found, including 16 that were upregulated and 31 that were downregulated. The screened differential miRNA showed a different expression in the cell resistant strains compared with the control group. Quantitative polymerase chain reaction analysis confirmed that the relative abundance of seven miRNAs and four target genes varied significantly. The encouraging issue is that we found Hsa-miR-206 significantly inhibited MFH proliferative activity.Conclusion: Hsa-miR-206 played a key role in regulating MFH CSC properties that might be a representative marker and target for the diagnosis and treatment of MFH in the future.

Molecular profiles of predictive biomarkers for platinum-based chemotherapy in Non-Small Cell Lung Cancer (NSCLC)

BackgroundNon-small cell lung cancer (NSCLC) is the principal subtype of lung cancer. Among all therapeutic options, platinum-based chemotherapy agents, especially Cisplatin, are still commonly used treatment for NSCLC patients. However, developing chemoresistance in NSCLC cells often gives rise to chemotherapy failure. Therefore, more studies are required to shed light on gene interaction and cellular pathways involved in initiating and developing resistance to platinum-based chemotherapy in NSCLC. Hence, it is urgent to find the key genes, microRNA (miRNAs), and potential molecular mechanisms implicated in chemoresistance and present markers to predict response to platinum-based chemotherapy in NSCLC patients. MethodsThe microarray datasets GSE6410, GSE7035, GSE14814, GSE26704, GSE73302 were downloaded from the Gene Expression Omnibus (GEO) database and were analyzed using R software. Functional and pathway enrichment analyses were performed using the Enrich R site. Then, the protein-protein interaction (PPI) network and hub genes were obtained using the Cytoscape software. Further, the miRSystem database was performed to predict the miRNAs regulating the hub genes. Moreover, Cytoscape software and the CytoHubba plugin were used to construct the miRNA-target interaction network and hub modules. Finally, the Kaplan–Meier curve was used to demonstrate the survival curves and assess the association of the genes signature with clinical outcomes. ResultsA total of 142 differentially expressed genes (DEGs) were found. The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses present the p53 signaling pathway as the most significant pathway.Besides, from the top ten terms obtained of Biological Process, Molecular Function, and Cellular Component, the first ones, including cholesterol biosynthetic process, the extrinsic component of external side of plasma membrane, cytokine activity, were selected respectively. Based on the PPI network, the ten nodes with the highest degree were screened as hub genes. In addition, from the miRNA–target regulatory network in Cytoscape, ten hub nodes were found. Ultimately, according to Kaplan–Meier curve, BTG2 and TP53I3 with p-value <0.05 were associated with a better prognosis. ConclusionsIn the present study, DEGs, candidate miRNAs, and underlying mechanisms involved in chemoresistance were identified to suggest potential biomarkers to provide new clues for the prediction of response to platinum-based chemotherapy.

A Rapid Evolving microRNA Cluster Rewires Its Target Regulatory Networks in Drosophila.

New miRNAs are evolutionarily important but their functional evolution remains unclear. Here we report that the evolution of a microRNA cluster, mir-972C rewires its downstream regulatory networks in Drosophila. Genomic analysis reveals that mir-972C originated in the common ancestor of Drosophila where it comprises six old miRNAs. It has subsequently recruited six new members in the melanogaster subgroup after evolving for at least 50 million years. Both the young and the old mir-972C members evolved rapidly in seed and non-seed regions. Combining target prediction and cell transfection experiments, we found that the seed and non-seed changes in individual mir-972C members cause extensive target divergence among D. melanogaster, D. simulans, and D. virilis, consistent with the functional evolution of mir-972C reported recently. Intriguingly, the target pool of the cluster as a whole remains relatively conserved. Our results suggest that clustering of young and old miRNAs broadens the target repertoires by acquiring new targets without losing many old ones. This may facilitate the establishment of new miRNAs in existing regulatory networks.

Genome-wide identification of miRNAs and target regulatory network in the invasive ectoparasitic mite Varroa destructor

Varroa destructor is an ectoparasite mite that attacks bees leading to colony disorders worldwide. microRNAs (miRNAs) are key molecules used by eukaryotes to post-transcriptional control of gene expression. Nevertheless, still lack information aboutV. destructor miRNAs and its regulatory networks. Here, we used an integrative strategy to characterize the miRNAs in the V. destructor mite. We identified 310 precursors that give rise to 500 mature miRNAs, which 257 are likely mite-specific elements. miRNAs showed canonical length ranging between 18 and 25 nucleotides and 5′ uracil preference. Top 10 elements concentrated over 80% of total miRNA expression, with bantam alone representing ~50%. We also detected non-templated bases in precursor-derived small RNAs, indicative of miRNA post-transcriptional regulatory mechanisms. Finally, we note that conserved miRNAs control similar processes in different organisms, suggesting a conservative role. Altogether, our findings contribute to the better understanding of the mite biology that can assist future studies on varroosis control.

Identification of Pivotal MicroRNAs and Target Genes Associated with Persistent Atrial Fibrillation Based on Bioinformatics Analysis.

Atrial fibrillation (AF) is one of the most common supraventricular arrhythmias worldwide. However, the specific molecular mechanism underlying AF remains unclear. Our study is aimed at identifying pivotal microRNAs (miRNAs) and targeting genes associated with persistent AF (pAF) using bioinformatics analysis. Three gene expression array datasets (GSE31821, GSE41177, and GSE79768) and an miRNA expression array dataset (GSE68475) associated with pAF were downloaded. Differentially expressed genes (DEGs) were identified using the LIMMA package, and differentially expressed miRNAs (DEMs) were screened from GSE68475. Target genes for DEMs were predicted using the miRTarBase database, and intersections between these target genes and DEGs were selected for further analysis, including the generation of protein–protein interaction (PPI) network, miRNA–transcription factor–target regulatory network, and drug–gene network. A total of 264 DEGs and 40 DEMs were identified between the pAF and control groups. Functional and pathway enrichment analyses of up- and downregulated DEGs were performed. The common genes (CGs) were primarily enriched in the phosphoinositide 3-kinase- (PI3K-) protein kinase B (Akt) signaling pathway, negative regulation of cell division, and response to hypoxia. The PPI network, miRNA–transcription factor–target regulatory network, and drug–gene network were constructed using Cytoscape. The present study revealed several novel miRNAs and genes involved in pAF. We speculated that miR-4298, miR-3125, miR-4306, and miR-671-5p could represent significant miRNAs that act on the target gene superoxide dismutase 2 (SOD2) during the development of pAF and may serve as essential biomarkers for pAF diagnosis and treatment. Moreover, MYC might function in pAF pathogenesis through the PI3K–Akt signaling pathway.

Abnormal expression of TGFBR2, EGF, LRP10, and IQGAP1 is involved in the pathogenesis of coronary artery disease.

Coronary artery disease (CAD) is the most common cardiovascular disease worldwide. In this study, we investigated the pathogenesis of CAD. We downloaded the GSE98583 dataset, including 12 CAD samples and 6 normal samples, from the Gene Expression Omnibus (GEO) database and screened differentially expressed genes (DEGs) in CAD versus normal samples. Next, we performed functional enrichment analysis, protein-protein interaction (PPI) network, and functional module analyses to explore potential functions and regulatory functions of identified DEGs. Next, transcription factors (TFs) and microRNAs (miRNAs) targeting DEGs were predicted. In total, 456 DEGs were identified in CAD and normal samples, including 175 upregulated and 281 downregulated genes. These genes were enriched in the intestinal immune network for immunoglobulin A production and the mitogen-activated protein kinase signaling pathway (e.g., TGFBR2 and EGF). The PPI network contained 212 genes, and HIST1H2BJ, HIST1H2AC, EGF, and EP300 were hub genes with degrees higher than 10. Four significant modules were identified from the PPI network, with genes in the modules mainly enriched in the inflammatory response, protein ubiquitination involved in ubiquitin-dependent protein catabolic processes, protein transport, and mitochondrial translational elongation, respectively. Two TFs (E2F1 and FOXK1) and five miRNAs (miR-122A, miR-516-5P, miR-507, miR-342, and miR-520F) were predicted to target 112 DEGs. miR-122A reportedly targets both LRP10 and IQGAP1 in the TF-miRNA target regulatory network. The abnormal expression of TGFBR2, EGF, LRP10, and IQGAP1 may be implicated in CAD pathogenesis. Our study provides targets and potential regulators for investigating CAD pathogenesis.

MiR‑29 mediates exercise‑induced skeletal muscle angiogenesis by targeting VEGFA, COL4A1 and COL4A2 via the PI3K/Akt signaling pathway.

The present study investigated the molecular changes and related regulatory mechanisms in the response of skeletal muscle to exercise. The microarray dataset ‘GSE109657’ of the skeletal muscle response to high-intensity intermittent exercise training (HIIT) was downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were screened and analyzed using weighted gene co-expression network analysis (WGCNA) to identify the significant functional co-expressed gene modules. Moreover, functional enrichment analysis was performed for the DEGs in the significant modules. In addition, protein-protein interaction (PPI) network and microRNA (miR)-transcription factor (TF)-target regulatory network were constructed. A total of 530 DEGs in the skeletal muscle were screened after HIIT, suggesting an effect of HIIT on the skeletal muscle. Moreover, three significant modules (brown, blue and red modules) were identified after WGCNA, and the genes Collagen Type IV α1 Chain (COL4A1) and COL4A2 in the brown module showed the strongest correlation with HIIT. The DEGs in the three modules were significantly enriched in focal adhesion, extracellular matrix organization and the PI3K/Akt signaling pathway. Furthermore, the PPI network contained 104 nodes and 211 interactions. Vascular endothelial growth factor A (VEGFA), COL4A1 and COL4A2 were the hub genes in the PPI network, and were all regulated by miR-29a/b/c. In addition, VEGFA, COL4A1 and COL4A2 were significantly upregulated in the skeletal muscle response to HIIT. Therefore, the present results suggested that the growth and migration of vascular endothelial cells, and skeletal muscle angiogenesis may be regulated by miR-29a/b/c targeting VEGFA, COL4A1 and COL4A2 via the PI3K/Akt signaling pathway. The present results may provide a theoretical basis to investigate the effect of exercise on skeletal muscle.

MAP‑1B, PACS‑2 and AHCYL1 are regulated by miR‑34A/B/C and miR‑449 in neuroplasticity following traumatic spinal cord injury in rats: Preliminary explorative results from microarray data.

Spinal cord injury (SCI) is a specific type of damage to the central nervous system causing temporary or permanent changes in its function. The present aimed to identify the genetic changes in neuroplasticity following SCI in rats. The GSE52763 microarray dataset, which included 15 samples [3 sham (1 week), 4 injury only (1 week), 4 injury only (3 weeks), 4 injury + treadmill (3 weeks)] was downloaded from the Gene Expression Omnibus database. An empirical Bayes linear regression model in limma package was used to identify the differentially expressed genes (DEGs) in injury vs. sham and treadmill vs. non-treadmill comparison groups. Subsequently, time series and enrichment analyses were performed using pheatmap and clusterProfile packages, respectively. Additionally, protein-protein interaction (PPI) and transcription factor (TF)-microRNA (miRNA)-target regulatory networks were constructed using Cytoscape software. In total, 159 and 105 DEGs were identified in injury vs. sham groups and treadmill vs. non-treadmill groups, respectively. There were 40 genes in cluster 1 that presented increased expression levels in the injury (1 week/3 weeks) groups compared with the sham group, and decreased expression levels in the injury + treadmill group compared with the injury only groups; conversely, 52 genes in cluster 2 exhibited decreased expression levels in the injury (1 week/3 weeks) groups compared with the sham group, and increased expression levels in the injury + treadmill group compared with the injury only groups. Enrichment analysis indicated that clusters 1 and 2 were associated with immune response and signal transduction, respectively. Furthermore, microtubule associated protein 1B, phosphofurin acidic cluster sorting protein 2 and adenosylhomocysteinase-like 1 exhibited the highest degrees in the regulatory network, and were regulated by miRNAs including miR-34A, miR-34B, miR-34C and miR-449. These miRNAs and their target genes may serve important roles in neuroplasticity following traumatic SCI in rats. Nevertheless, additional in-depth studies are required to confirm these data.

Integrated analysis of miRNA and mRNA expression profiles in 2-, 6-, and 12-month-old Small Tail Han Sheep ovaries reveals that oar-miR-432 downregulates RPS6KA1 expression

Small Tail Han Sheep are an excellent local sheep breed in China, and their outstanding reproductive performance is one of their very important biological characteristics. Clarifying the ovary development process of these ewes should provide a theoretical basis for improving their reproductive efficiency. In this study, we identified the differentially expressed (DE) microRNAs (miRNAs) in 2-, 6-, and 12-month-old small-tail Han sheep ovaries by constructing and analyzing the miRNA expression profiles. These findings clarify the molecular mechanisms regulating the excellent reproductive performance of small-tail Han ewes. We used RNA-Seq technology and bioinformatic to analyze these profiles. Eleven, 13, and 19 DE miRNAs were identified in the 2- vs 6-, 6- vs 12-, and 2- vs 12-month-old ovaries, respectively. In total, 54, 37, and 198 predicted target genes of these DE miRNAs were identified in these three groups, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that in the 2- vs 6-month-old ovaries, the target genes of DE known sheep miRNAs were involved in 102 GO terms and seven signaling pathways; in the 6- vs 12-month-old ovaries, the target genes of DE known sheep miRNAs were involved in 52 GO terms and three signaling pathways; and in the 2- vs 12-month-old ovaries, the target genes of DE known sheep miRNAs were involved in 88 GO terms and six signaling pathways. Three miRNA–target regulatory networks were constructed based on these DE miRNA–target interactions. Nine miRNAs were selected to confirm to the accuracy of the miRNA sequencing data with qRT-PCR. The site at which oar-miR-432 binds RPS6KA1 was determined with a dual-luciferase system. This is the first integrated analysis the expression profiles of miRNAs and their targets during ovarian development in small-tail Han sheep. These data clarify the molecular regulatory mechanisms underlying sheep ovarian development and identify biomarkers that influence the reproductive performance of small-tail Han ewes.

Integrated Analysis of miRNA and mRNA Expression Profiles Reveals Functional miRNA-Targets in Development Testes of Small Tail Han Sheep.

Small Tail Han Sheep is a highly valued local breed in China because of their precocity, perennial estrus, and high fecundity. The average annual lambing rate of ewes is as high as 180–270%, the semen of ram has characteristics of high yield, high density, and good motility. To reveal the key miRNAs and miRNA-targets underlying testis development and spermatogenesis in male Small Tail Han Sheep, integrated analysis of miRNA and mRNA expression profiles in 2-, 6-, and 12-month-old testes was performed by RNA-seq technology and bioinformatics methods. The results showed that total of 153 known sheep miRNAs and 2712 novel miRNAs were obtained in 2-,6 - and 12-month-old Small Tail Han Sheep testes; 5, 1, and 4 differentially expressed (DE) known sheep miRNAs, and 132, 105, and 24 DE novel miRNAs were identified in 2- vs. 6-, 6- vs. 12-, and 2- vs. 12-month-old testes, respectively. We combined miRNA results of this study and the mRNA results obtained in our previous study to predict the target mRNAs of DE known sheep miRNAs; 131, 10, and 15 target mRNAs of DE known sheep miRNAs and 76, 1, and 11 DE miRNA–targets were identified in the three groups, respectively. GO and KEGG analyses showed that: in 2- vs. 6-month-olds, the target genes of DE known sheep miRNAs were involved in 100 biological processes and 11 signaling pathways; in 6- vs. 12-month-olds, the target genes of DE known sheep miRNAs were involved in 4 biological processes; and in 2- vs. 12-month-olds, the target genes of DE known sheep miRNAs were involved in 17 biological processes and 4 signaling pathways. Three miR–target regulatory networks were constructed based on these DE miRNA–targets. The key miRNA-Targets involved in testis development and spermatogenesis were screened. 6 known sheep miRNAs and 6 novel miRNAs were selected to validate the accuracy of miRNA sequencing data by qRT-PCR. The binding sites of oar-miR-379-5p with WNT8A was validated by a dual luciferase reporter gene detection system.

E2F7, EREG, miR-451a and miR-106b-5p are associated with the cervical cancer development.

We aimed to seek the crucial genes or microRNAs (miRNA) correlated with the cervical cancer development. The miRNA profiling GSE30656 and gene expression profiling GSE63514 were obtained from Gene Expression Omnibus database. Differentially expressed microRNAs (DEMiRs) and differentially expressed genes (DEGs) were screened. Then target genes of DEMiRs were obtained and matched with DEGs to obtain interaction pairs between DEMiRs and DEGs. Gene Ontology-biological process and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted for DEGs and DEMiRs in the DEMiRs-DEGs pairs. The DEMiRs-DEGs regulatory network, protein-protein interaction network and transcription factor (TF)-target regulatory network were constructed. Ultimately, long non-coding RNAs (lncRNAs) associated with DEMiRs were obtained, and then lncRNA-miRNA-target ceRNA network was established. Total 18 DEMiRs and 620 DEGs were identified. DEMiRs were enriched in 35 KEGG pathways, such as PI3K-Akt signaling pathway (involving miR-451a). DEGs were enriched in various functions, such as DNA replication (involving E2F7) and angiogenesis (involving EREG). There were 120 nodes and 216 interaction pairs in the DEMIR-DEG regulatory network, and miR-106b-5p has the greatest degree. EREG and E2F7 were regulated by miR-451a and miR-148a-3p, respectively. Besides, E2F7 was identified in the TF-target regulatory network, regulating CDC6. There were 15 lncRNAs, 11 miRNAs and 90 DEGs in the ceRNA network. Specially, miR-148a-3p was interacted with lncRNA HOTAIR in the ceRNA network. E2F7, EREG, miR-451a and miR-106b-5p were likely to be related to the cervical cancer development.

Identification of a microRNA signature associated with survivability in cervical squamous cell carcinoma

BackgroundThe aim of this study is to find the potential miRNA expression signature capable of predicting survival time for cervical squamous cell carcinoma (CSCC) patients.MethodsThe expression of 332 miRNAs was measured in 131 (Training cohort) and 130 (Validation cohort) patients with CSCC in the Cancer Genome Atlas (TCGA) data portal. The miRNA expression signature was identified by Cox Proportion Hazard regression model to the Training data set, and subsequently validated in an independent Validation set. Kaplan-Meier curves and the receiver operating characteristic analyses of 5 years were used to access the overall survival of miRNA signature. MiRNA signature-gene target analysis was performed, followed by the construction of the regulatory network. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis were used to explore the function of target genes of miRNA signature.ResultsA 2-miRNA expression signature of hsa-mir-642a and hsa-mir-378c associated with survivability was identified in CSCC. Both of them had a significant diagnostic and prognostic value of patients with CSCC. A total of 345 miRNA signature-target pairs were obtained in the miRNA signature-gene target regulatory network, in which 316 genes were targets of has-mir-378c and has-mir-642a. Functional analysis of target genes showed that MAPK signaling pathway, VEGF signaling pathway and endocytosis were the significantly enriched signal pathways that covered most genes.ConclusionsThe 2-miRNA signature adds to the prognostic value of CSCC. In-depth interrogation of the 2-miRNAs will provide important biological insights that finding and developing novel molecularly prediction to improve prognosis for CSCC patients.

Integrated analysis of gene expression signatures associated with colon cancer from three datasets

PurposeThe present study aimed to elucidate the pathogenesis of colon cancer and identify genes associated with tumor development. MethodsThree datasets, two (GSE74602 and GSE44861) from the Gene Expression Omnibus database and RNA-Seq colon cancer data from The Cancer Genome Atlas data portal, were downloaded. These three datasets were grouped using a meta-analysis approach, and differentially expressed genes (DEGs) were identified between colon tumor samples and adjacent normal samples. Functional enrichment analysis and regulatory factor predication were performed for significant genes. Additionally, small-molecule drugs associated with colon cancer were predicted, and a prognostic risk model was constructed. ResultsThere were 251 overlapping DEGs (135 up- and 116 downregulated) between cancer samples and control samples in the three datasets. The DEGs were mainly involved in protein transport and apoptotic and neurotrophin signaling pathways. A total of 70 small-molecule drugs were predicated to be associated with colon cancer. Additionally, in the miRNA–target regulatory network, we found that SLC44A1 can be targeted by hsa-miR-183, hsa-miR-206, and hsa-miR-147, while KLF13 can be regulated by hsa-miR-182, hsa-miR-206, and hsa-miR-153. Moreover, the results of the prognostic risk model showed that four genes (VAMP1, P2RX5, CACNB1, and CRY2) could divide the samples into high and low risk groups. ConclusionSLC44A1 and KLF13 may be involved in tumorigenesis and the metastasis of colon cancer by miRNA regulation. In addition, a four-gene (VAMP1, P2RX5, CACNB1, and CRY2) expression signature may have prognostic and predictive value in colon cancer.

Microarray Analysis of the Molecular Mechanism Involved in Parkinson's Disease.

Purpose This study aimed to investigate the underlying molecular mechanisms of Parkinson's disease (PD) by bioinformatics. Methods Using the microarray dataset GSE72267 from the Gene Expression Omnibus database, which included 40 blood samples from PD patients and 19 matched controls, differentially expressed genes (DEGs) were identified after data preprocessing, followed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Protein-protein interaction (PPI) network, microRNA- (miRNA-) target regulatory network, and transcription factor- (TF-) target regulatory networks were constructed. Results Of 819 DEGs obtained, 359 were upregulated and 460 were downregulated. Two GO terms, “rRNA processing” and “cytoplasm,” and two KEGG pathways, “metabolic pathways” and “TNF signaling pathway,” played roles in PD development. Intercellular adhesion molecule 1 (ICAM1) was the hub node in the PPI network; hsa-miR-7-5p, hsa-miR-433-3p, and hsa-miR-133b participated in PD pathogenesis. Six TFs, including zinc finger and BTB domain-containing 7A, ovo-like transcriptional repressor 1, GATA-binding protein 3, transcription factor dp-1, SMAD family member 1, and quiescin sulfhydryl oxidase 1, were related to PD. Conclusions “rRNA processing,” “cytoplasm,” “metabolic pathways,” and “TNF signaling pathway” were key pathways involved in PD. ICAM1, hsa-miR-7-5p, hsa-miR-433-3p, hsa-miR-133b, and the abovementioned six TFs might play important roles in PD development.

MicroRNA and target gene responses to salt stress in grafted cucumber seedlings

Grafting is widely applied to enhance cucumber growth and tolerance to various biotic and abiotic adversities. Except for regulating plant growth and plant development, microRNAs (miRNAs) respond to various adversities through negative mediation of target gene expression at the post-transcriptional level. In the present study, cucumber (Cucumis sativus) scions were grafted onto pumpkin (Cucurbita moschata) rootstocks to characterize the molecular mechanisms of miRNAs and target regulatory networks that confer increased salt stress tolerance. The expression levels of 17 miRNAs and their predicted target genes were determined by quantitative real-time PCR using auto-grafted (self-grafted cucumber seedlings) and hetero-grafted (cucumber scion grafted onto pumpkin rootstock) cucumber seedlings after 24 h of 100 mmol l−1 NaCl treatment. Compared with that of auto-grafted cucumber seedlings, the expression level of most miRNAs in the leaves of hetero-grafted cucumber seedlings increased under normal conditions, decreased considerably after 2 h of salt stress, and increased after 6 and 24 h of salt stress. The expression level of most miRNAs in the roots of hetero-grafted cucumber seedlings increased under normal conditions and salt-stressed conditions. Compared with that of auto-grafted cucumber seedlings, the expression level of most target genes in the leaves of hetero-grafted cucumber seedlings decreased under normal and salt-stressed conditions. However, the expression level of most target genes in the roots of hetero-grafted cucumber seedlings increased significantly. These results provide useful information relevant to the functional analysis of miRNAs mediating graft-dependent salt tolerance.