MicroRNA-mediated feedback and feedforward loops are recurrent network motifs in mammals

Current updates regarding biogenesis, functions and dysregulation of microRNAs in cancer: Innovative approaches for detection using CRISPR/Cas13‑based platforms (Review)

MicroRNAs (miRNAs) are short non‑coding RNAs, which perform a key role in cellular differentiation and development. Most human diseases, particularly cancer, are linked to miRNA functional dysregulation implicated in the expression of tumor‑suppressive or oncogenic targets. Cancer hallmarks such as continued proliferative signaling, dodging growth suppressors, invasion and metastasis, triggering angiogenesis, and avoiding cell death have all been demonstrated to be affected by dysregulated miRNAs. Thus, for the treatment of different cancer types, the detection and quantification of this type of RNA is significant. The classical and current methods of RNA detection, including northern blotting, reverse transcription‑quantitative PCR, rolling circle amplification and next‑generation sequencing, may be effective but differ in efficiency and accuracy. Furthermore, these approaches are expensive, and require special instrumentation and expertise. Thus, researchers are constantly looking for more innovative approaches for miRNA detection, which can be advantageous in all aspects. In this regard, an RNA manipulation tool known as the CRISPR and CRISPR‑associated sequence 13 (CRISPR/Cas13) system has been found to be more advantageous in miRNA detection. The Cas13‑based miRNA detection approach is cost effective and requires no special instrumentation or expertise. However, more research and validation are required to confirm the growing body of CRISPR/Cas13‑based research that has identified miRNAs as possible cancer biomarkers for diagnosis and prognosis, and as targets for treatment. In the present review, current updates regarding miRNA biogenesis, structural and functional aspects, and miRNA dysregulation during cancer are described. In addition, novel approaches using the CRISPR/Cas13 system as a next‑generation tool for miRNA detection are discussed. Furthermore, challenges and prospects of CRISPR/Cas13‑based miRNA detection approaches are described.

MicroRNA-210 Enhances Cell Survival and Paracrine Potential for Cardiac Cell Therapy While Targeting Mitophagy

The therapeutic potential of presumed cardiac progenitor cells (CPCs) in heart regeneration has garnered significant interest, yet clinical trials have revealed limited efficacy due to challenges in cell survival, retention, and expansion. Priming CPCs to survive the hostile hypoxic environment may be key to enhancing their regenerative capacity. We demonstrate that microRNA-210 (miR-210), known for its role in hypoxic adaptation, significantly improves CPC survival by inhibiting apoptosis through the downregulation of Casp8ap2, a ~40% reduction in caspase activity, and a ~90% decrease in DNA fragmentation. Contrary to the expected induction of Bnip3-dependent mitophagy by hypoxia, miR-210 did not upregulate Bnip3, indicating a distinct anti-apoptotic mechanism. Instead, miR-210 reduced markers of mitophagy and increased mitochondrial biogenesis and oxidative metabolism, suggesting a role in metabolic reprogramming. Furthermore, miR-210 enhanced the secretion of paracrine growth factors from CPCs, with a ~1.6-fold increase in the release of stem cell factor and of insulin growth factor 1, which promoted in vitro endothelial cell proliferation and cardiomyocyte survival. These findings elucidate the multifaceted role of miR-210 in CPC biology and its potential to enhance cell-based therapies for myocardial repair by promoting cell survival, metabolic adaptation, and paracrine signalling.

The Abnormal ERα-miRNA Cross-Talk in AD-Affected Human Hippocampus: A Bioinformatics Perspective

Estrogen's impact on Alzheimer's disease (AD) is multifaceted, with its receptors potentially influencing AD pathology in both beneficial and detrimental ways. This study aims to dissect the intricate cross-talk between estrogen receptor alpha (ERα) and microRNAs (miRNAs) in AD-affected human hippocampus. Through a comprehensive literature review in the PubMed database, coupled with a GeneCards database search, we obtained AD-related key miRNAs and genes in the hippocampus. Using bioinformatics tools, we predicted target genes and miRNAs of ERα, and the targets of the identified miRNAs. The integration of these elements resulted in the construction of an ERα-related FFL network, which includes 13 miRNAs and 56 core genes. Gene ontology (GO) and pathway enrichment analyses were conducted, revealing significant enrichment in biological processes such as neuron death and response to metal ions, and cellular components like membrane microdomains. Notably, the AKT-associated signaling pathway was prominently highlighted, with key genes including GSK3A, CDKN1A, AKT2, and MDM2, and key miRNAs including miR-485 and let-7f, suggesting a potential role of ERα in modulating this pathway in AD. The findings of this study provide a novel perspective on the regulatory network of ERα in the hippocampal region of AD and may pave the way for future research into the therapeutic potential of targeting the ERα pathway in neurodegenerative diseases.

Investigating Gene Expression Noise Reduction by MicroRNAs and MiRISC Reinforcement by Self-Feedback Regulation of mRNA Degradation

MicroRNA (miRNA) induced silencing complex (miRISC) is the targeting apparatus and arguably rate-limiting step of the miRNA-mediated regulatory subsystem - the major noise reducing though metabolically wasteful mechanism. Recently, we reported that miRISC channels miRNA-mediated regulatory activity back onto their own mRNAs to form negative self-feedback loops, a noise-reduction technique in engineering and synthetic/systems biology. Here, we describe mathematical modeling that predicts mRNA expression noise to correlate negatively with degradation rate (K deg ) and noise reduction by self-feedback control ofK deg . We also calculatedK deg and expression noise of mRNAs detected in a cutting-edge total-RNA single-cell RNA-seq (scRNA-seq) dataset. As predicted, miRNA-targeted mRNAs exhibited higherK deg values in conjunction with lower inter-cell expression noise. Moreover, as predicted by our self-feedback loop model, miRISC mRNAs (AGO1/2/3 and TNRC6A/B/C) exhibited further reduced expression noise. In short, mathematical-modeling and total-RNA scRNA-seq data-analysis shed insight into operational trade-off between noise reduction and metabolic/energetic expenditure in producing miRNA-targeted mRNAs destined for enhanced degradation and translational inhibition, as well as negative self-feedback loop reinforcement of miRISC - the core of miRNA-mediated noise-reduction subsystem. To our knowledge, this is the first report of concurrent mRNA degradation and expression noise analyses and of noise reduction by self-feedback control of mRNA degradation.

Identification of miRNA-TF Regulatory Pathways Related to Diseases from a Neuroendocrine-Immune Perspective

The neuroendocrine-immune (NEI) network is fundamental for maintaining body's homeostasis and health. While the roles of microRNAs (miRNAs) and transcription factors (TFs) in disease processes are well-established, their synergistic regulation within the NEI network has yet to be elucidated. In this study, we constructed a background NEI-related miRNA-TF regulatory network (NEI-miRTF-N) by integrating NEI signaling molecules (including miRNAs, genes, and TFs) and identifying miRNA-TF feed-forward loops. Our analysis reveals that the number of immune signaling molecules is the highest and suggests potential directions for signal transduction, primarily from the nervous system to both the endocrine and immune systems, as well as from the endocrine system to the immune system. Furthermore, disease-specific NEI-miRTF-Ns for depression, Alzheimer's disease (AD) and dilated cardiomyopathy (DCM) were constructed based on the known disease molecules and significantly differentially expressed (SDE) molecules. Additionally, we proposed a novel method using depth-first-search algorithm for identifying significantly dysregulated NEI-related miRNA-TF regulatory pathways (NEI-miRTF-Ps) and verified their reliability from multiple perspectives. Our study provides an effective approach for identifying disease-specific NEI-miRTF-Ps and offers new insights into the synergistic regulation of miRNAs and TFs within the NEI network. Our findings provide information for new therapeutic strategies targeting these regulatory pathways.

Advantage looping: Gene regulatory circuits between microRNAs and their target transcription factors in plants

The 20 to 24 nucleotide microRNAs (miRNAs) and their target transcription factors (TF) have emerged as key regulators of diverse processes in plants, including organ development and environmental resilience. In several instances, the mature miRNAs degrade the TF-encoding transcripts, while their protein products in turn bind to the promoters of the respective miRNA-encoding genes and regulate their expression, thus forming feedback loops (FBLs) or feedforward loops (FFLs). Computational analysis suggested that such miRNA-TF loops are recurrent motifs in gene regulatory networks (GRNs) in plants as well as animals. In recent years, modeling and experimental studies have suggested that plant miRNA-TF loops in GRNs play critical roles in driving organ development and abiotic stress responses. Here, we discuss the miRNA-TF FBLs and FFLs that have been identified and studied in plants over the past decade. We then provide some insights into the possible roles of such motifs within GRNs. Lastly, we provide perspectives on future directions for dissecting the functions of miRNA-centric GRNs in plants.

MicroRNA-mediated network redundancy is constrained by purifying selection and contributes to expression robustness in Drosophila melanogaster

MicroRNAs (miRNAs) are post-transcriptional, non-coding regulatory RNAs that function coordinately with transcription factors (TFs) in gene regulatory networks. TFs and their targets are often co-regulated by miRNAs, forming composite feedforward circuits (cFFCs) with varying degrees of redundancy, primarily mediated by miRNAs. However, the maintenance of miRNA-mediated regulatory redundancy and its impact on gene expression evolution remain elusive. By integrating ChIP-seq data from ENCODE and miRNA targeting from TargetScanFly, we quantified miRNA-mediated cFFC redundancy in Drosophila melanogaster embryos and larvae, revealing more than three quarters of miRNA targets are involved in redundant cFFCs. Higher cFFC redundancy, where more miRNAs target the same gene within a cFFC, is correlated with stronger purifying selection, reduced expression divergence between species, and increased expression stability under heat shock stress. Redundant cFFCs primarily regulate older or broadly expressed young genes. These findings highlight the role of miRNA-mediated cFFC redundancy in enhancing gene expression robustness through natural selection.

MicroRNA-mediated Krüppel-like factor 4 upregulation induces alternatively activated macrophage-associated marker and chemokine transcription in 4,4'-methylene diphenyl diisocyanate exposed macrophages

1. Occupational exposure to 4,4'-methylene diphenyl diisocyanate (MDI) is associated with occupational asthma (OA) development. Alveolar macrophage-induced recruitment of immune cells to the lung microenvironment plays an important role during asthma pathogenesis. Previous studies identified that MDI/MDI-glutathione (GSH)-exposure downregulates endogenous hsa-miR-206-3p/hsa-miR-381-3p. Our prior report shows that alternatively activated (M2) macrophage-associated markers/chemokines are induced by MDI/MDI-GSH-mediated Krüppel-Like Factor 4 (KLF4) upregulation in macrophages and stimulates immune cell chemotaxis. However, the underlying molecular mechanism(s) by which MDI/MDI-GSH upregulates KLF4 remain unclear.

2. Following MDI-GSH exposure, microRNA(miR)-inhibitors/mimics or plasmid transfection, endogenous hsa-miR-206-3p/hsa-miR-381-3p, KLF4, or M2 macrophage-associated markers (CD206, TGM2), and chemokines (CCL17, CCL22, CCL24) were measured by either RT-qPCR, western blot, or luciferase assay.

3. MDI-GSH exposure downregulates hsa-miR-206-3p/hsa-miR-381-3p by 1.46- to 9.75-fold whereas upregulates KLF4 by 1.68- to 1.99-fold, respectively. In silico analysis predicts binding between hsa-miR-206-3p/hsa-miR-381-3p and KLF4. Gain- and loss-of-function, luciferase reporter assays and RNA-induced silencing complex-immunoprecipitation (RISC-IP) studies confirm the posttranscriptional regulatory roles of hsa-miR-206-3p/hsa-miR-381-3p and KLF4 in macrophages. Furthermore, hsa-miR-206-3p/hsa-miR-381-3p regulate the expression of M2 macrophage-associated markers and chemokines via KLF4.

4. In conclusion, hsa-miR-206-3p/hsa-miR-381-3p play a major role in regulation of MDI/MDI-GSH-induced M2 macrophage-associated markers and chemokines by targeting the KLF4 transcript, and KLF4-mediated regulation in macrophages.

Exploring miRNA‑target gene profiles associated with drug resistance in patients with breast cancer receiving neoadjuvant chemotherapy

Exosomal microRNAs (miRNAs) are closely related to drug resistance in patients with breast cancer (BC); however, only a few roles of the exosomal miRNA-target gene networks have been clinically implicated in drug resistance in BC. Therefore, the present study aimed to identify the differential expression of exosomal miRNAs associated with drug resistance and their target mRNAs. In vitro microarray analysis was used to verify differentially expressed miRNAs (DEMs) in drug-resistant BC. Next, tumor-derived exosomes (TDEs) were isolated. Furthermore, it was determined whether the candidate drug-resistant miRNAs were also significant in TDEs, and then putative miRNAs in TDEs were validated in plasma samples from 35 patients with BC (20 patients with BC showing no response and 15 patients with BC showing a complete response). It was confirmed that the combination of five exosomal miRNAs, including miR-125b-5p, miR-146a-5p, miR-484, miR-1246-5p and miR-1260b, was effective for predicting therapeutic response to neoadjuvant chemotherapy, with an area under the curve value of 0.95, sensitivity of 75%, and specificity of 95%. Public datasets were analyzed to identify differentially expressed genes (DEGs) related to drug resistance and it was revealed that BAK1, NOVA1, PTGER4, RTKN2, AGO1, CAP1, and ETS1 were the target genes of exosomal miRNAs. Networks between DEMs and DEGs were highly correlated with mitosis, metabolism, drug transport, and immune responses. Consequently, these targets could be used as predictive markers and therapeutic targets for clinical applications to enhance treatment outcomes for patients with BC.

Synthetic dosage-compensating miRNA circuits allow precision gene therapy for Rett syndrome

A longstanding challenge in gene therapy is expressing a dosage-sensitive gene within a tight therapeutic window. For example, loss of MECP2 function causes Rett syndrome, while its duplication causes MECP2 duplication syndrome. Viral gene delivery methods generate variable numbers of gene copies in individual cells, creating a need for gene dosage-invariant expression systems. Here, we introduce a compact miRNA-based, incoherent feed-forward loop circuit that achieves precise control of Mecp2 expression in cells and brains, and improves outcomes in an AAV-based mouse model of Rett syndrome gene therapy. Single molecule analysis of endogenous and ectopic Mecp2 mRNA revealed precise, sustained expression across a broad range of gene dosages. Delivered systemically in a brain-targeting AAV capsid, the circuit strongly suppressed Rett behavioral symptoms for over 24 weeks, outperforming an unregulated gene therapy. These results demonstrate that synthetic miRNA-based regulatory circuits can enable precise in vivo expression to improve the safety and efficacy of gene therapy.

Self-feedback loop-containing synthetic mRNA switches for controlled microRNA sensing

Synthetic mRNA switches are powerful cell fate manipulation tools that sense cellular input molecules to directly control protein expression at the translational level. The lack of available switch designs that can mimic the natural sophisticated protein regulation is a fundamental issue that limits the application of synthetic mRNA switches. Here we report a new set of synthetic mRNA switches by incorporating self-feedback loop machineries to dynamically control protein expression levels upon sensing cellular microRNAs. We redesigned the coding region of the switch to express output protein along with mRNA regulatory proteins. RNA-binding proteins (RBPs) and RBP-binding RNA motifs (aptamers) guide the regulatory proteins to act on their own mRNAs, enhancing or flattening the effect of microRNA sensing. Importantly, we demonstrated that the switches with the positive feedback feature can enlarge a high-or-low microRNA effect into a nearly all-or-none pattern, substantially boosting the use of synthetic mRNA switches as high-performance microRNA sensors or binary cell regulation tools. We believe these novel mRNA switch designs provide new strategies to construct complex mRNA-based genetic circuits for future molecular sensing and cell engineering.

Tiny Guides, Big Impact: Focus on the Opportunities and Challenges of miR-Based Treatments for ARDS

Acute Respiratory Distress Syndrome (ARDS) is characterized by lung inflammation and increased membrane permeability, which represents the leading cause of mortality in ICUs. Mechanical ventilation strategies are at the forefront of supportive approaches for ARDS. Recently, an increasing understanding of RNA biology, function, and regulation, as well as the success of RNA vaccines, has spurred enthusiasm for the emergence of novel RNA-based therapeutics. The most common types of RNA seen in development are silencing (si)RNAs, antisense oligonucleotide therapy (ASO), and messenger (m)RNAs that collectively account for 80% of the RNA therapeutics pipeline. These three RNA platforms are the most mature, with approved products and demonstrated commercial success. Most recently, miRNAs have emerged as pivotal regulators of gene expression. Their dysregulation in various clinical conditions offers insights into ARDS pathogenesis and offers the innovative possibility of using microRNAs as targeted therapy. This review synthesizes the current state of the literature to contextualize the therapeutic potential of miRNA modulation. It considers the potential for miR-based therapeutics as a nuanced approach that incorporates the complexity of ARDS pathophysiology and the multifaceted nature of miRNA interactions.

Bioenergetic costs and the evolution of noise regulation by microRNAs

Noise control, together with other regulatory functions facilitated by microRNAs (miRNAs), is believed to have played important roles in the evolution of multicellular eukaryotic organisms. miRNAs can dampen protein fluctuations via enhanced degradation of messenger RNA (mRNA), but this requires compensation by increased mRNA transcription to maintain the same expression levels. The overall mechanism is metabolically expensive, leading to questions about how it might have evolved in the first place. We develop a stochastic model of miRNA noise regulation, coupled with a detailed analysis of the associated metabolic costs. Additionally, we calculate binding free energies for a range of miRNA seeds, the short sequences which govern target recognition. We argue that natural selection may have fine-tuned the Michaelis-Menten constant [Formula: see text] describing miRNA-mRNA affinity and show supporting evidence from analysis of experimental data. [Formula: see text] is constrained by seed length, and optimal noise control (minimum protein variance at a given energy cost) is achievable for seeds of 6 to 7 nucleotides in length, the most commonly observed types. Moreover, at optimality, the degree of noise reduction approaches the theoretical bound set by the Wiener-Kolmogorov linear filter. The results illustrate how selective pressure toward energy efficiency has potentially shaped a crucial regulatory pathway in eukaryotes.

Orthogonal inducible control of Cas13 circuits enables programmable RNA regulation in mammalian cells

RNA plays an indispensable role in mammalian cell functions. Cas13, a class of RNA-guided ribonuclease, is a flexible tool for modifying and regulating coding and non-coding RNAs, with enormous potential for creating new cell functions. However, the lack of control over Cas13 activity has limited its cell engineering capability. Here, we present the CRISTAL (Control of RNA with Inducible SpliT CAs13 Orthologs and Exogenous Ligands) platform. CRISTAL is powered by a collection (10 total) of orthogonal split inducible Cas13 effectors that can be turned ON or OFF via small molecules in multiple cell types, providing precise temporal control. Also, we engineer Cas13 logic circuits that can respond to endogenous signaling and exogenous small molecule inputs. Furthermore, the orthogonality, low leakiness, and high dynamic range of our inducible Cas13d and Cas13b enable the design and construction of a robust incoherent feedforward loop, leading to near-perfect and tunable adaptation response. Finally, using our inducible Cas13 effectors, we achieve simultaneous multiplexed control of multiple genes in vitro and in mice. Together, our CRISTAL design represents a powerful platform for precisely regulating RNA dynamics to advance cell engineering and elucidate RNA biology.

Gene dosage compensation: Origins, criteria to identify compensated genes, and mechanisms including sensor loops as an emerging systems-level property in cancer

The gene dosage compensation hypothesis presents a mechanism through which the expression of certain genes is modulated to compensate for differences in the dose of genes when additional chromosomes are present. It is one of the means through which cancer cells actively cope with the potential damaging effects of aneuploidy, a hallmark of most cancers. Dosage compensation arises through several processes, including downregulation or overexpression of specific genes and the relocation of dosage-sensitive genes. In cancer, a majority of compensated genes are generally thought to be regulated at the translational or post-translational level, and include the basic components of a compensation loop, including sensors of gene dosage and modulators of gene expression. Post-translational regulation is mostly undertaken by a general degradation or aggregation of remaining protein subunits of macromolecular complexes. An increasingly important role has also been observed for transcriptional level regulation. This article reviews the process of targeted gene dosage compensation in cancer and other biological conditions, along with the mechanisms by which cells regulate specific genes to restore cellular homeostasis. These mechanisms represent potential targets for the inhibition of dosage compensation of specific genes in aneuploid cancers. This article critically examines the process of targeted gene dosage compensation in cancer and other biological contexts, alongside the criteria for identifying genes subject to dosage compensation and the intricate mechanisms by which cells orchestrate the regulation of specific genes to reinstate cellular homeostasis. Ultimately, our aim is to gain a comprehensive understanding of the intricate nature of a systems-level property. This property hinges upon the kinetic parameters of regulatory motifs, which we have termed "gene dosage sensor loops." These loops have the potential to operate at both the transcriptional and translational levels, thus emerging as promising candidates for the inhibition of dosage compensation in specific genes. Additionally, they represent novel and highly specific therapeutic targets in the context of aneuploid cancer.

Effects of microRNA-mediated negative feedback on gene expression noise

MicroRNAs (miRNAs) are small noncoding RNAs that regulate gene expression post-transcriptionally in eukaryotes by binding with target mRNAs and preventing translation. miRNA-mediated feedback motifs are ubiquitous in various genetic networks that control cellular decision making. A key question is how such a feedback mechanism may affect gene expression noise. To answer this, we have developed a mathematical model to study the effects of a miRNA-dependent negative-feedback loop on mean expression and noise in target mRNAs. Combining analytics and simulations, we show the existence of an expression threshold demarcating repressed and expressed regimes in agreement with earlier studies. The steady-state mRNA distributions are bimodal near the threshold, where copy numbers of mRNAs and miRNAs exhibit enhanced anticorrelated fluctuations. Moreover, variation of negative-feedback strength shifts the threshold locations and modulates the noise profiles. Notably, the miRNA-mRNA binding affinity and feedback strength collectively shape the bimodality. We also compare our model with a direct auto-repression motif, where a gene produces its own repressor. Auto-repression fails to produce bimodal mRNA distributions as found in miRNA-based indirect repression, suggesting the crucial role of miRNAs in creating phenotypic diversity. Together, we demonstrate how miRNA-dependent negative feedback modifies the expression threshold and leads to a broader parameter regime of bimodality compared to the no-feedback case.

Signaling-induced systematic repression of miRNAs uncovers cancer vulnerabilities and targeted therapy sensitivity

Targeted therapies are effective in treating cancer, but success depends on identifying cancer vulnerabilities. In our study, we utilize small RNA sequencing to examine the impact of pathway activation on microRNA (miRNA) expression patterns. Interestingly, we discover that miRNAs capable of inhibiting key members of activated pathways are frequently diminished. Building on this observation, we develop an approach that integrates a low-miRNA-expression signature to identify druggable target genes in cancer. We train and validate our approach in colorectal cancer cells and extend it to diverse cancer models using patient-derived in vitro and in vivo systems. Finally, we demonstrate its additional value to support genomic and transcriptomic-based drug prediction strategies in a pan-cancer patient cohort from the National Center for Tumor Diseases (NCT)/German Cancer Consortium (DKTK) Molecularly Aided Stratification for Tumor Eradication (MASTER) precision oncology trial. In conclusion, our strategy can predict cancer vulnerabilities with high sensitivity and accuracy and might be suitable for future therapy recommendations in a variety of cancer subtypes.

Role of microRNAs and their downstream target transcription factors in zebrafish thrombopoiesis

Previous studies have shown that human platelets and megakaryocytes carry microRNAs suggesting their role in platelet function and megakaryocyte development, respectively. However, a comprehensive study on the microRNAs and their targets has not been undertaken. Zebrafish thrombocytes could be used as a model to study their role in megakaryocyte maturation and platelet function because thrombocytes have both megakaryocyte features and platelet properties. In our laboratory, we identified 15 microRNAs in thrombocytes using single-cell RNA sequencing. We knocked down each of these 15 microRNAs by the piggyback method and found knockdown of three microRNAs, mir-7148, let-7b, and mir-223 in adult zebrafish led to an increase in the percentage of thrombocytes. Functional thrombocyte analysis using plate tilt assay showed no modulatory effect of the three microRNAs on thrombocyte aggregation/agglutination. We also found enhanced thrombosis using arterial laser thrombosis assay in a group of zebrafish larvae after mir-7148, let-7b, and mir-223 knockdowns. These results suggested mir-7148, let-7b, and mir-223 are repressors for thrombocyte production. We then explored miRWalk database for let-7b downstream targets and then selected those that are expressed in thrombocytes, and from this list based on their role in differentiation selected 14 genes, rorca, tgif1, rfx1a, deaf1, zbtb18, mafba, cebpa, spi1a, spi1b, fhl3b, ikzf1, irf5, irf8, and lbx1b that encode transcriptional regulators. The qRT-PCR analysis of expression levels of the above genes following let-7b knockdown showed changes in the expression of 13 targets. We then studied the effect of the 13 targets on thrombocyte production and identified 5 genes, irf5, tgif1, irf8, cebpa, and rorca that showed thrombocytosis and one gene, ikzf1 that showed thrombocytopenia. Furthermore, we tested whether mir-223 regulates any of the above 13 transcription factors after mir-223 knockdown using qRT-PCR. Six of the 13 genes showed similar gene expression as observed with let-7b knockdown and 7 genes showed opposing results. Thus, our results suggested a possible regulatory network in common with both let-7b and mir-223. We also identified that tgif1, cebpa, ikzf1, irf5, irf8, and ikzf1 play a role in thrombopoiesis. Since the ikzf1 gene showed a differential expression profile in let-7b and mir-223 knockdowns but resulted in thrombocytopenia in ikzf1 knockdown in both adults and larvae we also studied an ikzf1 mutant and showed the mutant had thrombocytopenia. Taken together, these studies showed that thrombopoiesis is controlled by a network of transcription regulators that are regulated by multiple microRNAs in both positive and negative manner resulting in overall inhibition of thrombopoiesis.

A double-negative feedback loop between miR319c and JAW-TCPs establishes growth pattern in incipient leaf primordia in Arabidopsis thaliana

The microRNA miR319 and its target JAW-TCP transcription factors regulate the proliferation-to-differentiation transition of leaf pavement cells in diverse plant species. In young Arabidopsis leaf primordia, JAW-TCPs are detected towards the distal region whereas the major mRNA319-encoding gene MIR319C, is expressed at the base. Little is known about how this complementary expression pattern of MIR319C and JAW-TCPs is generated. Here, we show that MIR319C is initially expressed uniformly throughout the incipient primordia and is later abruptly down-regulated at the distal region, with concomitant distal appearance of JAW-TCPs, when leaves grow to ~100 μm long. Loss of JAW-TCPs causes distal extension of the MIR319C expression domain, whereas ectopic TCP activity restricts MIR319C more proximally. JAW-TCPs are recruited to and are capable of depositing histone H3K27me3 repressive marks on the MIR319C chromatin. JAW-TCPs fail to repress MIR319C in transgenic seedlings where the TCP-binding cis-elements on MIR319C are mutated, causing miR319 gain-of-function-like phenotype in the embryonic leaves. Based on these results, we propose a model for growth patterning in leaf primordia wherein MIR319C and JAW-TCPs repress each other and divide the uniformly growing primordia into distal differentiation zone and proximal proliferation domain.

Role of MicroRNAs and their Downstream Target Transcription Factors in Zebrafish Thrombopoiesis

Previous studies have shown that human platelets and megakaryocytes carry microRNAs suggesting their role in platelet function and megakaryocyte development, respectively. However, a comprehensive study on the microRNAs and their targets has not been undertaken. Zebrafish thrombocytes could be used as a model to study their role in megakaryocyte maturation and platelet function because thrombocytes have both megakaryocyte features and platelet properties. In our laboratory, we identified 15 microRNAs in thrombocytes using single-cell RNA sequencing. We knocked down each of these 15 microRNAs by the piggyback method and found knockdown of three microRNAs, mir-7148, let-7b , and mir-223 in adult zebrafish led to an increase in the percentage of thrombocytes. Functional thrombocyte analysis using plate tilt assay showed no modulatory effect of the three microRNAs on thrombocyte aggregation/agglutination. We also found enhanced thrombosis using arterial laser thrombosis assay in a group of zebrafish larvae after mir-7148, let-7b , and mir-223 knockdowns. These results suggested mir-7148, let-7b , and mir-223 are repressors for thrombocyte production. We then explored miRWalk database for let-7b downstream targets and then selected those that are expressed in thrombocytes, and from this list based on their role in differentiation selected 14 genes, rorca, tgif1, rfx1a, deaf1, zbtb18, mafba, cebpa, spi1a, spi1b, fhl3b, ikzf1, irf5, irf8 , and lbx1b that encode transcriptional regulators. The qRT-PCR analysis of expression levels of the above genes following let-7b knockdown showed changes in the expression of 13 targets. We then studied the effect of the 13 targets on thrombocyte production and identified 5 genes, irf5, tgif1, irf8, cebpa , and rorca that showed thrombocytosis and one gene, ikzf1 that showed thrombocytopenia. Furthermore, we tested whether mir-223 regulates any of the above 13 transcription factors after mir-223 knockdown using qRT-PCR. Six of the 13 genes showed similar gene expression as observed with let-7b knockdown and 7 genes showed opposing results. Thus, our results suggested a possible regulatory network in common with both let-7b and mir-223 . We also identified that tgif1, cebpa, ikzf1, irf5 , irf8 , and ikzf1 play a role in thrombopoiesis. Since the ikzf1 gene showed a differential expression profile in let-7b and mir-223 knockdowns but resulted in thrombocytopenia in ikzf1 knockdown in both adults and larvae we also studied an ikzf1 mutant and showed the mutant had thrombocytopenia. Taken together, these studies showed that thrombopoiesis is controlled by a network of transcription regulators that are regulated by multiple microRNAs in both positive and negative manner resulting in overall inhibition of thrombopoiesis.

hsa-mir-(4328, 4422, 548z and -628-5p) in diabetic retinopathy: diagnosis, prediction and linking a new therapeutic target

AIMS

Growing evidence suggests that microRNAs (miRNAs) are crucial in controlling how diabetic retinopathy (DR) develops. We intend to mine miRNAs with diagnostic and predictive value for DR and to investigate new drug therapeutic targets.

METHODS

After performing a differential analysis on the miRNA and mRNA datasets for DR and neovascularization (NEO), miRNA-mRNA networks were created. Combine the results of enrichment analysis, Protein-Protein Interaction Networks (PPI), and Cytoscape to identify key miRNAs. DrugBank was used to find drugs that interacted with transcription factors (TF) predicted by TransmiR. Finally, whole blood and clinical data were collected from 58 patients with type 2 diabetes mellitus (T2DM), and RT-qPCR, logistic analysis, and ROC were used to verify the value of key miRNAs.

RESULTS

Differential analysis indicated the presence of genes and miRNAs that co-regulate DR and NEO. Enrichment analysis showed that key genes are inextricably linked to neovascularization. Combining the results of PPI and Cytoscape identified four key miRNAs, namely hsa-mir-(4328, 4422, 548z and -628-5p). RT-qPCR, logistic, and ROC results showed that decreased expression levels of hsa-mir-(4328, 4422, 548z and -628-5p) signal the risk of evolution to DR in T2DM patients. Finally, we constructed a TF-miRNA network to find the 15 TFs and the 35 drugs that interact with these TFs.

CONCLUSION

hsa-mir-(4328, 4422, 548z and -628-5p) in whole blood are protective factors for DR as novel biomarkers for diagnosis and prediction. In addition, our research provides new drug directions for the treatment of DR, such as Diosmin, Atorvastatin, and so on.

Analyzing the Association of Visceral Adipose Tissue Growth Differentiation Factor-15 and MicroRNA in Type 2 Diabetes Mellitus

Background

Growth differentiation factor-15 (GDF-15) is involved in insulin resistance and diabetes. In this study, we determine the associations of GDF-15 with miR-181b-5p, miR-330-3p, mothers against decapentaplegic homolog 7 (SMAD7), and insulin resistance in visceral adipose tissue (VAT) and peripheral blood mononuclear cells (PBMCs) in type 2 diabetes mellitus (T2DM) patients.

Methods

Sixty patients, equally divided into those with T2DM and non-diabetic controls, were recruited for gene expression analysis. Protein-protein interaction (STRING), target prediction (miRNet), and functional enrichment were conducted accordingly.

Results

Our study showed that VAT and PBMCs had similar expression profiles, where GDF-15 and miR-181b-5p were upregulated, whereas SMAD7 and miR-330-3p were downregulated. Serum GDF-15 could differentiate between T2DM and non-diabetic patients (P<0.001). Target prediction revealed a microRNA (miRNA)-messenger RNA regulatory network, transcription factors, and functional enrichment for the miRNA that suggested involvement in T2DM pathogenesis.

Conclusion

VAT GDF-15 is associated with insulin resistance and is possibly regulated by miR-181b-5p, miR-330-3p, and SMAD7 in T2DM.

Orthogonal inducible control of Cas13 circuits enables programmable RNA regulation in mammalian cells

RNA plays an indispensable role in mammalian cell functions. Cas13, a class of RNA-guided ribonuclease, is a flexible tool for modifying and regulating coding and non-coding RNAs, with enormous potential for creating new cell functions. However, the lack of control over Cas13 activity has limited its cell engineering capability. Here, we present the CRISTAL ( C ontrol of R NA with Inducible S pli T C A s13 Orthologs and Exogenous L igands) platform. CRISTAL is powered by a collection (10 total) of orthogonal split inducible Cas13s that can be turned ON or OFF via small molecules in multiple cell types, providing precise temporal control. Also, we engineered Cas13 logic circuits that can respond to endogenous signaling and exogenous small molecule inputs. Furthermore, the orthogonality, low leakiness, and high dynamic range of our inducible Cas13d and Cas13b enable the design and construction of a robust incoherent feedforward loop, leading to near-perfect and tunable adaptation response. Finally, using our inducible Cas13s, we achieve simultaneous multiplexed control of multiple genes in vitro and in mice. Together, our CRISTAL design represents a powerful platform for precisely regulating RNA dynamics to advance cell engineering and elucidate RNA biology.

Mechanisms of robustness in gene regulatory networks involved in neural development

The functions of living organisms are affected by different kinds of perturbation, both internal and external, which in many cases have functional effects and phenotypic impact. The effects of these perturbations become particularly relevant for multicellular organisms with complex body patterns and cell type heterogeneity, where transcriptional programs controlled by gene regulatory networks determine, for example, the cell fate during embryonic development. Therefore, an essential aspect of development in these organisms is the ability to maintain the functionality of their genetic developmental programs even in the presence of genetic variation, changing environmental conditions and biochemical noise, a property commonly termed robustness. We discuss the implication of different molecular mechanisms of robustness involved in neurodevelopment, which is characterized by the interplay of many developmental programs at a molecular, cellular and systemic level. We specifically focus on processes affecting the function of gene regulatory networks, encompassing transcriptional regulatory elements and post-transcriptional processes such as miRNA-based regulation, but also higher order regulatory organization, such as gene network topology. We also present cases where impairment of robustness mechanisms can be associated with neurodevelopmental disorders, as well as reasons why understanding these mechanisms should represent an important part of the study of gene regulatory networks driving neural development.

miRNA Molecules-Late Breaking Treatment for Inflammatory Bowel Diseases?

MicroRNAs (miRNAs) are a group of non-coding RNAs that play a critical role in regulating epigenetic mechanisms in inflammation-related diseases. Inflammatory bowel diseases (IBDs), which primarily include ulcerative colitis (UC) and Crohn's disease (CD), are characterized by chronic recurrent inflammation of intestinal tissues. Due to the multifactorial etiology of these diseases, the development of innovative treatment strategies that can effectively maintain remission and alleviate disease symptoms is a major challenge. In recent years, evidence for the regulatory role of miRNAs in the pathogenetic mechanisms of various diseases, including IBD, has been accumulating. In light of these findings, miRNAs represent potential innovative candidates for therapeutic application in IBD. In this review, we discuss recent findings on the role of miRNAs in regulating inflammatory responses, maintaining intestinal barrier integrity, and developing fibrosis in clinical and experimental IBD. The focus is on the existing literature, indicating potential therapeutic application of miRNAs in both preclinical experimental IBD models and translational data in the context of clinical IBD. To date, a large and diverse data set, which is growing rapidly, supports the potential use of miRNA-based therapies in clinical practice, although many questions remain unanswered.

MicroRNA 3' ends shorten during adolescent brain maturation

MicroRNA (miRNA) dysregulation is well-documented in psychiatric disease, but miRNA dynamics remain poorly understood during adolescent and early adult brain maturation, when symptoms often first appear. Here, we use RNA sequencing to examine miRNAs and their mRNA targets in cortex and hippocampus from early-, mid-, and late-adolescent and adult mice. Furthermore, we use quantitative proteomics by tandem mass tag mass spectrometry (TMT-MS) to examine protein dynamics in cortex from the same subjects. We found that ~25% of miRNAs' 3' ends shorten with age due to increased 3' trimming and decreased U tailing. Particularly, shorter but functionally competent isoforms (isomiRs) of miR-338-3p increase up to 10-fold during adolescence and only in brain. MiRNAs that undergo 3' shortening exhibit stronger negative correlations with targets that decrease with age and stronger positive correlations with targets that increase with age, than miRNAs with stable 3' ends. Increased 3' shortening with age was also observed in available mouse and human miRNA-seq data sets, and stronger correlations between miRNAs that undergo shortening and their mRNA targets were observed in two of the three available data sets. We conclude that age-associated miRNA 3' shortening is a well-conserved feature of postnatal brain maturation.

Dual regulation of miR-375 and CREM genes in pancreatic beta cells

MicroRNA-375 (miR-375) is upregulated in the islets of some diabetics and is correlated with poor outcome. Previous work in our laboratory showed that cyclic adenosine monophosphate (cAMP) reduces miR-375 expression and could provide a way to restore normal miR-375 levels, however the transcription repression mechanism is unknown. Using a chromatin immunoprecipitation assay we show that cAMP response element modulator (CREM) binds to the miR-375 promoter 3-fold above background and we find that CREM represses transcription from the miR-375 promoter 1.8-fold. While investigating miR-375 target genes we discovered that several microRNA:mRNA target prediction algorithms listed human CREM as a target gene of miR-375. The predicted binding site is conserved in primates but not in other species. We found that indeed miR-375 binds to the predicted site on human CREM and represses translation of a green fluorescent protein reporter gene by 30%. These findings suggest a primate-specific double-negative feedback loop, a mechanism that would keep these important β-cell regulators in check.

MicroRNA-181a regulates Treg functions via TGF-β1/Smad axis in the spleen of mice with acute gouty arthritis induced by MSU crystals

Regulatory T cells (Tregs) play critical roles in restricting inflammatory pathogenesis and limiting undesirable Th2 response to environmental allergens. However, the role of miR-181a in regulating acute gouty arthritis (AGA) and Treg function remains unclear. This study aimed to investigate the potential roles of miR-181a in Treg immunity and the associated signaling pathway in the AGA mouse model. A solution with monosodium urate (MSU) crystals was injected into the joint tissue of mice to induce AGA. ELISA was used to examine inflammatory factors in blood samples, and flow cytometry was used to analyze Treg profile in mice with MSU-induced AGA. Cell proliferation and viability were assessed by CCK-8 assay. TGF-β1/Smad signaling activation was detected by western blot. We found that miR-181a expression showed a positive correlation with the changes of splenic Tregs percentage in AGA mice. miR-181a regulated the TGF-β1/Smad axis, since the transfection of miR-181a mimic increased the level of TGF-β1 and the phosphorylation of Smad2/3 in Tregs in AGA mice. Additionally, miR-181a mimic also promoted responses of Tregs via TGF-β1 in vitro and in vivo. Our work uncovered a vital role of miR-181a in the immune function of Treg cells by mediating the activity of the TGF-β1/Smad pathway in the AGA mouse model induced by MSU.

Sequential and additive expression of miR-9 precursors control timing of neurogenesis

MicroRNAs (miRs) have an important role in tuning dynamic gene expression. However, the mechanism by which they are quantitatively controlled is unknown. We show that the amount of mature miR-9, a key regulator of neuronal development, increases during zebrafish neurogenesis in a sharp stepwise manner. We characterize the spatiotemporal profile of seven distinct microRNA primary transcripts (pri-mir)-9s that produce the same mature miR-9 and show that they are sequentially expressed during hindbrain neurogenesis. Expression of late-onset pri-mir-9-1 is added on to, rather than replacing, the expression of early onset pri-mir-9-4 and -9-5 in single cells. CRISPR/Cas9 mutation of the late-onset pri-mir-9-1 prevents the developmental increase of mature miR-9, reduces late neuronal differentiation and fails to downregulate Her6 at late stages. Mathematical modelling shows that an adaptive network containing Her6 is insensitive to linear increases in miR-9 but responds to stepwise increases of miR-9. We suggest that a sharp stepwise increase of mature miR-9 is created by sequential and additive temporal activation of distinct loci. This may be a strategy to overcome adaptation and facilitate a transition of Her6 to a new dynamic regime or steady state.

elk1/miR-462-731 Feedback Loop Regulates Macrophages Polarization and Phagocytosis in Grass Carp (Ctenopharyngodon idella)

MicroRNA clusters are microRNAs (miRNAs) that are distributed in close proximity on chromosomes. In this study, we report a miRNA cluster identified from grass carp (Ctenopharyngodon idella), miR-462-731, which plays a positive role in host antibacterial immunity. The expression of miR-462-731 was disrupted after infection by Aeromonas hydrophila. Transcription factor ETS transcription factor ELK1 was identified to bind to the promoter of the miR-462-731 cluster and suppress its expression. In addition, miR-731 negatively regulates the expression of elk1, forms an elk1/miR-462-731 double negative feedback loop. In addition, we found that miR-731 directly targets ezrin a (ezra), participates in inducing PI3K/AKT signaling in macrophage, to induce macrophage polarization to the M1 phenotype with stronger phagocytosis. Our results demonstrate a novel elk1/miR-462-731 feedback loop. The data deepen our understanding of the relationship between macrophage polarization and phagocytosis in teleost fish.

The hypoxia-induced changes in miRNA-mRNA in RNA-induced silencing complexes and HIF-2 induced miRNAs in human endothelial cells

The cellular adaptive response to hypoxia relies on the expression of hypoxia-inducible factors (HIFs), HIF-1 and HIF-2. HIFs regulate global gene expression changes during hypoxia that are necessary for restoring oxygen homeostasis and promoting cell survival. In the early stages of hypoxia, HIF-1 is elevated, whereas at the later stages, HIF-2 becomes the predominant form. What governs the transition between the two HIFs (the HIF switch) and the role of miRNAs in this regulation are not completely clear. Genome-wide expression studies on the miRNA content of RNA-induced silencing complexes (RISC) in HUVECs exposed to hypoxia compared to the global miRNA-Seq analysis revealed very specific differences between these two populations. We analyzed the miRNA and mRNA composition of RISC at 2 h (mainly HIF-1 driven), 8 h (HIF-1 and HIF-2 elevated), and 16 h (mainly HIF-2 driven) in a gene ontology context. This allowed for determining the direct impact of the miRNAs in modulating the cellular signaling pathways involved in the hypoxic adaptive response. Our results indicate that the miRNA-mRNA RISC components control the adaptive responses, and this does not always rely on the miRNA transcriptional elevations during hypoxia. Furthermore, we demonstrate that the hypoxic levels of the vast majority of HIF-1-dependent miRNAs (including miR-210-3p) are also HIF-2 dependent and that HIF-2 governs the expression of 11 specific miRNAs. In summary, the switch from HIF-1 to HIF-2 during hypoxia provides an important level of miRNA-driven control in the adaptive pathways in endothelial cells.

Detection of features predictive of microRNA targets by integration of network data

Gene activity is controlled by multiple molecular mechanisms, for instance through transcription factors or by microRNAs (miRNAs), among others. Established bioinformatics tools for the prediction of miRNA target genes face the challenge of ensuring accuracy, due to high false positive rates. Further, these tools present poor overlap. However, we demonstrated that it is possible to filter good predictions of miRNA targets from the bulk of all predictions by using information from the gene regulatory network. Here, we take advantage of this strategy that selects a large subset of predicted microRNA binding sites as more likely to possess less false-positives because of their over-representation in RE1 silencing transcription factor (REST)-regulated genes from the background of TargetScanHuman 7.2 predictions to identify useful features for the prediction of microRNA targets. These enriched miRNA families would have silencing activity for neural transcripts overlapping the repressive activity on neural genes of REST. We analyze properties of associated microRNA binding sites and contrast the outcome to the background. We found that the selected subset presents significant differences respect to the background: (i) lower GC-content in the vicinity of the predicted miRNA binding site, (ii) more target genes with multiple identical microRNA binding sites and (iii) a higher density of predicted microRNA binding sites close to the 3’ terminal end of the 3’-UTR. These results suggest that network selection of miRNA-mRNA pairs could provide useful features to improve microRNA target prediction.

Multi-omics analysis predicts fibronectin 1 as a prognostic biomarker in glioblastoma multiforme

Glioblastoma (GBM) is one of the most malignant and intractable central nervous system tumors with high recurrence, low survival rate, and poor prognosis. Despite the advances of aggressive, multimodal treatment, a successful treatment strategy is still elusive, often leading to therapeutic resistance and fatality. Thus, it is imperative to search for and identify novel markers critically associated with GBM pathogenesis to improve the existing trend of diagnosis, prognosis, and treatment. Seven publicly available GEO microarray datasets containing 409 GBM samples were integrated and further data mining was conducted using several bioinformatics tools. A total of 209 differentially expressed genes (DEGs) were identified in the GBM tissue samples compared to the normal brains. Gene Ontology (GO) enrichment analysis of the DEGs revealed association of the upregulates genes with extracellular matrix (ECM), conceivably contributing to the invasive nature of GBM while downregulated DEGs were found to be predominantly related to neuronal processes and structures. Alongside, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome pathway analyses described the involvement of the DEGs with various crucial contributing pathways (PI3K-Akt signaling pathway, p53 signaling pathway, insulin secretion, etc.) in GBM progression and pathogenesis. Protein-protein interaction (PPI) network containing 879 nodes and 1237 edges revealed 3 significant modules and consecutive KEGG pathway analysis of these modules showed a significant connection to gliomagenesis. Later, 10 hub genes were screened out based on degree and their expressions were externally validated. Surprisingly, only fibronectin 1 (FN1) high expression appeared to be related to poor prognosis. Subsequently, 109 transcription factors and 211 miRNAs were detected to be involved with the hub genes where FN1 demonstrated the highest number of interactions. Considering its high connectivity and potential prognostic value FN1 could be a novel biomarker providing new insights into the prognosis and treatment for GBM, although experimental validation is required.

An miRNA derived from amelogenin exon4 regulates expression of transcription factor Runx2 by directly targeting upstream activators Nfia and Prkch

Amel, the gene encoding the amelogenin protein involved in enamel formation, is highly alternatively spliced. When exon4 is excised, it can form a mature miRNA (miR-exon4) that has previously been suggested to indirectly regulate expression of the Runt-related transcription factor 2 (Runx2) involved in bone development in ameloblasts and osteoblasts. However, the precise mechanism of this regulation is unclear. In this study, we aimed to identify direct targets of miR-exon4. The transcription factor family nuclear factor I/A (NFI/A) is known to negatively regulate expression of Runx2 and is among the most highly predicted direct targets of miR-exon4 that link to Runx2. Immunostaining detected NFI/A in osteoblasts and ameloblasts in vivo, and reporter assays confirmed direct interaction of the Nfia 3'-UTR and miR-exon4. In addition, silencing of Nfia in MC3T3-E1-M14 osteoblasts resulted in subsequent downregulation of Runx2. In a monoclonal subclone (mi2) of MC3T3-E1 cells wherein mature miR-exon4 was functionally inhibited, we observed significantly downregulated Runx2 expression. We showed that NFI/A was significantly upregulated in mi2 cells at both mRNA and protein levels. Furthermore, quantitative proteomics and pathway analysis of gene expression in mi2 cells suggested that miR-exon4 could directly target Prkch (protein kinase C-eta), possibly leading to RUNX2 regulation through mechanistic target of rapamycin kinase activation. Reporter assays also confirmed the direct interaction of miR-exon4 and the 3'-UTR of Prkch, and Western blot analysis confirmed significantly upregulated mechanistic target of rapamycin kinase phosphorylation in mi2 cells. Taken together, we conclude that Nfia and Prkch expression negatively correlates with miR-exon4-mediated Runx2 regulation in vivo and in vitro, suggesting miR-exon4 directly targets Nfia and Prkch to regulate Runx2.

MiR-222-3p inhibits formation of medulloblastoma stem-like cells by targeting Notch2/c-myc signaling pathway

BACKGROUND

Medulloblastoma (MB) is an embryonal tumor of the cerebellum, which commonly occurs in childhood. Herein, we investigated the effects of miR-222-3p on the formation of MB stem-like cells via the Notch2/c-myc pathway.

METHODS

Quantitative real-time PCR (qRT-PCR) or western blotting was performed to determine the expression of miR-222-3p and Notch2, c-myc, proliferating cell nuclear antigen (PCNA), and caspase-3. Luciferase reporter gene, RNA immunoprecipitation (RIP), and RNA pull-down assay were applied to confirm the interaction between miR-222-3p and Notch2. Cell growth was examined by Cell Counting Kit-8. Cell cycle distribution and the number of stem cell marker CD133⁺ cells were examined using flow cytometry. The sphere formation assay was performed.

RESULTS

miR-222-3p expression was decreased and Notch2 expression was increased in human medulloblastoma cells. miR-222-3p overexpression inhibited cell viability and the sphere formation, induced cell cycle arrest, decreased the number of CD133⁺ cells, and up-regulated caspase-3 expression and down-regulated PCNA, Notch2, and c-myc expression. However, Notch2 overexpression counteracted these effects of miR-222-3p overexpression. Simultaneous overexpression of Notch2 and miR-222-3p increased the c-myc promoter luciferase activity which was decreased by miR-222-3p overexpression. Luciferase reporter gene, RIP, and RNA pull-down assay revealed that miR-222-3p targeted Notch2.

CONCLUSION

MiR-222-3p suppressed cell viability, altered cell cycle distribution, and inhibited the formation of MB stem-like cells via the Notch2/c-myc pathway.

Transcriptional and post-transcriptional control of epithelial-mesenchymal plasticity: why so many regulators?

The dynamic transition between epithelial-like and mesenchymal-like cell states has been a focus for extensive investigation for decades, reflective of the importance of Epithelial-Mesenchymal Transition (EMT) through development, in the adult, and the contributing role EMT has to pathologies including metastasis and fibrosis. Not surprisingly, regulation of the complex genetic networks that underlie EMT have been attributed to multiple transcription factors and microRNAs. What is surprising, however, are the sheer number of different regulators (hundreds of transcription factors and microRNAs) for which critical roles have been described. This review seeks not to collate these studies, but to provide a perspective on the fundamental question of whether it is really feasible that so many regulators play important roles and if so, what does this tell us about EMT and more generally, the genetic machinery that controls complex biological processes.

Network Approaches to Study Endogenous RNA Competition and Its Impact on Tissue-Specific microRNA Functions

microRNAs are small non-coding RNAs that play a key role in regulating gene expression. These molecules exert their function through sequence complementarity with microRNA responsive elements and are typically located in the 3' untranslated region of mRNAs, negatively regulating expression. Even though the relevant role of miRNA-dependent regulation is broadly recognized, the principles governing their ability to lead to specific functional outcomes in distinct cell types are still not well understood. In recent years, an intriguing hypothesis proposed that miRNA-responsive elements act as communication links between different RNA species, making the investigation of microRNA function even more complex than previously thought. The competing endogenous RNA hypothesis suggests the presence of a new level of regulation, whereby a specific RNA transcript can indirectly influence the abundance of other transcripts by limiting the availability of a common miRNA, acting as a "molecular sponge". Since this idea has been proposed, several studies have tried to pinpoint the interaction networks that have been established between different RNA species and whether they contribute to normal cell function and disease. The focus of this review is to highlight recent developments and achievements made towards the process of characterizing competing endogenous RNA networks and their role in cellular function.

Generation of Realistic Gene Regulatory Networks by Enriching for Feed-Forward Loops

The regulatory relationships between genes and proteins in a cell form a gene regulatory network (GRN) that controls the cellular response to changes in the environment. A number of inference methods to reverse engineer the original GRN from large-scale expression data have recently been developed. However, the absence of ground-truth GRNs when evaluating the performance makes realistic simulations of GRNs necessary. One aspect of this is that local network motif analysis of real GRNs indicates that the feed-forward loop (FFL) is significantly enriched. To simulate this properly, we developed a novel motif-based preferential attachment algorithm, FFLatt, which outperformed the popular GeneNetWeaver network generation tool in reproducing the FFL motif occurrence observed in literature-based biological GRNs. It also preserves important topological properties such as scale-free topology, sparsity, and average in/out-degree per node. We conclude that FFLatt is well-suited as a network generation module for a benchmarking framework with the aim to provide fair and robust performance evaluation of GRN inference methods.

The HDAC2/SP1/miR-205 feedback loop contributes to tubular epithelial cell extracellular matrix production in diabetic kidney disease

Extracellular matrix (ECM) accumulation is considered an important pathological feature of diabetic kidney disease (DKD). Histone deacetylase (HDAC) inhibitors protect against kidney injury. However, the potential mechanisms of HDACs in DKD are still largely unknown. Here, we describe a novel feedback loop composed of HDAC2 and miR-205 that regulates ECM production in tubular epithelial cells in individuals with DKD. We found that HDAC2 mRNA expression in peripheral blood was markedly higher in patients with DKD than in patients with diabetes. Nuclear HDAC2 protein expression was increased in TGFβ1-stimulated tubular epithelial cells and db/db mice. We also found that miR-205 was regulated by HDAC2 and downregulated in TGFβ1-treated HK2 cells and db/db mice. In addition, HDAC2 reduced histone H3K9 acetylation in the miR-205 promoter region to inhibit its promoter activity and subsequently suppressed miR-205 expression through an SP1-mediated pathway. Furthermore, miR-205 directly targeted HDAC2 and inhibited HDAC2 expression. Intriguingly, miR-205 also regulated its own transcription by inhibiting HDAC2 and increasing histone H3K9 acetylation in its promoter, forming a feedback regulatory loop. Additionally, the miR-205 agonist attenuated ECM production in HK2 cells and renal interstitial fibrosis in db/db mice. In conclusion, the HDAC2/SP1/miR-205 feedback loop may be crucial for the pathogenesis of DKD.

microRNA-Mediated Encoding and Decoding of Time-Dependent Signals in Tumorigenesis

microRNAs, pivotal post-transcriptional regulators of gene expression, in the past decades have caught the attention of researchers for their involvement in different biological processes, ranging from cell development to cancer. Although lots of effort has been devoted to elucidate the topological features and the equilibrium properties of microRNA-mediated motifs, little is known about how the information encoded in frequency, amplitude, duration, and other features of their regulatory signals can affect the resulting gene expression patterns. Here, we review the current knowledge about microRNA-mediated gene regulatory networks characterized by time-dependent input signals, such as pulses, transient inputs, and oscillations. First, we identify the general characteristic of the main motifs underlying temporal patterns. Then, we analyze their impact on two commonly studied oncogenic networks, showing how their dysfunction can lead to tumorigenesis.

Co-Regulation of Protein Coding Genes by Transcription Factor and Long Non-Coding RNA in SARS-CoV-2 Infected Cells: An In Silico Analysis

Altered expression of protein coding gene (PCG) and long non-coding RNA (lncRNA) have been identified in SARS-CoV-2 infected cells and tissues from COVID-19 patients. The functional role and mechanism (s) of transcriptional regulation of deregulated genes in COVID-19 remain largely unknown. In the present communication, reanalyzing publicly available gene expression data, we observed that 66 lncRNA and 5491 PCG were deregulated in more than one experimental condition. Combining our earlier published results and using different publicly available resources, it was observed that 72 deregulated lncRNA interacted with 3228 genes/proteins. Many targets of deregulated lncRNA could also interact with SARS-CoV-2 coded proteins, modulated by IFN treatment and identified in CRISPR screening to modulate SARS-CoV-2 infection. The majority of the deregulated lncRNA and PCG were targets of at least one of the transcription factors (TFs), interferon responsive factors (IRFs), signal transducer, and activator of transcription (STATs), NFκB, MYC, and RELA/p65. Deregulated 1069 PCG was joint targets of lncRNA and TF. These joint targets are significantly enriched with pathways relevant for SARS-CoV-2 infection indicating that joint regulation of PCG could be one of the mechanisms for deregulation. Over all this manuscript showed possible involvement of lncRNA and mechanisms of deregulation of PCG in the pathogenesis of COVID-19.

RelA/MicroRNA-30a/NLRP3 signal axis is involved in rheumatoid arthritis via regulating NLRP3 inflammasome in macrophages

NLRP3 inflammasome plays an important role in the pathogenesis of rheumatoid arthritis (RA). However, the post-transcriptional regulation of NLRP3 expression by miRNA in synovial macrophages is still not well understood. The aim of the study is to elucidate the mechanisms of RA with the focus on miRNAs mediated post-transcriptional regulation of the NLRP3 inflammasome. Here, we used NLRP3-deficient mice (NLRP3^KO) to cross with TNFα-transgenic mice (TNF^TG) to generate NLRP3^KO/TNF^TG mice, and compared their joint phenotypes with those of their TNF^TG and wild-type (WT) littermates at 5 months of age. In comparison to WT mice, articular bone volume and cartilage area are decreased, whereas inflammed area, eroded surface, ALP+ osteoblast number, TRAP+ osteoclast number, and the areas of RelA+F4/80+, Caspase-1+F4/80+, IL-1β+F4/80+ synoviocytes are increased in the TNF^TG mice. Knockout of NLRP3 ameliorates joint inflammation and bone damage in TNF^TG mice. Further, in TNFα-primed BMDMs, RelA positively regulates NLRP3 expression, but negatively regulates miR-30a. Additionally, miR-30a negatively mediates NLRP3 expression by directly binding to its 3' UTR, suggesting a miR-30a-mediated feedforward loop acting on NLRP3. Finally, intra-articular injection of AAV-miR-30a inhibits NLRP3 inflammasome activation, reduces joint inflammation, and attenuates bone damage in TNF^TG mice. Thus, RelA/miR-30a/NLRP3 signal axis is involved in RA through regulating NLRP3 Inflammasome in macrophages.

Heterogeneous information network and its application to human health and disease

The molecular components with the functional interdependencies in human cell form complicated biological network. Diseases are mostly caused by the perturbations of the composite of the interaction multi-biomolecules, rather than an abnormality of a single biomolecule. Furthermore, new biological functions and processes could be revealed by discovering novel biological entity relationships. Hence, more and more biologists focus on studying the complex biological system instead of the individual biological components. The emergence of heterogeneous information network (HIN) offers a promising way to systematically explore complicated and heterogeneous relationships between various molecules for apparently distinct phenotypes. In this review, we first present the basic definition of HIN and the biological system considered as a complex HIN. Then, we discuss the topological properties of HIN and how these can be applied to detect network motif and functional module. Afterwards, methodologies of discovering relationships between disease and biomolecule are presented. Useful insights on how HIN aids in drug development and explores human interactome are provided. Finally, we analyze the challenges and opportunities for uncovering combinatorial patterns among pharmacogenomics and cell-type detection based on single-cell genomic data.

Silent FOSL1 Enhances the Radiosensitivity of Glioma Stem Cells by Down-Regulating miR-27a-5p

Since few reports have mentioned the role of FOSL1 in the radiotherapy sensitivity of glioma, this study would dig deep into this aspect. Cancer stem cells (CSCs) isolated by magnetic bead assay were identified by microscopy, qRT-PCR and western blot. The number of apoptotic cells was counted 72 h after X-ray irradiation to evaluate the sensitivity of cancer cells to radiotherapy. The effects of radiotherapy, FOSL1 and miR-27a-5p on basic cell functions were detected by functional experiments. The expressions of FOSL1, apoptosis-related genes and miR-27a-5p were detected by qRT-PCR and western blot as needed. The differential expression of FOSL1 and the effect of miR-27a-5p on survival rate were analyzed using GEPIA and UALCAN, respectively. FOSL1 was found highly expressed in glioma cells and patients. The appearance of spherical cells and high expressions of CSC-related markers indicated the successful isolation of CSC-like cells. The increment of X-ray dose enhanced the sensitivity of cancer cells to radiotherapy. Radiotherapy down-regulated cell viability and the expressions of FOSL1 and Bcl-2, but up-regulated cell apoptosis and the expressions of cleaved caspase-3 and Bax, which could be partially reversed by overexpressed FOSL1 or further enhanced by shFOSL1. MiR-27a-5p was highly expressed in in patients with glioma, which was associated with poor prognosis, while shFOSL1-inhibited miR-27a-5p expression enhanced the sensitivity of glioma stem cells to radiotherapy. In vivo experiments further verified the results obtained from in vitro experiments. Silent FOSL1 strengthened the radiosensitivity of glioma by down-regulating miR-27a-5p.

miR-183 and miR-96 orchestrate both glucose and fat utilization in skeletal muscle

Our knowledge of the coordination of fuel usage in skeletal muscle is incomplete. Whether and how microRNAs are involved in the substrate selection for oxidation is largely unknown. Here we show that mice lacking miR-183 and miR-96 have enhanced muscle oxidative phenotype and altered glucose/lipid homeostasis. Moreover, loss of miR-183 and miR-96 results in a shift in substrate utilization toward fat relative to carbohydrates in mice. Mechanistically, loss of miR-183 and miR-96 suppresses glucose utilization in skeletal muscle by increasing PDHA1 phosphorylation via targeting FoxO1 and PDK4. On the other hand, loss of miR-183 and miR-96 promotes fat usage in skeletal muscle by enhancing intramuscular lipolysis via targeting FoxO1 and ATGL. Thus, our study establishes miR-183 and miR-96 as master coordinators of fuel selection and metabolic homeostasis owing to their capability of modulating both glucose utilization and fat catabolism. Lastly, we show that loss of miR-183 and miR-96 can alleviate obesity and improve glucose metabolism in high-fat diet-induced mice, suggesting that miR-183 and miR-96 may serve as therapeutic targets for metabolic diseases.

Exact Probability Landscapes of Stochastic Phenotype Switching in Feed-Forward Loops: Phase Diagrams of Multimodality

Feed-forward loops (FFLs) are among the most ubiquitously found motifs of reaction networks in nature. However, little is known about their stochastic behavior and the variety of network phenotypes they can exhibit. In this study, we provide full characterizations of the properties of stochastic multimodality of FFLs, and how switching between different network phenotypes are controlled. We have computed the exact steady-state probability landscapes of all eight types of coherent and incoherent FFLs using the finite-butter Accurate Chemical Master Equation (ACME) algorithm, and quantified the exact topological features of their high-dimensional probability landscapes using persistent homology. Through analysis of the degree of multimodality for each of a set of 10,812 probability landscapes, where each landscape resides over 10⁵–10⁶ microstates, we have constructed comprehensive phase diagrams of all relevant behavior of FFL multimodality over broad ranges of input and regulation intensities, as well as different regimes of promoter binding dynamics. In addition, we have quantified the topological sensitivity of the multimodality of the landscapes to regulation intensities. Our results show that with slow binding and unbinding dynamics of transcription factor to promoter, FFLs exhibit strong stochastic behavior that is very different from what would be inferred from deterministic models. In addition, input intensity play major roles in the phenotypes of FFLs: At weak input intensity, FFL exhibit monomodality, but strong input intensity may result in up to 6 stable phenotypes. Furthermore, we found that gene duplication can enlarge stable regions of specific multimodalities and enrich the phenotypic diversity of FFL networks, providing means for cells toward better adaptation to changing environment. Our results are directly applicable to analysis of behavior of FFLs in biological processes such as stem cell differentiation and for design of synthetic networks when certain phenotypic behavior is desired.

A MicroRNA-Based Network Provides Potential Predictive Signatures and Reveals the Crucial Role of PI3K/AKT Signaling for Hepatic Lineage Maturation

Background

Functions of miRNAs involved in tumorigenesis are well reported, yet, their roles in normal cell lineage commitment remain ambiguous. Here, we investigated a specific "transcription factor (TF)-miRNA-Target" regulatory network during the lineage maturation of biliary tree stem cells (BTSCs) into adult hepatocytes (hAHeps).

Method

Bioinformatic analysis was conducted based on our RNA-seq and microRNA-seq datasets with four human hepatic-lineage cell lines, including hBTSCs, hepatic stem cells (hHpSCs), hepatoblasts (hHBs), and hAHeps. Short time-series expression miner (STEM) analysis was performed to reveal the time-dependent dynamically changed miRNAs and mRNAs. GO and KEGG analyses were applied to reveal the potential function of key miRNAs and mRNAs. Then, the miRDB, miRTarBase, TargetScan, miRWalk, and DIANA-microT-CDS databases were adopted to predict the potential targets of miRNAs while the TransmiR v2.0 database was used to obtain the experimentally supported TFs that regulate miRNAs. The TCGA, Kaplan-Meier Plotter, and human protein atlas (HPA) databases and more than 10 sequencing data, including bulk RNA-seq, microRNA-seq, and scRNA-seq data related to hepatic development or lineage reprogramming, were obtained from both our or other published studies for validation.

Results

STEM analysis showed that during the maturation from hBTSCs to hAHeps, 52 miRNAs were downwardly expressed and 928 mRNA were upwardly expressed. Enrichment analyses revealed that those 52 miRNAs acted as pluripotency regulators for stem cells and participated in various novel signaling pathways, including PI3K/AKT, MAPK, and etc., while 928 mRNAs played important roles in liver-functional metabolism. With an extensive sorting of those key miRNAs and mRNAs based on the target prediction results, 23 genes were obtained which not only functioned as the targets of 17 miRNAs but were considered critical for the hepatic lineage commitment. A "TF-miRNA-Target" regulatory network for hepatic lineage commitment was therefore established and had been well validated by various datasets. The network revealed that the PI3K/AKT pathway was gradually suppressed during the hepatic commitment.

Conclusion

A total of 17 miRNAs act as suppressors during hepatic maturation mainly by regulating 23 targets and modulating the PI3K/AKT signaling pathway. The regulatory network uncovers possible signatures and guidelines enabling us to identify or obtain the functional hepatocytes for future study.