SWIM: a computational tool to unveiling crucial nodes in complex biological networks

Expression and function of NF-Y subunits in cancer

NF-Y is a Transcription Factor (TF) targeting the CCAAT box regulatory element. It consists of the NF-YB/NF-YC heterodimer, each containing an Histone Fold Domain (HFD), and the sequence-specific subunit NF-YA. NF-YA expression is associated with cell proliferation and absent in some post-mitotic cells. The review summarizes recent findings impacting on cancer development. The logic of the NF-Y regulome points to pro-growth, oncogenic genes in the cell-cycle, metabolism and transcriptional regulation routes. NF-YA is involved in growth/differentiation decisions upon cell-cycle re-entry after mitosis and it is widely overexpressed in tumors, the HFD subunits in some tumor types or subtypes. Overexpression of NF-Y -mostly NF-YA- is oncogenic and decreases sensitivity to anti-neoplastic drugs. The specific roles of NF-YA and NF-YC isoforms generated by alternative splicing -AS- are discussed, including the prognostic value of their levels, although the specific molecular mechanisms of activity are still to be deciphered.

Inferring Interaction Networks from Transcriptomic Data: Methods and Applications

Transcriptomic data is a treasure trove in modern molecular biology, as it offers a comprehensive viewpoint into the intricate nuances of gene expression dynamics underlying biological systems. This genetic information must be utilized to infer biomolecular interaction networks that can provide insights into the complex regulatory mechanisms underpinning the dynamic cellular processes. Gene regulatory networks and protein-protein interaction networks are two major classes of such networks. This chapter thoroughly investigates the wide range of methodologies used for distilling insightful revelations from transcriptomic data that include association-based methods (based on correlation among expression vectors), probabilistic models (using Bayesian and Gaussian models), and interologous methods. We reviewed different approaches for evaluating the significance of interactions based on the network topology and biological functions of the interacting molecules and discuss various strategies for the identification of functional modules. The chapter concludes with highlighting network-based techniques of prioritizing key genes, outlining the centrality-based, diffusion- based, and subgraph-based methods. The chapter provides a meticulous framework for investigating transcriptomic data to uncover assembly of complex molecular networks for their adaptable analyses across a broad spectrum of biological domains.

Network medicine for patients' stratification: From single-layer to multi-omics

Precision medicine research increasingly relies on the integrated analysis of multiple types of omics. In the era of big data, the large availability of different health-related information represents a great, but at the same time untapped, chance with a potentially fundamental role in the prevention, diagnosis and prognosis of diseases. Computational methods are needed to combine this data to create a comprehensive view of a given disease. Network science can model biomedical data in terms of relationships among molecular players of different nature and has been successfully proposed as a new paradigm for studying human diseases. Patient stratification is an open challenge aimed at identifying subtypes with different disease manifestations, severity, and expected survival time. Several stratification approaches based on high-throughput gene expression measurements have been successfully applied. However, few attempts have been proposed to exploit the integration of various genotypic and phenotypic data to discover novel sub-types or improve the detection of known groupings. This article is categorized under: Cancer > Biomedical Engineering Cancer > Computational Models Cancer > Genetics/Genomics/Epigenetics.

Transcriptomic analysis reveals the gene regulatory networks involved in leaf and root response to osmotic stress in tomato

Introduction

Tomato (Solanum lycopersicum L.) is a major horticultural crop that is cultivated worldwide and is characteristic of the Mediterranean agricultural system. It represents a key component of the diet of billion people and an important source of vitamins and carotenoids. Tomato cultivation in open field often experiences drought episodes, leading to severe yield losses, since most modern cultivars are sensitive to water deficit. Water stress leads to changes in the expression of stress-responsive genes in different plant tissues, and transcriptomics can support the identification of genes and pathways regulating this response.

Methods

Here, we performed a transcriptomic analysis of two tomato genotypes, M82 and Tondo, in response to a PEG-mediated osmotic treatment. The analysis was conducted separately on leaves and roots to characterize the specific response of these two organs.

Results

A total of 6,267 differentially expressed transcripts related to stress response was detected. The construction of gene co-expression networks defined the molecular pathways of the common and specific responses of leaf and root. The common response was characterized by ABA-dependent and ABA-independent signaling pathways, and by the interconnection between ABA and JA signaling. The root-specific response concerned genes involved in cell wall metabolism and remodeling, whereas the leaf-specific response was principally related to leaf senescence and ethylene signaling. The transcription factors representing the hubs of these regulatory networks were identified. Some of them have not yet been characterized and can represent novel candidates for tolerance.

Discussion

This work shed new light on the regulatory networks occurring in tomato leaf and root under osmotic stress and set the base for an in-depth characterization of novel stress-related genes that may represent potential candidates for improving tolerance to abiotic stress in tomato.

Co-Expression Network Analysis Identifies Molecular Determinants of Loneliness Associated with Neuropsychiatric and Neurodegenerative Diseases

Loneliness and social isolation are detrimental to mental health and may lead to cognitive impairment and neurodegeneration. Although several molecular signatures of loneliness have been identified, the molecular mechanisms by which loneliness impacts the brain remain elusive. Here, we performed a bioinformatics approach to untangle the molecular underpinnings associated with loneliness. Co-expression network analysis identified molecular 'switches' responsible for dramatic transcriptional changes in the nucleus accumbens of individuals with known loneliness. Loneliness-related switch genes were enriched in cell cycle, cancer, TGF-β, FOXO, and PI3K-AKT signaling pathways. Analysis stratified by sex identified switch genes in males with chronic loneliness. Male-specific switch genes were enriched in infection, innate immunity, and cancer-related pathways. Correlation analysis revealed that loneliness-related switch genes significantly overlapped with 82% and 68% of human studies on Alzheimer's (AD) and Parkinson's diseases (PD), respectively, in gene expression databases. Loneliness-related switch genes, BCAM, NECTIN2, NPAS3, RBM38, PELI1, DPP10, and ASGR2, have been identified as genetic risk factors for AD. Likewise, switch genes HLA-DRB5, ALDOA, and GPNMB are known genetic loci in PD. Similarly, loneliness-related switch genes overlapped in 70% and 64% of human studies on major depressive disorder and schizophrenia, respectively. Nine switch genes, HLA-DRB5, ARHGAP15, COL4A1, RBM38, DMD, LGALS3BP, WSCD2, CYTH4, and CNTRL, overlapped with known genetic variants in depression. Seven switch genes, NPAS3, ARHGAP15, LGALS3BP, DPP10, SMYD3, CPXCR1, and HLA-DRB5 were associated with known risk factors for schizophrenia. Collectively, we identified molecular determinants of loneliness and dysregulated pathways in the brain of non-demented adults. The association of switch genes with known risk factors for neuropsychiatric and neurodegenerative diseases provides a molecular explanation for the observed prevalence of these diseases among lonely individuals.

Networks as Biomarkers: Uses and Purposes

Networks-based approaches are often used to analyze gene expression data or protein-protein interactions but are not usually applied to study the relationships between different biomarkers. Given the clinical need for more comprehensive and integrative biomarkers that can help to identify personalized therapies, the integration of biomarkers of different natures is an emerging trend in the literature. Network analysis can be used to analyze the relationships between different features of a disease; nodes can be disease-related phenotypes, gene expression, mutational events, protein quantification, imaging-derived features and more. Since different biomarkers can exert causal effects between them, describing such interrelationships can be used to better understand the underlying mechanisms of complex diseases. Networks as biomarkers are not yet commonly used, despite being proven to lead to interesting results. Here, we discuss in which ways they have been used to provide novel insights into disease susceptibility, disease development and severity.

A Network of MicroRNAs and mRNAs Involved in Melanosome Maturation and Trafficking Defines the Lower Response of Pigmentable Melanoma Cells to Targeted Therapy

BACKGROUND

The ability to increase their degree of pigmentation is an adaptive response that confers pigmentable melanoma cells higher resistance to BRAF inhibitors (BRAFi) compared to non-pigmentable melanoma cells.

METHODS

Here, we compared the miRNome and the transcriptome profile of pigmentable 501Mel and SK-Mel-5 melanoma cells vs. non-pigmentable A375 melanoma cells, following treatment with the BRAFi vemurafenib (vem). In depth bioinformatic analyses (clusterProfiler, WGCNA and SWIMmeR) allowed us to identify the miRNAs, mRNAs and biological processes (BPs) that specifically characterize the response of pigmentable melanoma cells to the drug. Such BPs were studied using appropriate assays in vitro and in vivo (xenograft in zebrafish embryos).

RESULTS

Upon vem treatment, miR-192-5p, miR-211-5p, miR-374a-5p, miR-486-5p, miR-582-5p, miR-1260a and miR-7977, as well as GPR143, OCA2, RAB27A, RAB32 and TYRP1 mRNAs, are differentially expressed only in pigmentable cells. These miRNAs and mRNAs belong to BPs related to pigmentation, specifically melanosome maturation and trafficking. In fact, an increase in the number of intracellular melanosomes-due to increased maturation and/or trafficking-confers resistance to vem.

CONCLUSION

We demonstrated that the ability of pigmentable cells to increase the number of intracellular melanosomes fully accounts for their higher resistance to vem compared to non-pigmentable cells. In addition, we identified a network of miRNAs and mRNAs that are involved in melanosome maturation and/or trafficking. Finally, we provide the rationale for testing BRAFi in combination with inhibitors of these biological processes, so that pigmentable melanoma cells can be turned into more sensitive non-pigmentable cells.

Sex-specific transcriptional rewiring in the brain of Alzheimer’s disease patients

Sex-specific differences may contribute to Alzheimer’s disease (AD) development. AD is more prevalent in women worldwide, and female sex has been suggested as a disease risk factor. Nevertheless, the molecular mechanisms underlying sex-biased differences in AD remain poorly characterized. To this end, we analyzed the transcriptional changes in the entorhinal cortex of symptomatic and asymptomatic AD patients stratified by sex. Co-expression network analysis implemented by SWItchMiner software identified sex-specific signatures of switch genes responsible for drastic transcriptional changes in the brain of AD and asymptomatic AD individuals. Pathway analysis of the switch genes revealed that morphine addiction, retrograde endocannabinoid signaling, and autophagy are associated with both females with AD (F-AD) and males with (M-AD). In contrast, nicotine addiction, cell adhesion molecules, oxytocin signaling, adipocytokine signaling, prolactin signaling, and alcoholism are uniquely associated with M-AD. Similarly, some of the unique pathways associated with F-AD switch genes are viral myocarditis, Hippo signaling pathway, endometrial cancer, insulin signaling, and PI3K-AKT signaling. Together these results reveal that there are many sex-specific pathways that may lead to AD. Approximately 20–30% of the elderly have an accumulation of amyloid beta in the brain, but show no cognitive deficit. Asymptomatic females (F-asymAD) and males (M-asymAD) both shared dysregulation of endocytosis. In contrast, pathways uniquely associated with F-asymAD switch genes are insulin secretion, progesterone-mediated oocyte maturation, axon guidance, renal cell carcinoma, and ErbB signaling pathway. Similarly, pathways uniquely associated with M-asymAD switch genes are fluid shear stress and atherosclerosis, FcγR mediated phagocytosis, and proteoglycans in cancer. These results reveal for the first time unique pathways associated with either disease progression or cognitive resilience in asymptomatic individuals. Additionally, we identified numerous sex-specific transcription factors and potential neurotoxic chemicals that may be involved in the pathogenesis of AD. Together these results reveal likely molecular drivers of sex differences in the brain of AD patients. Future molecular studies dissecting the functional role of these switch genes in driving sex differences in AD are warranted.

Physical Activity Rewires the Human Brain against Neurodegeneration

Physical activity may offset cognitive decline and dementia, but the molecular mechanisms by which it promotes neuroprotection remain elusive. In the absence of disease-modifying therapies, understanding the molecular effects of physical activity in the brain may be useful for identifying novel targets for disease management. Here we employed several bioinformatic methods to dissect the molecular underpinnings of physical activity in brain health. Network analysis identified 'switch genes' associated with drastic hippocampal transcriptional changes in aged cognitively intact individuals. Switch genes are key genes associated with dramatic transcriptional changes and thus may play a fundamental role in disease pathogenesis. Switch genes are associated with protein processing pathways and the metabolic control of glucose, lipids, and fatty acids. Correlation analysis showed that transcriptional patterns associated with physical activity significantly overlapped and negatively correlated with those of neurodegenerative diseases. Functional analysis revealed that physical activity might confer neuroprotection in Alzheimer's (AD), Parkinson's (PD), and Huntington's (HD) diseases via the upregulation of synaptic signaling pathways. In contrast, in frontotemporal dementia (FTD) its effects are mediated by restoring mitochondrial function and energy precursors. Additionally, physical activity is associated with the downregulation of genes involved in inflammation in AD, neurogenesis in FTD, regulation of growth and transcriptional repression in PD, and glial cell differentiation in HD. Collectively, these findings suggest that physical activity directs transcriptional changes in the brain through different pathways across the broad spectrum of neurodegenerative diseases. These results provide new evidence on the unique and shared mechanisms between physical activity and neurodegenerative diseases.

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.

GCEN: An Easy-to-Use Toolkit for Gene Co-Expression Network Analysis and lncRNAs Annotation

Gene co-expression network analysis has been widely used in gene function annotation, especially for long noncoding RNAs (lncRNAs). However, there is a lack of effective cross-platform analysis tools. For biologists to easily build a gene co-expression network and to predict gene function, we developed GCEN, a cross-platform command-line toolkit developed with C++. It is an efficient and easy-to-use solution that will allow everyone to perform gene co-expression network analysis without the requirement of sophisticated programming skills, especially in cases of RNA-Seq research and lncRNAs function annotation. Because of its modular design, GCEN can be easily integrated into other pipelines.

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.

Gene Network Analysis of Hepatocellular Carcinoma Identifies Modules Associated with Disease Progression, Survival, and Chemo Drug Resistance

Background

Hepatocellular carcinoma (HCC) is the second leading cause of cancer-related mortality worldwide. HCC transcriptome has been extensively studied; however, the progress in disease mechanisms, prognosis, and treatment is still slow.

Methods

A rank-based module-centric workflow was introduced to analyze important modules associated with HCC development, prognosis, and drug resistance. The currently largest HCC cell line RNA-Seq dataset from the LIMORE database was used to construct the reference modules by weighted gene co-expression network analysis.

Results

Thirteen reference modules were identified with validated reproducibility. These modules were all associated with specific biological functions. Differentially expressed module analysis revealed the crucial modules during HCC development. Modules and hub genes are indicative of patient survival. Modules can differentiate patients in different HCC stages. Furthermore, drug resistance was revealed by drug-module association analysis. Based on differentially expressed modules and hub genes, six candidate drugs were screened. The hub genes of those modules merit further investigation.

Conclusion

We proposed a reference module-based analysis of the HCC transcriptome. The identified modules are associated with HCC development, survival, and drug resistance. M3 and M6 may play important roles during HCV to HCC development. M1, M3, M5, and M7 are associated with HCC survival. High M4, high M9, low M1, and low M3 may be associated with dasatinib, doxorubicin, CD532, and simvastatin resistance. Our analysis provides useful information for HCC diagnosis and treatment.

Towards a global investigation of transcriptomic signatures through co-expression networks and pathway knowledge for the identification of disease mechanisms

We attempt to address a key question in the joint analysis of transcriptomic data: can we correlate the patterns we observe in transcriptomic datasets to known interactions and pathway knowledge to broaden our understanding of disease pathophysiology? We present a systematic approach that sheds light on the patterns observed in hundreds of transcriptomic datasets from over sixty indications by using pathways and molecular interactions as a template. Our analysis employs transcriptomic datasets to construct dozens of disease specific co-expression networks, alongside a human protein-protein interactome network. Leveraging the interoperability between these two network templates, we explore patterns both common and particular to these diseases on three different levels. Firstly, at the node-level, we identify most and least common proteins across diseases and evaluate their consistency against the interactome as a proxy for their prevalence in the scientific literature. Secondly, we overlay both network templates to analyze common correlations and interactions across diseases at the edge-level. Thirdly, we explore the similarity between patterns observed at the disease-level and pathway knowledge to identify signatures associated with specific diseases and indication areas. Finally, we present a case scenario in schizophrenia, where we show how our approach can be used to investigate disease pathophysiology.

Reference Module-Based Analysis of Ovarian Cancer Transcriptome Identifies Important Modules and Potential Drugs

Ovarian cancer (OVC) is often diagnosed at the advanced stage resulting in a poor overall outcome for the patient. The disease mechanisms, prognosis, and treatment require imperative elucidation. A rank-based module-centric framework was proposed to analyze the key modules related to the development, prognosis, and treatment of OVC. The ovarian cancer cell line microarray dataset GSE43765 from the Gene Expression Omnibus database was used to construct the reference modules by weighted gene correlation network analysis. Twenty-three reference modules were tested for stability and functionally annotated. Furthermore, to demonstrate the utility of reference modules, two more OVC datasets were collected, and their gene expression profiles were projected to the reference modules to generate a module-level expression. An epithelial-mesenchymal transition module was activated in OVC compared to the normal epithelium, and a pluripotency module was activated in ovarian cancer stroma compared to ovarian cancer epithelium. Seven differentially expressed modules were identified in OVC compared to the normal ovarian epithelium, with five up-regulated, and two down-regulated. One module was identified to be predictive of patient overall survival. Four modules were enriched with SNP signals. Based on differentially expressed modules and hub genes, five candidate drugs were screened. The hub genes of those modules merit further investigation. We firstly propose the reference module-based analysis of OVC. The utility of the analysis framework can be extended to transcriptome data of other kinds of diseases.

Gene network analysis using SWIM reveals interplay between the transcription factor-encoding genes HMGA1, FOXM1, and MYBL2 in triple-negative breast cancer

Among breast cancer subtypes, triple-negative breast cancer (TNBC) is the most aggressive with the worst prognosis and the highest rates of metastatic disease. To identify TNBC gene signatures, we applied the network-based methodology implemented by the SWIM software to gene expression data of TNBC patients in The Cancer Genome Atlas (TCGA) database. SWIM enables to predict key (switch) genes within the co-expression network, whose perturbations in expression pattern and abundance may contribute to the (patho)biological phenotype. Here, SWIM analysis revealed an interesting interplay between the genes encoding the transcription factors HMGA1, FOXM1, and MYBL2, suggesting a potential cooperation among these three switch genes in TNBC development. The correlative nature of this interplay in TNBC was assessed by in vitro experiments, demonstrating how they may actually modulate the expression of each other.

Key Disease Mechanisms Linked to Alzheimer's Disease in the Entorhinal Cortex

Alzheimer’s disease (AD) is a chronic, neurodegenerative brain disorder affecting millions of Americans that is expected to increase in incidence with the expanding aging population. Symptomatic AD patients show cognitive decline and often develop neuropsychiatric symptoms due to the accumulation of insoluble proteins that produce plaques and tangles seen in the brain at autopsy. Unexpectedly, some clinically normal individuals also show AD pathology in the brain at autopsy (asymptomatic AD, AsymAD). In this study, SWItchMiner software was used to identify key switch genes in the brain’s entorhinal cortex that lead to the development of AD or disease resilience. Seventy-two switch genes were identified that are differentially expressed in AD patients compared to healthy controls. These genes are involved in inflammation, platelet activation, and phospholipase D and estrogen signaling. Peroxisome proliferator-activated receptor γ (PPARG), zinc-finger transcription factor (YY1), sterol regulatory element-binding transcription factor 2 (SREBF2), and early growth response 1 (EGR1) were identified as transcription factors that potentially regulate switch genes in AD. Comparing AD patients to AsymAD individuals revealed 51 switch genes; PPARG as a potential regulator of these genes, and platelet activation and phospholipase D as critical signaling pathways. Chemical–protein interaction analysis revealed that valproic acid is a therapeutic agent that could prevent AD from progressing.

Clinical epigenetics settings for cancer and cardiovascular diseases: real-life applications of network medicine at the bedside

Despite impressive efforts invested in epigenetic research in the last 50 years, clinical applications are still lacking. Only a few university hospital centers currently use epigenetic biomarkers at the bedside. Moreover, the overall concept of precision medicine is not widely recognized in routine medical practice and the reductionist approach remains predominant in treating patients affected by major diseases such as cancer and cardiovascular diseases. By its' very nature, epigenetics is integrative of genetic networks. The study of epigenetic biomarkers has led to the identification of numerous drugs with an increasingly significant role in clinical therapy especially of cancer patients. Here, we provide an overview of clinical epigenetics within the context of network analysis. We illustrate achievements to date and discuss how we can move from traditional medicine into the era of network medicine (NM), where pathway-informed molecular diagnostics will allow treatment selection following the paradigm of precision medicine.

Network diffusion with centrality measures to identify disease-related genes

Disease-related gene prioritization is one of the most well-established pharmaceutical techniques used to identify genes that are important to a biological process relevant to a disease. In identifying these essential genes, the network diffusion (ND) approach is a widely used technique applied in gene prioritization. However, there is still a large number of candidate genes that need to be evaluated experimentally. Therefore, it would be of great value to develop a new strategy to improve the precision of the prioritization. Given the efficiency and simplicity of centrality measures in capturing a gene that might be important to the network structure, herein, we propose a technique that extends the scope of ND through a centrality measure to identify new disease-related genes. Five common centrality measures with different aspects were examined for integration in the traditional ND model. A total of 40 diseases were used to test our developed approach and to find new genes that might be related to a disease. Results indicated that the best measure to combine with the diffusion is closeness centrality. The novel candidate genes identified by the model for all 40 diseases were provided along with supporting evidence. In conclusion, the integration of network centrality in ND is a simple but effective technique to discover more precise disease-related genes, which is extremely useful for biomedical science.

Gene co-expression in the interactome: moving from correlation toward causation via an integrated approach to disease module discovery

In this study, we integrate the outcomes of co-expression network analysis with the human interactome network to predict novel putative disease genes and modules. We first apply the SWItch Miner (SWIM) methodology, which predicts important (switch) genes within the co-expression network that regulate disease state transitions, then map them to the human protein-protein interaction network (PPI, or interactome) to predict novel disease-disease relationships (i.e., a SWIM-informed diseasome). Although the relevance of switch genes to an observed phenotype has been recently assessed, their performance at the system or network level constitutes a new, potentially fascinating territory yet to be explored. Quantifying the interplay between switch genes and human diseases in the interactome network, we found that switch genes associated with specific disorders are closer to each other than to other nodes in the network, and tend to form localized connected subnetworks. These subnetworks overlap between similar diseases and are situated in different neighborhoods for pathologically distinct phenotypes, consistent with the well-known topological proximity property of disease genes. These findings allow us to demonstrate how SWIM-based correlation network analysis can serve as a useful tool for efficient screening of potentially new disease gene associations. When integrated with an interactome-based network analysis, it not only identifies novel candidate disease genes, but also may offer testable hypotheses by which to elucidate the molecular underpinnings of human disease and reveal commonalities between seemingly unrelated diseases.

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.

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.

Personalized analysis of breast cancer using sample-specific networks

Breast cancer is a disease with high heterogeneity. Cancer is not usually caused by a single gene, but by multiple genes and their interactions with others and surroundings. Estimating breast cancer-specific gene–gene interaction networks is critical to elucidate the mechanisms of breast cancer from a biological network perspective. In this study, sample-specific gene–gene interaction networks of breast cancer samples were established by using a sample-specific network analysis method based on gene expression profiles. Then, gene–gene interaction networks and pathways related to breast cancer and its subtypes and stages were further identified. The similarity and difference among these subtype-related (and stage-related) networks and pathways were studied, which showed highly specific for subtype Basal-like and Stages IV and V. Finally, gene pairwise interactions associated with breast cancer prognosis were identified by a Cox proportional hazards regression model, and a risk prediction model based on the gene pairs was established, which also performed very well on an independent validation data set. This work will help us to better understand the mechanism underlying the occurrence of breast cancer from the sample-specific network perspective.

Bioinformatic Analysis Reveals Phosphodiesterase 4D-Interacting Protein as a Key Frontal Cortex Dementia Switch Gene

The mechanisms that initiate dementia are poorly understood and there are currently no treatments that can slow their progression. The identification of key genes and molecular pathways that may trigger dementia should help reveal potential therapeutic reagents. In this study, SWItch Miner software was used to identify phosphodiesterase 4D-interacting protein as a key factor that may lead to the development of Alzheimer’s disease, vascular dementia, and frontotemporal dementia. Inflammation, PI3K-AKT, and ubiquitin-mediated proteolysis were identified as the main pathways that are dysregulated in these dementias. All of these dementias are regulated by 12 shared transcription factors. Protein–chemical interaction network analysis of dementia switch genes revealed that valproic acid may be neuroprotective for these dementias. Collectively, we identified shared and unique dysregulated gene expression, pathways and regulatory factors among dementias. New key mechanisms that lead to the development of dementia were revealed and it is expected that these data will advance personalized medicine for patients.

Prostate cancer screening research can benefit from network medicine: an emerging awareness

Up to date, screening for prostate cancer (PCa) remains one of the most appealing but also a very controversial topics in the urological community. PCa is the second most common cancer in men worldwide and it is universally acknowledged as a complex disease, with a multi-factorial etiology. The pathway of PCa diagnosis has changed dramatically in the last few years, with the multiparametric magnetic resonance (mpMRI) playing a starring role with the introduction of the “MRI Pathway”. In this scenario the basic tenet of network medicine (NM) that sees the disease as perturbation of a network of interconnected molecules and pathways, seems to fit perfectly with the challenges that PCa early detection must face to advance towards a more reliable technique. Integration of tests on body fluids, tissue samples, grading/staging classification, physiological parameters, MR multiparametric imaging and molecular profiling technologies must be integrated in a broader vision of “disease” and its complexity with a focus on early signs. PCa screening research can greatly benefit from NM vision since it provides a sound interpretation of data and a common language, facilitating exchange of ideas between clinicians and data analysts for exploring new research pathways in a rational, highly reliable, and reproducible way.

Molecular networks in Network Medicine: Development and applications

Network Medicine applies network science approaches to investigate disease pathogenesis. Many different analytical methods have been used to infer relevant molecular networks, including protein-protein interaction networks, correlation-based networks, gene regulatory networks, and Bayesian networks. Network Medicine applies these integrated approaches to Omics Big Data (including genetics, epigenetics, transcriptomics, metabolomics, and proteomics) using computational biology tools and, thereby, has the potential to provide improvements in the diagnosis, prognosis, and treatment of complex diseases. We discuss briefly the types of molecular data that are used in molecular network analyses, survey the analytical methods for inferring molecular networks, and review efforts to validate and visualize molecular networks. Successful applications of molecular network analysis have been reported in pulmonary arterial hypertension, coronary heart disease, diabetes mellitus, chronic lung diseases, and drug development. Important knowledge gaps in Network Medicine include incompleteness of the molecular interactome, challenges in identifying key genes within genetic association regions, and limited applications to human diseases. This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Translational, Genomic, and Systems Medicine > Translational Medicine Analytical and Computational Methods > Analytical Methods Analytical and Computational Methods > Computational Methods.

Integrated transcriptomic correlation network analysis identifies COPD molecular determinants

Chronic obstructive pulmonary disease (COPD) is a complex and heterogeneous syndrome. Network-based analysis implemented by SWIM software can be exploited to identify key molecular switches - called "switch genes" - for the disease. Genes contributing to common biological processes or defining given cell types are usually co-regulated and co-expressed, forming expression network modules. Consistently, we found that the COPD correlation network built by SWIM consists of three well-characterized modules: one populated by switch genes, all up-regulated in COPD cases and related to the regulation of immune response, inflammatory response, and hypoxia (like TIMP1, HIF1A, SYK, LY96, BLNK and PRDX4); one populated by well-recognized immune signature genes, all up-regulated in COPD cases; one where the GWAS genes AGER and CAVIN1 are the most representative module genes, both down-regulated in COPD cases. Interestingly, 70% of AGER negative interactors are switch genes including PRDX4, whose activation strongly correlates with the activation of known COPD GWAS interactors SERPINE2, CD79A, and POUF2AF1. These results suggest that SWIM analysis can identify key network modules related to complex diseases like COPD.

A genome-wide long noncoding RNA CRISPRi screen identifies PRANCR as a novel regulator of epidermal homeostasis

Genome-wide association studies indicate that many disease susceptibility regions reside in non-protein-coding regions of the genome. Long noncoding RNAs (lncRNAs) are a major component of the noncoding genome, but their biological impacts are not fully understood. Here, we performed a CRISPR interference (CRISPRi) screen on 2263 epidermis-expressed lncRNAs and identified nine novel candidate lncRNAs regulating keratinocyte proliferation. We further characterized a top hit from the screen, progenitor renewal associated non-coding RNA (PRANCR), using RNA interference–mediated knockdown and phenotypic analysis in organotypic human tissue. PRANCR regulates keratinocyte proliferation, cell cycle progression, and clonogenicity. PRANCR-deficient epidermis displayed impaired stratification with reduced expression of differentiation genes that are altered in human skin diseases, including keratins 1 and 10, filaggrin, and loricrin. Transcriptome analysis showed that PRANCR controls the expression of 1136 genes, with strong enrichment for late cell cycle genes containing a CHR promoter element. In addition, PRANCR depletion led to increased levels of both total and nuclear CDKN1A (also known as p21), which is known to govern both keratinocyte proliferation and differentiation. Collectively, these data show that PRANCR is a novel lncRNA regulating epidermal homeostasis and identify other lncRNA candidates that may have roles in this process as well.

Computational identification of key genes that may regulate gene expression reprogramming in Alzheimer's patients

The dementia epidemic is likely to expand worldwide as the aging population continues to grow. A better understanding of the molecular mechanisms that lead to dementia is expected to reveal potentially modifiable risk factors that could contribute to the development of prevention strategies. Alzheimer's disease is the most prevalent form of dementia. Currently we only partially understand some of the pathophysiological mechanisms that lead to development of the disease in aging individuals. In this study, Switch Miner software was used to identify key switch genes in the brain whose expression may lead to the development of Alzheimer's disease. The results indicate that switch genes are enriched in pathways involved in the proteasome, oxidative phosphorylation, Parkinson's disease, Huntington's disease, Alzheimer's disease and metabolism in the hippocampus and posterior cingulate cortex. Network analysis identified the krupel like factor 9 (KLF9), potassium channel tetramerization domain 2 (KCTD2), Sp1 transcription factor (SP1) and chromodomain helicase DNA binding protein 1 (CHD1) as key transcriptional regulators of switch genes in the brain of AD patients. These transcriptions factors have been implicated in conditions associated with Alzheimer's disease, including diabetes, glucocorticoid signaling, stroke, and sleep disorders. The specific pathways affected reveal potential modifiable risk factors by lifestyle changes.

A paradigm shift in medicine: A comprehensive review of network-based approaches

Network medicine is a rapidly evolving new field of medical research, which combines principles and approaches of systems biology and network science, holding the promise to uncovering the causes and to revolutionize the diagnosis and treatments of human diseases. This new paradigm reflects the fact that human diseases are not caused by single molecular defects, but driven by complex interactions among a variety of molecular mediators. The complexity of these interactions embraces different types of information: from the cellular-molecular level of protein-protein interactions to correlational studies of gene expression and regulation, to metabolic and disease pathways up to drug-disease relationships. The analysis of these complex networks can reveal new disease genes and/or disease pathways and identify possible targets for new drug development, as well as new uses for existing drugs. In this review, we offer a comprehensive overview of network types and algorithms used in the framework of network medicine. This article is part of a Special Issue entitled: Transcriptional Profiles and Regulatory Gene Networks edited by Dr. Dr. Federico Manuel Giorgi and Dr. Shaun Mahony.

Inferring the Disease-Associated miRNAs Based on Network Representation Learning and Convolutional Neural Networks

Identification of disease-associated miRNAs (disease miRNAs) are critical for understanding etiology and pathogenesis. Most previous methods focus on integrating similarities and associating information contained in heterogeneous miRNA-disease networks. However, these methods establish only shallow prediction models that fail to capture complex relationships among miRNA similarities, disease similarities, and miRNA-disease associations. We propose a prediction method on the basis of network representation learning and convolutional neural networks to predict disease miRNAs, called CNNMDA. CNNMDA deeply integrates the similarity information of miRNAs and diseases, miRNA-disease associations, and representations of miRNAs and diseases in low-dimensional feature space. The new framework based on deep learning was built to learn the original and global representation of a miRNA-disease pair. First, diverse biological premises about miRNAs and diseases were combined to construct the embedding layer in the left part of the framework, from a biological perspective. Second, the various connection edges in the miRNA-disease network, such as similarity and association connections, were dependent on each other. Therefore, it was necessary to learn the low-dimensional representations of the miRNA and disease nodes based on the entire network. The right part of the framework learnt the low-dimensional representation of each miRNA and disease node based on non-negative matrix factorization, and these representations were used to establish the corresponding embedding layer. Finally, the left and right embedding layers went through convolutional modules to deeply learn the complex and non-linear relationships among the similarities and associations between miRNAs and diseases. Experimental results based on cross validation indicated that CNNMDA yields superior performance compared to several state-of-the-art methods. Furthermore, case studies on lung, breast, and pancreatic neoplasms demonstrated the powerful ability of CNNMDA to discover potential disease miRNAs.

Next Generation Networks: Featuring the Potential Role of Emerging Applications in Translational Oncology

Nowadays, networks are pervasively used as examples of models suitable to mathematically represent and visualize the complexity of systems associated with many diseases, including cancer. In the cancer context, the concept of network entropy has guided many studies focused on comparing equilibrium to disequilibrium (i.e., perturbed) conditions. Since these conditions reflect both structural and dynamic properties of network interaction maps, the derived topological characterizations offer precious support to conduct cancer inference. Recent innovative directions have emerged in network medicine addressing especially experimental omics approaches integrated with a variety of other data, from molecular to clinical and also electronic records, bioimaging etc. This work considers a few theoretically relevant concepts likely to impact the future of applications in personalized/precision/translational oncology. The focus goes to specific properties of networks that are still not commonly utilized or studied in the oncological domain, and they are: controllability, synchronization and symmetry. The examples here provided take inspiration from the consideration of metastatic processes, especially their progression through stages and their hallmark characteristics. Casting these processes into computational frameworks and identifying network states with specific modular configurations may be extremely useful to interpret or even understand dysregulation patterns underlying cancer, and associated events (onset, progression) and disease phenotypes.

BRAFV600E-mutant cancers display a variety of networks by SWIM analysis: prediction of vemurafenib clinical response

PURPOSE

Several studies have shown that different tumour types sharing a driver gene mutation do not respond uniformly to the same targeted agent. Our aim was to use an unbiased network-based approach to investigate this fundamental issue using BRAF^V600E mutant tumours and the BRAF inhibitor vemurafenib.

METHODS

We applied SWIM, a software able to identify putative regulatory (switch) genes involved in drastic changes to the cell phenotype, to gene expression profiles of different BRAF^V600E mutant cancers and their normal counterparts in order to identify the switch genes that could potentially explain the heterogeneity of these tumours' responses to vemurafenib.

RESULTS

We identified lung adenocarcinoma as the tumour with the highest number of switch genes (298) compared to its normal counterpart. By looking for switch genes encoding for kinases with homology sequences similar to known vemurafenib targets, we found that thyroid cancer and lung adenocarcinoma have a similar number of putative targetable switch gene kinases (5 and 6, respectively) whereas colorectal cancer has just one.

CONCLUSIONS

We are persuaded that our network analysis may aid in the comprehension of molecular mechanisms underlying the different responses to vemurafenib in BRAF^V600E mutant tumours.

Identification of Disease-miRNA Networks Across Different Cancer Types Using SWIM

MicroRNAs (miRNAs) are small noncoding RNAs (ncRNAs) involved in several biological processes and diseases. MiRNAs regulate gene expression at the posttranscriptional level, mostly downregulating their targets by binding specific regions of transcripts through imperfect sequence complementarity. Prediction of miRNA-binding sites is challenging, and target prediction algorithms are usually based on sequence complementarity. In the last years, it has been shown that by adding miRNA and protein coding gene expression, we are able to build tissue-, cell line-, or disease-specific networks improving our understanding of complex biological scenarios. In this chapter, we present an application of a recently published software named SWIM, that allows to identify key genes in a network of interactions by defining appropriate "roles" of genes according to their local/global positioning in the overall network. Furthermore, we show how the SWIM software can be used to build miRNA-disease networks, by applying the approach to tumor data obtained from The Cancer Genome Atlas (TCGA).

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.

Network-Based Approaches to Explore Complex Biological Systems towards Network Medicine

Network medicine relies on different types of networks: from the molecular level of protein–protein interactions to gene regulatory network and correlation studies of gene expression. Among network approaches based on the analysis of the topological properties of protein–protein interaction (PPI) networks, we discuss the widespread DIAMOnD (disease module detection) algorithm. Starting from the assumption that PPI networks can be viewed as maps where diseases can be identified with localized perturbation within a specific neighborhood (i.e., disease modules), DIAMOnD performs a systematic analysis of the human PPI network to uncover new disease-associated genes by exploiting the connectivity significance instead of connection density. The past few years have witnessed the increasing interest in understanding the molecular mechanism of post-transcriptional regulation with a special emphasis on non-coding RNAs since they are emerging as key regulators of many cellular processes in both physiological and pathological states. Recent findings show that coding genes are not the only targets that microRNAs interact with. In fact, there is a pool of different RNAs—including long non-coding RNAs (lncRNAs) —competing with each other to attract microRNAs for interactions, thus acting as competing endogenous RNAs (ceRNAs). The framework of regulatory networks provides a powerful tool to gather new insights into ceRNA regulatory mechanisms. Here, we describe a data-driven model recently developed to explore the lncRNA-associated ceRNA activity in breast invasive carcinoma. On the other hand, a very promising example of the co-expression network is the one implemented by the software SWIM (switch miner), which combines topological properties of correlation networks with gene expression data in order to identify a small pool of genes—called switch genes—critically associated with drastic changes in cell phenotype. Here, we describe SWIM tool along with its applications to cancer research and compare its predictions with DIAMOnD disease genes.

Large Scale Proteomic Data and Network-Based Systems Biology Approaches to Explore the Plant World

The investigation of plant organisms by means of data-derived systems biology approaches based on network modeling is mainly characterized by genomic data, while the potential of proteomics is largely unexplored. This delay is mainly caused by the paucity of plant genomic/proteomic sequences and annotations which are fundamental to perform mass-spectrometry (MS) data interpretation. However, Next Generation Sequencing (NGS) techniques are contributing to filling this gap and an increasing number of studies are focusing on plant proteome profiling and protein-protein interactions (PPIs) identification. Interesting results were obtained by evaluating the topology of PPI networks in the context of organ-associated biological processes as well as plant-pathogen relationships. These examples foreshadow well the benefits that these approaches may provide to plant research. Thus, in addition to providing an overview of the main-omic technologies recently used on plant organisms, we will focus on studies that rely on concepts of module, hub and shortest path, and how they can contribute to the plant discovery processes. In this scenario, we will also consider gene co-expression networks, and some examples of integration with metabolomic data and genome-wide association studies (GWAS) to select candidate genes will be mentioned.

Inverse Problems in Systems Biology: A Critical Review

Systems Biology may be assimilated to a symbiotic cyclic interplaying between the forward and inverse problems. Computational models need to be continuously refined through experiments and in turn they help us to make limited experimental resources more efficient. Every time one does an experiment we know that there will be some noise that can disrupt our measurements. Despite the noise certainly is a problem, the inverse problems already involve the inference of missing information, even if the data is entirely reliable. So the addition of a certain limited noise does not fundamentally change the situation but can be used to solve the so-called ill-posed problem, as defined by Hadamard. It can be seen as an extra source of information. Recent studies have shown that complex systems, among others the systems biology, are poorly constrained and ill-conditioned because it is difficult to use experimental data to fully estimate their parameters. For these reasons was born the concept of sloppy models, a sequence of models of increasing complexity that become sloppy in the limit of microscopic accuracy. Furthermore the concept of sloppy models contains also the concept of un-identifiability, because the models are characterized by many parameters that are poorly constrained by experimental data. Then a strategy needs to be designed to infer, analyze, and understand biological systems. The aim of this work is to provide a critical review to the inverse problems in systems biology defining a strategy to determine the minimal set of information needed to overcome the problems arising from dynamic biological models that generally may have many unknown, non-measurable parameters.

Cell cycle arrest through indirect transcriptional repression by p53: I have a DREAM

Activation of the p53 tumor suppressor can lead to cell cycle arrest. The key mechanism of p53-mediated arrest is transcriptional downregulation of many cell cycle genes. In recent years it has become evident that p53-dependent repression is controlled by the p53–p21–DREAM–E2F/CHR pathway (p53–DREAM pathway). DREAM is a transcriptional repressor that binds to E2F or CHR promoter sites. Gene regulation and deregulation by DREAM shares many mechanistic characteristics with the retinoblastoma pRB tumor suppressor that acts through E2F elements. However, because of its binding to E2F and CHR elements, DREAM regulates a larger set of target genes leading to regulatory functions distinct from pRB/E2F. The p53–DREAM pathway controls more than 250 mostly cell cycle-associated genes. The functional spectrum of these pathway targets spans from the G₁ phase to the end of mitosis. Consequently, through downregulating the expression of gene products which are essential for progression through the cell cycle, the p53–DREAM pathway participates in the control of all checkpoints from DNA synthesis to cytokinesis including G₁/S, G₂/M and spindle assembly checkpoints. Therefore, defects in the p53–DREAM pathway contribute to a general loss of checkpoint control. Furthermore, deregulation of DREAM target genes promotes chromosomal instability and aneuploidy of cancer cells. Also, DREAM regulation is abrogated by the human papilloma virus HPV E7 protein linking the p53–DREAM pathway to carcinogenesis by HPV. Another feature of the pathway is that it downregulates many genes involved in DNA repair and telomere maintenance as well as Fanconi anemia. Importantly, when DREAM function is lost, CDK inhibitor drugs employed in cancer treatment such as Palbociclib, Abemaciclib and Ribociclib can compensate for defects in early steps in the pathway upstream from cyclin/CDK complexes. In summary, the p53–p21–DREAM–E2F/CHR pathway controls a plethora of cell cycle genes, can contribute to cell cycle arrest and is a target for cancer therapy.

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.