Functional Regulatory Network Research Articles

Establishing gene regulatory networks from Parkinson's disease risk loci.

The latest meta-analysis of genome-wide association studies identified 90 independent variants across 78 genomic regions associated with Parkinson’s disease, yet the mechanisms by which these variants influence the development of the disease remains largely elusive.To establish the functional gene regulatory networks associated with Parkinson’s disease risk variants, we utilized an approach combining spatial (chromosomal conformation capture) and functional (expression quantitative trait loci) data.We identified 518 genes subject to regulation by 76 Parkinson’s variants across 49 tissues, whicih encompass 36 peripheral and 13 CNS tissues. Notably, one-third of these genes were regulated via trans-acting mechanisms (distal; risk locus-gene separated by >1 Mb, or on different chromosomes). Of particular interest is the identification of a novel trans-expression quantitative trait loci–gene connection between rs10847864 and SYNJ1 in the adult brain cortex, highlighting a convergence between familial studies and Parkinson’s disease genome-wide association studies loci for SYNJ1 (PARK20) for the first time. Furthermore, we identified 16 neurodevelopment-specific expression quantitative trait loci–gene regulatory connections within the foetal cortex, consistent with hypotheses suggesting a neurodevelopmental involvement in the pathogenesis of Parkinson’s disease. Through utilizing Louvain clustering we extracted nine significant and highly intraconnected clusters within the entire gene regulatory network. The nine clusters are enriched for specific biological processes and pathways, some of which have not previously been associated with Parkinson’s disease.Together, our results not only contribute to an overall understanding of the mechanisms and impact of specific combinations of Parkinson’s disease variants, but also highlight the potential impact gene regulatory networks may have when elucidating aetiological subtypes of Parkinson’s disease.

Comparative Physiological and Transcriptomic Analyses Reveal Altered Fe-Deficiency Responses in Tomato Epimutant Colorless Non-ripening.

The mechanisms associated with the regulation of iron (Fe) homeostasis have been extensively examined, however, epigenetic regulation of these processes remains largely unknown. Here, we report that a naturally occurring epigenetic mutant, Colorless non-ripening (Cnr), displayed increased Fe-deficiency responses compared to its wild-type Ailsa Craig (AC). RNA-sequencing revealed that a total of 947 and 1,432 genes were up-regulated by Fe deficiency in AC and Cnr roots, respectively, while 923 and 1,432 genes were, respectively, down-regulated. Gene ontology analysis of differentially expressed genes showed that genes encoding enzymes, transporters, and transcription factors were preferentially affected by Fe deficiency. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed differential metabolic responses to Fe deficiency between AC and Cnr. Based on comparative transcriptomic analyses, 24 genes were identified as potential targets of Cnr epimutation, and many of them were found to be implicated in Fe homeostasis. By developing CRISPR/Cas9 genome editing SlSPL-CNR knockout (KO) lines, we found that some Cnr-mediated Fe-deficiency responsive genes showed similar expression patterns between SlSPL-CNR KO plants and the Cnr epimutant. Moreover, both two KO lines displayed Fe-deficiency-induced chlorosis more severe than AC plants. Additionally, the Cnr mutant displayed hypermethylation in the 286-bp epi-mutated region on the SlSPL-CNR promoter, which contributes to repressed expression of SlSPL-CNR when compared with AC plants. However, Fe-deficiency induced no change in DNA methylation both at the 286-bp epi-allele region and the entire region of SlSPL-CNR gene. Taken together, using RNA-sequencing and genetic approaches, we identified Fe-deficiency responsive genes in tomato roots, and demonstrated that SlSPL-CNR is a novel regulator of Fe-deficiency responses in tomato, thereby, paving the way for further functional characterization and regulatory network dissection.

Manifold learning analysis suggests strategies to align single-cell multimodal data of neuronal electrophysiology and transcriptomics

Recent single-cell multimodal data reveal multi-scale characteristics of single cells, such as transcriptomics, morphology, and electrophysiology. However, integrating and analyzing such multimodal data to deeper understand functional genomics and gene regulation in various cellular characteristics remains elusive. To address this, we applied and benchmarked multiple machine learning methods to align gene expression and electrophysiological data of single neuronal cells in the mouse brain from the Brain Initiative. We found that nonlinear manifold learning outperforms other methods. After manifold alignment, the cells form clusters highly corresponding to transcriptomic and morphological cell types, suggesting a strong nonlinear relationship between gene expression and electrophysiology at the cell-type level. Also, the electrophysiological features are highly predictable by gene expression on the latent space from manifold alignment. The aligned cells further show continuous changes of electrophysiological features, implying cross-cluster gene expression transitions. Functional enrichment and gene regulatory network analyses for those cell clusters revealed potential genome functions and molecular mechanisms from gene expression to neuronal electrophysiology.

Transcriptome Analysis of Insulin Signaling-Associated Transcription Factors in C. elegans Reveal Their Genome-Wide Target Genes Specificity and Complexity.

Insulin/IGF-1-like signaling (IIS) plays a crucial, conserved role in development, growth, reproduction, stress tolerance, and longevity. In Caenorhabditis elegans, the enhanced longevity under reduced insulin signaling (rIIS) is primarily regulated by the transcription factors (TFs) DAF-16/FOXO, SKN-1/Nrf-1, and HSF1/HSF-1. The specific and coordinated regulation of gene expression by these TFs under rIIS has not been comprehensively elucidated. Here, using RNA-sequencing analysis, we report a systematic study of the complexity of TF-dependent target gene interactions during rIIS under analogous genetic and experimental conditions. We found that DAF-16 regulates only a fraction of the C. elegans transcriptome but controls a large set of genes under rIIS; SKN-1 and HSF-1 show the opposite trend. Both of the latter TFs function as activators and repressors to a similar extent, while DAF-16 is predominantly an activator. For expression of the genes commonly regulated by TFs under rIIS conditions, DAF-16 is the principal determining factor, dominating over the other two TFs, irrespective of whether they activate or repress these genes. The functional annotations and regulatory networks presented in this study provide novel insights into the complexity of the gene regulatory networks downstream of the IIS pathway that controls diverse phenotypes, including longevity.

Comprehensive analysis of the cancer driver genes in breast cancer demonstrates their roles in cancer prognosis and tumor microenvironment

BackgroundBreast cancer is the most common malignancy in women. Cancer driver gene-mediated alterations in the tumor microenvironment are critical factors affecting the biological behavior of breast cancer. The purpose of this study was to identify the expression characteristics and prognostic value of cancer driver genes in breast cancer.MethodsThe Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets are used as the training and test sets. Classified according to cancer and paracancerous tissues, we identified differentially expressed cancer driver genes. We further screened prognosis-associated genes, and candidate genes were submitted for the construction of a risk signature. Functional enrichment analysis and transcriptional regulatory networks were performed to search for possible mechanisms by which cancer driver genes affect breast cancer prognosis.ResultsWe identified more than 200 differentially expressed driver genes and 27 prognosis-related genes. High-risk group patients had a lower survival rate compared to the low-risk group (P<0.05), and risk signature showed high specificity and sensitivity in predicting the patient prognosis (AUC 0.790). Multivariate regression analysis suggested that risk scores can independently predict patient prognosis. Further, we found differences in PD-1 expression, immune score, and stromal score among different risk groups.ConclusionOur study confirms the critical prognosis role of cancer driver genes in breast cancer. The cancer driver gene risk signature may provide a novel biomarker for clinical treatment strategy and survival prediction of breast cancer.

Comparative Label-Free Quantitative Proteomics Analysis Reveals the Essential Roles of N-Glycans in Salt Tolerance by Modulating Protein Abundance in Arabidopsis.

Many pieces of evidence show that the adaptive response of plants to salt stress requires the maturation of N-glycan on associated proteins. However, it is still little known about the salt-responsive glycoproteins that function in this process. In the present study, we identified salt-responsive glycoproteins in wild-type (WT) Arabidopsis and two mutants defective in N-glycan maturation, mns1 mns2 and cgl1. A total of 97 proteins with abundance changes of >1.5‐ or <0.67-fold were identified against salt stress by label-free liquid chromatography coupled mass spectrometry (LC-MS/MS) quantitative analyses. A comparison of differentially abundant glycoproteins (DAGs) indicated the substrate preferences regulated by MNS1/MNS2 and CGL1. In addition, the DAGs in mns1 mns2 hardly form functional regulatory networks in STRING analysis. Comparably, the regulatory network in cgl1 was visible and shared overlapping with that in WT. Such difference may supply the evidence to partially explain the lower salt sensitivity of mutant cgl1 than mns1 mns2. We further confirmed that two N-glycosylation clients, peroxidases PRX32 and PRX34, were involved in the salt stress response since the double mutants showed enhanced salt sensitivity. Together, our study provided proteomic evidence that N-glycans are crucial for modulating stress-responsive protein levels, and several novel glycoproteins responsible for salt stress tolerance in Arabidopsis were listed. Data are available via ProteomeXchange with identifier PXD006893.

Long Noncoding RNA RP11-357H14.17 Plays an Oncogene Role in Gastric Cancer by Activating ATF2 Signaling and Enhancing Treg Cells.

Background Gastric cancer (GC) is one of the most common malignant tumors in the world. The potential functions and mechanisms of long noncoding RNAs (lncRNAs) in GC development are still unclear. It is of great significance to explore the prognostic value of LncRNA signatures for GC. Methods LncRNAs differently expressed in GC and their prognostic value were studied based on The Cancer Genome Atlas (TCGA) database. The functional regulatory network and immune infiltration of RP11-357H14.17 were further studied using a variety of bioinformatics tools and databases. Results We found that the high expression of RP11-357H14.17 was closely associated with shortened overall survival (OS) and poor prognosis in gastric cancer patients. We also found that its expression was related to clinical features including tumor volume, metastasis, and differentiation. Functional enrichment analysis revealed that RP11-357H14.17 is closely related to enhanced DNA replication and metabolism; ssGSEA analysis implied the oncogenic roles of RP11-357H14.17 was related to ATF2 signaling and Treg cell differentiation. Furthermore, we verified such link by using real-time PCR and IHC staining in human GC samples. Conclusion We demonstrate that RP11-357H14.17 may play a crucial role in the occurrence, development, and malignant biological behavior of gastric cancer as a potential prognostic marker for gastric cancer.

Identification of Primary and Metastatic Lung Cancer-Related lncRNAs and Potential Targeted Drugs Based on ceRNA Network.

Lung cancer metastasis is the leading cause of poor prognosis and death for patients. Long noncoding RNAs (lncRNAs) have been validated the close correlation with lung cancer metastasis, but few comprehensive analyses have reported the specific association between lncRNA and cancer metastasis, especially via both competing endogenous RNA (ceRNA) regulatory relationships and functional regulatory networks. Here, we constructed primary and metastatic ceRNA networks, identified 12 and 3 candidate lncRNAs for lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) respectively and excavated some drugs that might have potential therapeutic effects on lung cancer progression. In summary, this study systematically analyzed the competitive relationships and regulatory mechanism of the repeatedly dysregulated lncRNAs in lung cancer carcinogenesis and metastasis, and provided a new idea for screening potential therapeutic drugs for lung cancer.

Expression and gene regulatory network of SNHG1 in hepatocellular carcinoma

BackgroundSmall nucleolar RNA host gene 1 (SNHG1), a long noncoding RNA (lncRNA), is a transcript that negatively regulates tumour suppressor genes, such as p53. Abnormal SNHG1 expression is associated with cell proliferation and cancer. We used sequencing data downloaded from Genomic Data Commons to analyse the expression and interaction networks of SNHG1 in hepatocellular carcinoma (HCC).MethodsExpression was examined using the limma package of R and verified by Gene Expression Profiling Interactive Analysis. We also obtained miRNA expression data from StarBase to determine the lncRNA-miRNA-mRNA–related RNA regulatory network in HCC. Kaplan–Meier (KM) analysis was performed using the survival package of R. Gene Ontology annotation of genes was carried out using Metascape.ResultsWe found that SNHG1 was overexpressed and often amplified in HCC patients. In addition, SNHG1 upregulation was associated with the promotion of several primary biological functions, including cell proliferation, transcription and protein binding. Moreover, we found similar trends of small nucleolar RNA host gene 1 (SNHG1), E2F8 (E2F transcription factor 8), FANCE (FA complementation group E) and LMNB2 (encodes lamin B2) expression. In the SNHG1-associated network, high expression levels of SNHG1 (log-rank P value = 0.0643), E2F8 (log-rank P value = 0.000048), FANCE (log-rank P value = 0.00125) and LMNB2 (log-rank P value = 0.0392) were significantly associated with poor survival. Single-cell analysis showed that E2F8 may play an important role in tumorigenesis or cancer development.ConclusionsOur results highlight the benefit of utilizing multiple datasets to understand the functional potential regulatory networks of SNHG1 and the role of SNHG1 in tumours.

Targeting Nestin+ hepatic stellate cells ameliorates liver fibrosis by facilitating TβRI degradation

Liver fibrosis is a wound healing response that arises from various aetiologies. The intermediate filament protein Nestin has been reported to participate in maintaining tissue homeostasis during wound healing responses. However, little is known about the role Nestin plays in liver fibrosis. This study investigated the function and precise regulatory network of Nestin during liver fibrosis. Nestin expression was assessed via immunostaining and quantitative real-time PCR (qPCR) in fibrotic/cirrhotic samples. The induction of Nestin expression by transforming growth factor beta (TGFβ)-Smad2/3 signalling was investigated through luciferase reporter assays. The functional role of Nestin in hepatic stellate cells (HSCs) was investigated by examining the pathway activity of profibrogenic TGFβ-Smad2/3 signalling and degradation of TGFβ receptor I (TβRI) after interfering with Nestin. The invivo effects of knocking down Nestin were examined with an adeno-associated virus vector (serotype 6, AAV6) carrying short-hairpin RNA targeting Nestin in fibrotic mouse models. Nestin was mainly expressed in activated HSCs and increased with the progression of liver fibrosis. The profibrogenic pathway TGFβ-Smad2/3 induced Nestin expression directly. Knocking down Nestin promoted caveolin 1-mediated TβRI degradation, resulting in TGFβ-Smad2/3 pathway impairment and reduced fibrosis marker expression in HSCs. In AAV6-treated murine fibrotic models, knocking down Nestin resulted in decreased levels of inflammatory infiltration, hepatocellular damage, and a reduced degree of fibrosis. The expression of Nestin in HSCs was induced by TGFβ and positively correlated with the degree of liver fibrosis. Knockdown of Nestin decreased activation of the TGFβ pathway and alleviated liver fibrosis both invitro and invivo. Our data demonstrate a novel role of Nestin in controlling HSC activation in liver fibrosis. Liver fibrosis has various aetiologies but represents a common process in chronic liver diseases that is associated with high morbidity and mortality. Herein, we demonstrate that the intermediate filament protein Nestin plays an essential profibrogenic role in liver fibrosis by forming a positive feedback loop with the TGFβ-Smad2/3 pathway, providing a potential therapeutic target for the treatment of liver fibrosis.

Genome wide in-silico miRNA and target network prediction from stress responsive Horsegram (Macrotyloma uniflorum) accessions

Horsegram (Macrotyloma uniflorum (Lam.) Verdc.) is a drought hardy food and fodder legume of Indo-African continents with diverse germplasm sources demonstrating alternating mechanisms depicting contrasting adaptations to different climatic zones. Tissue specific expression of genes contributes substantially to location specific adaptations. Regulatory networks of such adaptive genes are elucidated for downstream translational research. MicroRNAs are small endogenous regulatory RNAs which alters the gene expression profiles at a particular time and type of tissue. Identification of such small regulatory RNAs in low moisture stress hardy crops can help in cross species transfer and validation confirming stress tolerance ability. This study outlined prediction of conserved miRNAs from transcriptome shotgun assembled sequences and EST sequences of horsegram. We could validate eight out of 15 of the identified miRNAs to demonstrate their role in deficit moisture stress tolerance mechanism of horsegram variety Paiyur1 with their target networks. The putative mumiRs were related to other food legumes indicating the presence of gene regulatory networks. Differential miRNA expression among drought specific tissues indicted the probable energy conservation mechanism. Targets were identified for functional characterization and regulatory network was constructed to find out the probable pathways of post-transcriptional regulation. The functional network revealed mechanism of biotic and abiotic stress tolerance, energy conservation and photoperiod responsiveness.

Crosstalk Between MYC and lncRNAs in Hematological Malignancies.

The human genome project revealed the existence of many thousands of long non-coding RNAs (lncRNAs). These transcripts that are over 200 nucleotides long were soon recognized for their importance in regulating gene expression. However, their poor conservation among species and their still controversial annotation has limited their study to some extent. Moreover, a generally lower expression of lncRNAs as compared to protein coding genes and their enigmatic biochemical mechanisms have impeded progress in the understanding of their biological roles. It is, however, known that lncRNAs engage in various kinds of interactions and can form complexes with other RNAs, with genomic DNA or proteins rendering their functional regulatory network quite complex. It has emerged from recent studies that lncRNAs exert important roles in gene expression that affect many cellular processes underlying development, cellular differentiation, but also the pathogenesis of blood cancers like leukemia and lymphoma. A number of lncRNAs have been found to be regulated by several well-known transcription factors including Myelocytomatosis viral oncogene homolog (MYC). The c-MYC gene is known to be one of the most frequently deregulated oncogenes and a driver for many human cancers. The c-MYC gene is very frequently activated by chromosomal translocations in hematopoietic cancers most prominently in B- or T-cell lymphoma or leukemia and much is already known about its role as a DNA binding transcriptional regulator. Although the understanding of MYC’s regulatory role controlling lncRNA expression and how MYC itself is controlled by lncRNA in blood cancers is still at the beginning, an intriguing picture emerges indicating that c-MYC may execute part of its oncogenic function through lncRNAs. Several studies have identified lncRNAs regulating c-MYC expression and c-MYC regulated lncRNAs in different blood cancers and have unveiled new mechanisms how these RNA molecules act. In this review, we give an overview of lncRNAs that have been recognized as critical in the context of activated c-MYC in leukemia and lymphoma, describe their mechanism of action and their effect on transcriptional reprogramming in cancer cells. Finally, we discuss possible ways how an interference with their molecular function could be exploited for new cancer therapies.

An integrated approach to unravel a putative crosstalk network in Alzheimer’s disease and Parkinson’s disease

Integration of multiple profiling data and construction of functional regulatory networks provide a powerful approach to uncover functional relationships and significant molecular entities from transcriptomic data, highlighting the molecular mechanisms of complex diseases. Despite having an overlap in the neuropathologies of AD and PD, the molecular entities overlapped and mechanisms behind them are less known. Here we used an integrated strategy to analyze miRNA and gene transcriptomic data to understand the role of miRNAs and genes in regulatory activities taking place in cells, and find transcriptomic signatures linking AD and PD. We preprocessed and analyzed publicly available microarray datasets and identified 97 DEGs and 21 DEmiRs that may be involved in the overlapped mechanisms between these two disorders. Among the DEGs, we found HSPA9, PGK1, SDHC, FH, DLD, YWHAZ and ACLY as the major protein-coding genes involved in the crosstalk for AD-PD pathogenesis. Further we integrated these DEGs and DEmiRs with regulatory TFs to construct an overlapped dysregulated network of AD and PD. In the network, miR-27a-3p, miR-148a-3p and miR-15a-5p were found to be the most relevant with maximum interactions, describing their significance in the potential crosstalk. We also looked into the dysregulated biological processes and pathways overlapped in AD and PD. In conclusion, we highlighted the DEGs, DEmiRs, their interactions and related pathways overlapped in AD and PD pathogenesis, also describing a potential crosstalk at molecular level. Besides, our findings can further be used for molecular studies to reveal an assured AD-PD crosstalk.

Integrated Bioinformatics Analysis Reveals Function and Regulatory Network of miR-200b-3p in Endometriosis.

Objective MicroRNAs play vital roles in the development of endometriosis. It is reported that miR-200b-3p is downregulated in endometriosis, although its mechanisms in this disease remain still unclear. Therefore, the purpose of this study was to explore the function and potential regulatory network of miR-200b-3p in endometriosis through database analysis. Methods The endometriosis gene expression profiles were downloaded from the GEO database to screen differentially expressed genes (DEGs). The predicted and validated target genes of miR-200b-3p were obtained from miRWalk and miRTarBase database. Then, a comparison was performed between miR-200b-3p target genes and DEGs. GO enrichment and KEGG pathway analysis of the target genes was performed using clusterProfiler package. STRING was used to predict the protein-protein interaction among the proteins encoded by the target genes. Then, TransmiR, LncBase, StarBase, PROMO, and AnimalTFDB were employed to identify interactive transcription factors and lncRNAs of miR-200b-3p. Results miR-200b-3p was associated with the transcription factors DNMT1, EZH2, HNF1B, JUN, MYB, ZEB1, and ZEB2 during the pathogenesis of endometriosis. The downstream 110 target genes were involved in the biological processes of positive regulation of MAPK cascade, muscle cell proliferation, organ growth, vasculogenesis, and axon development. KEGG analysis revealed that the main pathways related to miR-200b-3p were microRNAs in cancer, PI3K-Akt signaling pathway, colorectal cancer, and tight junction. In addition, four lncRNAs such as MALAT1, NEAT1, SNHG22, and XIST interacted with miR-200b-3p and were associated with transcription factors FOXP3 and YY1. Conclusion The predicted target genes and molecular regulatory network of miR-200b-3p in endometriosis not only revealed its biological function but also provided a valuable guideline for further research.

Ddx56 maintains proliferation of mouse embryonic stem cells via ribosome assembly and interaction with the Oct4/Sox2 complex

BackgroundEmbryonic stem cells (ESCs) are important source of clinical stem cells for therapy, so dissecting the functional gene regulatory network involved in their self-renewal and proliferation is an urgent task. We previously reported that Ddx56 interacts with the core transcriptional factor Oct4 by mass spectrometry analysis in ESCs. However, the exact function of Ddx56 in ESCs remains unclear.MethodsWe investigated the role of Ddx56 in mouse ESCs (mESCs) through both gain- and loss-of-function strategies. The effect of Ddx56 on mESCs was determined based on morphological changes, involvement in the network of pluripotency markers (Nanog, Oct4, Sox2), and altered lineage marker expression. In addition, the role of Ddx56 in mESCs was evaluated by polysome fractionation, qRT-PCR, and co-immunoprecipitation (co-IP). Finally, RNA sequencing was applied to explore potential network regulation by Ddx56 in mESCs.ResultWe found that Ddx56 participated in ribosome assembly, as knockout or RNAi knockdown of Ddx56 led to ribosome dysfunction and cell lethality. Surprisingly, exogenous expression of C-terminal domain truncated Ddx56 (Ddx56 ΔC-ter) did not affect ribosome assembly, but decreased mESC proliferation by downregulation of proliferation-related genes and cell cycle changing. In terms of mechanism, Ddx56 interacted with the Oct4 and Sox2 complex by binding to Sox2, whereas Ddx56 ΔC-ter showed weaker interaction with Sox2 and led to retardation of mESC proliferation.ConclusionsDdx56 maintains ESC proliferation by conventional regulation of ribosome assembly and interaction with the Oct4 and Sox2 complex.

Identification of long noncoding RNA RP11-89K21.1 and RP11-357H14.17 as prognostic signature of endometrial carcinoma via integrated bioinformatics analysis

BackgroundEndometrial carcinoma (EC) is one of the most common malignant tumors in gynecology. The potential functions and mechanisms of long noncoding RNAs (lncRNAs) in the occurrence and progression of EC remains unclear. It’s meaningful to explore lncRNAs signature for providing prognostic value of EC.MethodsThe differentially expressed lncRNAs and their prognostic values in EC were investigated based on The Cancer Genome Atlas (TCGA) database; the transcriptional factors (TFs), the competing endogenous RNA (ceRNA) mechanism, functional regulatory network and immune infiltration of RP11-89K21.1 and RP11-357H14.17 were further explored by various bioinformatics tools and databases.ResultsWe firstly identified high expression of RP11-89K21.1 and RP11-357H14.17 were closely associated with shorten overall survival (OS) and poor prognosis in patients with EC. We also elucidated the networks of transcription factor and co-expression genes associated with RP11-89K21.1 and RP11-357H14.17. Furthermore, the ceRNA network mechanism was successfully constructed through 2 lncRNAs (RP11-89K21.1 and RP11-357H14.17), 11 miRNAs and 183 mRNAs. Functional enrichment analysis revealed that the targeting genes of RP11-89K21.1 and RP11-357H14.17 were strongly associated with microRNAs in cancer, vessel development, growth regulation, growth factor and cell differentiation, and involved in pathways including pathways in cancer, microRNAs in cancer and apoptotic signaling pathway.ConclusionsWe demonstrated for the first time that RP11-89K21.1 and RP11-357H14.17 may play crucial roles in the occurrence, development and malignant biological behavior of EC, and can be regarded as potential prognostic biomarkers for EC.

Multi-level transcriptome sequencing identifies COL1A1 as a candidate marker in human heart failure progression

BackgroundHeart failure (HF) has been recognized as a global pandemic with a high rate of hospitalization, morbidity, and mortality. Although numerous advances have been made, its representative molecular signatures remain largely unknown, especially the role of genes in HF progression. The aim of the present prospective follow-up study was to reveal potential biomarkers associated with the progression of heart failure.MethodsWe generated multi-level transcriptomic data from a cohort of left ventricular heart tissue collected from 21 HF patients and 9 healthy donors. By using Masson staining to calculate the fibrosis percentage for each sample, we applied lasso regression model to identify the genes associated with fibrosis as well as progression. The genes were further validated by immunohistochemistry (IHC) staining in the same cohort and qRT-PCR using another independent cohort (20 HF and 9 healthy donors). Enzyme-linked immunosorbent assay (ELISA) was used to measure the plasma level in a validation cohort (139 HF patients) for predicting HF progression.ResultsBased on the multi-level transcriptomic data, we examined differentially expressed genes [mRNAs, microRNAs, and long non-coding RNAs (lncRNAs)] in the study cohort. The follow-up functional annotation and regulatory network analyses revealed their potential roles in regulating extracellular matrix. We further identified several genes that were associated with fibrosis. By using the survival time before transplantation, COL1A1 was identified as a potential biomarker for HF progression and its upregulation was confirmed by both IHC and qRT-PCR. Furthermore, COL1A1 content ≥ 256.5 ng/ml in plasma was found to be associated with poor survival within 1 year of heart transplantation from heart failure [hazard ratio (HR) 7.4, 95% confidence interval (CI) 3.5 to 15.8, Log-rank p value < 1.0 × 10− 4].ConclusionsOur results suggested that COL1A1 might be a plasma biomarker of HF and associated with HF progression, especially to predict the 1-year survival from HF onset to transplantation.

Dysregulated circulating SOCS3 and haptoglobin expression associated with stable coronary artery disease and acute coronary syndrome: An integrated study based on bioinformatics analysis and case-control validation.

Objective:To extensively use blood transcriptome analysis to identify potential diagnostic and therapeutic targets for cardiovascular diseases.Methods:Two gene expression datasets (GSE59867 and GSE62646) were downloaded from GEO DataSets to identify altered blood transcriptomes in patients with ST-segment elevation myocardial infarction (STEMI) compared to stable coronary artery disease (CAD). Thereafter, several computational approaches were taken to determine functional roles and regulatory networks of differentially expressed genes (DEGs). Finally, the expression of dysregulated two hub genes–suppressor of cytokine signaling 3 (SOCS3) and haptoglobin (HP)–were validated in a case–control study.Results:A total of 119 DEGs were identified in the discovery phase, consisting of 71 downregulated genes and 48 upregulated genes; two hub modules consisting of two hub genes–SOCS3 and HP–were identified. In the validation phase, both SOCS3 and HP were significantly downregulated in the stable CAD and acute coronary syndrome (ACS) patients when compared with healthy controls. Meanwhile, HP was significantly upregulated in STEMI patients when compared with stable CAD patients (p=0.041). Logistic regression analysis indicated that: downregulated expression of HP correlated with increased risk of CAD [odds ratio (OR)=0.52, 95% confidence interval (CI)=0.31~0.87, p=0.013]; and downregulated expression of SOCS3 correlated with increased risk of ACS (OR=0.66, 95% CI=0.46~0.94, p=0.023) when age, gender, history of hyperlipidemia, diabetes and hypertension were included as covariates.Conclusion:Future clarification of how SOCS3 and HP influence the pathogenesis of disease may pave the way for the development of novel diagnostic and therapeutic methods.

Rewired functional regulatory networks among miRNA isoforms (isomiRs) from let-7 and miR-10 gene families in cancer

Classical microRNA (miRNA) has been so far believed as a single sequence, but it indeed contains multiple miRNA isoforms (isomiR) with various sequences and expression patterns. It is not clear whether these diverse isomiRs have potential relationships and whether they contribute to miRNA:mRNA interactions. Here, we aimed to reveal the potential evolutionary and functional relationships of multiple isomiRs based on let-7 and miR-10 gene families that are prone to clustering together on chromosomes. Multiple isomiRs within gene families showed similar functions to their canonical miRNAs, indicating selection of the predominant sequence. IsomiRs containing novel seed regions showed increased/decreased biological function depending on whether they had more/less specific target mRNAs than their annotated seed. Few gene ontology(GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were shared among the target genes of the annotated seeds and the novel seeds. Various let-7 isomiRs with novel seed regions may cause opposing drug responses despite the fact that they are generated from the same miRNA locus and have highly similar sequences. IsomiRs, especially the dominant isomiRs with shifted seeds, may disturb the coding-non-coding RNA regulatory network. These findings provide insight into the multiple isomiRs and isomiR-mediated control of gene expression in the pathogenesis of cancer.

Exploration of the effects of a degS mutant on the growth of Vibrio cholerae and the global regulatory function of degS by RNA sequencing.

BackgroundDegS is a periplasmic serine protease that is considered to be the initiator of the σE stress response pathway, and this protein plays an important role in the regulation of the stress response in E. coli. However, knowledge of the biological function and global regulatory network of DegS in Vibrio cholerae remains limited. In this study, we aimed to characterize the molecular functions and further investigate the regulatory network of degS in V. cholerae.MethodsA deletion mutant of degS was constructed in the V. cholerae HN375 strain. Bacterial colony morphology was observed by a plate-based growth experiment, and bacterial growth ability was observed by a growth curve experiment. High-throughput RNA sequencing (RNA-Seq) technology was used to analyze the differential transcriptomic profiles between the wild-type and degS mutant strains. Gene ontology (GO), pathway analysis and Gene-Act-network analysis were performed to explore the main functions of the differentially expressed genes. Quantitative real-time PCR (qRT-PCR) was performed to validate the reliability and accuracy of the RNA-Seq analysis. The complementation experiments were used to test the roles of degS and ropS in the small colony degS mutant phenotype.ResultsWhen degS was deleted, the degS mutant exhibited smaller colonies on various media and slower growth than the wild-type strain. A total of 423 differentially expressed genes were identified, including 187 genes that were upregulated in the degS mutant compared to the wild-type strain and 236 genes that were relatively downregulated. GO categories and pathway analysis showed that many differentially expressed genes were associated with various cellular metabolic pathways and the cell cycle. Furthermore, Gene-Act network analysis showed that many differentially expressed genes were involved in cellular metabolic pathways and bacterial chemotaxis. The cAMP-CRP-RpoS signaling pathway and the LuxPQ signal transduction system were also affected by the degS mutant. The expression patterns of nine randomly selected differentially expressed genes were consistent between the qRT-PCR and RNA-seq results. The complementation experiments showed that the small colony degS mutant phenotype could be partially restored by complementation with the pBAD24-degS or pBAD24-rpoS plasmid.DiscussionThese results suggest that the degS gene is important for normal growth of V. cholerae. Some of the differentially expressed genes were involved in various cellular metabolic processes and the cell cycle, which may be associated with bacterial growth. Several new degS-related regulatory networks were identified. In addition, our results suggested that the cAMP-CRP-RpoS signaling pathway may be involved in the small colony degS mutant phenotype. Overall, we believe that these transcriptomic data will serve as useful genetic resources for research on the functions of degS in V. cholerae.