SWIMmeR: an R-based software to unveiling crucial nodes in complex biological networks.

SUMMARY: We present SWIMmeR, an open-source version of its predecessor SWIM (SWitchMiner) that is a network-based tool for mining key (switch) genes that are associated with intriguing patterns of molecular co-abundance and may play a crucial role in phenotypic transitions in various biological settings. SWIM was originally written in MATLAB®, a proprietary programming language that requires the purchase of a license to install, manipulate, operate and run the software. Over the last years, SWIM has sparked a widespread interest within the scientific community thanks to the promising results obtained through its application in a broad range of phenotype-specific scenarios, spanning from complex diseases to grapevine berry maturation. This success has created the call for it to be distributed in a freely accessible, open-source, runtime environment, such as R, aimed at a general audience of non-expert users that cannot afford the leading proprietary solution. SWIMmeR is provided as a comprehensive collection of R functions and it also includes several additional features that make it less intensive in terms of computer time and more efficient in terms of usability and further implementation and extension. AVAILABILITY AND IMPLEMENTATION: The SWIMmeR source code is freely available at https://github.com/sportingCode/SWIMmeR.git, along with a practical user guide, including a usage example of its application on breast cancer dataset. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Hepatitis B Virus and microRNAs: A Bioinformatics Approach.

In recent decades, microRNAs (miRNAs) have emerged as key regulators of gene expression, and the identification of viral miRNAs (v-miRNAs) within some viruses, including hepatitis B virus (HBV), has attracted significant attention. HBV infections often progress to chronic states (CHB) and may induce fibrosis/cirrhosis and hepatocellular carcinoma (HCC). The presence of HBV can dysregulate host miRNA expression, influencing several biological pathways, such as apoptosis, innate and immune response, viral replication, and pathogenesis. Consequently, miRNAs are considered a promising biomarker for diagnostic, prognostic, and treatment response. The dynamics of miRNAs during HBV infection are multifaceted, influenced by host variability and miRNA interactions. Given the ability of miRNAs to target multiple messenger RNA (mRNA), understanding the viral-host (human) interplay is complex but essential to develop novel clinical applications. Therefore, bioinformatics can help to analyze, identify, and interpret a vast amount of miRNA data. This review explores the bioinformatics tools available for viral and host miRNA research. Moreover, we introduce a brief overview focusing on the role of miRNAs during HBV infection. In this way, this review aims to help the selection of the most appropriate bioinformatics tools based on requirements and research goals.

SPINNAKER: an R-based tool to highlight key RNA interactions in complex biological networks.

BACKGROUND: Recently, we developed a mathematical model for identifying putative competing endogenous RNA (ceRNA) interactions. This methodology has aroused a broad acknowledgment within the scientific community thanks to the encouraging results achieved when applied to breast invasive carcinoma, leading to the identification of PVT1, a long non-coding RNA functioning as ceRNA for the miR-200 family. The main shortcoming of the model is that it is no freely available and implemented in MATLAB®, a proprietary programming platform requiring a paid license for installing, operating, manipulating, and running the software. RESULTS: Breaking through these model limitations demands to distribute it in an open-source, freely accessible environment, such as R, designed for an ordinary audience of users that are not able to afford a proprietary solution. Here, we present SPINNAKER (SPongeINteractionNetworkmAKER), the open-source version of our widely established mathematical model for predicting ceRNAs crosstalk, that is released as an exhaustive collection of R functions. SPINNAKER has been even designed for providing many additional features that facilitate its usability, make it more efficient in terms of further implementation and extension, and less intense in terms of computational execution time. CONCLUSIONS: SPINNAKER source code is freely available at https://github.com/sportingCode/SPINNAKER.git together with a thoroughgoing PPT-based guideline. In order to help users get the key points more conveniently, also a practical R-styled plain-text guideline is provided. Finally, a short movie is available to help the user to set the own directory, properly.

StaRTrEK:in silico estimation of RNA half-lives from genome-wide time-course experiments without transcriptional inhibition.

BACKGROUND: Gene expression is the result of the balance between transcription and degradation. Recent experimental findings have shown fine and specific regulation of RNA degradation and the presence of various molecular machinery purposely devoted to this task, such as RNA binding proteins, non-coding RNAs, etc. A biological process can be studied by measuring time-courses of RNA abundance in response of internal and/or external stimuli, using recent technologies, such as the microarrays or the Next Generation Sequencing devices. Unfortunately, the picture provided by looking only at the transcriptome abundance may not gain insight into its dynamic regulation. By contrast, independent simultaneous measurement of RNA expression and half-lives could provide such valuable additional insight. A computational approach to the estimation of RNAs half-lives from RNA expression time profiles data, can be a low-cost alternative to its experimental measurement which may be also affected by various artifacts. RESULTS: Here we present a computational methodology, called StaRTrEK (STAbility Rates ThRough Expression Kinetics), able to estimate half-life values basing only on genome-wide gene expression time series without transcriptional inhibition. The StaRTrEK algorithm makes use of a simple first order kinetic model and of a [Formula: see text]-norm regularized least square optimization approach to find its parameter values. Estimates provided by StaRTrEK are validated using simulated data and three independent experimental datasets of two short (6 samples) and one long (48 samples) time-courses. CONCLUSIONS: We believe that our algorithm can be used as a fast valuable computational complement to time-course experimental gene expression studies by adding a relevant kinetic property, i.e. the RNA half-life, with a strong biological interpretation, thus providing a dynamic picture of what is going in a cell during the biological process under study.

A gene dosage-dependent effect unveils NBS1 as both a haploinsufficient tumour suppressor and an essential gene for SHH-medulloblastoma.

AIMS: Inherited or somatic mutations in the MRE11, RAD50 and NBN genes increase the incidence of tumours, including medulloblastoma (MB). On the other hand, MRE11, RAD50 and NBS1 protein components of the MRN complex are often overexpressed and sometimes essential in cancer. In order to solve the apparent conundrum about the oncosuppressive or oncopromoting role of the MRN complex, we explored the functions of NBS1 in an MB-prone animal model. MATERIALS AND METHODS: We generated and analysed the monoallelic or biallelic deletion of the Nbn gene in the context of the SmoA1 transgenic mouse, a Sonic Hedgehog (SHH)-dependent MB-prone animal model. We used normal and tumour tissues from these animal models, primary granule cell progenitors (GCPs) from genetically modified animals and NBS1-depleted primary MB cells, to uncover the effects of NBS1 depletion by RNA-Seq, by biochemical characterisation of the SHH pathway and the DNA damage response (DDR) as well as on the growth and clonogenic properties of GCPs. RESULTS: We found that monoallelic Nbn deletion increases SmoA1-dependent MB incidence. In addition to a defective DDR, Nbn+/- GCPs show increased clonogenicity compared to Nbn+/+ GCPs, dependent on an enhanced Notch signalling. In contrast, full NbnKO impairs MB development both in SmoA1 mice and in an SHH-driven tumour allograft. CONCLUSIONS: Our study indicates that Nbn is haploinsufficient for SHH-MB development whereas full NbnKO is epistatic on SHH-driven MB development, thus revealing a gene dosage-dependent effect of Nbn inactivation on SHH-MB development.

Convolutional Neural Networks for Automated Classification of Prostate Multiparametric Magnetic Resonance Imaging Based on Image Quality.

BACKGROUND: Prostate magnetic resonance imaging (MRI) is technically demanding, requiring high image quality to reach its full diagnostic potential. An automated method to identify diagnostically inadequate images could help optimize image quality. PURPOSE: To develop a convolutional neural networks (CNNs) based analysis pipeline for the classification of prostate MRI image quality. STUDY TYPE: Retrospective. SUBJECTS: Three hundred sixteen prostate mpMRI scans and 312 men (median age 67). FIELD STRENGTH/SEQUENCE: A 3 T; fast spin echo T2WI, echo planar imaging DWI, ADC, gradient-echo dynamic contrast enhanced (DCE). ASSESSMENT: MRI scans were reviewed by three genitourinary radiologists (V.P., M.D.M., S.C.) with 21, 12, and 5 years of experience, respectively. Sequences were labeled as high quality (Q1) or low quality (Q0) and used as the reference standard for all analyses. STATISTICAL TESTS: Sequences were split into training, validation, and testing sets (869, 250, and 120 sequences, respectively). Inter-reader agreement was assessed with the Fleiss kappa. Following preprocessing and data augmentation, 28 CNNs were trained on MRI slices for each sequence. Model performance was assessed on both a per-slice and a per-sequence basis. A pairwise t-test was performed to compare performances of the classifiers. RESULTS: The number of sequences labeled as Q0 or Q1 was 38 vs. 278 for T2WI, 43 vs. 273 for DWI, 41 vs. 275 for ADC, and 38 vs. 253 for DCE. Inter-reader agreement was almost perfect for T2WI and DCE and substantial for DWI and ADC. On the per-slice analysis, accuracy was 89.95% ± 0.02% for T2WI, 79.83% ± 0.04% for DWI, 76.64% ± 0.04% for ADC, 96.62% ± 0.01% for DCE. On the per-sequence analysis, accuracy was 100% ± 0.00% for T2WI, DWI, and DCE, and 92.31% ± 0.00% for ADC. The three best algorithms performed significantly better than the remaining ones on every sequence (P-value < 0.05). DATA CONCLUSION: CNNs achieved high accuracy in classifying prostate MRI image quality on an individual-slice basis and almost perfect accuracy when classifying the entire sequences. EVIDENCE LEVEL: 4 TECHNICAL EFFICACY: Stage 1.

SAveRUNNER: A network-based algorithm for drug repurposing and its application to COVID-19.

The novelty of new human coronavirus COVID-19/SARS-CoV-2 and the lack of effective drugs and vaccines gave rise to a wide variety of strategies employed to fight this worldwide pandemic. Many of these strategies rely on the repositioning of existing drugs that could shorten the time and reduce the cost compared to de novo drug discovery. In this study, we presented a new network-based algorithm for drug repositioning, called SAveRUNNER (Searching off-lAbel dRUg aNd NEtwoRk), which predicts drug-disease associations by quantifying the interplay between the drug targets and the disease-specific proteins in the human interactome via a novel network-based similarity measure that prioritizes associations between drugs and diseases locating in the same network neighborhoods. Specifically, we applied SAveRUNNER on a panel of 14 selected diseases with a consolidated knowledge about their disease-causing genes and that have been found to be related to COVID-19 for genetic similarity (i.e., SARS), comorbidity (e.g., cardiovascular diseases), or for their association to drugs tentatively repurposed to treat COVID-19 (e.g., malaria, HIV, rheumatoid arthritis). Focusing specifically on SARS subnetwork, we identified 282 repurposable drugs, including some the most rumored off-label drugs for COVID-19 treatments (e.g., chloroquine, hydroxychloroquine, tocilizumab, heparin), as well as a new combination therapy of 5 drugs (hydroxychloroquine, chloroquine, lopinavir, ritonavir, remdesivir), actually used in clinical practice. Furthermore, to maximize the efficiency of putative downstream validation experiments, we prioritized 24 potential anti-SARS-CoV repurposable drugs based on their network-based similarity values. These top-ranked drugs include ACE-inhibitors, monoclonal antibodies (e.g., anti-IFNγ, anti-TNFα, anti-IL12, anti-IL1ß, anti-IL6), and thrombin inhibitors. Finally, our findings were in-silico validated by performing a gene set enrichment analysis, which confirmed that most of the network-predicted repurposable drugs may have a potential treatment effect against human coronavirus infections.

A Comparison of Network-Based Methods for Drug Repurposing along with an Application to Human Complex Diseases.

Drug repurposing strategy, proposing a therapeutic switching of already approved drugs with known medical indications to new therapeutic purposes, has been considered as an efficient approach to unveil novel drug candidates with new pharmacological activities, significantly reducing the cost and shortening the time of de novo drug discovery. Meaningful computational approaches for drug repurposing exploit the principles of the emerging field of Network Medicine, according to which human diseases can be interpreted as local perturbations of the human interactome network, where the molecular determinants of each disease (disease genes) are not randomly scattered, but co-localized in highly interconnected subnetworks (disease modules), whose perturbation is linked to the pathophenotype manifestation. By interpreting drug effects as local perturbations of the interactome, for a drug to be on-target effective against a specific disease or to cause off-target adverse effects, its targets should be in the nearby of disease-associated genes. Here, we used the network-based proximity measure to compute the distance between the drug module and the disease module in the human interactome by exploiting five different metrics (minimum, maximum, mean, median, mode), with the aim to compare different frameworks for highlighting putative repurposable drugs to treat complex human diseases, including malignant breast and prostate neoplasms, schizophrenia, and liver cirrhosis. Whilst the standard metric (that is the minimum) for the network-based proximity remained a valid tool for efficiently screening off-label drugs, we observed that the other implemented metrics specifically predicted further interesting drug candidates worthy of investigation for yielding a potentially significant clinical benefit.

Repurposing Histaminergic Drugs in Multiple Sclerosis.

Multiple sclerosis is an autoimmune disease with a strong neuroinflammatory component that contributes to severe demyelination, neurodegeneration and lesions formation in white and grey matter of the spinal cord and brain. Increasing attention is being paid to the signaling of the biogenic amine histamine in the context of several pathological conditions. In multiple sclerosis, histamine regulates the differentiation of oligodendrocyte precursors, reduces demyelination, and improves the remyelination process. However, the concomitant activation of histamine H1-H4 receptors can sustain either damaging or favorable effects, depending on the specifically activated receptor subtype/s, the timing of receptor engagement, and the central versus peripheral target district. Conventional drug development has failed so far to identify curative drugs for multiple sclerosis, thus causing a severe delay in therapeutic options available to patients. In this perspective, drug repurposing offers an exciting and complementary alternative for rapidly approving some medicines already approved for other indications. In the present work, we have adopted a new network-medicine-based algorithm for drug repurposing called SAveRUNNER, for quantifying the interplay between multiple sclerosis-associated genes and drug targets in the human interactome. We have identified new histamine drug-disease associations and predicted off-label novel use of the histaminergic drugs amodiaquine, rupatadine, and diphenhydramine among others, for multiple sclerosis. Our work suggests that selected histamine-related molecules might get to the root causes of multiple sclerosis and emerge as new potential therapeutic strategies for the disease.

SAveRUNNER: an R-based tool for drug repurposing.

BACKGROUND: Currently, no proven effective drugs for the novel coronavirus disease COVID-19 exist and despite widespread vaccination campaigns, we are far short from herd immunity. The number of people who are still vulnerable to the virus is too high to hamper new outbreaks, leading a compelling need to find new therapeutic options devoted to combat SARS-CoV-2 infection. Drug repurposing represents an effective drug discovery strategy from existing drugs that could shorten the time and reduce the cost compared to de novo drug discovery. RESULTS: We developed a network-based tool for drug repurposing provided as a freely available R-code, called SAveRUNNER (Searching off-lAbel dRUg aNd NEtwoRk), with the aim to offer a promising framework to efficiently detect putative novel indications for currently marketed drugs against diseases of interest. SAveRUNNER predicts drug-disease associations by quantifying the interplay between the drug targets and the disease-associated proteins in the human interactome through the computation of a novel network-based similarity measure, which prioritizes associations between drugs and diseases located in the same network neighborhoods. CONCLUSIONS: The algorithm was successfully applied to predict off-label drugs to be repositioned against the new human coronavirus (2019-nCoV/SARS-CoV-2), and it achieved a high accuracy in the identification of well-known drug indications, thus revealing itself as a powerful tool to rapidly detect potential novel medical indications for various drugs that are worth of further investigation. SAveRUNNER source code is freely available at https://github.com/giuliafiscon/SAveRUNNER.git , along with a comprehensive user guide.

Genetic, epigenetic and genomic effects on variation of gene expression among grape varieties.

The transcriptional regulatory structure of plant genomes is still relatively unexplored, and little is known about factors that influence expression variation in plants. We used a genetic system consisting of 10 heterozygous grape varieties with high consanguinity and high haplotypic diversity to: (i) identify regions of haplotype sharing through whole-genome resequencing and single-nucleotide polymorphism (SNP) genotyping; (ii) analyse gene expression through RNA-seq in four stages of berry development; and (iii) associate gene expression variation with genetic and epigenetic properties. We found that haplotype sharing in and around genes was positively correlated with similarity in expression and was negatively correlated with the fraction of differentially expressed genes. Genetic and epigenetic properties of the gene and the surrounding region showed significant effects on the extent of expression variation, with negative associations for the level of gene body methylation and mean expression level, and with positive associations for nucleotide diversity, structural diversity and ratio of non-synonymous to synonymous nucleotide diversity. We also observed a spatial dependency of covariation of gene expression among varieties. These results highlight relevant roles for cis-acting factors, selective constraints and epigenetic features of the gene, and the regional context in which the gene is located, in the determination of expression variation. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/bioproject/PRJNA385116; https://www.ncbi.nlm.nih.gov/bioproject/PRJNA392287; https://www.ncbi.nlm.nih.gov/bioproject/PRJNA373967 (released upon publication); https://www.ncbi.nlm.nih.gov/bioproject/PRJNA490160 (released upon publication); https://www.ncbi.nlm.nih.gov/bioproject/PRJNA265039; https://www.ncbi.nlm.nih.gov/bioproject/PRJNA265040.

The New Paradigm of Network Medicine to Analyze Breast Cancer Phenotypes.

Breast cancer (BC) is a heterogeneous and complex disease as witnessed by the existence of different subtypes and clinical characteristics that poses significant challenges in disease management. The complexity of this tumor may rely on the highly interconnected nature of the various biological processes as stated by the new paradigm of Network Medicine. We explored The Cancer Genome Atlas (TCGA)-BRCA data set, by applying the network-based algorithm named SWItch Miner, and mapping the findings on the human interactome to capture the molecular interconnections associated with the disease modules. To characterize BC phenotypes, we constructed protein-protein interaction modules based on "hub genes", called switch genes, both common and specific to the four tumor subtypes. Transcriptomic profiles of patients were stratified according to both clinical (immunohistochemistry) and genetic (PAM50) classifications. 266 and 372 switch genes were identified from immunohistochemistry and PAM50 classifications, respectively. Moreover, the identified switch genes were functionally characterized to select an interconnected pathway of disease genes. By intersecting the common switch genes of the two classifications, we selected a unique signature of 28 disease genes that were BC subtype-independent and classification subtype-independent. Data were validated both in vitro (10 BC cell lines) and ex vivo (66 BC tissues) experiments. Results showed that four of these hub proteins (AURKA, CDC45, ESPL1, and RAD54L) were over-expressed in all tumor subtypes. Moreover, the inhibition of one of the identified switch genes (AURKA) similarly affected all BC subtypes. In conclusion, using a network-based approach, we identified a common BC disease module which might reflect its pathological signature, suggesting a new vision to face with the disease heterogeneity.

Network Analysis of Gut Microbiome and Metabolome to Discover Microbiota-Linked Biomarkers in Patients Affected by Non-Small Cell Lung Cancer.

Several studies in recent times have linked gut microbiome (GM) diversity to the pathogenesis of cancer and its role in disease progression through immune response, inflammation and metabolism modulation. This study focused on the use of network analysis and weighted gene co-expression network analysis (WGCNA) to identify the biological interaction between the gut ecosystem and its metabolites that could impact the immunotherapy response in non-small cell lung cancer (NSCLC) patients undergoing second-line treatment with anti-PD1. Metabolomic data were merged with operational taxonomic units (OTUs) from 16S RNA-targeted metagenomics and classified by chemometric models. The traits considered for the analyses were: (i) condition: disease or control (CTRLs), and (ii) treatment: responder (R) or non-responder (NR). Network analysis indicated that indole and its derivatives, aldehydes and alcohols could play a signaling role in GM functionality. WGCNA generated, instead, strong correlations between short-chain fatty acids (SCFAs) and a healthy GM. Furthermore, commensal bacteria such as Akkermansia muciniphila, Rikenellaceae, Bacteroides, Peptostreptococcaceae, Mogibacteriaceae and Clostridiaceae were found to be more abundant in CTRLs than in NSCLC patients. Our preliminary study demonstrates that the discovery of microbiota-linked biomarkers could provide an indication on the road towards personalized management of NSCLC patients.

MIENTURNET: an interactive web tool for microRNA-target enrichment and network-based analysis.

BACKGROUND: miRNAs regulate the expression of several genes with one miRNA able to target multiple genes and with one gene able to be simultaneously targeted by more than one miRNA. Therefore, it has become indispensable to shorten the long list of miRNA-target interactions to put in the spotlight in order to gain insight into understanding the regulatory mechanism orchestrated by miRNAs in various cellular processes. A reasonable solution is certainly to prioritize miRNA-target interactions to maximize the effectiveness of the downstream analysis. RESULTS: We propose a new and easy-to-use web tool MIENTURNET (MicroRNA ENrichment TURned NETwork) that receives in input a list of miRNAs or mRNAs and tackles the problem of prioritizing miRNA-target interactions by performing a statistical analysis followed by a fully featured network-based visualization and analysis. The statistics is used to assess the significance of an over-representation of miRNA-target interactions and then MIENTURNET filters based on the statistical significance associated with each miRNA-target interaction. In addition, the holistic approach of the network theory is used to infer possible evidences of miRNA regulation by capturing emergent properties of the miRNA-target regulatory network that would be not evident through a pairwise analysis of the individual components. CONCLUSION: MIENTURNET offers the possibility to consistently perform both statistical and network-based analyses by using only a single tool leading to a more effective prioritization of the miRNA-target interactions. This has the potential to avoid researchers without computational and informatics skills to navigate multiple websites and thus to independently investigate miRNA activity in every cellular process of interest in an easy and at the same time exhaustive way thanks to the intuitive web interface. The web application along with a well-documented and comprehensive user guide are freely available at http://userver.bio.uniroma1.it/apps/mienturnet/ without any login requirement.

JARID1B expression and its function in DNA damage repair are tightly regulated by miRNAs in breast cancer.

JARID1B/KDM5B histone demethylase's mRNA is markedly overexpressed in breast cancer tissues and cell lines and the protein has been shown to have a prominent role in cancer cell proliferation and DNA repair. However, the mechanism of its post-transcriptional regulation in cancer cells remains elusive. We performed a computational analysis of transcriptomic data from a set of 103 breast cancer patients, which, along with JARID1B upregulation, showed a strong downregulation of 2 microRNAs (miRNAs), mir-381 and mir-486, potentially targeting its mRNA. We showed that both miRNAs can target JARID1B 3'UTR and reduce luciferase's activity in a complementarity-driven repression assay. Moreover, MCF7 breast cancer cells overexpressing JARID1B showed a strong protein reduction when transfected with mir-486. This protein's decrease is accompanied by accumulation of DNA damage, enhanced radiosensitivity and increase of BRCA1 mRNA, 3 features previously correlated with JARID1B silencing. These results enlighten an important role of a miRNA's circuit in regulating JARID1B's activity and suggest new perspectives for epigenetic therapies.

SWIM tool application to expression data of glioblastoma stem-like cell lines, corresponding primary tumors and conventional glioma cell lines.

BACKGROUND: It is well-known that glioblastoma contains self-renewing, stem-like subpopulation with the ability to sustain tumor growth. These cells - called cancer stem-like cells - share certain phenotypic characteristics with untransformed stem cells and are resistant to many conventional cancer therapies, which might explain the limitations in curing human malignancies. Thus, the identification of genes controlling the differentiation of these stem-like cells is becoming a successful therapeutic strategy, owing to the promise of novel targets for treating malignancies. METHODS: Recently, we developed SWIM, a software able to unveil a small pool of genes - called switch genes - critically associated with drastic changes in cell phenotype. Here, we applied SWIM to the expression profiling of glioblastoma stem-like cells and conventional glioma cell lines, in order to identify switch genes related to stem-like phenotype. RESULTS: SWIM identifies 171 switch genes that are all down-regulated in glioblastoma stem-like cells. This list encompasses genes like CAV1, COL5A1, COL6A3, FLNB, HMMR, ITGA3, ITGA5, MET, SDC1, THBS1, and VEGFC, involved in "ECM-receptor interaction" and "focal adhesion" pathways. The inhibition of switch genes highly correlates with the activation of genes related to neural development and differentiation, such as the 4-core OLIG2, POU3F2, SALL2, SOX2, whose induction has been shown to be sufficient to reprogram differentiated glioblastoma into stem-like cells. Among switch genes, the transcription factor FOSL1 appears as the brightest star since: it is down-regulated in stem-like cells; it highly negatively correlates with the 4-core genes that are all up-regulated in stem-like cells; the promoter regions of the 4-core genes harbor a consensus binding motif for FOSL1. CONCLUSIONS: We suggest that the inhibition of switch genes in stem-like cells could induce the deregulation of cell communication pathways, contributing to neoplastic progression and tumor invasiveness. Conversely, their activation could restore the physiological equilibrium between cell adhesion and migration, hampering the progression of cancer. Moreover, we posit FOSL1 as promising candidate to orchestrate the differentiation of cancer stem-like cells by repressing the 4-core genes' expression, which severely halts cancer growth and might affect the therapeutic outcome. We suggest FOSL1 as novel putative therapeutic and prognostic biomarker, worthy of further investigation.

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).

Ripening Transcriptomic Program in Red and White Grapevine Varieties Correlates with Berry Skin Anthocyanin Accumulation.

Grapevine (Vitis vinifera) berry development involves a succession of physiological and biochemical changes reflecting the transcriptional modulation of thousands of genes. Although recent studies have investigated the dynamic transcriptome during berry development, most have focused on a single grapevine variety, so there is a lack of comparative data representing different cultivars. Here, we report, to our knowledge, the first genome-wide transcriptional analysis of 120 RNA samples corresponding to 10 Italian grapevine varieties collected at four growth stages. The 10 varieties, representing five red-skinned and five white-skinned berries, were all cultivated in the same experimental vineyard to reduce environmental variability. The comparison of transcriptional changes during berry formation and ripening allowed us to determine the transcriptomic traits common to all varieties, thus defining the core transcriptome of berry development, as well as the transcriptional dynamics underlying differences between red and white berry varieties. A greater variation among the red cultivars than between red and white cultivars at the transcriptome level was revealed, suggesting that anthocyanin accumulation during berry maturation has a direct impact on the transcriptomic regulation of multiple biological processes. The expression of genes related to phenylpropanoid/flavonoid biosynthesis clearly distinguished the behavior of red and white berry genotypes during ripening but also reflected the differential accumulation of anthocyanins in the red berries, indicating some form of cross talk between the activation of stilbene biosynthesis and the accumulation of anthocyanins in ripening berries.

A feature-based integrated scoring scheme for cell cycle-regulated genes prioritization.

Prioritization of cell cycle-regulated genes from expression time-profiles is still an open problem. The point at issue is the surprisingly poor overlap among ranked lists obtained from different experimental protocols. Instead of developing a general-purpose computational methodology for detecting periodic signals, we focus on the budding yeast mitotic cell cycle. The reason being that the current availability of a total of 12 datasets, produced by 6 independent groups using 4 different synchronization methods, permits a re-analysis and re-consideration of this problem in a more reliable and extensive data domain. Notably, budding yeast is a model organism for studying cancer and testing new drugs. Here we propose a novel multi-feature score (called PERLA, PERiodicity, Regulation and Lag-Autocorrelation) that integrates different features of cell cycle-regulated gene expression time-profiles. We obtained increased performances on a wide range of benchmarks and, most importantly, a substantially increased overlap of the top ranking genes among different datasets, thus proving the effectiveness of the proposed prioritization algorithm. Examples on how to use PERLA to gain new insight into the biology of the cell cycle, are provided in a final dedicated section.

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.

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.