Functional characteristics of a double negative feedback loop mediated by microRNAs

Randomization based evaluation of distinct topological and cancer expression characteristics of mutually acting gene pairs

Small scale molecular network patterns and motifs are crucial for systems level understanding of cellular information transduction. Using randomizations, we statistically explored, previously overlooked basic patterns of mutually acting pairs, i.e. mutually positive (PP) or negative (NN) and positive-negative (PN) pairs, in two comprehensive and distinct large-scale molecular networks from literature; the human protein signaling network (PSN) and the human gene regulatory network (GRN). Only the positive and negative signs of all interacting pairs were randomized, while the gene pairs and the number of positive and negative signs in the original network were kept constant. While the numbers of NN and PN pairs were significantly higher, the number of PP pairs was significantly lower than randomly expected values. Genes participating in mutual pairs were more connected than other genes. NN genes were more connected than PP and PN in GRN for all types of degree values, including in, out, positive or negative connections, but less connected for in-degree and more connected for out-degree values in PSN. They also had significantly high number of intersections with each other and PN pairs than randomly expected values, indicating potential cooperative mechanisms. The three mutual interaction designs we examined had distinct RNA and protein expression correlation characteristics. NN protein pairs were uniquely over-represented across normal tissue samples, whose negative correlations were lost across cancer tissue samples. PP and PN pairs showed non-random positive RNA or protein expression correlation across normal or cancer tissue samples. Moreover, we developed an online tool, i.e. MGPNet, for further user specific analysis of mutual gene pairs. We identified SNCA with significantly enriched negatively correlated NN pairs. Unique non-random characteristics of mutual gene pairs identified in two different comprehensive molecular networks could provide valuable information for a better comparative understanding of molecular design principles between normal and cancer states. Insight Box/Paragraph Statement: This study provides a systems-level perspective on cellular information transduction by analyzing mutually acting pairs of genes. By examining mutually positive (PP), mutually negative (NN), and positive-negative (PN) pairs in the human protein signaling network (PSN) and the human gene regulatory network (GRN), we uncover significant variations in their connectivity and expression correlation. Our findings highlight the unique features of NN pairs across normal and cancer tissues and offer insights into molecular design principles. The development of the MGPNet tool further enhances user-specific analyses, enabling a deeper understanding of gene pair mechanisms and their potential cooperative roles in cellular processes.

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.

A regulatory network controlling developmental boundaries and meristem fates contributed to maize domestication

During domestication, early farmers selected different vegetative and reproductive traits, but identifying the causative loci has been hampered by their epistasis and functional redundancy. Using chromatin immunoprecipitation sequencing combined with genome-wide association analysis, we uncovered a developmental regulator that controls both types of trait while acting upstream of multiple domestication loci. tasselsheath4 (tsh4) is a new maize domestication gene that establishes developmental boundaries and specifies meristem fates despite not being expressed within them. TSH4 accomplishes this by using a double-negative feedback loop that targets and represses the very same microRNAs that negatively regulate it. TSH4 functions redundantly with a pair of homologs to positively regulate a suite of domestication loci while specifying the meristem that doubled seed yield in modern maize. TSH4 has a critical role in yield gain and helped generate ideal crop plant architecture, thus explaining why it was a major domestication target.

Activation of limbal epithelial proliferation is partly controlled by the ACE2-LCN2 pathway

In response to corneal injury, an activation of corneal epithelial stem cells and their direct progeny the early transit amplifying (eTA) cells to rapidly proliferate is critical for proper re-epithelialization. Thus, it is important to understand how such stem/eTA cell activation is regulated. Angiotensin-converting enzyme 2 (ACE2) is predominantly expressed in the stem/eTA-enriched limbal epithelium but its role in the limbal epithelium was unclear. Single cell RNA sequencing (scRNA-seq) suggested that Ace2 involved the proliferation of the stem/eTA cells. Ace2 was reduced following corneal injury. Such reduction enhanced limbal epithelial proliferation and downregulated LCN2, a negative regulator of proliferation in a variety of tissues, via upregulating TGFA and consequently activating epidermal growth factor receptor (EGFR). Inhibition of EGFR or overexpression of LCN2 reversed the increased proliferation in limbal epithelial cells lacking ACE2. Our findings demonstrate that after corneal injury, ACE2 is downregulated, which activates limbal epithelial cell proliferation via a TGFA/EGFR/LCN2 pathway.

Fetal Brain-Derived Exosomal miRNAs from Maternal Blood: Potential Diagnostic Biomarkers for Fetal Alcohol Spectrum Disorders (FASDs)

Fetal alcohol spectrum disorders (FASDs) are leading causes of neurodevelopmental disability but cannot be diagnosed early in utero. Because several microRNAs (miRNAs) are implicated in other neurological and neurodevelopmental disorders, the effects of EtOH exposure on the expression of these miRNAs and their target genes and pathways were assessed. In women who drank alcohol (EtOH) during pregnancy and non-drinking controls, matched individually for fetal sex and gestational age, the levels of miRNAs in fetal brain-derived exosomes (FB-Es) isolated from the mothers' serum correlated well with the contents of the corresponding fetal brain tissues obtained after voluntary pregnancy termination. In six EtOH-exposed cases and six matched controls, the levels of fetal brain and maternal serum miRNAs were quantified on the array by qRT-PCR. In FB-Es from 10 EtOH-exposed cases and 10 controls, selected miRNAs were quantified by ddPCR. Protein levels were quantified by ELISA. There were significant EtOH-associated reductions in the expression of several miRNAs, including miR-9 and its downstream neuronal targets BDNF, REST, Synapsin, and Sonic hedgehog. In 20 paired cases, reductions in FB-E miR-9 levels correlated strongly with reductions in fetal eye diameter, a prominent feature of FASDs. Thus, FB-E miR-9 levels might serve as a biomarker to predict FASDs in at-risk fetuses.

A double-negative feedback loop between NtrBC and a small RNA rewires nitrogen metabolism in legume symbionts

IMPORTANCE

Root nodule endosymbioses between diazotrophic rhizobia and legumes provide the largest input of combined N to the biosphere, thus representing an alternative to harmful chemical fertilizers for sustainable crop production. Rhizobia have evolved intricate strategies to coordinate N assimilation for their own benefit with N2 fixation to sustain plant growth. The rhizobial N status is transduced by the NtrBC two-component system, the seemingly ubiquitous form of N signal transduction in Proteobacteria. Here, we show that the regulatory sRNA NfeR1 (nodule formation efficiency RNA) of the alfalfa symbiont Sinorhizobium meliloti is transcribed from a complex promoter repressed by NtrC in a N-dependent manner and feedback silences ntrBC by complementary base-pairing. These findings unveil a more prominent role of NtrC as a transcriptional repressor than hitherto anticipated and a novel RNA-based mechanism for NtrBC regulation. The NtrBC-NfeR1 double-negative feedback loop accurately rewires symbiotic S. meliloti N metabolism and is likely conserved in α-rhizobia.

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.

Oscillation dynamic mechanism driven by time delays in the competent gene regulatory circuit of B. subtilis

Bacillus subtilis with competent states absorbs DNA and may improve the growth of bacteria by integrating new genetic material. Therefore, it is important to clarify how the genes interact in the circuit so that cells enter into a competent state or return to a vegetative state. The gene regulatory circuit consists of two positive feedback loops and one negative feedback loop. In this paper, a mathematical model is developed by considering transcription time delays to further study dynamic behavior of the B. subtilis competent gene regulatory network. Combined with theoretical calculation and numerical simulation, it is verified that the time delay in indirect transcription inhibition indeed has the effect of inducing the periodic oscillation of the B. subtilis competent system. In addition, some important chemical reaction rates can also regulate system dynamic behavior. However, under the control of time delay, the effects of the important chemical reaction rates have changed significantly. In particular, the time delay can advance critical value of the important chemical reaction rates where vibration occurs and can also weaken or even eliminate the effect of the important chemical reaction rates. These results will help us to analyze the competent state of B. subtilis.

The miR-200-Zeb1 axis regulates key aspects of β-cell function and survival in vivo

OBJECTIVE

The miR-200-Zeb1 axis regulates the epithelial-to-mesenchymal transition (EMT), differentiation, and resistance to apoptosis. A better understanding of these processes in diabetes is highly relevant, as β-cell dedifferentiation and apoptosis contribute to the loss of functional β-cell mass and diabetes progression. Furthermore, EMT promotes the loss of β-cell identity in the in vitro expansion of human islets. Though the miR-200 family has previously been identified as a regulator of β-cell apoptosis in vivo, studies focusing on Zeb1 are lacking. The aim of this study was thus to investigate the role of Zeb1 in β-cell function and survival in vivo.

METHODS

miR-200 and Zeb1 are involved in a double-negative feedback loop. We characterized a mouse model in which miR-200 binding sites in the Zeb1 3'UTR are mutated (Zeb1200), leading to a physiologically relevant upregulation of Zeb1 mRNA expression. The role of Zeb1 was investigated in this model via metabolic tests and analysis of isolated islets. Further insights into the distinct contributions of the miR-200 and Zeb1 branches of the feedback loop were obtained by crossing the Zeb1200 allele into a background of miR-141-200c overexpression.

RESULTS

Mild Zeb1 derepression in vivo led to broad transcriptional changes in islets affecting β-cell identity, EMT, insulin secretion, cell-cell junctions, the unfolded protein response (UPR), and the response to ER stress. The aggregation and insulin secretion of dissociated islets of mice homozygous for the Zeb1200 mutation (Zeb1200M) were impaired, and Zeb1200M islets were resistant to thapsigargin-induced ER stress ex vivo. Zeb1200M mice had increased circulating proinsulin levels but no overt metabolic phenotype, reflecting the strong compensatory ability of islets to maintain glucose homeostasis.

CONCLUSIONS

This study signifies the importance of the miR-200-Zeb1 axis in regulating key aspects of β-cell function and survival. A better understanding of this axis is highly relevant in developing therapeutic strategies for inducing β-cell redifferentiation and maintaining β-cell identity in in vitro islet expansion.

Dynamical Mechanisms for Gene Regulation Mediated by Two Noncoding RNAs in Long-Term Memory Formation

Noncoding RNAs such as miRNAs and piRNAs have long-lasting effects on the regulation of gene expression involved in long-term synaptic changes. To characterize gene regulation mediated by small noncoding RNAs associated with long-term memory in Aplysia, we consider two noncoding RNAs stimulated by 5-HT into a gene regulatory network motif model, including miR-124 that binds to and inhibits the mRNA of CREB1 and piR-F that facilitates serotonin-dependent DNA methylation to lead to repression of CREB2. Codimension-1 and -2 bifurcation analyses of 5-HT regulating both miR-124 and piR-F and a negative feedback strength for oscillation reveal rich dynamical properties of bistability and oscillations robust to variations in all other parameters. More importantly, we verify three stimulus protocols of 5-HT in experiments by our model and find that application of five pulses of 5-HT leads to a transient decrease of miR-124 but increase of piR-F concentrations, which matters sustained high level of CREB1 concentration associated with long-term memory. Furthermore, we perform bifurcation analyses for the concentrations of miR-124 and piR-F as two parameters to explore dynamical mechanisms underlying the epigenetic regulation in long-term memory formation. This study provides insights into revealing regulatory roles of epigenetic changes in gene expression involving noncoding RNAs associated with synaptic plasticity.

Dysregulation of EZH2/miR-138 axis contributes to drug resistance in multiple myeloma by downregulating RBPMS

EZH2 is highly expressed in multiple myeloma (MM). However, the molecular mechanisms underlying EZH2 overexpression and its role in drug resistance of MM remain undefined. Here we show that EZH2 is upregulated in drug-resistant MM cells and its aberrant overexpression is associated with poor prognosis of MM patients. Overexpression of EZH2 in parental MM cells renders them resistant to anti-myeloma drugs and suppression of EZH2 displays the opposite effects. Using miRNA target scan algorithms, we identify miR-138 as a regulator of EZH2, which is conversely repressed by EZH2-induced H3K27 trimethylation in MM-resistant cell lines and primary tumor cells. Analysis of ChIP-seq dataset and H3K27me3 ChIP reveals that RBPMS is a direct and functionally relevant target of EZH2. RBPMS silencing confers resistance to MM cells and restoration of RBPMS by miR-138 overexpression re-sensitizes the resistant cells to drug. Importantly, in vivo delivery of miR-138 mimics or pharmacological inhibitor of EZH2 in combination with a proteasome inhibitor, bortezomib, induces significant regression of tumors in xenograft model. This study establishes EZH2/miR-138 axis as a potential therapeutic target for MM.

Alteration in MicroRNA Expression Governs the Nature and Timing of Cellular Fate Commitment

In the central nervous system, the expression level of transcriptional repressor Hes1 (hairy and enhancer of split-1) tightly controls the alternative cell fate commitment during differentiation as well as the time required for such cellular transitions. A microRNA, miR-9, that interacts with Hes1 in a mutually antagonistic manner, influences both the process of lineage specification and timing of differentiation significantly, but the impact of the miR-9 in guiding these events still remains poorly understood. Here, we proposed a stochastic mathematical model of the miR-9/Hes1 double-negative feedback interaction network that at the outset shows how alternative cell fate such as quiescence, progenitor, and neuronal states can be accomplished through fine-tuning the Hes1 dynamics by altering the expression level of miR-9. The model simulations further foretell a correlated variation of the period of oscillation of Hes1, and the time delay observed between Hes1 mRNA and protein as the transcription rate of miR-9 increases during the neural progenitor state attainment. Importantly, the model simulations aided by the systematic sensitivity analysis predict that the timing of differentiation to the neuronal state crucially depends on the negative regulators (miR-9 and Hes6) of the Hes1. Our results indicate that miR-9/Hes1 interaction network can be effectively exploited for an efficient and well-timed neuronal transformation.

LINC01016 promotes the malignant phenotype of endometrial cancer cells by regulating the miR-302a-3p/miR-3130-3p/NFYA/SATB1 axis

Long noncoding RNAs (lncRNAs) have been implicated in tumorigenesis and cancer progression and are tightly associated with the phenotypes of numerous cancers. However, the functional roles underlying these effects are unknown. The expression levels of LINC01016, miR-302a-3p, miR-3130-3p, NFYA, and SATB1 were evaluated by quantitative real-time polymerase chain reaction (qRT-PCR) in 33 endometrial cancer tissues and 20 normal tissues. Bioinformatics analyses, luciferase reporter analyses, chromatin immunoprecipitation (ChIP) assays, and qRT-PCR assays were performed to verify potential binding sites. The qRT-PCR and western blot were used to identify the regulatory mechanisms of LINC01016 in cell biological behavior, which were also examined by cell counting kit -8 (CCK-8), 5-ethynyl-2′-deoxyuridine (EdU) assays, flow cytometry, wound healing assays, and transwell assays. LINC01016 was substantially upregulated in endometrial cancer tissues, and LINC01016 silencing abolished the malignant behavior of endometrial cancer cells. LINC01016 positively rescued the downstream gene nuclear factor YA (NFYA) by competitively “sponging” miR-302a-3p and miR-3130-3p. In turn, these two miRNAs could inhibit LINC01016 transcription, thus forming two reciprocal repression cycles, which influenced the biological behavior of endometrial cancer cells. MiR-302a-3p and miR-3130-3p could specifically bind with the 3′-UTR regions of NFYA, and NFYA could upregulate the expression of special AT-rich sequence-binding protein 1 (SATB1) as a transcriptional factor. This study was the first to show that the LINC01016–miR-302a-3p/miR-3130-3p/NFYA/SATB1 axis played a crucial role in the occurrence of endometrial cancer. These findings may provide relevant insights into the diagnosis and therapy of endometrial cancer.

ANXA1 inhibits miRNA-196a in a negative feedback loop through NF-kB and c-Myc to reduce breast cancer proliferation

MiRNAs are endogenous ~22 nt RNAs which play critical regulatory roles in a wide range of biological and pathological processes, which can act as oncogenes or tumor suppressor genes depending on their target genes. We have recently shown that ANXA1 inhibits the expression of miRNAs including miR196a. Here, we show that miR196a was highly expressed in ER+ MCF-7 breast cancer cells when compared to normal mammary gland cells, with expression levels negatively correlating to ANXA1. ANXA1 inhibits the biogenesis of oncogenic miR-196a by suppressing primary-miR196a indirectly through the stimulation of c-myc and NFkB expression and activity in breast cancer cells. In a negative feedback loop, miR-196a directly inhibits ANXA1 and enhances breast cancer cell proliferation in vitro. Finally, miR196a promotes breast tumor growth in vivo. This study reports a novel regulatory circuit between ANXA1, NF-kB, c-myc and miR-196a which regulates breast cancer cell proliferation and tumor growth.

Computational modeling of cytokine signaling in microglia

Neuroinflammation due to glial activation has been linked to many CNS diseases. We developed a computational model of a microglial cytokine interaction network to study the regulatory mechanisms of microglia-mediated neuroinflammation. We established a literature-based cytokine network, including TNFα, TGFβ, and IL-10, and fitted a mathematical model to published data from LPS-treated microglia. The addition of a previously unreported TGFβ autoregulation loop to our model was required to account for experimental data. Global sensitivity analysis revealed that TGFβ- and IL-10-mediated inhibition of TNFα was critical for regulating network behavior. We assessed the sensitivity of the LPS-induced TNFα response profile to the initial TGFβ and IL-10 levels. The analysis showed two relatively shifted TNFα response profiles within separate domains of initial condition space. Further analysis revealed that TNFα exhibited adaptation to sustained LPS stimulation. We simulated the effects of functionally inhibiting TGFβ and IL-10 on TNFα adaptation. Our analysis showed that TGFβ and IL-10 knockouts (TGFβ KO and IL-10 KO) exert divergent effects on adaptation. TFGβ KO attenuated TNFα adaptation whereas IL-10 KO enhanced TNFα adaptation. We experimentally tested the hypothesis that IL-10 KO enhances TNFα adaptation in murine macrophages and found supporting evidence. These opposing effects could be explained by differential kinetics of negative feedback. Inhibition of IL-10 reduced early negative feedback that results in enhanced TNFα-mediated TGFβ expression. We propose that differential kinetics in parallel negative feedback loops constitute a novel mechanism underlying the complex and non-intuitive pro- versus anti-inflammatory effects of individual cytokine perturbations.

Understanding microRNA-mediated gene regulatory networks through mathematical modelling

The discovery of microRNAs (miRNAs) has added a new player to the regulation of gene expression. With the increasing number of molecular species involved in gene regulatory networks, it is hard to obtain an intuitive understanding of network dynamics. Mathematical modelling can help dissecting the role of miRNAs in gene regulatory networks, and we shall here review the most recent developments that utilise different mathematical modelling approaches to provide quantitative insights into the function of miRNAs in the regulation of gene expression. Key miRNA regulation features that have been elucidated via modelling include: (i) the role of miRNA-mediated feedback and feedforward loops in fine-tuning of gene expression; (ii) the miRNA–target interaction properties determining the effectiveness of miRNA-mediated gene repression; and (iii) the competition for shared miRNAs leading to the cross-regulation of genes. However, there is still lack of mechanistic understanding of many other properties of miRNA regulation like unconventional miRNA–target interactions, miRNA regulation at different sub-cellular locations and functional miRNA variant, which will need future modelling efforts to deal with. This review provides an overview of recent developments and challenges in this field.

Current View of microRNA Processing

Small evolutionarily conserved noncoding RNAs, microRNAs (miRNAs), regulate gene expression either by translational repression or by mRNA degradation in mammals. miRNAs play functional roles in diverse physiological and pathological processes. miRNA processing is accurately regulated through multifarious factors. The canonical miRNA processing pathway consists of four sequential steps: (a) miRNA gene is transcribed into primary miRNA (pri-miRNA) mainly by RNA polymerase II; (b) pri-miRNA is processed into precursor miRNA (pre-miRNA) through microprocessor complex; (c) pre-miRNA is exported from the nucleus to the cytoplasm with the assistance of Exportin 5 (EXP5/XP05) protein; and (d) pre-miRNA is further processed into mature miRNA via Dicer. Emerging evidence has also demonstrated that some miRNAs undergo alternative processing pathways. Dysregulation of miRNA processing is closely related to tumorigenesis. Here, we review the current advances in the knowledge of miRNA processing and briefly discuss its impact on human cancers.

Bistable switch in let-7 miRNA biogenesis pathway involving Lin28

miRNAs are small noncoding RNAs capable of regulating gene expression at the post-transcriptional level. A growing body of evidence demonstrated that let-7 family of miRNAs, as one of the highly conserved miRNAs, plays an important role in cell differentiation and development, as well as tumor suppressor function depending on their levels of expression. To explore the physiological significance of let-7 in regulating cell fate decisions, we present a coarse grained model of let-7 biogenesis network, in which let-7 and its regulator Lin28 inhibit mutually. The dynamics of this minimal network architecture indicates that, as the concentration of Lin28 increases, the system undergoes a transition from monostability to a bistability and then to a one-way switch with increasing strength of positive feedback of let-7, while in the absence of Lin28 inhibition, the system loses bistability. Moreover, the ratio of degradation rates of let-7 and Lin28 is critical for the switching sensitivity and resistance to stimulus fluctuations. These findings may highlight why let-7 is required for normal gene expression in the context of embryonic development and oncogenesis, which will facilitate the development of approaches to exploit this regulatory pathway by manipulating Lin28/let-7 axis for novel treatments of human diseases.

microRNA input into a neural ultradian oscillator controls emergence and timing of alternative cell states

Progenitor maintenance, timed differentiation and the potential to enter quiescence are three fundamental processes that underlie the development of any organ system. In the nervous system, progenitor cells show short-period oscillations in the expression of the transcriptional repressor Hes1, while neurons and quiescent progenitors show stable low and high levels of Hes1, respectively. Here we use experimental data to develop a mathematical model of the double-negative interaction between Hes1 and a microRNA, miR-9, with the aim of understanding how cells transition from one state to another. We show that the input of miR-9 into the Hes1 oscillator tunes its oscillatory dynamics, and endows the system with bistability and the ability to measure time to differentiation. Our results suggest that a relatively simple and widespread network of cross-repressive interactions provides a unifying framework for progenitor maintenance, the timing of differentiation and the emergence of alternative cell states.

Hes1 is an important regulator of progenitor maintenance and timed differentiation, which shows oscillatory expression. Here, the authors combine experimental data and mathematical modelling to show that the interaction between miRNA-9 and Hes1 can predict progenitor transition from one cell state to another, as well as the timing of this transition.