Paola Paci

Computational analysis identifies a sponge interaction network between long non-coding RNAs and messenger RNAs in human breast cancer

BackgroundNon-coding RNAs (ncRNAs) are emerging as key regulators of many cellular processes in both physiological and pathological states. Moreover, the constant discovery of new non-coding RNA species suggests that the study of their complex functions is still in its very early stages. This variegated class of RNA species encompasses the well-known microRNAs (miRNAs) and the most recently acknowledged long non-coding RNAs (lncRNAs). Interestingly, in the last couple of years, a few studies have shown that some lncRNAs can act as miRNA sponges, i.e. as competing endogenous RNAs (ceRNAs), able to reduce the amount of miRNAs available to target messenger RNAs (mRNAs).ResultsWe propose a computational approach to explore the ability of lncRNAs to act as ceRNAs by protecting mRNAs from miRNA repression. A seed match analysis was performed to validate the underlying regression model. We built normal and cancer networks of miRNA-mediated sponge interactions (MMI-networks) using breast cancer expression data provided by The Cancer Genome Atlas.ConclusionsOur study highlights a marked rewiring in the ceRNA program between normal and pathological breast tissue, documented by its “on/off” switch from normal to cancer, and vice-versa. This mutually exclusive activation confers an interesting character to ceRNAs as potential oncosuppressive, or oncogenic, protagonists in cancer. At the heart of this phenomenon is the lncRNA PVT1, as illustrated by both the width of its antagonist mRNAs in normal-MMI-network, and the relevance of the latter in breast cancer. Interestingly, PVT1 revealed a net binding preference towards the mir-200 family as the bone of contention with its rival mRNAs.

Gut microbiota profiling of pediatric nonalcoholic fatty liver disease and obese patients unveiled by an integrated meta-omics-based approach

There is evidence that nonalcoholic fatty liver disease (NAFLD) is affected by gut microbiota. Therefore, we investigated its modifications in pediatric NAFLD patients using targeted metagenomics and metabolomics. Stools were collected from 61 consecutive patients diagnosed with nonalcoholic fatty liver (NAFL), nonalcoholic steatohepatitis (NASH), or obesity and 54 healthy controls (CTRLs), matched in a case‐control fashion. Operational taxonomic units were pyrosequenced targeting 16S ribosomal RNA and volatile organic compounds determined by solid‐phase microextraction gas chromatography‐mass spectrometry. The α‐diversity was highest in CTRLs, followed by obese, NASH, and NAFL patients; and β‐diversity distinguished between patients and CTRLs but not NAFL and NASH. Compared to CTRLs, in NAFLD patients Actinobacteria were significantly increased and Bacteroidetes reduced. There were no significant differences among the NAFL, NASH, and obese groups. Overall NAFLD patients had increased levels of Bradyrhizobium, Anaerococcus, Peptoniphilus, Propionibacterium acnes, Dorea, and Ruminococcus and reduced proportions of Oscillospira and Rikenellaceae compared to CTRLs. After reducing metagenomics and metabolomics data dimensionality, multivariate analyses indicated a decrease of Oscillospira in NAFL and NASH groups and increases of Ruminococcus, Blautia, and Dorea in NASH patients compared to CTRLs. Of the 292 volatile organic compounds, 26 were up‐regulated and 2 down‐regulated in NAFLD patients. Multivariate analyses found that combination of Oscillospira, Rickenellaceae, Parabacteroides, Bacteroides fragilis, Sutterella, Lachnospiraceae, 4‐methyl‐2‐pentanone, 1‐butanol, and 2‐butanone could discriminate NAFLD patients from CTRLs. Univariate analyses found significantly lower levels of Oscillospira and higher levels of 1‐pentanol and 2‐butanone in NAFL patients compared to CTRLs. In NASH, lower levels of Oscillospira were associated with higher abundance of Dorea and Ruminococcus and higher levels of 2‐butanone and 4‐methyl‐2‐pentanone compared to CTRLs. Conclusion: An Oscillospira decrease coupled to a 2‐butanone up‐regulation and increases in Ruminococcus and Dorea were identified as gut microbiota signatures of NAFL onset and NAFL‐NASH progression, respectively. (Hepatology 2017;65:451‐464)

PVT1: A Rising Star among Oncogenic Long Noncoding RNAs

It is becoming increasingly clear that short and long noncoding RNAs critically participate in the regulation of cell growth, differentiation, and (mis)function. However, while the functional characterization of short non-coding RNAs has been reaching maturity, there is still a paucity of well characterized long noncoding RNAs, even though large studies in recent years are rapidly increasing the number of annotated ones. The long noncoding RNA PVT1 is encoded by a gene that has been long known since it resides in the well-known cancer risk region 8q24. However, a couple of accidental concurrent conditions have slowed down the study of this gene, that is, a preconception on the primacy of the protein-coding over noncoding RNAs and the prevalent interest in its neighbor MYC oncogene. Recent studies have brought PVT1 under the spotlight suggesting interesting models of functioning, such as competing endogenous RNA activity and regulation of protein stability of important oncogenes, primarily of the MYC oncogene. Despite some advancements in modelling the PVT1 role in cancer, there are many questions that remain unanswered concerning the precise molecular mechanisms underlying its functioning.

ImmunoGrid, an integrative environment for large-scale simulation of the immune system for vaccine discovery, design and optimization

Vaccine research is a combinatorial science requiring computational analysis of vaccine components, formulations and optimization. We have developed a framework that combines computational tools for the study of immune function and vaccine development. This framework, named ImmunoGrid combines conceptual models of the immune system, models of antigen processing and presentation, system-level models of the immune system, Grid computing, and database technology to facilitate discovery, formulation and optimization of vaccines. ImmunoGrid modules share common conceptual models and ontologies. The ImmunoGrid portal offers access to educational simulators where previously defined cases can be displayed, and to research simulators that allow the development of new, or tuning of existing, computational models. The portal is accessible at .

Integrated Network Analysis Identifies Fight-Club Nodes as a Class of Hubs Encompassing Key Putative Switch Genes That Induce Major Transcriptome Reprogramming during Grapevine Development

A small pool of putative switch genes, whose expression is negatively correlated with that of genes outside their community in the network, appear to mark the shift from immature to mature growth. We developed an approach that integrates different network-based methods to analyze the correlation network arising from large-scale gene expression data. By studying grapevine (Vitis vinifera) and tomato (Solanum lycopersicum) gene expression atlases and a grapevine berry transcriptomic data set during the transition from immature to mature growth, we identified a category named “fight-club hubs” characterized by a marked negative correlation with the expression profiles of neighboring genes in the network. A special subset named “switch genes” was identified, with the additional property of many significant negative correlations outside their own group in the network. Switch genes are involved in multiple processes and include transcription factors that may be considered master regulators of the previously reported transcriptome remodeling that marks the developmental shift from immature to mature growth. All switch genes, expressed at low levels in vegetative/green tissues, showed a significant increase in mature/woody organs, suggesting a potential regulatory role during the developmental transition. Finally, our analysis of tomato gene expression data sets showed that wild-type switch genes are downregulated in ripening-deficient mutants. The identification of known master regulators of tomato fruit maturation suggests our method is suitable for the detection of key regulators of organ development in different fleshy fruit crops.

Shedding light on protein–ligand binding by graph theory: The topological nature of allostery

Allostery is a very important feature of proteins; we propose a mesoscopic approach to allosteric mechanisms elucidation, based on protein contact matrices. The application of graph theory methods to the characterization of the allosteric process and, more broadly, to obtain the conformational changes upon binding, reveals key features of the protein function. The proposed method highlights the leading role played by topological over geometrical changes in allosteric transitions. Topological invariants were able to discriminate between true allosteric motions and generic protein motions upon binding.

Modeling lymphocyte homing and encounters in lymph nodes

BackgroundThe efficiency of lymph nodes depends on tissue structure and organization, which allow the coordination of lymphocyte traffic. Despite their essential role, our understanding of lymph node specific mechanisms is still incomplete and currently a topic of intense research.ResultsIn this paper, we present a hybrid discrete/continuous model of the lymph node, accounting for differences in cell velocity and chemotactic response, influenced by the spatial compartmentalization of the lymph node and the regulation of cells migration, encounter, and antigen presentation during the inflammation process.ConclusionOur model reproduces the correct timing of an immune response, including the observed time delay between duplication of T helper cells and duplication of B cells in response to antigen exposure. Furthermore, we investigate the consequences of the absence of dendritic cells at different times during infection, and the dependence of system dynamics on the regulation of lymphocyte exit from lymph nodes. In both cases, the model predicts the emergence of an impaired immune response, i.e., the response is significantly reduced in magnitude. Dendritic cell removal is also shown to delay the response time with respect to normal conditions.

Phylogenetic and Metabolic Tracking of Gut Microbiota during Perinatal Development

The colonization and development of gut microbiota immediately after birth is highly variable and depends on several factors, such as delivery mode and modality of feeding during the first months of life. A cohort of 31 mother and neonate pairs, including 25 at-term caesarean (CS) and 6 vaginally (V) delivered neonates (DNs), were included in this study and 121 meconium/faecal samples were collected at days 1 through 30 following birth. Operational taxonomic units (OTUs) were assessed in 69 stool samples by phylogenetic microarray HITChip and inter- and intra-individual distributions were established by inter-OTUs correlation matrices and OTUs co-occurrence or co-exclusion networks. 1H-NMR metabolites were determined in 70 stool samples, PCA analysis was performed on 55 CS DNs samples, and metabolome/OTUs co-correlations were assessed in 45 CS samples, providing an integrated map of the early microbiota OTUs-metabolome. A microbiota “core” of OTUs was identified that was independent of delivery mode and lactation stage, suggesting highly specialized communities that act as seminal colonizers of microbial networks. Correlations among OTUs, metabolites, and OTUs-metabolites revealed metabolic profiles associated with early microbial ecological dynamics, maturation of milk components, and host physiology.

Role of the long non-coding RNA PVT1 in the dysregulation of the ceRNA-ceRNA network in human breast cancer.

Recent findings have identified competing endogenous RNAs (ceRNAs) as the drivers in many disease conditions, including cancers. The ceRNAs indirectly regulate each other by reducing the amount of microRNAs (miRNAs) available to target messenger RNAs (mRNAs). The ceRNA interactions mediated by miRNAs are modulated by a titration mechanism, i.e. large changes in the ceRNA expression levels either overcome, or relieve, the miRNA repression on competing RNAs; similarly, a very large miRNA overexpression may abolish competition. The ceRNAs are also called miRNA “decoys” or miRNA “sponges” and encompass different RNAs competing with each other to attract miRNAs for interactions: mRNA, long non-coding RNAs (lncRNAs), pseudogenes, or circular RNAs. Recently, we developed a computational method for identifying ceRNA-ceRNA interactions in breast invasive carcinoma. We were interested in unveiling which lncRNAs could exert the ceRNA activity. We found a drastic rewiring in the cross-talks between ceRNAs from the physiological to the pathological condition. The main actor of this dysregulated lncRNA-associated ceRNA network was the lncRNA PVT1, which revealed a net biding preference towards the miR-200 family members in normal breast tissues. Despite its up-regulation in breast cancer tissues, mimicked by the miR-200 family members, PVT1 stops working as ceRNA in the cancerous state. The specific conditions required for a ceRNA landscape to occur are still far from being determined. Here, we emphasized the importance of the relative concentration of the ceRNAs, and their related miRNAs. In particular, we focused on the withdrawal in breast cancer tissues of the PVT1 ceRNA activity and performed a gene expression and sequence analysis of its multiple isoforms. We found that the PVT1 isoform harbouring the binding site for a representative miRNA of the miR-200 family shows a drastic decrease in its relative concentration with respect to the miRNA abundance in breast cancer tissues, providing a plausibility argument to the breakdown of the sponge program orchestrated by the oncogene PVT1.

Immune control of HIV-1 infection after therapy interruption: immediate versus deferred antiretroviral therapy

BackgroundThe optimal stage for initiating antiretroviral therapies in HIV-1 bearing patients is still a matter of debate.MethodsWe present computer simulations of HIV-1 infection aimed at identifying the pro et contra of immediate as compared to deferred Highly Active Antiretroviral Therapy (HAART).ResultsOur simulations highlight that a prompt specific CD8+ cytotoxic T lymphocytes response is detected when therapy is delayed. Compared to very early initiation of HAART, in deferred treated patients CD8+ T cells manage to mediate the decline of viremia in a shorter time and, at interruption of therapy, the virus experiences a stronger immune pressure. We also observe, however, that the immunological effects of the therapy fade with time in both therapeutic regimens. Thus, within one year from discontinuation, viral burden recovers to the value at which it would level off in the absence of therapy.In summary, simulations show that immediate therapy does not prolong the disease-free period and does not confer a survival benefit when compared to treatment started during the chronic infection phase.ConclusionOur conclusion is that, since there is no therapy to date that guarantees life-long protection, deferral of therapy should be preferred in order to minimize the risk of adverse effects, the occurrence of drug resistances and the costs of treatment.