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

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.

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.

RNA-binding proteins regulate aldosterone homeostasis in human steroidogenic cells

Angiotensin II (AngII) stimulates adrenocortical cells to produce aldosterone, a master regulator of blood pressure. Despite extensive characterization of the transcriptional and enzymatic control of adrenocortical steroidogenesis, there are still major gaps in the precise regulation of AII-induced gene expression kinetics. Specifically, we do not know the regulatory contribution of RNA-binding proteins (RBPs) and RNA decay, which can control the timing of stimulus-induced gene expression. To investigate this question, we performed a high-resolution RNA-seq time course of the AngII stimulation response and 4-thiouridine pulse labeling in a steroidogenic human cell line (H295R). We identified twelve temporally distinct gene expression responses that contained mRNA encoding proteins known to be important for various steps of aldosterone production, such as cAMP signaling components and steroidogenic enzymes. AngII response kinetics for many of these mRNAs revealed a coordinated increase in both synthesis and decay. These findings were validated in primary human adrenocortical cells stimulated ex vivo with AngII. Using a candidate screen, we identified a subset of RNA-binding protein and RNA decay factors that activate or repress AngII-stimulated aldosterone production. Among the repressors of aldosterone were BTG2, which promotes deadenylation and global RNA decay. BTG2 was induced in response to AngII stimulation and promoted the repression of mRNAs encoding pro-steroidogenic factors indicating the existence of an incoherent feedforward loop controlling aldosterone homeostasis. These data support a model in which coordinated increases in transcription and decay facilitate the major transcriptomic changes required to implement a pro-steroidogenic expression program that actively resolved to prevent aldosterone overproduction.

Perfect adaptation in biology

A distinctive feature of many biological systems is their ability to adapt to persistent stimuli or disturbances that would otherwise drive them away from a desirable steady state. The resulting stasis enables organisms to function reliably while being subjected to very different external environments. This perspective concerns a stringent type of biological adaptation, robust perfect adaptation (RPA), that is resilient to certain network and parameter perturbations. As in engineered control systems, RPA requires that the regulating network satisfy certain structural constraints that cannot be avoided. We elucidate these ideas using biological examples from systems and synthetic biology. We then argue that understanding the structural constraints underlying RPA allows us to look past implementation details and offers a compelling means to unravel regulatory biological complexity.

Acute 4,4'-Methylene Diphenyl Diisocyanate Exposure-Mediated Downregulation of miR-206-3p and miR-381-3p Activates Inducible Nitric Oxide Synthase Transcription by Targeting Calcineurin/NFAT Signaling in Macrophages

Exposure to 4,4'-methylene diphenyl diisocyanate (MDI) in the occupational setting may lead to development of occupational asthma (OA), and the underlying molecular mechanisms of MDI-induced disease pathogenesis remain an active area of research. Using a nose-only mouse inhalation model, we find that circulating microRNA (miR)-206-3p and miR-381-3p are downregulated after MDI exposure; however, cellular miR-206-3p and miR-381-3p responses after MDI aerosol exposure and their pathophysiological roles in MDI-OA are unknown. We hypothesize that miR-206-3p and miR-381-3p-regulated mechanisms cause increased expression of the inducible nitric oxide synthase (iNOS) after MDI aerosol exposure. We examined cellular miR-206-3p and miR-381-3p, calcineurins, nuclear factors of activated T cells (NFATs), and iNOS levels from both nose-only exposed murine bronchoalveolar lavage cells (BALCs) and differentiated THP-1 macrophages treated with MDI-glutathione (GSH) conjugates. Both in vivo murine MDI aerosol exposure and in vitro MDI-GSH exposures in THP-1 macrophages result in downregulation of endogenous miR-206-3p and miR-381-3p and upregulation of PPP3CA and iNOS expression. Transfection of THP-1 macrophages with miR-inhibitor-206-3p and miR-inhibitor-381-3p resulted in the upregulation of PPP3CA and iNOS. Using RNA-induced silencing complex immunoprecipitation and translational reporter assays, we verified that PPP3CA, but not iNOS, is directly targeted by both miR-206-3p and miR-381-3p. Downregulation of miR-206-3p and miR-381-3p following by MDI exposure induces calcineurin/NFAT signaling-mediated iNOS transcription in macrophages and BALCs.

Investigating the effects of chronic perinatal alcohol and combined nicotine and alcohol exposure on dopaminergic and non-dopaminergic neurons in the VTA

The ventral tegmental area (VTA) is the origin of dopaminergic neurons and the dopamine (DA) reward pathway. This pathway has been widely studied in addiction and drug reinforcement studies and is believed to be the central processing component of the reward circuit. In this study, we used a well-established rat model to expose mother dams to alcohol, nicotine-alcohol, and saline perinatally. DA and non-DA neurons collected from the VTA of the rat pups were used to study expression profiles of miRNAs and mRNAs. miRNA pathway interactions, putative miRNA-mRNA target pairs, and downstream modulated biological pathways were analyzed. In the DA neurons, 4607 genes were differentially upregulated and 4682 were differentially downregulated following nicotine-alcohol exposure. However, in the non-DA neurons, only 543 genes were differentially upregulated and 506 were differentially downregulated. Cell proliferation, differentiation, and survival pathways were enriched after the treatments. Specifically, in the PI3K/AKT signaling pathway, there were 41 miRNAs and 136 mRNAs differentially expressed in the DA neurons while only 16 miRNAs and 20 mRNAs were differentially expressed in the non-DA neurons after the nicotine-alcohol exposure. These results depicted that chronic nicotine and alcohol exposures during pregnancy differentially affect both miRNA and gene expression profiles more in DA than the non-DA neurons in the VTA. Understanding how the expression signatures representing specific neuronal subpopulations become enriched in the VTA after addictive substance administration helps us to identify how neuronal functions may be altered in the brain.

MicroRNA governs bistable cell differentiation and lineage segregation via a noncanonical feedback

Positive feedback driven by transcriptional regulation has long been considered a key mechanism underlying cell lineage segregation during embryogenesis. Using the developing spinal cord as a paradigm, we found that canonical, transcription-driven feedback cannot explain robust lineage segregation of motor neuron subtypes marked by two cardinal factors, Hoxa5 and Hoxc8. We propose a feedback mechanism involving elementary microRNA-mRNA reaction circuits that differ from known feedback loop-like structures. Strikingly, we show that a wide range of biologically plausible post-transcriptional regulatory parameters are sufficient to generate bistable switches, a hallmark of positive feedback. Through mathematical analysis, we explain intuitively the hidden source of this feedback. Using embryonic stem cell differentiation and mouse genetics, we corroborate that microRNA-mRNA circuits govern tissue boundaries and hysteresis upon motor neuron differentiation with respect to transient morphogen signals. Our findings reveal a previously underappreciated feedback mechanism that may have widespread functions in cell fate decisions and tissue patterning.

Distinct miRNA Signatures and Networks Discern Fetal from Adult Erythroid Differentiation and Primary from Immortalized Erythroid Cells

MicroRNAs (miRNAs) are small non-coding RNAs crucial for post-transcriptional and translational regulation of cellular and developmental pathways. The study of miRNAs in erythropoiesis elucidates underlying regulatory mechanisms and facilitates related diagnostic and therapy development. Here, we used DNA Nanoball (DNB) small RNA sequencing to comprehensively characterize miRNAs in human erythroid cell cultures. Based on primary human peripheral-blood-derived CD34+ (hCD34+) cells and two influential erythroid cell lines with adult and fetal hemoglobin expression patterns, HUDEP-2 and HUDEP-1, respectively, our study links differential miRNA expression to erythroid differentiation, cell type, and hemoglobin expression profile. Sequencing results validated by reverse-transcription quantitative PCR (RT-qPCR) of selected miRNAs indicate shared differentiation signatures in primary and immortalized cells, characterized by reduced overall miRNA expression and reciprocal expression increases for individual lineage-specific miRNAs in late-stage erythropoiesis. Despite the high similarity of same-stage hCD34+ and HUDEP-2 cells, differential expression of several miRNAs highlighted informative discrepancies between both cell types. Moreover, a comparison between HUDEP-2 and HUDEP-1 cells displayed changes in miRNAs, transcription factors (TFs), target genes, and pathways associated with globin switching. In resulting TF-miRNA co-regulatory networks, major therapeutically relevant regulators of globin expression were targeted by many co-expressed miRNAs, outlining intricate combinatorial miRNA regulation of globin expression in erythroid cells.

A critical approach for successful use of circulating microRNAs as biomarkers in cardiovascular diseases: the case of hypertrophic cardiomyopathy

MicroRNAs (miRNAs) are small noncoding RNA molecules that act as major regulators of gene expression at the post-transcriptional level. As the potential applications of miRNAs in the diagnosis and treatment of human diseases have become more evident, many studies of hypertrophic cardiomyopathy (HCM) have focused on the systemic identification and quantification of miRNAs in biofluids and myocardial tissues. HCM is a hereditary cardiomyopathy caused by mutations in genes encoding proteins of the sarcomere. Despite overall improvements in survival, progression to heart failure, stroke, and sudden cardiac death remain prominent features of living with HCM. Several miRNAs have been shown to be promising biomarkers of HCM; however, there are many challenges to ensuring the validity, consistency, and reproducibility of these biomarkers for clinical use. In particular, miRNA testing may be limited by pre-analytical and analytical caveats, making our interpretation of results challenging. Such factors that may affect miRNA testing include sample type selection, hemolysis, platelet activation, and renal dysfunction. Therefore, researchers should be careful when developing appropriate standards for the design of miRNA profiling studies in order to ensure that all results provided are both accurate and reliable. In this review, we discuss the application of miRNAs as biomarkers for HCM.

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

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