It is well known that, under suitable conditions, microRNAs are able to fine tune the relative concentration of their targets to any desired value. We show that this function is particularly effective when one of the targets is a Transcription Factor (TF) which regulates the other targets. This combination defines a new class of feed-forward loops (FFLs) in which the microRNA plays the role of master regulator. Using both deterministic and stochastic equations, we show that these FFLs are indeed able not only to fine-tune the TF/target ratio to any desired value as a function of the miRNA concentration but also, thanks to the peculiar topology of the circuit, to ensure the stability of this ratio against stochastic fluctuations. These two effects are due to the interplay between the direct transcriptional regulation and the indirect TF/Target interaction due to competition of TF and target for miRNA binding (the so called “sponge effect”). We then perform a genome wide search of these FFLs in the human regulatory network and show that they are characterized by a very peculiar enrichment pattern. In particular, they are strongly enriched in all the situations in which the TF and its target have to be precisely kept at the same concentration notwithstanding the environmental noise. As an example we discuss the FFL involving E2F1 as Transcription Factor, RB1 as target and miR-17 family as master regulator. These FFLs ensure a tight control of the E2F/RB ratio which in turns ensures the stability of the transition from the G0/G1 to the S phase in quiescent cells.

BackgroundMicroRNAs, post-transcriptional repressors of gene expression, play a pivotal role in gene regulatory networks. They are involved in core cellular processes and their dysregulation is associated to a broad range of human diseases. This paper focus on a minimal microRNA-mediated regulatory circuit, in which a protein-coding gene (host gene) is targeted by a microRNA located inside one of its introns.ResultsAutoregulation via intronic microRNAs is widespread in the human regulatory network, as confirmed by our bioinformatic analysis, and can perform several regulatory tasks despite its simple topology. Our analysis, based on analytical calculations and simulations, indicates that this circuitry alters the dynamics of the host gene expression, can induce complex responses implementing adaptation and Weber’s law, and efficiently filters fluctuations propagating from the upstream network to the host gene. A fine-tuning of the circuit parameters can optimize each of these functions. Interestingly, they are all related to gene expression homeostasis, in agreement with the increasing evidence suggesting a role of microRNA regulation in conferring robustness to biological processes. In addition to model analysis, we present a list of bioinformatically predicted candidate circuits in human for future experimental tests.ConclusionsThe results presented here suggest a potentially relevant functional role for negative self-regulation via intronic microRNAs, in particular as a homeostatic control mechanism of gene expression. Moreover, the map of circuit functions in terms of experimentally measurable parameters, resulting from our analysis, can be a useful guideline for possible applications in synthetic biology.

BackgroundThe MYC transcription factors are known to be involved in the biology of many human cancer types. But little is known about the Myc/microRNAs cooperation in the regulation of genes at the transcriptional and post-transcriptional level.Methodology/Principal FindingsEmploying independent databases with experimentally validated data, we identified several mixed microRNA/Transcription Factor Feed-Forward Loops regulated by Myc and characterized completely by experimentally supported regulatory interactions, in human. We then studied the statistical and functional properties of these circuits and discussed in more detail a few interesting examples involving E2F1, PTEN, RB1 and VEGF.Conclusions/SignificanceWe have assembled and characterized a catalogue of human mixed Transcription Factor/microRNA Feed-Forward Loops, having Myc as master regulator and completely defined by experimentally verified regulatory interactions.

Mechanical unloading by left ventricular assist devices (LVADs) in advanced heart failure (HF), in addition to improving symptoms and end-organ perfusion, is supposed to stimulate cellular and molecular responses which can reverse maladaptive cardiac remodeling. As microRNAs (miRNAs) are key regulators in remodeling processes, a comparative miRNA profiling in transplanted hearts of HF patients with/without LVAD assistance could aid to comprehend underlying molecular mechanisms. Next generation sequencing (NGS) was used to analyze miRNA differential expression in left ventricles of HF patients who underwent heart transplantation directly (n = 9) or following a period of LVAD support (n = 8). After data validation by quantitative real-time PCR, association with functional clinical parameters was investigated. Bioinformatics' tools were then used for prediction of putative targets of modulated miRNAs and relative pathway enrichment. The analysis revealed 13 upregulated and 10 downregulated miRNAs in failing hearts subjected to LVAD assistance. In particular, the expression level of some of them (miR-338-3p, miR-142-5p and -3p, miR-216a-5p, miR-223-3p, miR-27a-5p, and miR-378g) showed correlation with off-pump cardiac index values. Predicted targets of these miRNAs were involved in focal adhesion/integrin pathway and in actin cytoskeleton regulation. The identified miRNAs might contribute to molecular regulation of reverse remodeling and heart recovery mechanisms.

The tumor suppressor gene TP53 is the most frequently mutated gene in human papillomavirus (HPV)-negative head and neck squamous cell carcinoma (HNSCC). It represents a known transcription factor that controls different microRNAs (miRNA) and target genes involved in the regulation of cellular stress, apoptosis and response to DNA damage. We used The Cancer Genome Atlas database to investigate the difference in transcriptome and proteome levels between mutated and wild-type TP53 HPV-negative HNSCC. Using different databases and an extensive literature review, we built the transcriptional and post-transcriptional network regulated by TP53. TP53 mutation was associated with poor overall survival in 203 HPV-negative patients compared to 40 patients with TP53 wild-type tumors. Using the enrichment analysis, we found that UHRF1BP1 and SESN1 mRNA were linked to prognosis in the TP53 mutated group. This is also the case for miR-377-3p, an important miRNA regulator of SESN1. Our study shows that SESN1 mRNA, UHRF1BP11 mRNA and miRNA-377-3p levels are prognostically relevant in HPV-negative HNSCC patients. This finding may help with patient stratification and the development of potential new therapeutic targets to treat patients with HNSCC.

Cell-cycle pathway impairments resulting in CDK4 and 6 activation are frequently observed in human papillomavirus (HPV)-negative squamous cell carcinoma of the head and neck (SCCHN). We investigated the activity of ribociclib, a CDK4/6 inhibitor, in SCCHN models with the aim of identifying predictive biomarkers of response. HPV-negative or HPV-positive SCCHN cell lines (n ¼ 8) and patient-derived tumor xenograft (PDTX) models (n ¼ 6) were used. The models were classified according to their sensitivity to ribociclib to investigate potential predictive biomarkers. Ribociclib had a cytostatic effect in some HPV-negative SCCHN models but had no effect in HPV-positive models. In SCCHN cell lines and PDTXs, the retinoblastoma (Rb) protein expression level correlated with ribociclib activity. Rb knockdown was, however, not sufficient to block G 0 -G 1 arrest induced by ribociclib in Detroit-562 where p107, p130, and Forkhead BOX M1 (FOXM1) were also implicated in ribociclib activity. Cell lines harboring epithelial-to-mesenchymal transition (EMT) features were less sensitive to ribociclib than those with an epithelial phenotype. Rb downregulation induced EMT in our Rb-expressing SCCHN cell lines. However, ribociclib still had significant activity in one PDTX model with high Rb and vimentin expression, suggesting that the presence of vimentin alone is not enough to induce ribociclib resistance. These findings suggest that CDK4/6 inhibitors should be investigated in patients with HPV-negative SCCHN with high Rb expression and an epithelial phenotype. Although these biomarkers are not predictive in all cases, they may enrich the population that could benefit from CDK4/6 inhibitors.

BackgroundEpigenetic variation is a main regulation mechanism of gene expression in various cancer histotypes, and due to its reversibility, the potential impact in therapy can be very relevant.MethodsBased on a selected pair, breast cancer (BC) and melanoma, we conducted inference analysis in parallel on a few cell lines (MCF-7 for BC and A375 for melanoma). Starting from differential expression after treatment with a demethylating agent, the 5-Aza-2'-deoxycytidine (DAC), we provided pathway enrichment analysis and gene regulatory maps with cross-linked microRNAs and transcription factors.ResultsSeveral oncogenic signaling pathways altered upon DAC treatment were detected with significant enrichment. We represented the association between these cancers by depicting the landscape of common and specific variation affecting them.

In this study, we report that immortal mouse embryonic fibroblasts (I-MEFs) have a baseline level of cells positive for alkaline phosphatase (AP+) staining. Environmental stresses, including long-lasting growth in the absence of expansion and treatment with drugs, enhance the frequency of AP+ I-MEFs. By adapting fast red AP staining to the sorting procedure, we separated AP+ and AP− I-MEFs and demonstrated that the differentially expressed genes are consistent with a reprogrammed phenotype. In particular, we found that sestrin 1 is upregulated in AP+ I-MEFs. We focused on this gene and demonstrated that increased sestrin 1 expression is accompanied by the growth of I-MEFs in the absence of expansion and occurs before the formation of AP+ I-MEFs. Together with sestrin 1 upregulation, we found that AP+ I-MEFs accumulated in the G1 phase of the cell cycle, suggesting that the two events are causally related. Accordingly, we found that silencing sestrin 1 expression reduced the frequency and G1 accumulation of AP+ I-MEFs. Taken together, our data suggested that I-MEFs stressed by environmental changes acquire the AP+ phenotype and achieve a quiescent state characterized by a new transcriptional network.