Search

CD8+T cells are a major prognostic determinant in solid tumors, including colorectal cancer (CRC). However, understanding how the interplay between different immune cells impacts on clinical outcome is still in its infancy. Here, we describe that the interaction of tumor infiltrating neutrophils expressing high levels of CD15 with CD8+T effector memory cells (TEM) correlates with tumor progression. Mechanistically, stromal cell-derived factor-1 (CXCL12/SDF-1) promotes the retention of neutrophils within tumors, increasing the crosstalk with CD8+T cells. As a consequence of the contact-mediated interaction with neutrophils, CD8+T cells are skewed to produce high levels of GZMK, which in turn decreases E-cadherin on the intestinal epithelium and favors tumor progression. Overall, our results highlight the emergence of GZMKhighCD8+TEMin non-metastatic CRC tumors as a hallmark driven by the interaction with neutrophils, which could implement current patient stratification and be targeted by novel therapeutics.

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

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.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.We found that monoallelic Nbn deletion increases SmoA1-dependent MB incidence. In addition to a defective DDR, NbnOur study indicates that Nbn is haploinsufficient for SHH-MB development whereas full Nbn

Alzheimer’s disease (AD) is the most common neurodegenerative disease that currently lacks available effective therapy. Thus, identifying novel molecular biomarkers for diagnosis and treatment of AD is urgently demanded. In this study, we exploited tools and concepts of the emerging research area of Network Medicine to unveil a novel putative disease gene signature associated with AD. We proposed a new pipeline, which combines the strengths of two consolidated algorithms of the Network Medicine: DIseAse MOdule Detection (DIAMOnD), designed to predict new disease-associated genes within the human interactome network; and SWItch Miner (SWIM), designed to predict important (switch) genes within the co-expression network. Our integrated computational analysis allowed us to enlarge the set of the known disease genes associated to AD with additional 14 genes that may be proposed as new potential diagnostic biomarkers and therapeutic targets for AD phenotype.

Autism spectrum disorders (ASDs) is a multifactorial neurodevelopmental disorder. The communication between the gastrointestinal (GI) tract and the central nervous system seems driven by gut microbiota (GM). Herein, we provide GM profiling, considering GI functional symptoms, neurological impairment, and dietary habits. Forty-one and 35 fecal samples collected from ASD and neurotypical children (CTRLs), respectively, (age range, 3–15 years) were analyzed by 16S targeted-metagenomics (the V3–V4 region) and inflammation and permeability markers (i.e., sIgA, zonulin lysozyme), and then correlated with subjects’ metadata. Our ASD cohort was characterized as follows: 30/41 (73%) with GI functional symptoms; 24/41 (58%) picky eaters (PEs), with one or more dietary needs, including 10/41 (24%) with food selectivity (FS); 36/41 (88%) presenting high and medium autism severity symptoms (HMASSs). Among the cohort with GI symptoms, 28/30 (93%) showed HMASSs, 17/30 (57%) were picky eaters and only 8/30 (27%) with food selectivity. The remaining 11/41 (27%) ASDs without GI symptoms that were characterized by HMASS for 8/11 (72%) and 7/11 (63%) were picky eaters. GM ecology was investigated for the overall ASD cohort versus CTRLs; ASDs with GI and without GI, respectively, versus CTRLs; ASD with GI versus ASD without GI; ASDs with HMASS versus low ASSs; PEs versus no-PEs; and FS versus absence of FS. In particular, the GM of ASDs, compared to CTRLs, was characterized by the increase of Proteobacteria, Bacteroidetes, Rikenellaceae, Pasteurellaceae, Klebsiella, Bacteroides, Roseburia, Lactobacillus, Prevotella, Sutterella, Staphylococcus, and Haemophilus. Moreover, Sutterella, Roseburia and Fusobacterium were associated to ASD with GI symptoms compared to CTRLs. Interestingly, ASD with GI symptoms showed higher value of zonulin and lower levels of lysozyme, which were also characterized by differentially expressed predicted functional pathways. Multiple machine learning models classified correctly 80% overall ASDs, compared with CTRLs, based on Bacteroides, Lactobacillus, Prevotella, Staphylococcus, Sutterella, and Haemophilus features. In conclusion, in our patient cohort, regardless of the evaluation of many factors potentially modulating the GM profile, the major phenotypic determinant affecting the GM was represented by GI hallmarks and patients’ age.

Tumor contexture has emerged as a major prognostic determinant and tumor infiltrating CD8+ T cells have been associated with a better prognosis in several solid tumors, including early-stage colorectal cancer (CRC). However, the tumor immune infiltrate is highly heterogeneous and understanding how the interplay between different immune cell compartments impacts on the clinical outcome is still in its infancy.Here, we describe in a prospective cohort a novel CD8+ T effector memory population, which is characterized by high levels of Granzyme K (GZMKhigh CD8+ TEM) and is correlated with CD15high tumor infiltrating neutrophils. We provide both in vitro and in vivo evidence of the role of stromal cell-derived factor 1 (CXCL12/SDF-1) in driving functional changes on neutrophils at the tumor site, promoting their retention and increasing the crosstalk with CD8+ T cells. Mechanistically, as a consequence of the interaction with neutrophils, CD8+ T cells are skewed towards a CD8+ TEM phenotype and produce high levels of GZMK, which in turn decreases E-cadherin pathway. The correlations of GZMKhigh CD8+ TEM and neutrophils with both tumor progression in mice and early relapse in CRC patients demonstrate the role of GZMKhigh CD8+ TEM in promoting malignancy. Indeed, a gene signature defining GZMKhigh CD8+ TEM was associated with worse prognosis on a larger independent cohort of CRC patients and a similar analysis was extended to lung cancer (TCGA).Overall, our results highlight the emergence of GZMKhigh CD8+ TEM in early-stage CRC tumors as a hallmark driven by the interaction with neutrophils, which could implement current patient stratification and be targeted by novel therapeutics.

CD8

Despite advances in modern medicine that led to improvements in cardiovascular outcomes, cardiovascular disease (CVD) remains the leading cause of mortality and morbidity globally. Thus, there is an urgent need for new approaches to improve CVD drug treatments. As the development time and cost of drug discovery to clinical application are excessive, alternate strategies for drug development are warranted. Among these are included computational approaches based on omics data for drug repositioning, which have attracted increasing attention. In this work, we developed an adjusted similarity measure implemented by the algorithm SAveRUNNER to reposition drugs for cardiovascular diseases while, at the same time, considering the side effects of drug candidates. We analyzed nine cardiovascular disorders and two side effects. We formulated both disease disorders and side effects as network modules in the human interactome, and considered those drug candidates that are proximal to disease modules but far from side-effects modules as ideal. Our method provides a list of drug candidates for cardiovascular diseases that are unlikely to produce common, adverse side-effects. This approach incorporating side effects is applicable to other diseases, as well.

Breast cancer (BC) is a heterogeneous and complex disease characterized by different subtypes with distinct morphologies and clinical implications and for which new and effective treatment options are urgently demanded. The computational approaches recently developed for drug repurposing provide a very promising opportunity to offer tools that efficiently screen potential novel medical indications for various drugs that are already approved and used in clinical practice. Here, we started with disease-associated genes that were identified through a transcriptome-based analysis, which we used to predict potential repurposable drugs for various breast cancer subtypes by using an algorithm that we developed for drug repurposing called SAveRUNNER. Our findings were also in silico validated by performing a gene set enrichment analysis, which confirmed that most of the predicted repurposable drugs may have a potential treatment effect against breast cancer pathophenotypes.

Alzheimer’s disease (AD) is a neurologic disorder causing brain atrophy and the death of brain cells. It is a progressive condition marked by cognitive and behavioral impairment that significantly interferes with daily activities. AD symptoms develop gradually over many years and eventually become more severe, and no cure has been found yet to arrest this process. The present study is directed towards suggesting putative novel solutions and paradigms for fighting AD pathogenesis by exploiting new insights from network medicine and drug repurposing strategies. To identify new drug–AD associations, we exploited SAveRUNNER, a recently developed network-based algorithm for drug repurposing, which quantifies the vicinity of disease-associated genes to drug targets in the human interactome. We complemented the analysis with an in silico validation of the candidate compounds through a gene set enrichment analysis, aiming to determine if the modulation of the gene expression induced by the predicted drugs could be counteracted by the modulation elicited by the disease. We identified some interesting compounds belonging to the beta-blocker family, originally approved for treating hypertension, such as betaxolol, bisoprolol, and metoprolol, whose connection with a lower risk to develop Alzheimer’s disease has already been observed. Moreover, our algorithm predicted multi-kinase inhibitors such as regorafenib, whose beneficial effects were recently investigated for neuroinflammation and AD pathology, and mTOR inhibitors such as sirolimus, whose modulation has been associated with AD.

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.

Pools of RNA molecules can act as competing endogenous RNAs (ceRNAs) and indirectly alter their expression levels by competitively binding shared microRNAs. This ceRNA cross talk yields an additional posttranscriptional regulatory layer, which plays key roles in both physiological and pathological processes. MicroRNAs can act as decoys by binding multiple RNAs, as well as RNAs can act as ceRNAs by competing for binding multiple microRNAs, leading to many cross talk interactions that could favor significant large-scale effects in spite of the weakness of single interactions. Identifying and studying these extended ceRNA interaction networks could provide a global view of the fine-tuning gene regulation in a wide range of biological processes and tumor progressions. In this chapter, we review current progress of predicting ceRNA cross talk, by summarizing the most up-to-date databases, which collect computationally predicted and/or experimentally validated miRNA-target and ceRNA-ceRNA interactions, as well as the widespread computational methods for discovering and modeling possible evidences of ceRNA-ceRNA interaction networks. These methods can be grouped in two categories: statistics-based methods exploit multivariate analysis to build ceRNA networks, by considering the miRNA expression levels when evaluating miRNA sponging relationships; mathematical methods build deterministic or stochastic models to analyze and predict the behavior of ceRNA cross talk.

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 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 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 mutant tumours.

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.

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 $$l_1$$ l 1 -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.

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.

Non-nonavalent vaccine (9v) Human papillomavirus (HPV) types have been shown to have high prevalence among HIV-positive women. Here, 1444 cervical samples were tested for HPV DNA positivity. Co-infections of the 9v HPV types with other HPV types were evaluated. The HPV81 L1 and L2 genes were used to investigate the genetic variability of antigenic epitopes. HPV-positive samples were genotyped using the HPVCLART2 assay. The L1 and L2 protein sequences were analyzed using a self-optimized prediction method to predict their secondary structure. Co-occurrence probabilities of the 9v HPV types were calculated. Non9v types represented 49% of the HPV infections; 31.2% of the non9v HPV types were among the low-grade squamous intraepithelial lesion samples, and 27.3% among the high-grade squamous intraepithelial lesion samples, and several genotypes were low risk. The co-occurrence of 9v HPV types with the other genotypes was not correlated with the filogenetic distance. HPV81 showed an amino-acid substitution within the BC loop (N75Q) and the FGb loop (T315N). In the L2 protein, all of the mutations were located outside antigenic sites. The weak cross-protection of the 9v types suggests the relevance of a sustainable and effective screening program, which should be implemented by HPV DNA testing that does not include only high-risk types.

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.

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.

Background Triple-negative breast cancers (TNBCs) display poor prognosis, have a high risk of tumour recurrence, and exhibit high resistance to drug treatments. Based on their gene expression profiles, the majority of TNBCs are classified as basal-like breast cancers. Currently, there are not available widely-accepted prognostic markers to predict outcomes in basal-like subtype, so the selection of new prognostic indicators for this BC phenotype represents an unmet clinical challenge. Results Here, we attempted to address this challenging issue by exploiting a bioinformatics pipeline able to integrate transcriptomic, genomic, epigenomic, and clinical data freely accessible from public repositories. This pipeline starts from the application of the well-established network-based SWIM methodology on the transcriptomic data to unveil important (switch) genes in relation with a complex disease of interest. Then, survival and linear regression analyses are performed to associate the gene expression profiles of the switch genes with both the patients’ clinical outcome and the disease aggressiveness. This allows us to identify a prognostic gene signature that in turn is fed to the last step of the pipeline consisting of an analysis at DNA level, to investigate whether variations in the expression of identified prognostic switch genes could be related to genetic (copy number variations) or epigenetic (DNA methylation differences) alterations in their gene loci, or to the activities of transcription factors binding to their promoter regions. Finally, changes in the protein expression levels corresponding to the so far identified prognostic switch genes are evaluated by immunohistochemical staining results taking advantage of the Human Protein Atlas. Conclusion The application of the proposed pipeline on the dataset of The Cancer Genome Atlas (TCGA)-Breast Invasive Carcinoma (BRCA) patients affected by basal-like subtype led to an in silico recognition of a basal-like specific gene signature composed of 11 potential prognostic biomarkers to be further investigated.

The COVID-19 pandemic has overwhelmed the life and security of most of the world countries, and especially of the Western countries, without similar experiences in the recent past. In a first phase, the response of health systems and governments was disorganized, but then incisive, also driven by the fear of a new and dramatic phenomenon. In the second phase, several governments, including Italy, accepted the doctrine of “coexistence with the virus” by putting into practice a series of containment measures aimed at limiting the dramatic sanitary consequences while not jeopardizing the economic and social stability of the country. Here, we present a new mathematical approach to modeling the COVID-19 dynamics that accounts for typical evolution parameters (i.e., virus variants, vaccinations, containment measurements). Reproducing the COVID-19 epidemic spread is an extremely challenging task due to the low reliability of the available data, the lack of recurrent patterns, and the considerable amount and variability of the involved parameters. However, the adoption of fairly uniform criteria among the Italian regions enabled to test and optimize the model in various conditions leading to robust and interesting results. Although the regional variability is quite large and difficult to predict, we have retrospectively obtained reliable indications on which measures were the most appropriate to limit the transmissibility coefficients within detectable ranges for all the regions. To complicate matters further, the rapid spread of the English variant has upset contexts where the propagation of contagion was close to equilibrium conditions, decreeing success or failure of a certain measure. Finally, we assessed the effectiveness of the zone assignment criteria, highlighting how the reactivity of the measures plays a fundamental role in limiting the spread of the infection and thus the total number of deaths, the most important factor in assessing the success of epidemic management.

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.

Extensive molecular characterization of human colorectal cancer (CRC) via Next Generation Sequencing (NGS) indicated that genetic or epigenetic dysregulation of a relevant, but limited, number of molecular pathways typically occurs in this tumor. The molecular picture of the disease is significantly complicated by the frequent occurrence of individually rare genetic aberrations, which expand tumor heterogeneity. Inter- and intratumor molecular heterogeneity is very likely responsible for the remarkable individual variability in the response to conventional and target-driven first-line therapies, in metastatic CRC (mCRC) patients, whose median overall survival remains unsatisfactory. Implementation of an extensive molecular characterization of mCRC in the clinical routine does not yet appear feasible on a large scale, while multigene panel sequencing of most commonly mutated oncogene/oncosuppressor hotspots is more easily achievable. Here, we report that clinical multigene panel sequencing performed for anti-EGFR therapy predictive purposes in 639 formalin-fixed paraffin-embedded (FFPE) mCRC specimens revealed previously unknown pairwise mutation associations and a high proportion of cases carrying actionable gene mutations. Most importantly, a simple principal component analysis directed the delineation of a new molecular stratification of mCRC patients in eight groups characterized by non-random, specific mutational association patterns (MAPs), aggregating samples with similar biology. These data were validated on a The Cancer Genome Atlas (TCGA) CRC dataset. The proposed stratification may provide great opportunities to direct more informed therapeutic decisions in the majority of mCRC cases.

Chronic obstructive pulmonary disease (COPD) is a heterogeneous and complex syndrome. Network-based analysis implemented by SWIM software can be exploited to identify key molecular switches - called “switch genes” - for disease. Genes contributing to common biological processes or define given cell types are frequently co-regulated and co-expressed, giving rise to 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.

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

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 BRAFWe 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 BRAFWe 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.We are persuaded that our network analysis may aid in the comprehension of molecular mechanisms underlying the different responses to vemurafenib in BRAF

Systems biology focuses on the integration of experimental, mathematical and computational techniques to develop systemic views and predictive models of biological systems. In this perspective, the concept of network has been a powerful tool for the representation and the analysis of complex systems: during the last two decades, the so-called network biology approach has been fruitfully applied in many different biological areas, from gene regulation, to protein-protein interactions, to neural signals. Here, making no claim to completeness, we briefly account for the processes of reconstructing several among the most significant types of biological networks in molecular biology and neuroscience, as well as for some of the most promising methodologies applied in the recent field of network medicine.

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., Author summary The global pandemic caused by the new coronavirus SARS-CoV-2 (COVID-19) leads a compelling need to find new therapeutic options devoted to fight the disease progression in the short term and to prevent it from happening in the future. The strategy of reusing an ‘old drug’ for new therapeutic purposes (known as drug repurposing) appears as a powerful solution for emerging diseases, such as COVID-19, since it allows to shorten the time and reduce the cost compared to de novo drug discovery. In this context, we propose SAveRUNNER (Searching off-lAbel dRUg aNd NEtwoRk), a new network-medicine-based algorithm for drug repurposing, with the aim to offer a promising framework to efficiently screen potential novel medical indications for various drugs, which are already approved and used in clinical practice, against the new human coronavirus (2019-nCoV/SARS-CoV-2). Our computational analysis predicts several repurposable drugs, including some of the most rumored off-label drugs for COVID-19 treatments as well as new promising candidates worthy of further exploration.

Downregulation of microRNAs (miRNAs) is commonly observed in cancers and promotes tumorigenesis suggesting that miRNAs may function as tumor suppressors. However, the mechanism through which miRNAs are regulated in cancer, and the connection between oncogenes and miRNA biogenesis remain poorly understood. The TP53 tumor-suppressor gene is mutated in half of human cancers resulting in an oncogene with gain-of-function activities. Here we demonstrate that mutant p53 (mutp53) oncoproteins modulate the biogenesis of a subset of miRNAs in cancer cells inhibiting their post-transcriptional maturation. Interestingly, among these miRNAs several are also downregulated in human tumors. By confocal, co-immunoprecipitation and RNA-chromatin immunoprecipitation experiments, we show that endogenous mutp53 binds and sequesters RNA helicases p72/82 from the microprocessor complex, interfering with Drosha-pri-miRNAs association. In agreement with this, the overexpression of p72 leads to an increase of mature miRNAs levels. Moreover, functional experiments demonstrate the oncosuppressive role of mutp53-dependent miRNAs (miR-517a, - 519a, - 218, - 105). Our study highlights a previously undescribed mechanism by which mutp53 interferes with Drosha-p72/82 association leading, at least in part, to miRNA deregulation observed in cancer.

In the last decade noncoding RNAs (ncRNAs) have been extensively studied in several biological processes and human diseases including cancer. microRNAs (miRNAs) are the best-known class of ncRNAs. miRNAs are small ncRNAs of around 20-22 nucleotides (nt) and are crucial posttranscriptional regulators of protein coding genes. Recently, new classes of ncRNAs, longer than miRNAs have been discovered. Those include intergenic noncoding RNAs (lincRNAs) and circular RNAs (circRNAs). These novel types of ncRNAs opened a very exciting field in biology, leading researchers to discover new relationships between miRNAs and long noncoding RNAs (lncRNAs), which act together to control protein coding gene expression. One of these new discoveries led to the formulation of the "competing endogenous RNA (ceRNA) hypothesis." This hypothesis suggests that an lncRNA acts as a sponge for miRNAs reducing their expression and causing the upregulation of miRNA targets. In this chapter we first discuss some recent discoveries in this field showing the mutual regulation of miRNAs, lncRNAs, and protein-coding genes in cancer. We then discuss the general approaches for the study of ceRNAs and present in more detail a recent computational approach to explore the ability of lncRNAs to act as ceRNAs in human breast cancer that has been shown to be, among the others, the most precise and promising.

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. Electronic supplementary material The online version of this article (10.1186/s12859-018-2421-x) contains supplementary material, which is available to authorized users.

Glioblastoma, the most malignant brain cancer, contains self-renewing, stem-like cells that sustain tumor growth and therapeutic resistance. Identifying genes promoting stem-like cell differentiation might unveil targets for novel treatments. To detect them, here we apply SWIM - a software able to unveil genes (named switch genes) involved in drastic changes of cell phenotype - to public datasets of gene expression profiles from human glioblastoma cells. By analyzing matched pairs of stem-like and differentiated glioblastoma cells, SWIM identified 336 switch genes, potentially involved in the transition from stem-like to differentiated state. A subset of them was significantly related to focal adhesion and extracellular matrix and strongly down-regulated in stem-like cells, suggesting that they may promote differentiation and restrain tumor growth. Their expression in differentiated cells strongly correlated with the down-regulation of transcription factors like OLIG2, POU3F2, SALL2, SOX2, capable of reprogramming differentiated glioblastoma cells into stem-like cells. These findings were corroborated by the analysis of expression profiles from glioblastoma stem-like cell lines, the corresponding primary tumors, and conventional glioma cell lines. Switch genes represent a distinguishing feature of stem-like cells and we are persuaded that they may reveal novel potential therapeutic targets worthy of further investigation.

Glioblastoma, the most malignant brain cancer, contains self-renewing, stem-like cells that sustain tumor growth and therapeutic resistance. Identifying genes promoting stem-like cell differentiation might unveil targets for novel treatments. To detect them, here we apply SWIM – a software able to unveil genes (named switch genes) involved in drastic changes of cell phenotype – to public datasets of gene expression profiles from human glioblastoma cells. By analyzing matched pairs of stem-like and differentiated glioblastoma cells, SWIM identified 336 switch genes, potentially involved in the transition from stem-like to differentiated state. A subset of them was significantly related to focal adhesion and extracellular matrix and strongly down-regulated in stem-like cells, suggesting that they may promote differentiation and restrain tumor growth. Their expression in differentiated cells strongly correlated with the down-regulation of transcription factors like OLIG2, POU3F2, SALL2, SOX2, capable of reprogramming differentiated glioblastoma cells into stem-like cells. These findings were corroborated by the analysis of expression profiles from glioblastoma stem-like cell lines, the corresponding primary tumors, and conventional glioma cell lines. Switch genes represent a distinguishing feature of stem-like cells and we are persuaded that they may reveal novel potential therapeutic targets worthy of further investigation.

Background Circular RNAs are a class of endogenous RNAs with various functions in eukaryotic cells. Worthy of note, circular RNAs play a critical role in cancer. Currently, nothing is known about their role in head and neck squamous cell carcinoma (HNSCC). The identification of circular RNAs in HNSCC might become useful for diagnostic and therapeutic strategies in HNSCC. Results Using samples from 115 HNSCC patients, we find that circPVT1 is over-expressed in tumors compared to matched non-tumoral tissues, with particular enrichment in patients with TP53 mutations. circPVT1 up- and down-regulation determine, respectively, an increase and a reduction of the malignant phenotype in HNSCC cell lines. We show that circPVT1 expression is transcriptionally enhanced by the mut-p53/YAP/TEAD complex. circPVT1 acts as an oncogene modulating the expression of miR-497-5p and genes involved in the control of cell proliferation. Conclusions This study shows the oncogenic role of circPVT1 in HNSCC, extending current knowledge about the role of circular RNAs in cancer. Electronic supplementary material The online version of this article (doi:10.1186/s13059-017-1368-y) contains supplementary material, which is available to authorized users.

Broader comprehension of gene expression regulatory mechanisms can be gained from a global analysis of how transcription and degradation are coordinated to orchestrate complex cell responses. The role of messenger RNA (mRNA) turnover modulation in gene expression levels has become increasingly recognized. From such perspective, in this review we briefly illustrate how a simple but effective mathematical model of mRNA turnover and some experimental findings, may together shed light on the molecular mechanisms underpinning the major role of mRNA decay rates in shaping the kinetics of gene activation and repression.

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.

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.

Scientific Reports 7: Article number: 44797; published online: 20 March 2017; updated: 16 June 2017. In the HTML version of this Article, Figure 6 was incorrect. The correct Figure 6 appears below as Figure 1. This error has been corrected in the HTML version of the Article; the PDF version was correct at the time of publication.

SWItchMiner (SWIM) is a wizard-like software implementation of a procedure, previously described, able to extract information contained in complex networks. Specifically, SWIM allows unearthing the existence of a new class of hubs, called “fight-club hubs”, characterized by a marked negative correlation with their first nearest neighbors. Among them, a special subset of genes, called “switch genes”, appears to be characterized by an unusual pattern of intra- and inter-module connections that confers them a crucial topological role, interestingly mirrored by the evidence of their clinic-biological relevance. Here, we applied SWIM to a large panel of cancer datasets from The Cancer Genome Atlas, in order to highlight switch genes that could be critically associated with the drastic changes in the physiological state of cells or tissues induced by the cancer development. We discovered that switch genes are found in all cancers we studied and they encompass protein coding genes and non-coding RNAs, recovering many known key cancer players but also many new potential biomarkers not yet characterized in cancer context. Furthermore, SWIM is amenable to detect switch genes in different organisms and cell conditions, with the potential to uncover important players in biologically relevant scenarios, including but not limited to human cancer.

Direct-acting antiviral drugs (DAA) regimen improve the SVR rate. However, adverse effects often lead to therapy interruption. This underlines the importance to find some predictive parameters of response in order to consider the possibility of a shorter time of antiviral treatment in the appearance of adverse effects without affecting the success of the therapy.We aimed to examine the HCVAg kinetics in the early phase of treatment and its predictive value of SVR in patients undergoing TPV/Peg-IFN/RBV treatment.Twenty-three patients infected by HCV genotype 1 (1a n=11; 1b n=12) were included in this prospective study.At baseline the median Log of HCVAg concentration in RVR and EVR patients were 3.15 fmol/L and 3.45 fmol/L, respectively with no significant differences. The baseline median HCV-RNA to HCVAg ratio was 233.77, this ratio was significantly lower when measured on day 1 (27.52) and on day 6 (24.84) (p0.001). The two-tailed Fisher's exact test indicated that the SVR response is statistically significantly different in patients with detected HCVAg at week1 compared to patients with no detectable HCVAg (p=0.05). The sensitivity, specificity, and negative and positive predictive values (NPV, PPV) were 53.8, 87.5, 53.8 and 87.5%, respectively. The area under the ROC curve was 0.71 at day T6, the best cut-off of 3 fmol/L when evaluated with the HCVAg plasma concentration at day T6.Undetectable HCVAg in the early phase of TPV/Peg-IFN/RBV treatment could represent an important parameter for predicting SVR.

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

Topic: Psychological and stressful social situations, as well as factors of an unhealthy lifestyle (such as insufficient exercise, weight increase and/or social isolation) are potential causes of MUS (Medically Unexplained Symptoms) and obesity, usually overlooked by general practitioners. Obesity is often underestimated using BMI, since it does not reflect the real loss of fat mass after nutritional strategies. Scope: We analyzed the changes in fat mass (FM) and fat-free mass (FFM) in the presence of decreasing BMI, using a novel advanced bio-impedance device, BIA-ACC (Bioelectric Impedance Analyzer for Analisi Composizione Corporea) that measures numerous parameters of body composition. Methods: 109 patients were enrolled at the center for primary care for a routine check-up, all apparently free of disease but with problems of excess weight and MUS to varying degrees. The evaluation of the parameters BIA-ACC was performed at the first visit T0, and at follow-up T1, after about 2 months of a standard nutritional strategy. Results: BMI decreased in all patients. FFM, particularly skeletal muscle, was the 1st principal component (PCA1: 58%) that determines the well-being of the patients and a decrease of MUS. The 2nd component (PCA2: 22%) is the FM. We found 3 different patterns of patient response to BMI reduction: 75.23% showed loss of FM and minimally FFM; 17.43% showed loss of FM with preservation or increase FFM; 7.34% showed loss of FFM alone. Conclusions: Body composition analysis is a valuable non-invasive tool to monitor patients in the early stages of immune-metabolic dysregulation when clinical symptoms are not yet evident. Our results demonstrate the importance of the maintenance of FFM rather than the loss of FM alone during weight loss therapy. Preservation of skeletal muscle is essential to facilitate the stabilization of loss of only fat and thus to eliminate the MUS correlated to chronic stress.

Objectives Identifying the predictive factors of Sustained Virological Response (SVR) represents an important challenge in new interferon-based DAA therapies. Here, we analyzed the kinetics of antiviral response associated with a triple drug regimen, and the association between negative residual viral load at different time points during treatment. Methods Twenty-three HCV genotype 1 (GT 1a n = 11; GT1b n = 12) infected patients were included in the study. Linear Discriminant Analysis (LDA) was used to establish possible association between HCV RNA values at days 1 and 4 from start of therapy and SVR. Principal compo- nent analysis (PCA) was applied to analyze the correlation between HCV RNA slope and SVR. A ultrasensitive (US) method was established to measure the residual HCV viral load in those samples which resulted "detected

The scope of this review is, by briefly discussing some applications in the rapidly emerging field of protein contact network, to sketch an at least initial answer to the quest for a new "structural formula" language for proteins. This quest will be pursued by presenting side-by-side the different complex network invariants developed by graph theory and their protein counterparts.

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

The time of initiation of antiretroviral therapy in HIV-1 infected patients has a determinant effect on the viral dynamics. The question is, how far can the therapy be delayed? Is sooner always better? We resort to clinical data and to microsimulations to forecast the dynamics of the viral load at therapy interruption after prolonged antiretroviral treatment. A computational model previously evaluated, produces results that are statistically adherent to clinical data. In addition, it allows a finer grain analysis of the impact of the therapy initiation point to the disease course. We find a swift increase of the viral density as a function of the time of initiation of the therapy measured when the therapy is stopped. In particular there is a critical time delay with respect to the infection instant beyond which the therapy does not affect the viral rebound. Initiation of the treatment is beneficial because it can down-regulate the immune activation, hence limiting viral replication and spread.