Identification of a novel circRNA, hsa_circ_0065898, that regulates tumor growth in cervical squamous cell carcinoma

Background

Circular RNAs (circRNAs) were reported to play an important role in regulating tumor pathogenesis. The molecular mechanism of circRNAs in cervical squamous cell carcinoma (CSCC) remains poorly understood. We aimed to identify the circRNAs differentially expressed, and to investigate the role of a novel circRNA, hsa_circ_0065898, in regulating proliferation, migration, and invasion in CSCC.

Methods

The online Kaplan-Meier Plotter was used to analyze the relationship between miRNA expression and overall survival. Bioinformatics tools, such as R, Cytoscape, and Perl, were used to analyze the Gene Ontology (GO) enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, protein-protein interaction (PPI) network, and regulatory network. The expression level of hsa_circ_0065898 in CSCC cell lines was evaluated using quantitative polymerase chain reaction in vitro. The cell counting kit-8 (CCK-8) and transwell assays were used to assess cell proliferation, migration, and invasion.

Results

circRNA expression data (GSE102686) was downloaded from the Gene Expression Omnibus database, and this included data from 5 CSCC patients and 5 normal tissues. 13 differentially expressed circRNAs were identified, which included 9 upregulated circRNAs and 4 downregulated circRNAs. GO enrichment analysis showed that the target genes of miRNAs associated with hsa_circ_0065898 were enriched in ubiquitin-protein transferase activity, ubiquitin-like protein transferase activity, core promoter sequence-specific DNA binding, mRNA 3’-UTR AU-rich region binding, core promoter binding, and so on. KEGG showed that the Hippo and p53 signaling pathways played significant role in the pathway network. Hsa_circ_0065898 was significantly overexpressed in the CSCC cell lines. Hsa_circ_0065898 facilitated cell proliferation, migration, and invasion in CSCC.

Conclusions

This study identified differentially expressed circRNAs and constructed the regulatory network of hsa_circ_0065898 targeting microRNAs and mRNAs. We demonstrated that hsa_circ_0065898 promoted CSCC cell proliferation, migration, and invasion. Hence, hsa_circ_0065898 might be useful as a biomarker for CSCC diagnosis and targeted therapy.

Bioinformatic analysis of a microRNA regulatory network in Huntington's disease

Huntington's disease is an autosomal dominant hereditary neurodegenerative disease characterized by progressive dystonia, chorea and cognitive or psychiatric disturbances. The leading cause is the Huntington gene mutation on the patient's chromosome 4 that produces a mutated protein. Recently, attention has focused on the relationship between microRNAs and Huntington's disease's pathogenesis. In Huntington's disease, microRNAs can interact with various transcription factors; dysregulated microRNAs may be associated with the Cytosine deoxynucleotide-Adenine ribonucleotides-Guanine ribonucleotide length and Huntington's disease's progression and severity. This study explores the role of microRNAs in the pathogenesis of Huntington's disease through bioinformatics analysis. By analyzing data from the Gene Expression Omnibus database, we identified a total of 9 differentially expressed microRNA. Subsequently, target genes and long non-coding RNAs were predicted, and a comprehensive regulatory network centered on microRNA was constructed. The microRNA integrated regulatory network, Homo sapiens (hsa)-miR-144-3p, interacted with the largest number of long non-coding RNAs, including X-inactive specific transcript and taurine upregulated gene 1. The miRNAs, hsa-miR-10b-5p and hsa-miR-196a-5p, regulated most of the target genes, including class I homeobox and brain-derived neurotrophic factor genes. Additionally, 59 Gene Ontology terms and eight enrichment pathways were identified by analyzing the target genes of hsa-miR-196a-5p and hsa-microRNA-10b-5p. In conclusion, hsa-miR-10b-5p and hsa-miR-196a-5p were significantly and differentially expressed in Huntington's disease, the long non-coding RNAs X-inactive specific transcript, taurine upregulated gene 1, and target genes such as homeobox or brain-derived neurotrophic factor may play critical roles in the pathogenesis of Huntington's disease.

Identification and Validation of Key miRNAs and a microRNA-mRNA Regulatory Network Associated with Ulcerative Colitis

Ulcerative colitis (UC) is a chronic, nonspecific, intestinal inflammatory disease that involves various genes and pathways in its pathogenesis. The current study revealed the key miRNAs and potential target gene regulatory network as a model for predicting the molecular mechanism of UC. This may provide novel insights for unraveling the pathogenesis of UC. Differentially expressed miRNAs (DEMIs) and mRNAs (differentially expressed genes [DEGs]) between UC patients and normal controls were screened using the Gene Expression Omnibus database. DEMI target genes were predicted using the miRDB, miRWalk, starBase, TarBase, and TargetScan databases, and an miRNA-mRNA network was established using DEGs that altered in opposition to DEMIs. We verified the expression of key DEMIs in a rodent UC model. The miRNA-mRNA network contained 31 DEMIs and 199 DEGs, which showed enrichment in inflammatory bowel disease. We selected 2 key miRNAs and 4 hub genes. In addition, we identified six DEMIs and genes from the preliminary validation analysis in model tissues. In the pathophysiological process of UC, various genes and proteins display expression differences and complex interactions with each other. These findings provide new insights into the potential key mechanisms associated with UC development.

Identification of Candidate Circular RNAs Underlying Intramuscular Fat Content in the Donkey

Intramuscular fat (IMF) content is a crucial indicator of meat quality. Circular RNAs (circRNAs) are a large class of endogenous RNAs that are involved in many physiological processes. However, the expression and function of circRNA in IMF in the donkey remains unresolved. Here we performed an expression profiling of circRNAs in the donkey longissimus dorsi muscle and identified 12,727 candidate circRNAs. Among these, 70% were derived from the exons of protein genes. Furthermore, a total of 127 differentially expressed (DE) circRNAs were identified in high (H) and low (L) IMF content groups, including 63 upregulated and 64 downregulated circRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of the host genes of the DE circRNAs showed that the host genes were enriched in lipid metabolism related GO terms (e.g., fatty acid beta-oxidation using acyl-CoA dehydrogenase and MLL3/4 complex), and signaling pathways (e.g., TGF-beta and lysine degradation signaling pathway). Further analyses indicated that 127 DE circRNAs were predicted to potentially interact with miRNAs, leading to the construction of circRNA-miRNA regulatory network. Multiple circRNAs can potentially function as sponges of miRNAs that regulate the differentiation of adipocytes. Our results provide valuable expression profile information for circRNA in the donkey and new insight into the regulatory mechanisms of circRNAs in the regulation of IMF content.

A Network-Based Analysis Reveals the Mechanism Underlying Vitamin D in Suppressing Cytokine Storm and Virus in SARS-CoV-2 Infection

Background

SARS-CoV-2 causes ongoing pandemic coronavirus disease of 2019 (COVID-19), infects the cells of the lower respiratory tract that leads to a cytokine storm in a significant number of patients resulting in severe pneumonia, shortness of breathing, respiratory and organ failure. Extensive studies suggested the role of Vitamin D in suppressing cytokine storm in COVID-19 and reducing viral infection; however, the precise molecular mechanism is not clearly known. In this work, bioinformatics and systems biology approaches were used to understand SARS-CoV-2 induced cytokine pathways and the potential mechanism of Vitamin D in suppressing cytokine storm and enhancing antiviral response.

Results

This study used transcriptome data and identified 108 differentially expressed host genes (DEHGs) in SARS-CoV-2 infected normal human bronchial epithelial (NHBE) cells compared to control. Then, the DEHGs was integrated with the human protein-protein interaction data to generate a SARS-CoV-2 induced host gene regulatory network (SiHgrn). Analysis of SiHgrn identified a sub-network "Cluster 1" with the highest MCODE score, 31 up-regulated genes, and predominantly associated immune and inflammatory response. Interestingly, the iRegulone tool identified that "Cluster 1" is under the regulation of transcription factors STAT1, STAT2, STAT3, POU2F2, and NFkB1, collectively referred to as "host response signature network". Functional enrichment analysis with NDEx revealed that the "host response signature network" is predominantly associated with critical pathways, including "cytokines and inflammatory response", "non-genomic action of Vitamin D", "the human immune response to tuberculosis", and "lung fibrosis". Finally, in-depth analysis and literature mining revealed that Vitamin D binds with its receptor and could work through two different pathways: (i) it inhibits the expression of pro-inflammatory cytokines through blocking the TNF induced NFkB1 signaling pathway; and (ii) it initiates the expression of interferon-stimulating genes (ISGs) for antiviral defense program through activating the IFN-α induced Jak-STAT signaling pathway.

Conclusion

This comprehensive study identified the pathways associated with cytokine storm in SARS-CoV-2 infection. The proposed underlying mechanism of Vitamin D could be promising in suppressing the cytokine storm and inducing a robust antiviral response in severe COVID-19 patients. The finding in this study urgently needs further experimental validations for the suitability of Vitamin D in combination with IFN-α to control severe COVID-19.

Networks of Physiological Adjustments and Defenses, and Their Synergy With Sodium (Na+) Homeostasis Explain the Hidden Variation for Salinity Tolerance Across the Cultivated Gossypium hirsutum Germplasm

The abilities to mobilize and/or sequester excess ions within and outside the plant cell are important components of salt-tolerance mechanisms. Mobilization and sequestration of Na+ involves three transport systems facilitated by the plasma membrane H+/Na+ antiporter (SOS1), vacuolar H+/Na+ antiporter (NHX1), and Na+/K+ transporter in vascular tissues (HKT1). Many of these mechanisms are conserved across the plant kingdom. While Gossypium hirsutum (upland cotton) is significantly more salt-tolerant relative to other crops, the critical factors contributing to the phenotypic variation hidden across the germplasm have not been fully unraveled. In this study, the spatio-temporal patterns of Na+ accumulation along with other physiological and biochemical interactions were investigated at different severities of salinity across a meaningful genetic diversity panel across cultivated upland Gossypium. The aim was to define the importance of holistic or integrated effects relative to the direct effects of Na+ homeostasis mechanisms mediated by GhHKT1, GhSOS1, and GhNHX1. Multi-dimensional physio-morphometric attributes were investigated in a systems-level context using univariate and multivariate statistics, randomForest, and path analysis. Results showed that mobilized or sequestered Na+ contributes significantly to the baseline tolerance mechanisms. However, the observed variance in overall tolerance potential across a meaningful diversity panel were more significantly attributed to antioxidant capacity, maintenance of stomatal conductance, chlorophyll content, and divalent cation (Mg2+) contents other than Ca2+ through a complex interaction with Na+ homeostasis. The multi-tier macro-physiological, biochemical and molecular data generated in this study, and the networks of interactions uncovered strongly suggest that a complex physiological and biochemical synergy beyond the first-line-of defense (Na+ sequestration and mobilization) accounts for the total phenotypic variance across the primary germplasm of Gossypium hirsutum. These findings are consistent with the recently proposed Omnigenic Theory for quantitative traits and should contribute to a modern look at phenotypic selection for salt tolerance in cotton breeding.

Identification of functional circRNA/miRNA/mRNA regulatory network for exploring prospective therapy strategy of colorectal cancer

Circular RNA (circRNA) has been reported to have great scientific significance and clinical value in multiple cancers including colorectal cancer (CRC). However, the biological function of most circRNAs in CRC is still in its infancy. Herein, we discovered the differential expressed circRNAs (DECs) between CRC tissues and matched adjacent using deep RNA sequencing and further confirmed the DECs expression by combining with another Gene Expression Omnibus dataset. Furthermore, we validated the expression of the top four upregulated circRNAs (hsa_circ_0030632, hsa_circ_0004887, hsa_circ_0001550, and hsa_circ_0001681) in both of paired CRC tissues and CRC cell lines. Then, a circRNA/microRNA/messenger RNA regulatory network was established and the Gene Ontology analysis and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis showed these four circRNAs participated in various biological processed including apoptotic process and multiple metabolic processes. Moreover, based on the regulatory network, three bioactive compounds (pergolide, pivampicillin, and methylergometrine) for the treatment of CRC were also found. In conclusion, this study improved our understanding of circRNAs and may also facilitate the finding of promising targets and biomarkers in CRC.

Integrated Analysis of Seed microRNA and mRNA Transcriptome Reveals Important Functional Genes and microRNA-Targets in the Process of Walnut (Juglans regia) Seed Oil Accumulation

Walnut (Juglans regia) is known as a promising woody oil crop with abundant polyunsaturated fatty acids in its kernel. However, the regulation mechanism of walnut oil accumulation and fatty acid metabolism is still poorly understood, which restricted the breeding and genetic improvement of high-quality oil-bearing walnuts. To reveal the molecular mechanism of walnut oil accumulation, considering the potential regulation of microRNA (miRNA) in seed development, in this study, the oil content of walnut kernel on the 80th, 100th and 120th day after flowering (DAF) was tested and the corresponding proportions are 11.51%, 40.40% and 53.20%. Between DAF of 80th~120th, the content of stearic acid and oleic acid tended to increase, but the proportion of other fatty acids tended to decrease. Meanwhile, comparative transcriptome and sRNA-seq analysis on three stages (80th, 100th and 120th DAF), found 204 conserved miRNAs and 554 novel miRNAs in walnut kernels, among which 104 key genes related to walnut oil accumulation were screened. The phospholipid:diacylglycerol acyltransferase metabolic pathway may contribute more to oil accumulation in walnut. 16 miRNA-mRNA regulatory modules related to walnut oil accumulation and fatty acid synthesis were constructed. 8 known miRNAs and 9 novel miRNAs regulate 28 genes involved in fatty acid (FA) metabolism and lipid synthesis. Among them, jre-miRn105, jre-miRn434, jre-miR477d and jre-miR156a.2 are key miRNAs that regulate walnut FA synthesis. Jre-miRn411 and jre-miR399a.1 are closely related to oil accumulation. These data provide new insights and lay the foundation for subsequent studies on walnut FA synthesis and oil accumulation.

Expression Profiles of Circular RNA in Human Atrial Fibrillation With Valvular Heart Diseases

Circular RNAs (circRNA) are involved in a variety of human heart diseases, however, circRNA expression profiles and circRNA-miRNA-mRNA regulatory network in human atrial fibrillation (AF) especially with valvular heart diseases (VHD) remain poorly understood. A high-throughput RNA sequencing was used to investigate the differentially expressed circRNAs in left atrial appendage from VHD patients with or without persistent AF. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed to predict the potential functions of the host genes of differentially expressed circRNA and their downstream targets. CircRNA-miRNA-mRNA regulatory network was constructed to identify mechanisms underlying circRNAs. qRT-PCR and sanger sequencing were further performed to validate the results. Compared with sinus rhythm (SR) patients, there were 3094 upregulated and 4472 downregulated circRNAs in AF patients respectively. The expression of 10 most differentially expressed circRNAs (circ 255-ITGA7, circ 418-KCNN2, circ 13913-MIB1, circ 44670-BARD1, circ 44782-LAMA2, circ 81906-RYR2, circ 35880-ANO5, circ 22249-TNNI3K, circ 3136-TNNI3K, circ 56186-TNNI3K) between SR and persistent AF patients were verified by qRT-PCR. In addition, specific back-splicing sites of these circRNAs was confirmed by sanger sequencing. GO and KEGG pathway analysis indicated that cAMP signal pathway and Wnt signal pathway might play important role in the development of AF in VHD patients, which might be affected by circRNAs. This study provided a preliminary landscape of circRNAs expression profiles which are involved in persistent AF due to VHD, and established the possibility for future related researches in this field.

Elucidation of the miR164c-Guided Gene/Protein Interaction Network Controlling Seed Vigor in Rice

MicroRNAs (miRNAs) play important roles in various aspects of plant physiology and metabolism. The expression level of miR164c is negatively correlated with seed vigor in rice (Oryza sativa L.); however, the mechanism of seed vigor regulation by miR164c remains unknown. Anti-aging capacity is an important indicator of seed vigor. Here, we report an miR164c-guided gene/protein interaction network that regulates the anti-aging ability of rice seeds. Seeds of the wild-type (WT) rice cultivar "Kasalath" and its transgenic derivatives, miR164c-silenced line (MIM164c) and miR164c overexpression line (OE164c), with significant differences in anti-aging capacity, showed significant differences in gene and protein expression levels. The differentially expressed genes (DEGs) or proteins were significantly enriched in six metabolic functional categories related to seed vigor, including "stress response," "protein processing in endoplasmic reticulum (ER)," "embryo development," "serine-type endopeptidase inhibitor," "energy metabolism," and "other." Differences in the expression levels of genes or proteins related to energy metabolism, serine endopeptidase, and stress response in seeds under normal storage conditions may be associated with anti-aging capacity. The results of gene/protein interaction analyses suggest that miR164c first targets PSK5, and the PSK5 protein then interacts with the ubiquitin-associated gene RPS27AA, which simultaneously impacts the genes/proteins in the six above-mentioned functional categories. Expression levels of some of the key genes and proteins in the interaction network were verified by real-time fluorescence quantitative PCR (RT-qPCR) and multiple reaction monitoring mass spectrometry (MRM-MS), respectively. Thus, the present study provides new insights into the miRNA-mediated gene and protein interaction network that regulates seed vigor.

Integrated Analysis of the Functions and Prognostic Values of RNA-Binding Proteins in Colorectal Cancer

Colorectal cancer (CRC) is one of the most common malignant tumors. Selecting effective treatment for CRC patients, especially in the early stages, remains a challenge because of the lack of adequate biomarkers. Recent evidence suggests that RNA-binding proteins (RBPs) play a vital role in development and progression of carcinogenesis. However, their mechanisms in cancer progression are still limited. The role of RBPs in CRC has been poorly understood. There were 1,542 reported RBPs analyzed between CRC tissues and normal tissues using the Wilcoxon test to identify differentially expressed RBPs (DE RBPs). Then, the potential functions and the prognostic value of these DE RBPs were explored through systematic bioinformatics analysis. There were 177 DE RBPs identified between CRC tissues and normal tissues. A protein–protein interaction network was constructed based on DE RBPs, and critical modules were screened. A regulatory network between prognostic DE RBPs and differentially expressed transcription factors was constructed. Besides, a risk signature was built based on prognostic DE RBPs, which is able to predict overall survival of CRC patients with high accuracy. In conclusion, the results provided a comprehensive understanding of the functions of RBPs in CRC, as well as an RBP-related prognostic signature.

Regulatory network of fruit ripening: current understanding and future challenges

Fruit ripening is a developmental process that is spatio-temporally tuned at multiple levels. Molecular dissections of the mechanisms underlying the ripening process have revealed a network encompassed by hormones, transcriptional regulators, epigenomic modifications and other regulatory elements that directly determine fruit quality and the postharvest commodity of fresh produce. Many studies have addressed the important roles of ethylene, abscisic acid (ABA) and other hormones in regulating fruit ripening. Recent studies have shown that some spontaneous mutants for tomato transcription factors (TFs) have resulted from loss-of-function or dominant-negative mutations. Unlike in DNA methylation variation, the histone mark H3K27me3 may be conserved and prevents the transcriptional feedback circuit from generating autocatalytic ethylene. These observations of a network of partially redundant component indicate the need to improve our current understanding. Here, we focussed on the recent advances and future challenges in investigations of the molecular mechanisms of fruit ripening. We also identified several issues that still need to be addressed in future studies.

Gene Regulatory Networks of Penicillium echinulatum 2HH and Penicillium oxalicum 114-2 Inferred by a Computational Biology Approach

Penicillium echinulatum 2HH and Penicillium oxalicum 114-2 are well-known cellulase fungal producers. However, few studies addressing global mechanisms for gene regulation of these two important organisms are available so far. A recent finding that the 2HH wild-type is closely related to P. oxalicum leads to a combined study of these two species. Firstly, we provide a global gene regulatory network for P. echinulatum 2HH and P. oxalicum 114-2, based on TF-TG orthology relationships, considering three related species with well-known regulatory interactions combined with TFBSs prediction. The network was then analyzed in terms of topology, identifying TFs as hubs, and modules. Based on this approach, we explore numerous identified modules, such as the expression of cellulolytic and xylanolytic systems, where XlnR plays a key role in positive regulation of the xylanolytic system. It also regulates positively the cellulolytic system by acting indirectly through the cellodextrin induction system. This remarkable finding suggests that the XlnR-dependent cellulolytic and xylanolytic regulatory systems are probably conserved in both P. echinulatum and P. oxalicum. Finally, we explore the functional congruency on the genes clustered in terms of communities, where the genes related to cellular nitrogen, compound metabolic process and macromolecule metabolic process were the most abundant. Therefore, our approach allows us to confer a degree of accuracy regarding the existence of each inferred interaction.

Identification of Novel SARS-CoV-2 Drug Targets by Host MicroRNAs and Transcription Factors Co-regulatory Interaction Network Analysis

Understanding the host regulatory mechanisms opposing virus infection and virulence can provide actionable insights to identify novel therapeutics against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We have used a network biology approach to elucidate the crucial factors involved in host responses involving host-microRNA (miRNA) interactions with host and virus genes using recently published experimentally verified protein-protein interaction data. We were able to identify 311 host genes to be potentially targetable by 2,197 human miRNAs. These miRNAs are known to be involved in various biological processes, such as T-cell differentiation and activation, virus replication, and immune system. Among these, the anti-viral activity of 38 miRNAs to target 148 host genes is experimentally validated. Six anti-viral miRNAs, namely, hsa-miR-1-3p, hsa-miR-17-5p, hsa-miR-199a-3p, hsa-miR-429, hsa-miR-15a-5p, and hsa-miR-20a-5p, are previously reported to be anti-viral in respiratory diseases and were found to be downregulated. The interaction network of the 2,197 human miRNAs and interacting transcription factors (TFs) enabled the identification of 51 miRNAs to interact with 77 TFs inducing activation or repression and affecting gene expression of linked genes. Further, from the gene regulatory network analysis, the top five hub genes HMOX1, DNMT1, PLAT, GDF1, and ITGB1 are found to be involved in interferon (IFN)-α2b induction, epigenetic modification, and modulation of anti-viral activity. The comparative miRNAs target identification analysis in other respiratory viruses revealed the presence of 98 unique host miRNAs targeting SARS-CoV-2 genome. Our findings identify prioritized key regulatory interactions that include miRNAs and TFs that provide opportunities for the identification of novel drug targets and development of anti-viral drugs.

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.

A semi-automatic cell type annotation method for single-cell RNA sequencing dataset

Single-cell RNA sequencing (scRNA-seq) has been widely applied to provide insights into the cell-by-cell expression difference in a given bulk sample. Accordingly, numerous analysis methods have been developed. As it involves simultaneous analyses of many cell and genes, efficiency of the methods is crucial. The conventional cell type annotation method is laborious and subjective. Here we propose a semi-automatic method that calculates a normalized score for each cell type based on user-supplied cell type–specific marker gene list. The method was applied to a publicly available scRNA-seq data of mouse cardiac non-myocyte cell pool. Annotating the 35 t-stochastic neighbor embedding clusters into 12 cell types was straightforward, and its accuracy was evaluated by constructing co-expression network for each cell type. Gene Ontology analysis was congruent with the annotated cell type and the corollary regulatory network analysis showed upstream transcription factors that have well supported literature evidences. The source code is available as an R script upon request.

Integrative microRNA and mRNA expression profiling in acute aristolochic acid nephropathy in mice

In acute aristolochic acid nephropathy (AAN), aristolochic acid (AA) induces renal injury and tubulointerstitial fibrosis. However, the roles of microRNAs (miRNAs/miRs) and mRNAs involved in AAN are not clearly understood. The aim of the present study was to examine AA‑induced genome‑wide differentially expressed (DE) miRNAs and DE mRNAs using deep sequencing in mouse kidneys, and to analyze their regulatory networks. In the present self‑controlled study, mice were treated with 5 mg/kg/day AA for 5 days, following unilateral nephrectomy. AA‑induced renal injury and tubulointerstitial fibrosis were detected using hematoxylin and eosin staining and Masson's trichrome staining in the mouse kidneys. A total of 82 DE miRNAs and 4,605 DE mRNAs were identified between the AA‑treated group and the self‑control group. Of these DE miRNAs and mRNAs, some were validated using reverse transcription‑quantitative PCR. Expression levels of the profibrotic miR‑21, miR‑433 and miR‑132 families were significantly increased, whereas expression levels of the anti‑fibrotic miR‑122‑5p and let‑7a‑1‑3p were significantly decreased. Functions and signaling pathways associated with the DE miRNAs and mRNAs were analyzed using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG). A total of 767 DE pairs (in opposing directions) of miRNAs and their mRNA targets were identified. Among these, regulatory networks of miRNAs and mRNAs were analyzed using KEGG to identify enriched signaling pathways and extracellular matrix‑associated pathways. In conclusion, the present study identified genome‑wide DE miRNAs and mRNAs in the kidneys of AA‑treated mice, as well as their regulatory pairs and signaling networks. The present results may improve the understanding of the role of DE miRNAs and their mRNA targets in the pathophysiology of acute AAN.

Systematic Analysis of Competing Endogenous RNA Networks in Diffuse Large B-Cell Lymphoma and Hodgkin's Lymphoma

Lymphoma is a systemic malignancy, originating from the lymphatic system, which accounts for 3 to 4% of all tumors. There are two major subtypes of lymphoma, namely, diffuse large B-cell lymphoma (DLBCL) and Hodgkin’s lymphoma (HL). Elucidation of the pathogenesis of these two lymphoma types is crucial for the identification of potential therapeutic targets. Compared with the corresponding knowledge of other diseases, the understanding of the regulatory networks involved in DLBCL and HL is relatively deficient. To address this, we comprehensively analyzed the mRNAs, lncRNAs, and miRNAs that were differentially expressed between normal and tumor samples of DLBCL and HL. In addition, functional enrichment analysis of the differentially expressed mRNAs was performed. We constructed two specific ceRNA networks of DLBCL and HL. The pathways enriched by dysregulated mRNAs in DLBCL and HL were mainly involved in immune responses, transcription process, and metabolism process. The ceRNA network analysis revealed that 45 ceRNAs were shared between the two ceRNA networks, including five pivotal lncRNAs (MALAT1, CTBP1-AS, THUMPD3-AS, PSMA3-AS1, and NUTM2A-AS1). In addition, we proposed a DLBCL survival risk model based on a DLBCL-specific network constructed by Lasso regression analysis. The model, which is based on eight mRNAs, exhibited excellent performance in regard to predicting outcomes in DLBCL patients, with a p value of 0.0017 and AUC of 0.9783. In summary, although the molecular mechanisms underlying tumorigenesis in DLBCL and HL were quite different, the same pivotal lncRNAs acted as key regulators. Our findings identify novel potential prognostic and therapeutic targets for DLBCL and HL.

Multiple Isolated Transcription Factors Act as Switches and Contribute to Species Uniqueness

Mammals have variable numbers (1300-2000) of transcription factors (TFs), but the reasons for this large variation are unclear. To investigate general TF patterns, we de novo identified 156,906 TFs from 96 mammalian species. We identified more than 500 human isolated TFs that are rarely reported in human TF-to-TF networks. Mutations in the genes of these TFs were less lethal than those of connected TFs. Consequently, these isolated TFs are more tolerant of changes and have become unique during speciation. They may also serve as a source of variation for TF evolution. Reconciliation of TF-family phylogenetic trees with a mammalian species tree revealed an average of 37.8% TF gains and 15.0% TF losses over 177 million years, which implies that isolated TFs are pervasive in mammals. Compared with non-TF interacting genes, TF-interacting genes have unique TF profiles and have higher expression levels in mice than in humans. Different expression levels of the same TF-interacting gene contribute to species-specific phenotypes. Formation and loss of isolated TFs enabling unique TF profiles may provide variable switches that adjust divergent expression profiles of target genes to generate species-specific phenotypes, thereby making species unique.

Develop a circular RNA-related regulatory network associated with prognosis of gastric cancer

Background

In gastric cancer (GC), circular RNAs (circRNAs) mainly play an important role in miRNA sponge, which not only indicate long‐term survival and prognosis but also increase resistance to the apoptosis. The purpose of the study is to explore new circRNAs and their underlying mechanisms in GC.

Method

Through rigorous retrieval strategies, we used the sva package to analyze and identify differentially expressed circRNAs (DECs) from three Gene Expression Omnibus microarray datasets (GSE83521, GSE89143, and GSE78092). Online website CSCD and CircInteractome were used to reveal the binding sites between miRNAs and DECs. The possible target miRNAs of the DECs identified based on miRNAs, and Cytoscape was used to create a regulatory network of circRNA‐miRNA‐mRNA and identified the hub genes which were further validated using The Cancer Genome Atlas database and Human Protein Atlas.

Results

Twenty‐eight DECs were obtained using the sva package. A regulatory network of circRNA‐miRNA‐mRNA (competing endogenous RNA) containing 15 circRNAs, 24 miRNAs, and 158 genes was identified. A protein‐protein interaction network based on the 158 genes was established, and further determined that 10 hub genes (SKA1, ANLN, CHEK1, SKA3, TOP2A, BIRC5, RRM2, NCAPG2, FANCI, and RAD51) were associated with some cancer‐related pathways based on the functional enrichment analysis. Finally, six hub genes (BIRC5, TOP2A, FANCI, NCAPG2, RAD51, and RRM2) were proven to influence the overall survival of GC.

Conclusion

Our study established a circRNA‐miRNA‐mRNA regulatory network and defined six circRNA‐related hub genes in GC, which could serve as potential therapeutic targets or prognostic biomarker for GC treatment.

Construction of a circRNA-miRNA-mRNA Network Related to Macrophage Infiltration in Hepatocellular Carcinoma

Immune cells in the tumor microenvironment play a crucial role in regulating tumor progression. The circular RNA (circRNA) regulatory network involved in immune cell infiltration in hepatocellular carcinoma (HCC) remains largely unknown. In this study, the “estimate the proportion of immune and cancer cells” (EPIC) application is used to evaluate the fractions of immune cells, cancer-associated fibroblasts, and endothelial cells in HCC from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Patients with a high macrophage fraction have better overall survival, and macrophage fraction is an independent prognostic factor for HCC. Next, the common differentially expressed mRNAs (DEmRNAs), miRNAs (DEmiRNAs), and circRNAs (DEcircRNAs) between paired tumor and non-tumor tissues are screened out from the TCGA and/or GEO databases. Through spearman correlation analysis, the macrophage-related DEmRNAs are identified to construct a circRNA-miRNA-mRNA regulatory network, which includes 6 DEcircRNAs, 7 DEmiRNAs, and 45 DEmRNAs. Functional enrichment analysis reveals that these DEmRNAs are mainly involved in immune-related processes. Furthermore, six hub DEmRNAs are identified to establish a hub circRNA regulatory network. Among the DEmRNAs in the network, PRC1 is identified as the most influential node. PRC1 high expression is correlated with poor prognosis and low macrophage infiltration in HCC. Taken together, we identify a certain circRNA regulatory network related to macrophage infiltration and provide novel insight into the mechanism of study and therapeutic targets for HCC.

NetExtractor: Extracting a Cerebellar Tissue Gene Regulatory Network Using Differentially Expressed High Mutual Information Binary RNA Profiles

Bigenic expression relationships are conventionally defined based on metrics such as Pearson or Spearman correlation that cannot typically detect latent, non-linear dependencies or require the relationship to be monotonic. Further, the combination of intrinsic and extrinsic noise as well as embedded relationships between sample sub-populations reduces the probability of extracting biologically relevant edges during the construction of gene co-expression networks (GCNs). In this report, we address these problems via our NetExtractor algorithm. NetExtractor examines all pairwise gene expression profiles first with Gaussian mixture models (GMMs) to identify sample sub-populations followed by mutual information (MI) analysis that is capable of detecting non-linear differential bigenic expression relationships. We applied NetExtractor to brain tissue RNA profiles from the Genotype-Tissue Expression (GTEx) project to obtain a brain tissue specific gene expression relationship network centered on cerebellar and cerebellar hemisphere enriched edges. We leveraged the PsychENCODE pre-frontal cortex (PFC) gene regulatory network (GRN) to construct a cerebellar cortex (cerebellar) GRN associated with transcriptionally active regions in cerebellar tissue. Thus, we demonstrate the utility of our NetExtractor approach to detect biologically relevant and novel non-linear binary gene relationships.

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.

The differences of regulatory networks between papillary and anaplastic thyroid carcinoma: an integrative transcriptomics study

BACKGROUND

Unlike papillary thyroid cancer (PTC), anaplastic thyroid carcinoma (ATC) is extremely aggressive and rapidly lethal without effective therapies. However, the differences of master regulators and regulatory networks between PTC and ATC remain unclear. Methods: Three representative datasets comprising 32 ATC, 69 PTC, and 78 normal thyroid tissue samples were combined to form a large dataset. Differentially expressed genes (DEGs) were identified and enriched by limma package and gene set enrichment analysis, respectively. Subsequently, protein-protein interaction network and transcription factors (TFs) regulatory network were constructed to identify gene modules and master regulators. Further, master regulators were validated by RT-PCR and western blot. Finally, Kaplan-Meier plotter was applied to evaluate their prognostic values. Results: A total of 560 DEGs were identified as ATC-specific malignant signature. The regulatory network analysis showed that nine master regulators were significantly correlated with three gene modules and potentially regulated the expression of DEGs in three gene modules, respectively. Furthermore, CREB3L1, FOSL2, E2F1 and CAT were significantly associated with overall survival of thyroid cancer patients. FOXM1, FOSL2, MYBL2, AVEN and E2F1 were unfavorable factors of recurrence-free survival (RFS), while CAT was a favorable factor of RFS. RT-PCR and western blot confirmed that six TFs were obviously up-regulated in ATC tissues/cell line as compared with PTC and normal thyroid tissues/cell lines, respectively. In addition, 19 ATC-specific kinases were identified to illustrate the potential post-translational modification. Conclusions: Our findings provide a comprehensive insight into malignant mechanism of ATC, which may indicate their value in the future investigation of ATC.

Melatonin confers cadmium tolerance by modulating critical heavy metal chelators and transporters in radish plants

Cadmium (Cd) is an environmental pollutant that causes health hazard to living organisms. Melatonin (MT) has emerged as a ubiquitous pleiotropic molecule capable of coordinating heavy metal (HM) stresses in plants. However, it remains unclear how melatonin mediates Cd homeostasis and detoxification at transcriptional and/or post-transcriptional levels in radish. Herein, the activities of five key antioxidant enzymes were increased, while root and shoot Cd contents were dramatically decreased by melatonin. A combined small RNA and transcriptome sequencing analysis showed that 14 differentially expressed microRNAs (DEMs) and 966 differentially expressed genes (DEGs) were shared between the Cd and Cd + MT conditions. In all, 23 and ten correlated miRNA-DEG pairs were identified in Con vs. Cd and Con vs. Cd + MT comparisons, respectively. Several DEGs encoding yellow stripe 1-like (YSL), heavy metal ATPases (HMA), and ATP-binding cassette (ABC) transporters were involved in Cd transportation and sequestration in radish. Root exposure to Cd²⁺ induced several specific signaling molecules, which consequently trigger some HM chelators, transporters, and antioxidants to achieve reactive oxygen species (ROS) scavenging and detoxification and eliminate Cd toxicity in radish plants. Notably, transgenic analysis revealed that overexpression of the RsMT1 (Metallothionein 1) gene could enhance Cd tolerance of tobacco plants, indicating that the exogenous melatonin confers Cd tolerance, which might be attributable to melatonin-mediated upregulation of RsMT1 gene in radish plants. These results could contribute to dissecting the molecular basis governing melatonin-mediated Cd stress response in plants and pave the way for high-efficient genetically engineering low-Cd-content cultivars in radish breeding programs.

Multi-Locus GWAS of Quality Traits in Bread Wheat: Mining More Candidate Genes and Possible Regulatory Network

In wheat breeding, improved quality traits, including grain quality and dough rheological properties, have long been a critical goal. To understand the genetic basis of key quality traits of wheat, two single-locus and five multi-locus GWAS models were performed for six grain quality traits and three dough rheological properties based on 19, 254 SNPs in 267 bread wheat accessions. As a result, 299 quantitative trait nucleotides (QTNs) within 105 regions were identified to be associated with these quality traits in four environments. Of which, 40 core QTN regions were stably detected in at least three environments, 19 of which were novel. Compared with the previous studies, these novel QTN regions explained smaller phenotypic variation, which verified the advantages of the multi-locus GWAS models in detecting important small effect QTNs associated with complex traits. After characterization of the function and expression in-depth, 67 core candidate genes involved in protein/sugar synthesis, histone modification and the regulation of transcription factor were observed to be associated with the formation of grain quality, which showed that multi-level regulations influenced wheat grain quality. Finally, a preliminary network of gene regulation that may affect wheat quality formation was inferred. This study verified the power and reliability of multi-locus GWAS methods in wheat quality trait research, and increased the understanding of wheat quality formation mechanisms. The detected QTN regions and candidate genes in this study could be further used for gene cloning and marker-assisted selection in high-quality breeding of bread wheat.

A MARTX Toxin rtxA Gene Is Controlled by Host Environmental Signals through a CRP-Coordinated Regulatory Network in Vibrio vulnificus

A MARTX toxin, RtxA, is an essential virulence factor of many pathogens, including Vibrio species. H-NS and HlyU repress and derepress, respectively, rtxA expression of a life-threatening pathogen, Vibrio vulnificus. We found that Lrp directly activates rtxA independently of H-NS and HlyU, and leucine inhibits the Lrp-mediated activation of rtxA. Furthermore, we demonstrated that CRP represses rtxA but derepresses in the presence of exogenous glucose. CRP represses rtxA not only directly by binding to upstream of rtxA but also indirectly by repressing lrp and hlyU. This is the first report of a regulatory network comprising CRP, Lrp, H-NS, and HlyU, which coordinates the rtxA expression in response to environmental signals such as leucine and glucose during infection. This elaborate regulatory network will enhance the fitness of V. vulnificus and contribute to its successful infection within the host.

A multifunctional autoprocessing repeats-in-toxin (MARTX) toxin plays an essential role in the virulence of many pathogens, including a fulminating human pathogen Vibrio vulnificus. H-NS and HlyU repress and derepress expression of the MARTX toxin gene rtxA in V. vulnificus, respectively. However, little is known about other regulatory proteins and environmental signals involved in rtxA regulation. In this study, we found that a leucine-responsive regulatory protein (Lrp) activates rtxA by binding directly and specifically to the rtxA promoter, P_rtxA. Phased hypersensitivity resulting from DNase I cleavage of the P_rtxA regulatory region suggests that Lrp probably induces DNA bending in P_rtxA. Lrp activates P_rtxA independently of H-NS and HlyU, and leucine inhibits Lrp binding to P_rtxA and reduces the Lrp-mediated activation. Furthermore, a cyclic AMP receptor protein (CRP) represses P_rtxA, and exogenous glucose relieves the CRP-mediated repression. Biochemical and mutational analyses demonstrated that CRP binds directly and specifically to the upstream region of P_rtxA, which presumably alters the DNA conformation in P_rtxA and thus represses rtxA. Moreover, CRP represses expression of lrp and hlyU by binding directly to their upstream regions, forming coherent feed-forward loops with Lrp and HlyU. In conclusion, expression of rtxA is controlled by a regulatory network comprising CRP, Lrp, H-NS, and HlyU in response to changes in host environmental signals such as leucine and glucose. This collaborative regulation enables the elaborate expression of rtxA, thereby enhancing the fitness and pathogenesis of V. vulnificus during the course of infection.

ATAC-pipe: general analysis of genome-wide chromatin accessibility

Assay of Transposase-Accessible Chromatin by deep sequencing (ATAC-seq) has been widely used to profile the chromatin accessibility genome-wide. For the absence of an integrated scheme for deep data mining of specific biological issues, here we present ATAC-pipe, an efficient pipeline for general analysis of chromatin accessibility data obtained from ATAC-seq experiments. ATAC-pipe captures information includes not only the quality of original data and genome-wide chromatin accessibility but also signatures of significant differential peaks, transcription factor (TF) occupancy and nucleosome positions around regulatory sites. In addition, ATAC-pipe automatically converts statistic results into intuitive plots at publication quality, such as the read length distribution, heatmaps of sample clustering and cell-type-specific regulatory elements, enriched TF occupancy with motifs footprints and TF-driven regulatory networks. ATAC-pipe provides convenient workflow for researchers to study chromatin accessibility and gene regulation. Availability https://github.com/QuKunLab/ATAC-pipe.

hsa-miR-106b-5p participates in the development of chronic thromboembolic pulmonary hypertension via targeting matrix metalloproteinase 2

Background

Chronic thromboembolic pulmonary hypertension (CTEPH) is characterized by elevated pressure in pulmonary arteries. This study was performed to explore the critical miRNAs and genes affecting the pathogenesis of CTEPH.

Methods

GSE56914 dataset (10 CTEPH whole blood samples and 10 control samples) was downloaded from the Gene Expression Omnibus database. Using limma package, the differentially expressed miRNAs (DE-miRNAs) were acquired. After miRNA-target pairs were obtained using miRWalk2.0 tool, a miRNA-target regulatory network was built by Cytoscape software. Using DAVID tool, significantly enriched pathways involving the target genes were identified. Moreover, the protein–protein interaction network and transcription factor-target regulatory network were built by the Cytoscape software. Additionally, quantitative real-time PCR (qRT-PCR) experiments and luciferase assay were conducted to validate miRNA/gene expression and miRNA–target regulatory relationship, respectively.

Results

There were 25 DE-miRNAs (8 up-regulated and 17 down-regulated) between CTEPH and control groups. The target genes of has-let-7b-3p, has-miR-17-5p, has-miR-3202, has-miR-106b-5p, and has-miR-665 were enriched in multiple pathways such as “Insulin secretion”. qRT-PCR analysis confirmed upregulation of hsa-miR-3202, hsa-miR-665, and matrix metalloproteinase 2 (MMP2) as well as downregulation of hsa-let-7b-3p, hsa-miR-17-5p, and hsa-miR-106b-5p. Luciferase assay indicated that MMP2 was negatively mediated by hsa-miR-106b-5p.

Conclusions

These miRNAs and genes were associated with the pathogenesis of CTEPH. Besides, hsa-miR-106b-5p was involved in the development of CTEPH via targeting MMP2.

Unraveling the Complex Interplay of Fis and IHF Through Synthetic Promoter Engineering

Bacterial promoters are usually formed by multiple cis-regulatory elements recognized by a plethora of transcriptional factors (TFs). From those, global regulators are key elements since these TFs are responsible for the regulation of hundreds of genes in the bacterial genome. For instance, Fis and IHF are global regulators that play a major role in gene expression control in Escherichia coli, and usually, multiple cis-regulatory elements for these proteins are present at target promoters. Here, we investigated the relationship between the architecture of the cis-regulatory elements for Fis and IHF in E. coli. For this, we analyze 42 synthetic promoter variants harboring consensus cis-elements for Fis and IHF at different distances from the core -35/-10 region and in various numbers and combinations. We first demonstrated that although Fis preferentially recognizes its consensus cis-element, it can also recognize, to some extent, the consensus-binding site for IHF, and the same was true for IHF, which was also able to recognize Fis binding sites. However, changing the arrangement of the cis-elements (i.e., the position or number of sites) can completely abolish the non-specific binding of both TFs. More remarkably, we demonstrated that combining cis-elements for both TFs could result in Fis and IHF repressed or activated promoters depending on the final architecture of the promoters in an unpredictable way. Taken together, the data presented here demonstrate how small changes in the architecture of bacterial promoters could result in drastic changes in the final regulatory logic of the system, with important implications for the understanding of natural complex promoters in bacteria and their engineering for novel applications.

Effects of CORO2A on Cell Migration and Proliferation and Its Potential Regulatory Network in Breast Cancer

Coronin 2A (CORO2A) is a novel component of the N-CoR (nuclear receptor co-repressor) complex. Abnormal CORO2A expression is associated with carcinogenesis. We used databases from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO), and analyzed CORO2A expression and gene regulation networks in breast cancer. Expression was analyzed using GEO and TCGA database and further validated in breast cancer samples collected in our clinic. The prognostic value of CORO2A was explored by using the Kaplan–Meier survival analysis and Cox proportional hazards regression analysis. LinkedOmics was used to identify coexpressed genes associated with CORO2A. After analyzing the intersection of coexpressed genes correlated with CORO2A and differentially expressed genes after CORO2A silencing, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of the intersecting genes were conducted by using FunRich software. Transwell assays were performed in breast cancer cells to determine the effect of CORO2A on cell migration. MTS, colony formation, and cell cycle distribution assays were performed in breast cancer cells to determine the effect of CORO2A on cell proliferation. Gene enrichment analysis was employed to explore the target networks of transcription factors and miRNAs. We found that CORO2A was upregulated and that the elevated expression of CORO2A was associated with poor overall survival (OS) and relapse-free survival (RFS) in TNBC patients. Further bioinformatics analysis of public sequencing data and our own RNA-Seq data revealed that CORO2A was probably involved in the epithelial-to-mesenchymal transition process and might have a significant effect on the migration of breast cancer cells, which might be mediated via pathways involving several miRNAs and MYC transcription factors. Functionally, the knockdown of CORO2A inhibited cell migration, decreased viability, and colony formation and induced cell cycle arrest in the G0/G1 phase in breast cancer cells. These results demonstrate that bioinformatics-based analysis efficiently reveals information about CORO2A expression and its potential regulatory networks in breast cancer, laying a foundation for further mechanistic research on the role of CORO2A in carcinogenesis. Moreover, CORO2A promotes the migration and proliferation of breast cancer cells and may have an important function in breast cancer progression. CORO2A is a potential prognostic predictor for TNBC patients. Targeting CORO2A may provide promising therapy strategies for breast cancer treatment.

Transcriptional regulatory network controlling the ontogeny of hematopoietic stem cells

In this study from Gao et al., the authors performed RNA-seq and histone mark ChIP-seq to define the transcriptomes and epigenomes of cells representing key developmental stages of HSC ontogeny in mice. Using a novel computational algorithm, target inference via physical connection (TIPC), they constructed developmental stage-specific transcriptional regulatory networks by linking enhancers and predicted bound transcription factors to their target promoters, thus providing a useful resource for uncovering regulators of HSC formation.

Hematopoietic stem cell (HSC) ontogeny is accompanied by dynamic changes in gene regulatory networks. We performed RNA-seq and histone mark ChIP-seq to define the transcriptomes and epigenomes of cells representing key developmental stages of HSC ontogeny in mice. The five populations analyzed were embryonic day 10.5 (E10.5) endothelium and hemogenic endothelium from the major arteries, an enriched population of prehematopoietic stem cells (pre-HSCs), fetal liver HSCs, and adult bone marrow HSCs. Using epigenetic signatures, we identified enhancers for each developmental stage. Only 12% of enhancers are primed, and 78% are active, suggesting the vast majority of enhancers are established de novo without prior priming in earlier stages. We constructed developmental stage-specific transcriptional regulatory networks by linking enhancers and predicted bound transcription factors to their target promoters using a novel computational algorithm, target inference via physical connection (TIPC). TIPC predicted known transcriptional regulators for the endothelial-to-hematopoietic transition, validating our overall approach, and identified putative novel transcription factors, including the broadly expressed transcription factors SP3 and MAZ. Finally, we validated a role for SP3 and MAZ in the formation of hemogenic endothelium. Our data and computational analyses provide a useful resource for uncovering regulators of HSC formation.

Weighted gene co-expression network analysis to investigate the key genes implicated in global brain ischemia/reperfusion injury in rats

BACKGROUND

Ischemia/reperfusion (I/R) refers to situations where blood is perfused into ischemic or hypoxic tissues, potentially resulting in an inflammatory response and oxidative injury.

OBJECTIVES

This study was conducted to explore the pathogenesis of I/R injury.

MATERIAL AND METHODS

GSE82146 was extracted from the Gene Expression Omnibus, consisting of 15 complete global brain ischemia (CGBI) reperfusion hippocampus samples and 12 non-ischemic control (NIC) hippocampus samples. The differentially expressed genes (DEGs) between the CGBI and NIC samples were selected using LIMMA package, and were then analyzed with weighted gene co-expression network analysis (WGCNA). Using DAVID software, the DEGs in significant modules were run through enrichment analysis. The DEGs in significant modules were merged, and then a protein-protein interaction (PPI) network was built for them using Cytoscape software. After miRNAs and transcription factors (TFs) were predicted for the DEGs using the WebGestalt tool, a TF-miRNA-target regulatory network was built using Cytoscape software. Furthermore, quantitative real-time polymerase chain reaction (qRT-PCR) analysis was conducted to detect the levels of key genes.

RESULTS

There were 390 DEGs in the CGBI samples. Based on WGCNA, brown and turquoise modules were screened as CGBI-associated modules. In the PPI network, key nodes HSP90AA1 and HSPA5 were able to interact with each other. In the regulatory network, MYC, HSF1 and miR-22 had higher degree values. Moreover, HSPA5 was targeted by MYC in the regulatory network. In addition, upregulated HSPB1 and HMOX1, as well as downregulated NR4A2, were confirmed with qRT-PCR analysis.

CONCLUSIONS

HSPB1, HMOX1 and NR4A2 were the key genes correlated with I/R injury. Additionally, HSP90AA1, HSPA5, MYC, HSF1, and miR-22 might be related to the pathogenesis of I/R injury.

ZmIBH1-1 regulates plant architecture in maize

Leaf angle (LA) is a critical agronomic trait in maize, with more upright leaves allowing higher planting density, leading to more efficient light capture and higher yields. A few genes responsible for variation in LA have been identified by map-based cloning. In this study, we cloned maize ZmIBH1-1, which encodes a bHLH transcription factor with both a basic binding region and a helix-loop-helix domain, and the results of qRT-PCR showed that it is a negative regulator of LA. Histological analysis indicated that changes in LA were mainly caused by differential cell wall lignification and cell elongation in the ligular region. To determine the regulatory framework of ZmIBH1-1, we conducted RNA-seq and DNA affinity purification (DAP)-seq analyses. The combined results revealed 59 ZmIBH1-1-modulated target genes with annotations, and they were mainly related to the cell wall, cell development, and hormones. Based on the data, we propose a regulatory model for the control of plant architecture by ZmIBH1-1 in maize.

In silico analysis of the molecular regulatory networks in peripheral arterial occlusive disease

BACKGROUND

Peripheral arterial occlusive disease (PAOD) is a global public health concern that decreases the quality of life of the patients and can lead to disabilities and death. The aim of this study was to identify the genes and pathways associated with PAOD pathogenesis, and the potential therapeutic targets.

METHODS

Differentially expressed genes (DEGs) and miRNAs related to PAOD were extracted from the GSE57691 dataset and through text mining. Additionally, bioinformatics analysis was applied to explore gene ontology, pathways and protein-protein interaction of those DEGs. The potential miRNAs targeting the DEGs and the transcription factors (TFs) regulating miRNAs were predicted by multiple different databases.

RESULTS

A total of 59 DEGs were identified, which were significantly enriched in the inflammatory response, immune response, chemokine-mediated signaling pathway and JAK-STAT signaling pathway. Thirteen genes including IL6, CXCL12, IL1B, and STAT3 were hub genes in protein-protein interaction network. In addition, 513 miRNA-target gene pairs were identified, of which CXCL12 and PTPN11 were the potential targets of miRNA-143, and IL1B of miRNA-21. STAT3 was differentially expressed and regulated 27 potential target miRNAs including miRNA-143 and miRNA-21 in TF-miRNA regulatory network.

CONCLUSION

In summary, inflammation, immune response and STAT3-mediated miRNA-target genes axis play an important role in PAOD development and progression.

Identification of MicroRNA-Target Gene-Transcription Factor Regulatory Networks in Colorectal Adenoma Using Microarray Expression Data

Objective

The aim of the study was to find the key genes, microRNAs (miRNAs) and transcription factors (TFs) and construct miRNA-target gene-TF regulatory networks to investigate the underlying molecular mechanism in colorectal adenoma (CRA).

Methods

Four mRNA expression datasets and one miRNA expression dataset were downloaded from Gene Expression Omnibus (GEO) database. Differentially expressed miRNAs (DEMs) and differentially expressed genes (DEGs) were identified between CRA and normal samples. Moreover, functional enrichment analysis for DEGs was carried out utilizing the Cytoscape-plugin, known as ClueGO. These DEGs were mapped to STRING database to construct a protein-protein interaction (PPI) network. Then, a miRNA-target gene regulatory network was established to screen key DEMs. In addition, similar workflow of the analyses were also performed comparing the CRC samples with CRA ones to screen key DEMs. Finally, miRNA-target gene-TF regulatory networks were constructed for these key DEMs using iRegulon plug-in in Cytoscape.

Results

We identified 514 DEGs and 167 DEMs in CRA samples compared to healthy samples. Functional enrichment analysis revealed that these DEGs were significantly enriched in several terms and pathways, such as regulation of cell migration and bile secretion pathway. A PPI network was constructed including 325 nodes as well as 890 edges. A total of 59 DEGs and 65 DEMs were identified in CRC samples compared to CRA ones. In addition, Two key DEMs in CRA samples compared to healthy samples were identified, such as hsa-miR-34a and hsa-miR-96. One key DEM, hsa-miR-29c, which was identified when we compared the differentially expressed molecules found in the comparison CRA versus normal samples to the ones obtained in the comparison CRC versus CRA, was also identified in CRC samples compared to CRA ones. The miRNA-target gene-TF regulatory networks for these key miRNAs included two TFs, one TF and five TFs, respectively.

Conclusion

These identified key genes, miRNA, TFs and miRNA-target gene-TF regulatory networks associated with CRA, to a certain degree, may provide some hints to enable us to better understand the underlying pathogenesis of CRA.

Global control of bacterial nitrogen and carbon metabolism by a PTSNtr-regulated switch

The nitrogen-related phosphotransferase system (PTS^Ntr) of Rhizobium leguminosarum bv. viciae 3841 transfers phosphate from PEP via PtsP and NPr to two output regulators, ManX and PtsN. ManX controls central carbon metabolism via the tricarboxylic acid (TCA) cycle, while PtsN controls nitrogen uptake, exopolysaccharide production, and potassium homeostasis, each of which is critical for cellular adaptation and survival. Cellular nitrogen status modulates phosphorylation when glutamine, an abundant amino acid when nitrogen is available, binds to the GAF sensory domain of PtsP, preventing PtsP phosphorylation and subsequent modification of ManX and PtsN. Under nitrogen-rich, carbon-limiting conditions, unphosphorylated ManX stimulates the TCA cycle and carbon oxidation, while unphosphorylated PtsN stimulates potassium uptake. The effects are reversed with the phosphorylation of ManX and PtsN, occurring under nitrogen-limiting, carbon-rich conditions; phosphorylated PtsN triggers uptake and nitrogen metabolism, the TCA cycle and carbon oxidation are decreased, while carbon-storage polymers such as surface polysaccharide are increased. Deleting the GAF domain from PtsP makes cells "blind" to the cellular nitrogen status. PTS^Ntr constitutes a switch through which carbon and nitrogen metabolism are rapidly, and reversibly, regulated by protein:protein interactions. PTS^Ntr is widely conserved in proteobacteria, highlighting its global importance.

Identification of Autophagy-Associated Biomarkers and Corresponding Regulatory Factors in the Progression of Colorectal Cancer

Autophagy is a self-degradation process that maintains homeostasis against stress in cells. Autophagy dysfunction plays a central role in the development of tumors, such as colorectal cancer (CRC). In this study, autophagy-related differentially expressed genes, their downstream functions, and upstream regulatory factors including RNA-binding proteins (RBP) involved in programmed cell death in the CRC were investigated. Transcription factors (TFs) and miRNAs have been shown to mainly regulate autophagy genes. Interestingly, we found that some of the RBP in the CRC, such as DDX17, SETDB1, and POLR3A, play an important regulatory role in maintaining autophagy at a basal level during growth by acting as TFs that regulate autophagy. Promoter methylations showed negative regulations on differentially expressed autophagy gene (DEAG), while copy number variations revealed a positive role in them. A proportional hazards regression analysis indicated that using autophagy-related prognostic signature can divide patients into high-risk and low-risk groups. Autophagy associated FDA-approved drugs were studied by a prognostic network. This would contribute to the identifications of new potential molecular therapeutic targets for CRC.

Identification of Differentially Expressed MicroRNAs and Their Potential Target Genes in Adipose Tissue from Pigs with Highly Divergent Backfat Thickness

Simple Summary

The role of microRNA in fat deposition is very important and not clearly understood. We detected 318 pig microRNAs (miRNAs), among high and low backfat tissue samples, by high throughput sequencing. Among them, 18 miRNAs were differentially expressed between the high and low backfat groups. Some of the differentially expressed miRNAs were involved mainly in lipid and carbohydrate metabolism, and glycan biosynthesis and metabolism. In addition, in silico analysis of the mRNA and miRNA transcriptomes, revealed possible regulatory relationships for fat deposition. In particular, three miRNA–mRNA pairs, miR-137–PPARGC1A, miR-141–FASN, and miR-122-5p–PKM, were identified as candidate key regulators of fat deposition. Our findings provide an important insight into miRNA expression patterns in backfat tissue of pig and new insights into the regulatory mechanisms of fat deposition in pig.

Abstract

Fatty traits are very important in pig production. However, the role of microRNAs (miRNAs) in fat deposition is not clearly understood. In this study, we compared adipose miRNAs from three full-sibling pairs of female Landrace pigs, with high and low backfat thickness, to investigate the associated regulatory network. We obtained an average of 17.29 million raw reads from six libraries, 62.27% of which mapped to the pig reference genome. A total of 318 pig miRNAs were detected among the samples. Among them, 18 miRNAs were differentially expressed (p-value < 0.05, |log₂fold change| ≥ 1) between the high and low backfat groups; 6 were up-regulated and 12 were down-regulated. Functional enrichment of the predicted target genes of the differentially expressed miRNAs, indicated that these miRNAs were involved mainly in lipid and carbohydrate metabolism, and glycan biosynthesis and metabolism. Comprehensive analysis of the mRNA and miRNA transcriptomes revealed possible regulatory relationships for fat deposition. Negatively correlated mRNA–miRNA pairs included miR-137–PPARGC1A, miR-141–FASN, and miR-122-5p–PKM, indicating these interactions may be key regulators of fat deposition. Our findings provide important insights into miRNA expression patterns in the backfat tissue of pig and new insights into the regulatory mechanisms of fat deposition in pig.

ETS2 and microRNA-155 regulate the pathogenesis of heart failure through targeting and regulating GPR18 expression

Heart failure (HF) is a global pandemic cardiovascular disease with increasing prevalence, but the pathogenesis remains to be elucidated. The present study aimed to investigate the underlying mechanism in heart failure (HF) using bioinformatics and experimental validation. A HF-associated dataset GSE84796 was downloaded from the Gene Expression Omnibus database and differentially expressed genes (DEGs) were screened for using Bayes method in the Limma package. Kyoto Encyclopedia of Genes and Genomes pathway analysis was used to perform pathway enrichment analysis of these DEGs using The Database for Annotation, Visualization and Integrated Discovery. A protein-protein interaction (PPI) network of DEG-encoded proteins was subsequently constructed using the Search tool for the Retrieval of Interacting Genes/Proteins, and a transcription factor (TF)/miRNA-target network was constructed according to the WEB-based Gene SeT AnaLysis Tookit. The expression levels of microRNA (miRNA/miR)-155, G-protein coupled receptor 18 (GRP18) and E26 transformation-specific transcription factor 2 (ETS2) were analyzed in clinical HF samples, and functional validations were performed in H9c2 (2-1) cells. A total of 419 DEGs were identified, including 366 upregulated genes and 53 downregulated genes. The upregulated DEGs were significantly enriched in the pathways of ‘cytokine-cytokine receptor interaction’, ‘natural killer cell mediated cytotoxicity’ and ‘primary immunodeficiency’. A total of two functional modules were identified in the PPI network: Module A was enriched in 3 KEGG pathways and module B was enriched in 15 KEGG pathways. Furthermore, a total of three miRNAs and eight TFs were identified in the TF/miRNA-target network. Specifically, GPR18 was discovered to be targeted by both ETS2 and miR-155. Clinical validation revealed that the expression levels of miR-155 were significantly decreased in the HF samples, whereas the expression levels of ETS2 and GPR18 were significantly increased in HF samples. In conclusion, the present study suggested that GPR18 may be a target of ETS2 and miR-155, and miR-155 may regulate cell viability and apoptosis in H9c2 (2-1) cells through targeting and regulating GPR18.

Genome-wide analysis of circular RNAs involved in Marek's disease tumourigenesis in chickens

Marek's disease (MD), induced by Marek's disease virus (MDV), is a lymphotropic neoplastic disease and causes huge economic losses to the poultry industry. Non-coding RNAs (ncRNAs) play important regulatory roles in disease pathogenesis. To investigate host circular RNA (circRNA) and microRNA (miRNA) expression profile, RNA sequencing was performed in tumourous spleens (TS), spleens from the survivors (SS) without any lesion after MDV infection, and non-infected chicken spleens (NS). A total of 2,169 circRNAs were identified and more than 80% of circRNAs were derived from exon. The flanking introns of 1,744 exonic circRNAs possessed 579 reverse complementary matches (RCMs), which mainly overlapped with chicken repeat 1 family (CR1F). It suggested that CR1F mediated the cyclization of exons by intron pairing. Out of 2,169 circRNAs, 113 were differentially expressed circRNAs (DECs). The Q-PCR and Rnase R digestion experiments showed circRNA possessed high stability compared with their linear RNAs. Integrated with previous transcriptome data, we profiled regulatory networks of circRNA/long non-coding RNA (lncRNA)-miRNA-mRNA. Extensive competing endogenous RNA (ceRNA) networks were predicted to be involved in MD tumourigenesis. Interestingly, circZMYM3, an intronic circRNA, interacted with seven miRNAs which targeted some immune genes, such as SWAP70 and CCL4. Gga-miR-155 not only interacted with circGTDC1 and circMYO1B, but also targeted immune-related genes, such as GATA4, which indicated the roles of non-coding RNAs played to mediate immune responsive genes. Collectively, this is the first study that integrated RNA expression profiles in MD model. Our results provided comprehensive interactions of ncRNAs and mRNA in MD tumourigenesis.

Identification and molecular analysis of cashmere goat lncRNAs reveal their integrated regulatory network and potential roles in secondary hair follicle

Long non-coding RNAs (lncRNAs) is a class of eukaryotic transcripts with length of more than 200 bp. They contribute to the regulation of gene expressions involved in multiple processes including the skin cell proliferation, differentiation, and reconstruction of the secondary hair follicles (SHFs). In this study, firstly, we identified 16 putative lncRNAs from SHFs of cashmere goat based on the EST sequences from NCBI database. Secondly, we investigated their transcriptional pattern in SHFs of cashmere goat, and constructed their ceRNA regulatory networks. The RT-qPCR results showed four lncRNAs (lncRNA-475074, -052149, -052140, and -051789) were significantly up-regulated, and nine lncRNAs (lncRNA-711032, -475083, -475070, -052139, -052127, -052037, -051903, -051847, and -051804) were significantly down-regulatd in anagen SHFs of cashmere goat. CeRNA networks analysis revealed complex interactional relationship among lncRNAs, miRNAs and mRNAs. Further, the KEGG pathway enrichment was performed for the potential target genes of the identified lncRNAs based on bioinformatics technique, and the results indicated that differentially expressed lncRNAs directly or indirectly might regulate potential target genes. Our results from this study will provide a significant information for further exploring the functions and possible mechanisms of the identified lncRNAs in SHFs of cashmere goat.

Characterization of Circular RNA and microRNA Profiles in Septic Myocardial Depression: a Lipopolysaccharide-Induced Rat Septic Shock Model

Septic shock with heart dysfunction is common in intensive care units. However, the mechanism underlying myocardial depression is still unclear. Whether circular RNA (circRNA) or microRNA (miRNA) profiles differ between patients with and without myocardial depression is unknown. We generated a hypodynamic septic shock model induced by lipopolysaccharide (LPS) in adolescent rats. A total of 12 rats were utilized and heart tissue from each was collected. RNA sequencing was performed on left ventricular tissue. We focused on features of circRNAs and miRNAs, predicting their function by bioinformatic analysis and constructing circRNA-associated and miRNA-associated regulatory networks in heart tissue. We detected 851 circRNAs in heart samples, and 11 showed differential expression. A total of 639 annotated miRNAs and 91 novel miRNAs were explored including 78 showing differential expression between the two groups. We then constructed the most comprehensive circRNA-associated and miRNA-associated networks to explore their regulatory relationship in septic heart tissue, and demonstrated that different networks could potentially participate in and regulate the pathological process of sepsis. Furthermore, gene ontology term enrichment indicated miRNAs, and miRNA-mRNA networks could be associated with regulation and metabolic process, or influence cellular functions. The construction of regulator networks could improve the understanding of the basic molecular mechanisms underlying myocardial depression. It will be important for future investigations to ascertain the biological mechanisms present during the development of sepsis-induced myocardial depression to influence approaches to treatment.

[Identification of potential traumatic spinal cord injury related circular RNA-microRNA networks by sequence analysis]

Objective

To systematically profile and characterize the circular RNA (circRNA) and microRNA (miRNA) expression pattern in the lesion epicenter of spinal tissues after traumatic spinal cord injury (TSCI) and predict the structure and potential functions of the regulatory network.

Methods

Forty-eight adult male C57BL/6 mice (weighing, 18-22 g) were randomly divided into the TSCI ( n=24) and sham ( n=24) groups. Mice in the TSCI group underwent T _8-10 vertebral laminectomy and Allen's weight-drop spinal cord injury. Mice in the sham group underwent the same laminectomy without TSCI. The spinal tissues were harvested after 3 days. Some tissues were stained with HE staining to observe the structure. The others were used for sequencing. The RNA-Seq, gene ontology (GO) analysis, and circRNA-miRNA network analyses (TargetScan and miRanda) were used to profile the expression and regulation patterns of network of mice models after TSCI.

Results

HE staining showed the severe damage to the spinal cord in TSCI group compared with sham group. A total of 17 440 circRNAs and 1 228 miRNAs were identified. The host gene of significant differentially expressed circRNA enriched in the cytoplasm, associated with positive regulation of transcription and protein phosphorylation. mmu-miR-21-5p was the most significant differentially expressed miRNA after TSCI, and circRNA6730 was predicted to be its targeted circRNA. Then a potential regulatory circRNA-miRNA network was constructed.

Conclusion

The significant differentially expressed circRNAs and miRNAs may play important roles after TSCI. A targeted interaction network with mmu-miR-21-5p at the core of circRNA6730 could provide basis of pathophysiological mechanism, as well as help guide therapeutic strategies for TSCI.

Inferring Regulatory Networks From Mixed Observational Data Using Directed Acyclic Graphs

Construction of regulatory networks using cross-sectional expression profiling of genes is desired, but challenging. The Directed Acyclic Graph (DAG) provides a general framework to infer causal effects from observational data. However, most existing DAG methods assume that all nodes follow the same type of distribution, which prohibit a joint modeling of continuous gene expression and categorical variables. We present a new mixed DAG (mDAG) algorithm to infer the regulatory pathway from mixed observational data containing both continuous variables (e.g. expression of genes) and categorical variables (e.g. categorical phenotypes or single nucleotide polymorphisms). Our method can identify upstream causal factors and downstream effectors closely linked to a variable and generate hypotheses for causal direction of regulatory pathways. We propose a new permutation method to test the conditional independence of variables of mixed types, which is the key for mDAG. We also utilize an L₁ regularization in mDAG to ensure it can recover a large sparse DAG with limited sample size. We demonstrate through extensive simulations that mDAG outperforms two well-known methods in recovering the true underlying DAG. We apply mDAG to a cross-sectional immunological study of Chlamydia trachomatis infection and successfully infer the regularity network of cytokines. We also apply mDAG to a large cohort study, generating sensible mechanistic hypotheses underlying plasma adiponectin level. The R package mDAG is publicly available from CRAN at https://CRAN.R-project.org/package=mDAG.

From Modules to Networks: a Systems-Level Analysis of the Bacitracin Stress Response in Bacillus subtilis

Antibiotic resistance poses a major threat to global health, and systematic studies to understand the underlying resistance mechanisms are urgently needed. Although significant progress has been made in deciphering the mechanistic basis of individual resistance determinants, many bacterial species rely on the induction of a whole battery of resistance modules, and the complex regulatory networks controlling these modules in response to antibiotic stress are often poorly understood. In this work we combined experiments and theoretical modeling to decipher the resistance network of Bacillus subtilis against bacitracin, which inhibits cell wall biosynthesis in Gram-positive bacteria. We found a high level of cross-regulation between the two major resistance modules in response to bacitracin stress and quantified their effects on bacterial resistance. To rationalize our experimental data, we expanded a previously established computational model for the lipid II cycle through incorporating the quantitative action of the resistance modules. This led us to a systems-level description of the bacitracin stress response network that captures the complex interplay between resistance modules and the essential lipid II cycle of cell wall biosynthesis and accurately predicts the minimal inhibitory bacitracin concentration in all the studied mutants. With this, our study highlights how bacterial resistance emerges from an interlaced network of redundant homeostasis and stress response modules.

Bacterial resistance against antibiotics often involves multiple mechanisms that are interconnected to ensure robust protection. So far, the knowledge about underlying regulatory features of those resistance networks is sparse, since they can hardly be determined by experimentation alone. Here, we present the first computational approach to elucidate the interplay between multiple resistance modules against a single antibiotic and how regulatory network structure allows the cell to respond to and compensate for perturbations of resistance. Based on the response of Bacillus subtilis toward the cell wall synthesis-inhibiting antibiotic bacitracin, we developed a mathematical model that comprehensively describes the protective effect of two well-studied resistance modules (BceAB and BcrC) on the progression of the lipid II cycle. By integrating experimental measurements of expression levels, the model accurately predicts the efficacy of bacitracin against the B. subtilis wild type as well as mutant strains lacking one or both of the resistance modules. Our study reveals that bacitracin-induced changes in the properties of the lipid II cycle itself control the interplay between the two resistance modules. In particular, variations in the concentrations of UPP, the lipid II cycle intermediate that is targeted by bacitracin, connect the effect of the BceAB transporter and the homeostatic response via BcrC to an overall resistance response. We propose that monitoring changes in pathway properties caused by a stressor allows the cell to fine-tune deployment of multiple resistance systems and may serve as a cost-beneficial strategy to control the overall response toward this stressor.

IMPORTANCE Antibiotic resistance poses a major threat to global health, and systematic studies to understand the underlying resistance mechanisms are urgently needed. Although significant progress has been made in deciphering the mechanistic basis of individual resistance determinants, many bacterial species rely on the induction of a whole battery of resistance modules, and the complex regulatory networks controlling these modules in response to antibiotic stress are often poorly understood. In this work we combined experiments and theoretical modeling to decipher the resistance network of Bacillus subtilis against bacitracin, which inhibits cell wall biosynthesis in Gram-positive bacteria. We found a high level of cross-regulation between the two major resistance modules in response to bacitracin stress and quantified their effects on bacterial resistance. To rationalize our experimental data, we expanded a previously established computational model for the lipid II cycle through incorporating the quantitative action of the resistance modules. This led us to a systems-level description of the bacitracin stress response network that captures the complex interplay between resistance modules and the essential lipid II cycle of cell wall biosynthesis and accurately predicts the minimal inhibitory bacitracin concentration in all the studied mutants. With this, our study highlights how bacterial resistance emerges from an interlaced network of redundant homeostasis and stress response modules.

Function and regulation annotation of up-regulated long non-coding RNA LINC01234 in gastric cancer

Background

Accumulated evidences indicate that long non‐coding RNAs (lncRNAs) participate in many biological mechanisms. Moreover, it acts as an essential regulator in various human diseases such as gastric cancer (GC). Nevertheless, the comprehensive regulatory roles and clinical significance of most lncRNAs in GC are not fully understood.

Methods

In this research, our aim was to investigate the underlying mechanism of lncRNA LINC01234 in GC. Firstly, the usage of qRT‐PCR helped to establish expression pattern of LINC01234 in GC tissues. Following this, appropriate statistical tests were applied to analyze the relation between expression level and clinicopathological factors. Ultimately, potential functions and regulatory network of LINC01234 were concluded via GSEA and a series of bioinformatics tools or databases, respectively.

Results

Consequently, at the end of research we found LINC01234 is up‐regulated in GC tissues in comparison with adjacent normal tissues. Furthermore, its expression level is correlated with differentiation of patients with GC. It is also important to highlight bioinformatics analysis revealed that LINC01234 is involved in cancer‐associated pathways such as cell cycle and mismatch repair. Also, regulatory network of LINC01234 presented a probability in the involvement of tumorigenesis through regulating cancer‐associated genes.

Conclusion

Overall, our results suggested that LINC01234 may play a crucial role in GC.

Systematically Dissecting the Function of RNA-Binding Proteins During Glioma Progression

RNA-binding proteins (RBPs) play important roles in regulating gene expression and dysregulation of RBPs have been observed in various types of cancer. However, the role of RBPs during glioma progression, and particular in Chinese patients, is only starting to be unveiled. Here, we systematically analyzed the somatic mutation, gene expression patterns of 2949 RBPs during glioma progression. Our comprehensive study reveals several of highly mutated genes (such as ATRX, TTN and SETD2) and differentially expressed genes (such as KIF4A, TTK and CEP55). Integration of the expression of RBPs and genes, we constructed a regulatory network in glioma and revealed the functional links between RBPs and cancer-related genes. Moreover, we identified the prognosis spectrum of RBPs during glioma progression. The expression of a number of RBPs, such as SNRPN and IGF2BP3, are significantly associated with overall survival of patients in all grades. Taken together, our analyses provided a valuable RBP resource during glioma progression, and revealed several candidates that potentially contribute to development of therapeutic targets for glioma.

Identification of biomarkers and construction of a microRNA-mRNA regulatory network for ependymoma using integrated bioinformatics analysis

Ependymomas (EPNs) are one of the most common types of malignant neuroepithelial tumors. In an effort to identify potential biomarkers involved in the pathogenesis of EPN, the mRNA expression profiles of the GSE25604, GSE50161, GSE66354, GSE74195 and GSE86574 datasets, in addition to the microRNA (miRNA/miR) expression profiles of GSE42657 were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs) between EPN and normal brain tissue samples were identified using the Limma package in R and GEO2R, respectively. Functional and pathway enrichment analyses were conducted using the Database for Annotation, Visualization and Integrated Discovery. A protein-protein interaction network was constructed using the Search Tool for Retrieval of Interacting Genes database, which was visualized using Cytoscape. The targeted genes of DEMs were predicted using miRWalk2.0 and a miRNA-mRNA regulatory network was constructed. Following analysis, a total of 948 DEGs and 129 DEMs were identified. Functional enrichment analysis revealed that 609 upregulated DEGs were significantly enriched in ‘PI3K-Akt signaling pathway’, while 339 downregulated DEGs were primarily involved in ‘cell junction’ and ‘retrograde endocannabinoid signaling’. In addition, 6 hub genes [cyclin dependent kinase 1, CD44 molecule (Indian blood group) (CD44), proliferating cell nuclear antigen (PCNA), MYC, synaptotagmin 1 (SYT1) and kinesin family member 4A] and 6 crucial miRNAs [homo sapiens (hsa)-miR-34a-5p, hsa-miR-449a, hsa-miR-106a-5p, hsa-miR-124-3p, hsa-miR-128-3p and hsa-miR-330-3p] were identified as biomarkers and potential therapeutic targets for EPN. Furthermore, a microRNA-mRNA regulatory network was constructed to highlight the interactions between DEMs and their target DEGs; this included the hsa-miR-449a-SYT1, hsa-miR-34a-5p-SYT1, hsa-miR-330-3p-CD44 and hsa-miR-124-3p-PCNA pairs, whose expression levels were confirmed using reverse transcription-quantitative polymerase chain reaction. In conclusion, the present study may provide important data for the investigation of the molecular mechanisms of EPN pathogenesis.

Global Analysis of Enhancer Targets Reveals Convergent Enhancer-Driven Regulatory Modules

Single-cell screens enable high-throughput functional assessment of enhancers in their endogenous genomic context. However, the design of current studies limits their application to identifying the primary gene targets of enhancers. Here, we improve the experimental and computational parameters of single-cell enhancer screens to identify the secondary gene targets of enhancers. Our analysis of >500 putative enhancers in K562 cells reveals an interwoven enhancer-driven gene regulatory network. We find that enhancers from distinct genomic loci converge to modulate the expression of common sub-modules, including the α- and β-globin loci, by directly regulating transcription factors. Our analysis suggests that several genetic variants associated with myeloid blood cell traits alter the activity of a distal enhancer of MYB (~140 kb away), with downstream consequences on hemoglobin genes expression and cell state. These data have implications for the understanding of enhancer-associated traits and emphasize the flexibility of controlling transcriptional systems by modifying enhancer activity.

Xie et al. apply improved strategies for single-cell screens to identify an enhancer-driven transcriptional regulatory network in K562 cells. They demonstrate that the same group of genes can be indirectly regulated by enhancers from distinct genomic loci. These data have implications for the understanding of enhancer-associated traits.