Transcriptome analysis reveals genes associated with stem cell activation by physical exercise in the dentate gyrus of aged p16Ink4a knockout mice.

Throughout adulthood neural stem cells divide in neurogenic niches-the dentate gyrus of the hippocampus and the subventricular zone-producing progenitor cells and new neurons. Stem cells self-renew, thus preserving their pool. Furthermore, the number of stem/progenitor cells in the neurogenic niches decreases with age. We have previously demonstrated that the cyclin-dependent kinase inhibitor p16Ink4a maintains, in aged mice, the pool of dentate gyrus stem cells by preventing their activation after a neurogenic stimulus such as exercise (running). We showed that, although p16Ink4a ablation by itself does not activate stem/progenitor cells, exercise strongly induced stem cell proliferation in p16Ink4a knockout dentate gyrus, but not in wild-type. As p16Ink4a regulates stem cell self-renewal during aging, we sought to profile the dentate gyrus transcriptome from p16Ink4a wild-type and knockout aged mice, either sedentary or running for 12 days. By pairwise comparisons of differentially expressed genes and by correlative analyses through the DESeq2 software, we identified genes regulated by p16Ink4a deletion, either without stimulus (running) added, or following running. The p16Ink4a knockout basic gene signature, i.e., in sedentary mice, involves upregulation of apoptotic, neuroinflammation- and synaptic activity-associated genes, suggesting a reactive cellular state. Conversely, another set of 106 genes we identified, whose differential expression specifically reflects the pattern of proliferative response of p16 knockout stem cells to running, are involved in processes that regulate stem cell activation, such as synaptic function, neurotransmitter metabolism, stem cell proliferation control, and reactive oxygen species level regulation. Moreover, we analyzed the regulation of these stem cell-specific genes after a second running stimulus. Surprisingly, the second running neither activated stem cell proliferation in the p16Ink4a knockout dentate gyrus nor changed the expression of these genes, confirming that they are correlated to the stem cell reactivity to stimulus, a process where they may play a role regulating stem cell activation.

DeLA-Drug: A Deep Learning Algorithm for Automated Design of Druglike Analogues.

In this paper, we present a deep learning algorithm for automated design of druglike analogues (DeLA-Drug), a recurrent neural network (RNN) model composed of two long short-term memory (LSTM) layers and conceived for data-driven generation of similar-to-bioactive compounds. DeLA-Drug captures the syntax of SMILES strings of more than 1 million compounds belonging to the ChEMBL28 database and, by employing a new strategy called sampling with substitutions (SWS), generates molecules starting from a single user-defined query compound. Remarkably, the algorithm preserves druglikeness and synthetic accessibility of the known bioactive compounds present in the ChEMBL28 repository. The absence of any time-demanding fine-tuning procedure enables DeLA-Drug to perform a fast generation of focused libraries for further high-throughput screening and makes it a suitable tool for performing de novo design even in low-data regimes. To provide a concrete idea of its applicability, DeLA-Drug was applied to the cannabinoid receptor subtype 2 (CB2R), a known target involved in different pathological conditions such as cancer and neurodegeneration. DeLA-Drug, available as a free web platform (http://www.ba.ic.cnr.it/softwareic/deladrugportal/), can help medicinal chemists interested in generating analogues of compounds already available in their laboratories and, for this reason, good candidates for an easy and low-cost synthesis.

Transcriptome Analysis in a Mouse Model of Premature Aging of Dentate Gyrus: Rescue of Alpha-Synuclein Deficit by Virus-Driven Expression or by Running Restores the Defective Neurogenesis.

The dentate gyrus of the hippocampus and the subventricular zone are neurogenic niches where neural stem and progenitor cells replicate throughout life to generate new neurons. The Btg1 gene maintains the stem cells of the neurogenic niches in quiescence. The deletion of Btg1 leads to an early transient increase of stem/progenitor cells division, followed, however, by a decrease during adulthood of their proliferative capability, accompanied by apoptosis. Since a physiological decrease of neurogenesis occurs during aging, the Btg1 knockout mouse may represent a model of neural aging. We have previously observed that the defective neurogenesis of the Btg1 knockout model is rescued by the powerful neurogenic stimulus of physical exercise (running). To identify genes responsible for stem and progenitor cells maintenance, we sought here to find genes underlying this premature neural aging, and whose deregulated expression could be rescued by running. Through RNA sequencing we analyzed the transcriptomic profiles of the dentate gyrus isolated from Btg1 wild-type or Btg1 knockout adult (2-month-old) mice submitted to physical exercise or sedentary. In Btg1 knockout mice, 545 genes were deregulated, relative to wild-type, while 2081 genes were deregulated by running. We identified 42 genes whose expression was not only down-regulated in the dentate gyrus of Btg1 knockout, but was also counter-regulated to control levels by running in Btg1 knockout mice, vs. sedentary. Among these 42 counter-regulated genes, alpha-synuclein (Snca), Fos, Arc and Npas4 showed significantly greater differential regulation. These genes control neural proliferation, apoptosis, plasticity and memory and are involved in aging. In particular, Snca expression decreases during aging. We tested, therefore, whether an Snca-expressing lentivirus, by rescuing the defective Snca levels in the dentate gyrus of Btg1 knockout mice, could also reverse the aging phenotype, in particular the defective neurogenesis. We found that the exogenous expression of Snca reversed the Btg1 knockout-dependent decrease of stem cell proliferation as well as the increase of progenitor cell apoptosis. This indicates that Snca has a functional role in the process of neural aging observed in this model, and also suggests that Snca acts as a positive regulator of stem cell maintenance.

Structure-Based Prediction of hERG-Related Cardiotoxicity: A Benchmark Study.

Drug-induced blockade of the human ether-à-go-go-related gene (hERG) channel is today considered the main cause of cardiotoxicity in postmarketing surveillance. Hence, several ligand-based approaches were developed in the last years and are currently employed in the early stages of a drug discovery process for in silico cardiac safety assessment of drug candidates. Herein, we present the first structure-based classifiers able to discern hERG binders from nonbinders. LASSO regularized support vector machines were applied to integrate docking scores and protein-ligand interaction fingerprints. A total of 396 models were trained and validated based on: (i) high-quality experimental bioactivity information returned by 8337 curated compounds extracted from ChEMBL (version 25) and (ii) structural predictor data. Molecular docking simulations were performed using GLIDE and GOLD software programs and four different hERG structural models, namely, the recently published structures obtained by cryoelectron microscopy (PDB codes: 5VA1 and 7CN1) and two published homology models selected for comparison. Interestingly, some classifiers return performances comparable to ligand-based models in terms of area under the ROC curve (AUCMAX = 0.86 ± 0.01) and negative predictive values (NPVMAX = 0.81 ± 0.01), thus putting forward the herein proposed computational workflow as a valuable tool for predicting hERG-related cardiotoxicity without the limitations of ligand-based models, typically affected by low interpretability and a limited applicability domain. From a methodological point of view, our study represents the first example of a successful integration of docking scores and protein-ligand interaction fingerprints (IFs) through a support vector machine (SVM) LASSO regularized strategy. Finally, the study highlights the importance of using hERG structural models accounting for ligand-induced fit effects and allowed us to select the best-performing protein conformation (made available in the Supporting Information, SI) to be employed for a reliable structure-based prediction of hERG-related cardiotoxicity.

Cannabinoid Receptor Subtype 2 (CB2R) in a Multitarget Approach: Perspective of an Innovative Strategy in Cancer and Neurodegeneration.

The cannabinoid receptor subtype 2 (CB2R) represents an interesting and new therapeutic target for its involvement in the first steps of neurodegeneration as well as in cancer onset and progression. Several studies, focused on different types of tumors, report a promising anticancer activity induced by CB2R agonists due to their ability to reduce inflammation and cell proliferation. Moreover, in neuroinflammation, the stimulation of CB2R, overexpressed in microglial cells, exerts beneficial effects in neurodegenerative disorders. With the aim to overcome current treatment limitations, new drugs can be developed by specifically modulating, together with CB2R, other targets involved in such multifactorial disorders. Building on successful case studies of already developed multitarget strategies involving CB2R, in this Perspective we aim at prompting the scientific community to consider new promising target associations involving HDACs (histone deacetylases) and σ receptors by employing modern approaches based on molecular hybridization, computational polypharmacology, and machine learning algorithms.

Effect of different light-emitting diode (LED) irradiation on the shelf life and phytonutrient content of broccoli (Brassica oleracea L. var. italica).

Broccoli (Brassica oleracea L. var. italica) is largely cultivated in southern Italy. It is an important source of phytonutrients, which are partially lost during postharvest storage. The aim of this work was to evaluate the overall effect of five different low-intensity light-emitting diodes (LEDs) on the quality parameters of broccoli florets over 20 d of cold storage. The level of ascorbic acid, chlorophylls, carotenoids, phenolic compounds and soluble proteins, as well as colour analysis, were evaluated. Green LED increased the chlorophyll and ascorbic acid content; white, red and yellow LEDs had a positive effect on the redox status of broccoli. Globally, only green LED had a statistically significant positive effect when considering all analysed parameters and could be proposed to prolong the shelf life of broccoli during cold storage.

Proteome Response of Staphylococcus xylosus DSM 20266T to Anaerobiosis and Nitrite Exposure.

The viability and competitiveness of Staphylococcus xylosus in meat mostly depend on the ability to adapt itself to rapid oxygen and nutrients depletion during meat fermentation. The utilization of nitrite instead of oxygen becomes a successful strategy for this strain to improve its performance in anaerobiosis; however, metabolic pathways of this strain underlying this adaptation, are partially known. The aim of this study was to provide an overview on proteomic changes of S. xylosus DSM 20266T cultured under anaerobiosis and nitrite exposure. Thus, two different cultures of this strain, supplemented or not with nitrite, were in vitro incubated in aerobiosis and anaerobiosis monitoring cell viability, pH, oxidation reduction potential and nitrite content. Protein extracts, obtained from cells, collected as nitrite content was depleted, were analyzed by 2DE/MALDI-TOF/TOF-MS. Results showed that DSM 20266T growth was significantly sustained by nitrite in anaerobiosis, whereas no differences were found in aerobiosis. Accordingly, nitrite content was depleted after 13 h only in anaerobiosis. At this time of sampling, a comparative proteomic analysis showed 45 differentially expressed proteins. Most differences were found between aerobic and anaerobic cultures without nitrite; the induction of glycolytic enzymes and glyoxylate cycle, the reduction of TCA enzymes, and acetate fermentation were found in anaerobiosis to produce ATP and maintain the cell redox balance. In anaerobic cultures the nitrite supplementation partially restored TCA cycle, and reduced the amount of glycolytic enzymes. These results were confirmed by phenotypic microarray that, for the first time, was carried out on cell previously adapted at the different growth conditions. Overall, metabolic changes were similar between aerobiosis and anaerobiosis NO2-adapted cells, whilst cells grown under anaerobiosis showed different assimilation profiles by confirming proteomic data; indeed, these latter extensively assimilated substrates addressed at both supplying glucose for glycolysis or fueling alternative pathways to TCA cycle. In conclusion, metabolic pathways underlying the ability of S. xylosus to adapt itself to oxygen starvation were revealed; the addition of nitrite allowed S. xylosus to take advantage of nitrite to this condition, restoring some metabolic pathway underlying aerobic behavior of the strain.

Combined microRNA and mRNA expression analysis in pediatric multiple sclerosis: an integrated approach to uncover novel pathogenic mechanisms of the disease.

Multiple sclerosis (MS) is a complex disease of the CNS that usually affects young adults, although 3-5% of cases are diagnosed in childhood and adolescence (hence called pediatric MS, PedMS). Genetic predisposition, among other factors, seems to contribute to the risk of the onset, in pediatric as in adult ages, but few studies have investigated the genetic 'environmentally naïve' load of PedMS. The main goal of this study was to identify circulating markers (miRNAs), target genes (mRNAs) and functional pathways associated with PedMS; we also verified the impact of miRNAs on clinical features, i.e. disability and cognitive performances. The investigation was performed in 19 PedMS and 20 pediatric controls (PCs) using a High-Throughput Next-generation Sequencing (HT-NGS) approach followed by an integrated bioinformatics/biostatistics analysis. Twelve miRNAs were significantly upregulated (let-7a-5p, let-7b-5p, miR-25-3p, miR-125a-5p, miR-942-5p, miR-221-3p, miR-652-3p, miR-182-5p, miR-185-5p, miR-181a-5p, miR-320a, miR-99b-5p) and 1 miRNA was downregulated (miR-148b-3p) in PedMS compared with PCs. The interactions between the significant miRNAs and their targets uncovered predicted genes (i.e. TNFSF13B, TLR2, BACH2, KLF4) related to immunological functions, as well as genes involved in autophagy-related processes (i.e. ATG16L1, SORT1, LAMP2) and ATPase activity (i.e. ABCA1, GPX3). No significant molecular profiles were associated with any PedMS demographic/clinical features. Both miRNAs and mRNA expressions predicted the phenotypes (PedMS-PC) with an accuracy of 92% and 91%, respectively. In our view, this original strategy of contemporary miRNA/mRNA analysis may help to shed light in the genetic background of the disease, suggesting further molecular investigations in novel pathogenic mechanisms.

Functional Implications of MicroRNAs in Crohn's Disease Revealed by Integrating MicroRNA and Messenger RNA Expression Profiling.

Crohn's disease (CD) is a debilitating inflammatory bowel disease (IBD) that emerges due to the influence of genetic and environmental factors. microRNAs (miRNAs) have been identified in the tissue and sera of IBD patients and may play an important role in the induction of IBD. Our study aimed to identify differentially expressed miRNAs and miRNAs with the ability to alter transcriptome activity by comparing inflamed tissue samples with their non-inflamed counterparts. We studied changes in miRNA-mRNA interactions associated with CD by examining their differential co-expression relative to normal mucosa from the same patients. Correlation changes between the two conditions were incorporated into scores of predefined gene sets to identify biological processes with altered miRNA-mediated control. Our study identified 28 miRNAs differentially expressed (p-values < 0.01), of which 14 are up-regulated. Notably, our differential co-expression analysis highlights microRNAs (i.e., miR-4284, miR-3194 and miR-21) that have known functional interactions with key mechanisms implicated in IBD. Most of these miRNAs cannot be detected by differential expression analysis that do not take into account miRNA-mRNA interactions. The identification of differential miRNA-mRNA co-expression patterns will facilitate the investigation of the miRNA-mediated molecular mechanisms underlying CD pathogenesis and could suggest novel drug targets for validation.

The Immune Landscapes of Polypoid and Nonpolypoid Precancerous Colorectal Lesions.

Little is known about the immunoediting process in precancerous lesions. We explored this aspect of benign colorectal adenomas with a descriptive analysis of the immune pathways and immune cells whose regulation is linked to the morphology and size of these lesions. Two series of polypoid and nonpolypoid colorectal adenomas were used in this study: 1) 84 samples (42 lesions, each with matched samples of normal mucosa) whose gene expression data were used to quantify the tumor morphology- and size-related dysregulation of immune pathways collected in the Molecular Signature Database, using Gene Set Enrichment Analysis; 2) 40 other lesions examined with immunohistochemistry to quantify the presence of immune cells in the stromal compartment. In the analysis of transcriptomic data, 429 immune pathways displayed significant differential regulation in neoplasms of different morphology and size. Most pathways were significantly upregulated or downregulated in polypoid lesions versus nonpolypoid lesions (regardless of size). Differential pathway regulation associated with lesion size was observed only in polypoid neoplasms. These findings were mirrored by tissue immunostaining with CD4, CD8, FOXP3, MHC-I, CD68, and CD163 antibodies: stromal immune cell counts (mainly T lymphocytes and macrophages) were significantly higher in polypoid lesions. Certain markers displayed significant size-related differences regardless of lesion morphology. Multivariate analysis of variance showed that the marker panel clearly discriminated between precancerous lesions of different morphologies and sizes. Statistical analysis of immunostained cell counts fully support the results of the transcriptomic data analysis: the density of infiltration of most immune cells in the stroma of polypoid precancerous lesions was significantly higher than that observed in nonpolypoid lesions. Large neoplasms also have more immune cells in their stroma than small lesions. Immunoediting in precancerous colorectal tumors may vary with lesion morphology and stage of development, and this variability could influence a given lesion's trajectory to cancer.

Meta-Analysis of Differential Connectivity in Gene Co-Expression Networks in Multiple Sclerosis.

Differential gene expression analyses to investigate multiple sclerosis (MS) molecular pathogenesis cannot detect genes harboring genetic and/or epigenetic modifications that change the gene functions without affecting their expression. Differential co-expression network approaches may capture changes in functional interactions resulting from these alterations. We re-analyzed 595 mRNA arrays from publicly available datasets by studying changes in gene co-expression networks in MS and in response to interferon (IFN)-ß treatment. Interestingly, MS networks show a reduced connectivity relative to the healthy condition, and the treatment activates the transcription of genes and increases their connectivity in MS patients. Importantly, the analysis of changes in gene connectivity in MS patients provides new evidence of association for genes already implicated in MS by single-nucleotide polymorphism studies and that do not show differential expression. This is the case of amiloride-sensitive cation channel 1 neuronal (ACCN1) that shows a reduced number of interacting partners in MS networks, and it is known for its role in synaptic transmission and central nervous system (CNS) development. Furthermore, our study confirms a deregulation of the vitamin D system: among the transcription factors that potentially regulate the deregulated genes, we find TCF3 and SP1 that are both involved in vitamin D3-induced p27Kip1 expression. Unveiling differential network properties allows us to gain systems-level insights into disease mechanisms and may suggest putative targets for the treatment.

Systematic analysis of circadian genes using genome-wide cDNA microarrays in the inflammatory bowel disease transcriptome.

Simultaneous analysis of the transcripts of thousands of genes by cDNA microarrays allows the identification of genetic regulatory mechanisms involved in disease pathophysiology. The circadian clock circuitry controls essential cell processes and the functioning of organ systems, which are characterized by rhythmic variations with 24-hour periodicity. The derangement of these processes is involved in the basic mechanisms of inflammatory, metabolic, degenerative and neoplastic diseases. We evaluated by genome-wide cDNA microarray analysis the transcriptome of endoscopic mucosal biopsies of patients with inflammatory bowel diseases (IBD) focusing on the expression of circadian genes in Crohn's disease (CD) and ulcerative colitis (UC). Twenty-nine IBD patients (15 with CD and 14 with UC) were enrolled and mucosal biopsies were sampled at either inflamed or adjacent non-inflamed areas of the colon. A total of 150 circadian genes involved in pathways controlling crucial cell processes and tissue functions were investigated. In CD specimens 50 genes were differentially expressed, and 21 genes resulted up-regulated when compared to healthy colonic mucosa. In UC specimens 50 genes were differentially expressed, and 27 genes resulted up-regulated when compared to healthy colonic mucosa. Among the core clock genes ARNTL2 and RORA were up-regulated, while CSNK2B, NPAS2, PER1 and PER3 were down-regulated in CD specimens. Conversely, ARNTL2, CRY1, CSNK1E, RORA and TIPIN were up-regulated, while NR1D2 and PER3 were down-regulated in UC. In conclusion, in CD and UC patients there are differences in the expression of circadian genes between normal and diseased intestinal mucosa. The deregulated genes evidenced by transcriptome analysis in the major IBDs may play a crucial role in the pathophysiological mechanisms and may suggest novel therapeutic approaches.

Genome-wide Pathway Analysis Using Gene Expression Data of Colonic Mucosa in Patients with Inflammatory Bowel Disease.

BACKGROUND: Ulcerative colitis (UC) and Crohn's disease (CD) share some pathogenetic features. To provide new steps on the role of altered gene expression, and the involvement of gene networks, in the pathogenesis of these diseases, we performed a genome-wide analysis in 15 patients with CD and 14 patients with UC by comparing the RNA from inflamed and noninflamed colonic mucosa. METHODS: Two hundred ninety-eight differentially expressed genes in CD and 520 genes in UC were identified. By bioinformatic analyses, 34 pathways for CD, 6 of them enriched in noninflamed and 28 in inflamed tissues, and 19 pathways for UC, 17 in noninflamed and 2 in inflamed tissues, were also highlighted. RESULTS: In CD, the pathways included genes associated with cytokines and cytokine receptors connection, response to external stimuli, activation of cell proliferation or differentiation, cell migration, apoptosis, and immune regulation. In UC, the pathways were associated with genes related to metabolic and catabolic processes, biosynthesis and interconversion processes, leukocyte migration, regulation of cell proliferation, and epithelial-to-mesenchymal transition. CONCLUSIONS: In UC, the pattern of inflammation of colonic mucosa is due to a complex interaction network between host, gut microbiome, and diet, suggesting that bacterial products or endogenous synthetic/catabolic molecules contribute to impairment of the immune response, to breakdown of epithelial barrier, and to enhance the inflammatory process. In patients with CD, genes encoding a large variety of proteins, growth factors, cytokines, chemokines, and adhesion molecules may lead to uncontrolled inflammation with ensuing destruction of epithelial cells, inappropriate stimulation of antimicrobial and T cells differentiation, and inflammasome events.

Karyopherins: potential biological elements involved in the delayed graft function in renal transplant recipients.

BACKGROUND: Immediately after renal transplantation, patients experience rapid and significant improvement of their clinical conditions and undergo considerable systemic and cellular modifications. However, some patients present a slow recovery of the renal function commonly defined as delayed graft function (DGF). Although clinically well characterized, the molecular mechanisms underlying this condition are not totally defined, thus, we are currently missing specific clinical markers to predict and to make early diagnosis of this event. METHODS: We investigated, using a pathway analysis approach, the transcriptomic profile of peripheral blood mononuclear cells (PBMC) from renal transplant recipients with DGF and with early graft function (EGF), before (T0) and 24 hours (T24) after transplantation. RESULTS: Bioinformatics/statistical analysis showed that 15 pathways (8 up-regulated and 7 down-regulated) and 11 pathways (5 up-regulated and 6 down-regulated) were able to identify DGF patients at T0 and T24, respectively. Interestingly, the most up-regulated pathway at both time points was NLS-bearing substrate import into nucleus, which includes genes encoding for several subtypes of karyopherins, a group of proteins involved in nucleocytoplasmic transport. Signal transducers and activators of transcription (STAT) utilize karyopherins-alpha (KPNA) for their passage from cytoplasm into the nucleus. In vitro functional analysis demonstrated that in PBMCs of DGF patients, there was a significant KPNA-mediated nuclear translocation of the phosphorylated form of STAT3 (pSTAT3) after short-time stimulation (2 and 5 minutes) with interleukin-6. CONCLUSIONS: Our study suggests the involvement, immediately before transplantation, of karyopherin-mediated nuclear transport in the onset and development of DGF. Additionally, it reveals that karyopherins could be good candidates as potential DGF predictive clinical biomarkers and targets for pharmacological interventions in renal transplantation. However, because of the low number of patients analyzed and some methodological limitations, additional studies are needed to validate and to better address these points.

Loss of connectivity in cancer co-expression networks.

Differential gene expression profiling studies have lead to the identification of several disease biomarkers. However, the oncogenic alterations in coding regions can modify the gene functions without affecting their own expression profiles. Moreover, post-translational modifications can modify the activity of the coded protein without altering the expression levels of the coding gene, but eliciting variations to the expression levels of the regulated genes. These considerations motivate the study of the rewiring of networks co-expressed genes as a consequence of the aforementioned alterations in order to complement the informative content of differential expression. We analyzed 339 mRNAomes of five distinct cancer types to find single genes that presented co-expression patterns strongly differentiated between normal and tumor phenotypes. Our analysis of differentially connected genes indicates the loss of connectivity as a common topological trait of cancer networks, and unveils novel candidate cancer genes. Moreover, our integrated approach that combines the differential expression together with the differential connectivity improves the classic enrichment pathway analysis providing novel insights on putative cancer gene biosystems not still fully investigated.

A comparative study of covariance selection models for the inference of gene regulatory networks.

MOTIVATION: The inference, or 'reverse-engineering', of gene regulatory networks from expression data and the description of the complex dependency structures among genes are open issues in modern molecular biology. RESULTS: In this paper we compared three regularized methods of covariance selection for the inference of gene regulatory networks, developed to circumvent the problems raising when the number of observations n is smaller than the number of genes p. The examined approaches provided three alternative estimates of the inverse covariance matrix: (a) the 'PINV' method is based on the Moore-Penrose pseudoinverse, (b) the 'RCM' method performs correlation between regression residuals and (c) 'ℓ(2C)' method maximizes a properly regularized log-likelihood function. Our extensive simulation studies showed that ℓ(2C) outperformed the other two methods having the most predictive partial correlation estimates and the highest values of sensitivity to infer conditional dependencies between genes even when a few number of observations was available. The application of this method for inferring gene networks of the isoprenoid biosynthesis pathways in Arabidopsis thaliana allowed to enlighten a negative partial correlation coefficient between the two hubs in the two isoprenoid pathways and, more importantly, provided an evidence of cross-talk between genes in the plastidial and the cytosolic pathways. When applied to gene expression data relative to a signature of HRAS oncogene in human cell cultures, the method revealed 9 genes (p-value<0.0005) directly interacting with HRAS, sharing the same Ras-responsive binding site for the transcription factor RREB1. This result suggests that the transcriptional activation of these genes is mediated by a common transcription factor downstream of Ras signaling. AVAILABILITY: Software implementing the methods in the form of Matlab scripts are available at: http://users.ba.cnr.it/issia/iesina18/CovSelModelsCodes.zip.

Molecular pathways undergoing dramatic transcriptomic changes during tumor development in the human colon.

BACKGROUND: The malignant transformation of precancerous colorectal lesions involves progressive alterations at both the molecular and morphologic levels, the latter consisting of increases in size and in the degree of cellular atypia. Analyzing preinvasive tumors of different sizes can therefore shed light on the sequence of these alterations. METHODS: We used a molecular pathway-based approach to analyze transcriptomic profiles of 59 colorectal tumors representing early and late preinvasive stages and the invasive stage of tumorigenesis. Random set analysis was used to identify biological pathways enriched for genes differentially regulated in tumors (compared with 59 samples of normal mucosa). RESULTS: Of the 880 canonical pathways we investigated, 112 displayed significant tumor-related upregulation or downregulation at one or more stages of tumorigenesis. This allowed us to distinguish between pathways whose dysregulation is probably necessary throughout tumorigenesis and those whose involvement specifically drives progression from one stage to the next. We were also able to pinpoint specific changes within each gene set that seem to play key roles at each transition. The early preinvasive stage was characterized by cell-cycle checkpoint activation triggered by DNA replication stress and dramatic downregulation of basic transmembrane signaling processes that maintain epithelial/stromal homeostasis in the normal mucosa. In late preinvasive lesions, there was also downregulation of signal transduction pathways (e.g., those mediated by G proteins and nuclear hormone receptors) involved in cell differentiation and upregulation of pathways governing nuclear envelope dynamics and the G2>M transition in the cell cycle. The main features of the invasive stage were activation of the G1>S transition in the cell cycle, upregulated expression of tumor-promoting microenvironmental factors, and profound dysregulation of metabolic pathways (e.g., increased aerobic glycolysis, downregulation of pathways that metabolize drugs and xenobiotics). CONCLUSIONS: Our analysis revealed specific pathways whose dysregulation might play a role in each transition of the transformation process. This is the first study in which such an approach has been used to gain further insights into colorectal tumorigenesis. Therefore, these data provide a launchpad for further exploration of the molecular characterization of colorectal tumorigenesis using systems biology approaches.

The expression of leucine-rich repeat gene family members in colorectal cancer.

This study was conducted to evaluate the association of the leucine-rich repeat (LRR) gene family with colorectal cancer (CRC). The expression of members of the LRR gene family were analyzed in 17 CRC specimens and in 59 healthy colorectal tissues by using Human Exon1.0ST microarray, and in 25 CRC specimens and 32 healthy colorectal tissues by U133Plus2.0 microarray. An association was found for 25 genes belonging to the plant-specific (PS) class of LRR genes (P = 0.05 for Exon1.0 ST and P = 0.04 for U133Plus2.0). In both data-sets, in CRC, we found down-regulation of SHOC2 (P < 0.00003) and LRRC28 (P < 0.01) and up-regulation of LRSAM1 (P < 0.000001), while up-regulation of MFHAS1 (P = 0.0005) and down-regulation of WDFY3 (P = 0.026) were found only in the Exon1.0 ST data-set. The PS LLR gene class encodes proteins that activate immune cells and might play a key role in programmed cell death and autophagy. SHOC2 and LRRC28 genes involved in RAS-mediated signaling, which hinders nutrient deprivation-induced autophagy, might be a possible link between the negative control of autophagy and tumorigenesis.

Activated innate immunity and the involvement of CX3CR1-fractalkine in promoting hematuria in patients with IgA nephropathy.

A hallmark of immunoglobulin A nephropathy (IgAN) is episodes of gross hematuria coinciding with mucosal infections that can represent the disease-triggering event. Here we performed a whole genomic screen of IgAN patients during gross hematuria to clarify the link between mucosal antigens and glomerular hematuria. Modulated genes showed a clear involvement of the intracellular interferon signaling, antigen-presenting pathway, and the immunoproteasome. The mRNA and protein level of the chemokine receptor characterizing cytotoxic effector lymphocytes, CX3CR1, was upregulated. In vitro antigenic stimulation of peripheral blood mononuclear cells from IgAN patients, healthy blood donors, and other nephropathies with microscopic hematuria showed that only in IgAN patients was CX3CR1 enhanced in a dose-dependent manner. A significantly higher amount of glomerular and urinary fractalkine, the only ligand of CX3CR1, was also found in IgAN patients with recurrent episodes of gross hematuria compared with other patients with microscopic or no hematuria. This suggests a predisposition for cytotoxic cell extravasation only in patients with recurrent gross hematuria. Thus, we found a defect in antigen handling in peripheral blood mononuclear cells of IgAN patients with a specific increase of CX3CR1. This constitutive upregulation of glomerular and urinary fractalkine suggests an involvement of the CX3CR1-fractalkine axis in the exacerbation of gross hematuria.

Integrating genetic and gene expression evidence into genome-wide association analysis of gene sets.

Single variant or single gene analyses generally account for only a small proportion of the phenotypic variation in complex traits. Alternatively, gene set or pathway association analyses are playing an increasingly important role in uncovering genetic architectures of complex traits through the identification of systematic genetic interactions. Two dominant paradigms for gene set analyses are association analyses based on SNP genotypes and those based on gene expression profiles. However, gene-disease association can manifest in many ways, such as alterations of gene expression, genotype, and copy number; thus, an integrative approach combining multiple forms of evidence can more accurately and comprehensively capture pathway associations. We have developed a single statistical framework, Gene Set Association Analysis (GSAA), that simultaneously measures genome-wide patterns of genetic variation and gene expression variation to identify sets of genes enriched for differential expression and/or trait-associated genetic markers. Simulation studies illustrate that joint analyses of genomic data increase the power to detect real associations when compared with gene set methods that use only one genomic data type. The analysis of two human diseases, glioblastoma and Crohn's disease, detected abnormalities in previously identified disease-associated pathways, such as pathways related to PI3K signaling, DNA damage response, and the activation of NFKB. In addition, GSAA predicted novel pathway associations, for example, differential genetic and expression characteristics in genes from the ABC transporter family in glioblastoma and from the HLA system in Crohn's disease. These demonstrate that GSAA can help uncover biological pathways underlying human diseases and complex traits.