Comprehensive analysis of the genome transcriptome and proteome landscapes of three tumor cell lines

The use of fluid-phase 3D printing to pattern alginate-gelatin hydrogel properties to guide cell growth and behaviourin vitro

Three-dimensional (3D) (bio)printing technology has boosted the advancement of the biomedical field. However, tissue engineering is an evolving field and (bio)printing biomimetic constructions for tissue formation is still a challenge. As a new methodology to facilitate the construction of more complex structures, we suggest the use of the fluid-phase 3D printing to pattern the scaffold's properties. The methodology consists of an exchangeable fluid-phase printing medium in which the constructions are fabricated and patterned during the printing process. Using the fluid-phase methodology, the biological and mechanical properties can be tailored promoting cell behaviour guidance and compartmentalization. In this study, we first assessed different formulations of alginate/gelatin to create a stable substrate capable to promote massive cell colonizationin vitroover time. Overall, formulations with lower gelatin content and 2-(N-morpholino)ethanesulfonic acid (MES) buffer as a solvent showed better stability under cell culture conditions and enhanced U2OS cell growth. Next, the fluid-phase showed better printing fidelity and resolution in comparison to air printing as it diminished the collapsing and the spread of the hydrogel strand. In sequence, the fluid-phase methodology was used to create functionalized alginate-gelatin-arginylglycylaspartic acid peptide (RGD) hydrogels via carbodiimides chemistry. The alginate-gelatin-RGD hydrogels showed an increase of 2.97-fold in cell growth and more spread substrate colonization in comparison to alginate-gelatin hydrogel. Moreover, the fluid-phase methodology was used to add RGD molecules to pre-determined parts of the alginate-gelatin substrate during the printing process promoting U2OS cell compartmentalization. In addition, different substrate stiffnesses were also created via fluid-phase by crosslinking the hydrogel with different concentrations of CaCl2during the printing process. As a result, the U2OS cells were also compartmentalized on the stiffer parts of the printings. Finally, our results showed that by combining stiffer hydrogel with RGD increasing concentrations we can create a synergetic effect and boost cell metabolism by up to 3.17-fold. This work presents an idea of a new printing process for tailoring multiple parameters in hydrogel substrates by using fluid-phase to generate more faithful replication of thein vivoenvironment.

Evaluation of the antitumor effect of neoantigen peptide vaccines derived from the translatome of lung cancer

Emerging evidence suggests that tumor-specific neoantigens are ideal targets for cancer immunotherapy. However, how to predict tumor neoantigens based on translatome data remains obscure. Through the extraction of ribosome-nascent chain complexes (RNCs) from LLC cells, followed by RNC-mRNA extraction, RNC-mRNA sequencing, and comprehensive bioinformatic analysis, we successfully identified proteins undergoing translatome and exhibiting mutations in the cells. Subsequently, novel antigens identification was analyzed by the interaction between their high affinity and the Major Histocompatibility Complex (MHC). Neoantigens immunogenicity was analyzed by enzyme-linked immunospot assay (ELISpot). Finally, in vivo experiments in mice were conducted to evaluate the antitumor effects of translatome-derived neoantigen peptides on lung cancer. The results showed that ten neoantigen peptides were identified and synthesized by translatome data from LLC cells; 8 out of the 10 neoantigens had strong immunogenicity. The neoantigen peptide vaccine group exhibited significant tumor growth inhibition effect. In conclusion, neoantigen peptide vaccine derived from the translatome of lung cancer exhibited significant tumor growth inhibition effect.

Proteomics: Unraveling the Cross Talk Between Innate Immunity and Disease Pathophysiology, Diagnostics, and Treatment Options

Inflammation is crucial in diseases, and proteins play a key role in the interplay between innate immunity and pathology. This review explores how proteomics helps understanding this relationship, focusing on diagnosis and treatment. We explore the dynamic innate response and the significance of proteomic techniques in deciphering the complex network of proteins involved in prevalent diseases, including infections, cancer, autoimmune and neurodegenerative disorders. Proteomics identifies key proteins in host-pathogen interactions, shedding light on infection mechanisms and inflammation. These discoveries hold promise for diagnostic tools, therapies, and vaccines. In cancer research, proteomics reveals innate signatures associated with tumor development, immune evasion, and therapeutic response. Additionally, proteomic analysis has unveiled autoantigens and dysregulation of the innate immune system in autoimmunity, offering opportunities for early diagnosis, disease monitoring, and new therapeutic targets. Moreover, proteomic analysis has identified altered protein expression patterns in neurodegenerative diseases like Alzheimer's and Parkinson's, providing insights into potential therapeutic strategies. Proteomics of the innate immune system provides a comprehensive understanding of disease mechanisms, identifies biomarkers, and enables effective interventions in various diseases. Despite still in its early stages, this approach holds great promise to revolutionize innate immunity research and significantly improve patient outcomes across a wide range of diseases.

Feasibility of mass cytometry proteomic characterisation of circulating tumour cells in head and neck squamous cell carcinoma for deep phenotyping

BACKGROUND

Circulating tumour cells (CTCs) are a potential cancer biomarker, but current methods of CTC analysis at single-cell resolution are limited. Here, we describe high-dimensional single-cell mass cytometry proteomic analysis of CTCs in HNSCC.

METHODS

Parsortix microfluidic-enriched CTCs from 14 treatment-naïve HNSCC patients were analysed by mass cytometry analysis using 41 antibodies. Immune cell lineage, epithelial-mesenchymal transition (EMT), stemness, proliferation and immune checkpoint expression was assessed alongside phosphorylation status of multiple signalling proteins. Patient-matched tumour gene expression and CTC EMT profiles were compared. Standard bulk CTC RNAseq was performed as a baseline comparator to assess mass cytometry data.

RESULTS

CTCs were detected in 13/14 patients with CTC counts of 2-24 CTCs/ml blood. Unsupervised clustering separated CTCs into epithelial, early EMT and advanced EMT groups that differed in signalling pathway activation state. Patient-specific CTC cluster patterns separated into immune checkpoint low and high groups. Patient tumour and CTC EMT profiles differed. Mass cytometry outperformed bulk RNAseq to detect CTCs and characterise cell phenotype.

DISCUSSION

We demonstrate mass cytometry allows high-plex proteomic characterisation of CTCs at single-cell resolution and identify common CTC sub-groups with potential for novel biomarker development and immune checkpoint inhibitor treatment stratification.

The mutational spectrum of ATRX aberrations in neuroblastoma and the associated patient and tumor characteristics

Neuroblastoma is the most common extracranial solid tumor in children. The chromatin remodeler ATRX is frequently mutated in high‐risk patients with a poor prognosis. Although many studies have reported ATRX aberrations and the associated clinical characteristics in neuroblastoma, a comprehensive overview is currently lacking. In this study, we extensively characterize the mutational spectrum of ATRX aberrations in neuroblastoma tumors reported in previous studies and present an overview of patient and tumor characteristics. We collected the data of a total of 127 neuroblastoma patients and three cell lines with ATRX aberrations originating from 20 papers. We subdivide the ATRX aberrations into nonsense, missense, and multiexon deletions (MEDs) and show that 68% of them are MEDs. Of these MEDs, 75% are predicted to be in‐frame. Furthermore, we identify a missense mutational hotspot region in the helicase domain. We also confirm that all three ATRX mutation types are more often identified in patients diagnosed at an older age, but still approximately 40% of the patients are aged 5 years or younger at diagnosis. Surprisingly, we found that 11q deletions are enriched in neuroblastomas with ATRX deletions compared to a reference cohort, but not in neuroblastomas with ATRX point mutations. Taken together, our data emphasizes a distinct ATRX mutation spectrum in neuroblastoma, which should be considered when studying molecular phenotypes and therapeutic strategies.

Experimental investigation of a combinational iron chelating protoporphyrin IX prodrug for fluorescence detection and photodynamic therapy

Photodynamic therapy (PDT) is an oxygen-dependent, light-activated, and locally destructive drug treatment of cancer. Protoporphyrin IX (PpIX)-induced PDT exploits cancer cells' own innate heme biosynthesis to hyper-accumulate the naturally fluorescent and photoactive precursor to heme, PpIX. This occurs as a result of administering heme precursors (e.g., aminolevulinic acid; ALA) because the final step of the pathway (the insertion of ferrous iron into PpIX by ferrochelatase to form heme) is relatively slow. Separate administration of an iron chelating agent has previously been demonstrated to significantly improve dermatological PpIX-PDT by further limiting heme production. A newly synthesized combinational iron chelating PpIX prodrug (AP2-18) has been assessed experimentally in cultured primary human cells of bladder and dermatological origin, as an alternative photosensitizing agent to ALA or its methyl or hexyl esters (MAL and HAL respectively) for photodetection/PDT. Findings indicated that the technique of iron chelation (either through the separate administration of the established hydroxypyridinone iron chelator CP94 or the just as effective combined AP2-18) did not enhance either PpIX fluorescence or PDT-induced (neutral red assessed) cell death in human primary normal and malignant bladder cells. However, 500 µM AP2-18 significantly increased PpIX accumulation and produced a trend of increased cell death within epithelial squamous carcinoma cells. PpIX accumulation destabilized the actin cytoskeleton in bladder cancer cells prior to PDT and resulted in caspase-3 cleavage/early apoptosis afterwards. AP2-18 iron chelation should continue to be investigated for the enhancement of dermatological PpIX-PDT applications but not bladder photodetection/PDT.

Folate Deficiency Triggers the Abnormal Segregation of a Region With Large Cluster of CG-Rich Trinucleotide Repeats on Human Chromosome 2

Folate deficiency is associated with a broad range of human disorders, including anemia, fetal neural tube defects, age-associated dementia and several types of cancer. It is well established that a subgroup of rare fragile sites (RFSs) containing expanded CGG trinucleotide repeat (TNR) sequences display instability when cells are deprived of folate. However, given that folate sensitive RFSs exist in a very small percentage of the population, they are unlikely to be the cause of the widespread health problems associated with folate deficiency. We hypothesized that folate deficiency could specifically affect DNA replication at regions containing CG-rich repeat sequences. For this, we identified a region on human chromosome 2 (Chr2) comprising more than 300 CG-rich TNRs (termed “FOLD1”) by examining the human genome database. Via the analysis of chromosome shape and segregation in mitosis, we demonstrate that, when human cells are cultured under folate stress conditions, Chr2 is prone to display a “kink” or “bending” at FOLD1 in metaphase and nondisjunction in anaphase. Furthermore, long-term folate deprivation causes Chr2 aneuploidy. Our results provide new evidence on the abnormalities folate deficiency could cause in human cells. This could facilitate future studies on the deleterious health conditions associated with folate deficiency.

Computational Analysis Concerning the Impact of DNA Accessibility on CRISPR-Cas9 Cleavage Efficiency

Defining the variables that impact the specificity of CRISPR/Cas9 has been a major research focus. Whereas sequence complementarity between guide RNA and target DNA substantially dictates cleavage efficiency, DNA accessibility of the targeted loci has also been hypothesized to be an important factor. In this study, functional data from two genome-wide assays, genome-wide, unbiased identification of DSBs enabled by sequencing (GUIDE-seq) and circularization for in vitro reporting of cleavage effects by sequencing (CIRCLE-seq), have been computationally analyzed in conjunction with DNA accessibility determined via DNase I-hypersensitive sequencing from the Encyclopedia of DNA Elements (ENCODE) Database and transcriptome from the Sequence Read Archive to determine whether cellular factors influence CRISPR-induced cleavage efficiency. CIRCLE-seq and GUIDE-seq datasets were selected to represent the absence and presence of cellular factors, respectively. Data analysis revealed that correlations between sequence similarity and CRISPR-induced cleavage frequency were altered by the presence of cellular factors that modulated the level of DNA accessibility. The above-mentioned correlation was abolished when cleavage sites were located in less accessible regions. Furthermore, CRISPR-mediated edits were permissive even at regions that were insufficient for most endogenous genes to be expressed. These results provide a strong basis to dissect the contribution of local chromatin modulation markers on CRISPR-induced cleavage efficiency.

Multiscale, multimodal analysis of tumor heterogeneity in IDH1 mutant vs wild-type diffuse gliomas

Glioma is recognized to be a highly heterogeneous CNS malignancy, whose diverse cellular composition and cellular interactions have not been well characterized. To gain new clinical- and biological-insights into the genetically-bifurcated IDH1 mutant (mt) vs wildtype (wt) forms of glioma, we integrated data from protein, genomic and MR imaging from 20 treatment-naïve glioma cases and 16 recurrent GBM cases. Multiplexed immunofluorescence (MxIF) was used to generate single cell data for 43 protein markers representing all cancer hallmarks, Genomic sequencing (exome and RNA (normal and tumor) and magnetic resonance imaging (MRI) quantitative features (protocols were T1-post, FLAIR and ADC) from whole tumor, peritumoral edema and enhancing core vs equivalent normal region were also collected from patients. Based on MxIF analysis, 85,767 cells (glioma cases) and 56,304 cells (GBM cases) were used to generate cell-level data for 24 biomarkers. K-means clustering was used to generate 7 distinct groups of cells with divergent biomarker profiles and deconvolution was used to assign RNA data into three classes. Spatial and molecular heterogeneity metrics were generated for the cell data. All features were compared between IDH mt and IDHwt patients and were finally combined to provide a holistic/integrated comparison. Protein expression by hallmark was generally lower in the IDHmt vs wt patients. Molecular and spatial heterogeneity scores for angiogenesis and cell invasion also differed between IDHmt and wt gliomas irrespective of prior treatment and tumor grade; these differences also persisted in the MR imaging features of peritumoral edema and contrast enhancement volumes. A coherent picture of enhanced angiogenesis in IDHwt tumors was derived from multiple platforms (genomic, proteomic and imaging) and scales from individual proteins to cell clusters and heterogeneity, as well as bulk tumor RNA and imaging features. Longer overall survival for IDH1mt glioma patients may reflect mutation-driven alterations in cellular, molecular, and spatial heterogeneity which manifest in discernable radiological manifestations.

Heterogeneous glutamatergic receptor mRNA expression across phrenic motor neurons in rats

The diaphragm muscle comprises various types of motor units that are recruited in an orderly fashion governed by the intrinsic electrophysiological properties (membrane capacitance as a function of somal surface area) of phrenic motor neurons (PhMNs). Glutamate is the main excitatory neurotransmitter at PhMNs and acts primarily via fast acting AMPA and N-methyl-D-aspartic acid (NMDA) receptors. Differences in receptor expression may also contribute to motor unit recruitment order. We used single cell, multiplex fluorescence in situ hybridization to determine glutamatergic receptor mRNA expression across PhMNs based on their somal surface area. In adult male and female rats (n = 9) PhMNs were retrogradely labeled for analyses (n = 453 neurons). Differences in the total number and density of mRNA transcripts were evident across PhMNs grouped into tertiles according to somal surface area. A ~ 25% higher density of AMPA (Gria2) and NMDA (Grin1) mRNA expression was evident in PhMNs in the lower tertile compared to the upper tertile. These smaller PhMNs likely comprise type S motor units that are recruited first to accomplish lower force, ventilatory behaviors. In contrast, larger PhMNs with lower volume densities of AMPA and NMDA mRNA expression presumably comprise type FInt and FF motor units that are recruited during higher force, expulsive behaviors. Furthermore, there was a significantly higher cytosolic NMDA mRNA expression in small PhMNs suggesting a more important role for NMDA-mediated glutamatergic neurotransmission at smaller PhMNs. These results are consistent with the observed order of motor unit recruitment and suggest a role for glutamatergic receptors in support of this orderly recruitment. Cover Image for this issue: doi: 10.1111/jnc.14747.

Epigallocatechin-3-gallate induces telomere shortening and clastogenic damage in glioblastoma cells

Epigallocatechingallate (EGCG) is the major polyphenol in green tea, to which many anticancer features, such as antioxidative, antigenotoxic, and antiangiogenetic properties, are attributed. Moreover, it is also well known as a telomerase inhibitor. In this work, we have chronically treated U251 glioblastoma cells with low, physiologically realistic concentrations, of EGCG, in order to investigate its effects both on telomeres and on genome integrity. Inhibition of telomerase activity caused telomere shortening, ultimately leading to senescence and telomere dysfunction at 98 days. Remarkably, we have observed DNA damage through an increase of phosphorylation of γ-H2AX histone and micronuclei also with doses and at timepoints when telomere shortening was not present. Therefore, we concluded that this DNA damage was not correlated with telomere shortening and that EGCG treatment induced not only an increase of telomere-shortening-induced senescence but also telomere-independent genotoxicity. This study questions the common knowledge about EGCG properties, but confirms the few works that indicated the clastogenic properties of this molecule, probably due to DNA reductive damage and topoisomerase II poisoning. Environ. Mol. Mutagen., 60:683-692, 2019. © 2019 Wiley Periodicals, Inc.

MUDENG Expression Profiling in Cohorts and Brain Tumor Biospecimens to Evaluate Its Role in Cancer

Mu-2-related death-inducing gene (MUDENG, MuD) has been reported to be involved in the tumor necrosis factor-related apoptosis-inducing ligand (TRAIL)-associated apoptotic pathway of glioblastoma multiforme (GBM) cells; however, its expression level, interactors, and role in tumors are yet to be discovered. To investigate whether MuD expression correlates with cancer progression, we analyzed The Cancer Genome Atlas (TCGA) database using UALCAN and Gene Expression Profiling Interactive Analysis (GEPIA). Differential expression of MuD was detected in 6 and 10 cancer types, respectively. Validation performed using data from the Gene Expression Omnibus database showed that MuD expression is downregulated in KIRC tumor and correlate with higher chance of survival. Upregulation of MuD expression in GBM tumors was detected through GEPIA and high MuD expression correlated with higher survival in proneural GBM, whereas the opposite was observed in classical GBM subtype. GBM biospecimens analysis shows that MuD protein level was upregulated in three of six specimens, whereas mRNA level remained relatively unaltered. Therefore, MuD may exert differential effects according to subtypes, and/or be subjected to post-translational regulation in GBM. Correlation analysis between GBM cohort database and experiments using GBM cell lines revealed its positive effect on regulation of protein phosphatase 2 regulatory subunit B’Epsilon (PPP2R5E) and son of sevenless homolog 2 (SOS2). STRING database analysis indicated that the components of adaptor protein complexes putatively interacted with MuD but showed no correlation in terms of survival of patients with different GBM subtypes. In summary, we analyzed the expression of MuD in publicly available cancer patient data sets, GBM cell lines, and biospecimens to demonstrate its potential role as a biomarker for cancer prognosis and identified its candidate interacting molecules.

EviNet: a web platform for network enrichment analysis with flexible definition of gene sets

The new web resource EviNet provides an easily run interface to network enrichment analysis for exploration of novel, experimentally defined gene sets. The major advantages of this analysis are (i) applicability to any genes found in the global network rather than only to those with pathway/ontology term annotations, (ii) ability to connect genes via different molecular mechanisms rather than within one high-throughput platform, and (iii) statistical power sufficient to detect enrichment of very small sets, down to individual genes. The users’ gene sets are either defined prior to upload or derived interactively from an uploaded file by differential expression criteria. The pathways and networks used in the analysis can be chosen from the collection menu. The calculation is typically done within seconds or minutes and the stable URL is provided immediately. The results are presented in both visual (network graphs) and tabular formats using jQuery libraries. Uploaded data and analysis results are kept in separated project directories not accessible by other users. EviNet is available at https://www.evinet.org/.

A family of PIKFYVE inhibitors with therapeutic potential against autophagy-dependent cancer cells disrupt multiple events in lysosome homeostasis

High-throughput screening identified 5 chemical analogs (termed the WX8-family) that disrupted 3 events in lysosome homeostasis: (1) lysosome fission via tubulation without preventing homotypic lysosome fusion; (2) trafficking of molecules into lysosomes without altering lysosomal acidity, and (3) heterotypic fusion between lysosomes and autophagosomes. Remarkably, these compounds did not prevent homotypic fusion between lysosomes, despite the fact that homotypic fusion required some of the same machinery essential for heterotypic fusion. These effects varied 400-fold among WX8-family members, were time and concentration dependent, reversible, and resulted primarily from their ability to bind specifically to the PIKFYVE phosphoinositide kinase. The ability of the WX8-family to prevent lysosomes from participating in macroautophagy/autophagy suggested they have therapeutic potential in treating autophagy-dependent diseases. In fact, the most potent family member (WX8) was 100-times more lethal to 'autophagy-addicted' melanoma A375 cells than the lysosomal inhibitors hydroxychloroquine and chloroquine. In contrast, cells that were insensitive to hydroxychloroquine and chloroquine were also insensitive to WX8. Therefore, the WX8-family of PIKFYVE inhibitors provides a basis for developing drugs that could selectively kill autophagy-dependent cancer cells, as well as increasing the effectiveness of established anti-cancer therapies through combinatorial treatments. Abbreviations: ACTB: actin beta; Baf: bafilomycin A₁; BECN1: beclin 1; BODIPY: boron-dipyrromethene; BORC: BLOC-1 related complex; BRAF: B-Raf proto-oncogene, serine/threonine kinase; BSA: bovine serum albumin; CTSD: cathepsin D; CQ: chloroquine; DNA: deoxyribonucleic acid; EC₅₀: half maximal effective concentration; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; GFP: green fluorescent protein; HCQ: hydroxychloroquine; HOPS complex: homotypic fusion and protein sorting complex; Kd: equilibrium binding constant; IC50: half maximal inhibitory concentration; KO: knockout; LAMP1: lysosomal associated membrane protein 1; MAP1LC3A: microtubule associated protein 1 light chain 3 alpha; MES: 2-(N-morpholino)ethanesulphonic acid; MTOR: mechanistic target of rapamycin kinase; μM: micromolar; NDF: 3-methylbenzaldehyde (2,6-dimorpholin-4-ylpyrimidin-4-yl)hydrazine;NEM: N-ethylmaleimide; NSF: N-ethylmaleimide sensitive factor; PBS: phosphate-buffered saline; PIKFYVE: phosphoinositide kinase, FYVE-type zinc finger containing; PIP4K2C: phosphatidylinositol-5-phosphate 4-kinase type 2 gamma; PtdIns3P: phosphatidylinositol 3-phosphate; PtdIns(3,5)P₂: phosphatidylinositol 3,5-biphosphate; RFP: red fluorescent protein; RPS6: ribosomal protein S6; RPS6KB1: ribosomal protein S6 kinase B1; SQSTM1: sequestosome 1; TWEEN 20: polysorbate 20; V-ATPase: vacuolar-type H⁺-translocating ATPase; VPS39: VPS39 subunit of HOPS complex; VPS41: VPS41 subunit of HOPS complex; WWL: benzaldehyde [2,6-di(4-morpholinyl)-4-pyrimidinyl]hydrazone; WX8: 1H-indole-3-carbaldehyde [4-anilino-6-(4-morpholinyl)-1,3,5-triazin-2-yl]hydrazine; XBA: N-(3-chloro-4-fluorophenyl)-4,6-dimorpholino-1,3,5-triazin-2-amine hydrochloride; XB6: N-(4-ethylphenyl)-4,6-dimorpholino-1,3,5-triazin-2-amine hydrochloride.

Novel protein signatures suggest progression to muscular invasiveness in bladder cancer

Patients with bladder cancer need frequent controls over long follow-up time due to high recurrence rate and risk of conversion to muscle invasive cancer with poor prognosis. We identified cancer-related molecular signatures in apparently healthy bladder in patients with subsequent muscular invasiveness during follow-up. Global proteomics of the normal tissue biopsies revealed specific proteome fingerprints in these patients prior to subsequent muscular invasiveness. In these presumed normal samples, we detected modulations of proteins previously associated with different cancer types. This study indicates that analyzing apparently healthy tissue of a cancer-invaded organ may suggest disease progression.

CRISPR-Sirius: RNA scaffolds for signal amplification in genome imaging

Clustered regularly interspaced short palindromic repeats (CRISPR) guide RNA scaffolds have been adapted to carry multiple binding sites for fluorescent proteins to enhance brightness for live cell imaging of genomic loci. However, many of these modifications result in guide RNA instability and thus produce lower genome-labeling efficiency than anticipated. Here we introduce CRISPR-Sirius, based on octet arrays of aptamers conferring both enhanced guide RNA stability and brightness, and provide initial biological applications of this platform.

Discovery of coding regions in the human genome by integrated proteogenomics analysis workflow

Proteogenomics enable the discovery of novel peptides (from unannotated genomic protein-coding loci) and single amino acid variant peptides (derived from single-nucleotide polymorphisms and mutations). Increasing the reliability of these identifications is crucial to ensure their usefulness for genome annotation and potential application as neoantigens in cancer immunotherapy. We here present integrated proteogenomics analysis workflow (IPAW), which combines peptide discovery, curation, and validation. IPAW includes the SpectrumAI tool for automated inspection of MS/MS spectra, eliminating false identifications of single-residue substitution peptides. We employ IPAW to analyze two proteomics data sets acquired from A431 cells and five normal human tissues using extended (pH range, 3–10) high-resolution isoelectric focusing (HiRIEF) pre-fractionation and TMT-based peptide quantitation. The IPAW results provide evidence for the translation of pseudogenes, lncRNAs, short ORFs, alternative ORFs, N-terminal extensions, and intronic sequences. Moreover, our quantitative analysis indicates that protein production from certain pseudogenes and lncRNAs is tissue specific.

Proteogenomics enables the discovery of protein coding regions and disease-relevant mutations but their verification remains challenging. Here, the authors combine peptide discovery, curation and validation in an integrated proteogenomics workflow, robustly identifying unknown coding regions and mutations.

A game of substrates: replication fork remodeling and its roles in genome stability and chemo-resistance

During the hours that human cells spend in the DNA synthesis (S) phase of the cell cycle, they may encounter adversities such as DNA damage or shortage of nucleotides. Under these stresses, replication forks in DNA may experience slowing, stalling, and breakage. Fork remodeling mechanisms, which stabilize slow or stalled replication forks and ensure their ability to continue or resume replication, protect cells from genomic instability and carcinogenesis. Fork remodeling includes DNA strand exchanges that result in annealing of newly synthesized strands (fork reversal), controlled DNA resection, and cleavage of DNA strands. Defects in major tumor suppressor genes BRCA1 and BRCA2, and a subset of the Fanconi Anemia genes have been shown to result in deregulation in fork remodeling, and most prominently, loss of kilobases of nascent DNA from stalled replication forks. This phenomenon has recently gained spotlight as a potential marker and mediator of chemo-sensitivity in cancer cells and, conversely, its suppression - as a hallmark of acquired chemo-resistance. Moreover, nascent strand degradation at forks is now known to also trigger innate immune response to self-DNA. An increasingly sophisticated molecular description of these events now points at a combination of unbalanced fork reversal and end-resection as a root cause, yet also reveals the multi-layered complexity and heterogeneity of the underlying processes in normal and cancer cells.

Orchestration of DNA Damage Checkpoint Dynamics across the Human Cell Cycle

Although molecular mechanisms that prompt cell-cycle arrest in response to DNA damage have been elucidated, the systems-level properties of DNA damage checkpoints are not understood. Here, using time-lapse microscopy and simulations that model the cell cycle as a series of Poisson processes, we characterize DNA damage checkpoints in individual, asynchronously proliferating cells. We demonstrate that, within early G1 and G2, checkpoints are stringent: DNA damage triggers an abrupt, all-or-none cell-cycle arrest. The duration of this arrest correlates with the severity of DNA damage. After the cell passes commitment points within G1 and G2, checkpoint stringency is relaxed. By contrast, all of S phase is comparatively insensitive to DNA damage. This checkpoint is graded: instead of halting the cell cycle, increasing DNA damage leads to slower S phase progression. In sum, we show that a cell's response to DNA damage depends on its exact cell-cycle position and that checkpoints are phase-dependent, stringent or relaxed, and graded or all-or-none.

Inherited DNA lesions determine G1 duration in the next cell cycle

Replication stress is a major source of DNA damage and an important driver of cancer development. Replication intermediates that occur upon mild forms of replication stress frequently escape cell cycle checkpoints and can be transmitted through mitosis into the next cell cycle. The consequences of such inherited DNA lesions for cell fate and survival are poorly understood. By using time-lapse microscopy and quantitative image-based cytometry to simultaneously monitor inherited DNA lesions marked by the genome caretaker protein 53BP1 and cell cycle progression, we show that inheritance of 53BP1-marked lesions from the previous S-phase is associated with a prolonged G1 duration in the next cell cycle. These results suggest that cell-to-cell variation in S-phase commitment is determined, at least partially, by the amount of replication-born inherited DNA damage in individual cells. We further show that loss of the tumor suppressor protein p53 overrides replication stress-induced G1 prolongation and allows S-phase entry with excessive amounts of inherited DNA lesions. Thus, replication stress and p53 loss may synergize during cancer development by promoting cell cycle re-entry with unrepaired mutagenic DNA lesions originating from the previous cell cycle.

Bioinformatics Methods and Tools to Advance Clinical Care. Findings from the Yearbook 2015 Section on Bioinformatics and Translational Informatics

OBJECTIVES

To summarize excellent current research in the field of Bioinformatics and Translational Informatics with application in the health domain and clinical care.

METHOD

We provide a synopsis of the articles selected for the IMIA Yearbook 2015, from which we attempt to derive a synthetic overview of current and future activities in the field. As last year, a first step of selection was performed by querying MEDLINE with a list of MeSH descriptors completed by a list of terms adapted to the section. Each section editor has evaluated separately the set of 1,594 articles and the evaluation results were merged for retaining 15 articles for peer-review.

RESULTS

The selection and evaluation process of this Yearbook's section on Bioinformatics and Translational Informatics yielded four excellent articles regarding data management and genome medicine that are mainly tool-based papers. In the first article, the authors present PPISURV a tool for uncovering the role of specific genes in cancer survival outcome. The second article describes the classifier PredictSNP which combines six performing tools for predicting disease-related mutations. In the third article, by presenting a high-coverage map of the human proteome using high resolution mass spectrometry, the authors highlight the need for using mass spectrometry to complement genome annotation. The fourth article is also related to patient survival and decision support. The authors present datamining methods of large-scale datasets of past transplants. The objective is to identify chances of survival.

CONCLUSIONS

The current research activities still attest the continuous convergence of Bioinformatics and Medical Informatics, with a focus this year on dedicated tools and methods to advance clinical care. Indeed, there is a need for powerful tools for managing and interpreting complex, large-scale genomic and biological datasets, but also a need for user-friendly tools developed for the clinicians in their daily practice. All the recent research and development efforts contribute to the challenge of impacting clinically the obtained results towards a personalized medicine.

Impaired STING Pathway in Human Osteosarcoma U2OS Cells Contributes to the Growth of ICP0-Null Mutant Herpes Simplex Virus

Human herpes simplex virus 1 (HSV-1) is a widespread pathogen, with 80% of the population being latently infected. To successfully evade the host, the virus has evolved strategies to counteract antiviral responses, including the gene-silencing and innate immunity machineries. The immediately early protein of the virus, infected cell protein 0 (ICP0), plays a central role in these processes. ICP0 blocks innate immunity, and one mechanism is by degrading hostile factors with its intrinsic E3 ligase activity. ICP0 also functions as a promiscuous transactivator, and it blocks repressor complexes to enable viral gene transcription. For these reasons, the growth of a ΔICP0 virus is impaired in most cells, except cells of the human osteosarcoma cell line U2OS, and it is only partially impaired in cells of the human osteosarcoma cell line Saos-2. We found that the two human osteosarcoma cell lines that supported the growth of the ΔICP0 virus failed to activate innate immune responses upon treatment with 2'3'-cyclic GAMP (2'3'-cGAMP), the natural agonist of STING (i.e., stimulator of interferon genes) or after infection with the ΔICP0 mutant virus. Innate immune responses were restored in these cells by transient expression of the STING protein but not after overexpression of interferon-inducible protein 16 (IFI16). Restoration of STING expression resulted in suppression of ΔICP0 virus gene expression and a decrease in viral yields. Overexpression of IFI16 also suppressed ΔICP0 virus gene expression, albeit to a lesser extent than STING. These data suggest that the susceptibility of U2OS and Saos-2 cells to the ΔICP0 HSV-1 is in part due to an impaired STING pathway.IMPORTANCE The DNA sensor STING plays pivotal role in controlling HSV-1 infection both in cell culture and in mice. The HSV-1 genome encodes numerous proteins that are dedicated to combat host antiviral responses. The immediate early protein of the virus ICP0 plays major role in this process as it targets hostile host proteins for degradation with its E3 ligase activity, and it disrupts repressor complexes via protein-protein interaction to enable viral gene transcription. Therefore, the ΔICP0 HSV-1 virus is defective for growth in most cells, except the human osteosarcoma cell lines U2OS and Saos-2. We found that both cell lines that support ΔICP0 virus infection have defects in the STING DNA-sensing pathway, which partially accounts for the rescue of the ΔICP0 virus growth. Restoration of STING expression in these cells rescued innate immunity and suppressed ΔICP0 virus infection. This study underscores the importance of STING in the control of HSV-1.

Profiling DNA damage response following mitotic perturbations

Genome integrity relies on precise coordination between DNA replication and chromosome segregation. Whereas replication stress attracted much attention, the consequences of mitotic perturbations for genome integrity are less understood. Here, we knockdown 47 validated mitotic regulators to show that a broad spectrum of mitotic errors correlates with increased DNA breakage in daughter cells. Unexpectedly, we find that only a subset of these correlations are functionally linked. We identify the genuine mitosis-born DNA damage events and sub-classify them according to penetrance of the observed phenotypes. To demonstrate the potential of this resource, we show that DNA breakage after cytokinesis failure is preceded by replication stress, which mounts during consecutive cell cycles and coincides with decreased proliferation. Together, our results provide a resource to gauge the magnitude and dynamics of DNA breakage associated with mitotic aberrations and suggest that replication stress might limit propagation of cells with abnormal karyotypes.

DNA damage arising from replication stress is well studied, but the effect of mitotic errors on genome integrity is less understood. Here the authors knock down 47 mitotic regulators and record how they impact on DNA breakage events, providing a resource for future studies on the relation between cell division and genome integrity.

Cell line donor genotype and its influence on experimental phenotype: Toll-like receptor SNPs and potential variability in innate immunity

Cell lines are used to model a disease and provide valuable information regarding phenotype, mechanism, and response to novel therapies. Derived from individuals of diverse genetic backgrounds, the cell's genetic complement predicts the phenotype, and although some lines have been sequenced, little emphasis has been placed on genotyping. Toll-like receptors (TLRs) are essential in initiating the inflammatory cascade in response to infection. TLR single nucleotide polymorphism (SNP) alleles may predict an altered innate immune response: a SNP can affect TLR-dependent pathways and may alter experimental results. Thus, genotype variation may have far-reaching implications when using cell lines to model phenotypes. We recommend that cell lines be genotyped and cataloged in a fashion similar to that used for bacteria, with cumulative information being archived in an accessible central database to facilitate the understanding of SNP cell phenotypes reported in the literature.

Maternal age effect on mouse oocytes: new biological insight from proteomic analysis

The long-standing view of ‘immortal germline vs mortal soma’ poses a fundamental question in biology concerning how oocytes age in molecular terms. A mainstream hypothesis is that maternal ageing of oocytes has its roots in gene transcription. Investigating the proteins resulting from mRNA translation would reveal how far the levels of functionally available proteins correlate with mRNAs and would offer novel insights into the changes oocytes undergo during maternal ageing. Gene ontology (GO) semantic analysis revealed a high similarity of the detected proteome (2324 proteins) to the transcriptome (22 334 mRNAs), although not all proteins had a cognate mRNA. Concerning their dynamics, fourfold changes of abundance were more frequent in the proteome (3%) than the transcriptome (0.05%), with no correlation. Whereas proteins associated with the nucleus (e.g. structural maintenance of chromosomes and spindle-assembly checkpoints) were largely represented among those that change in oocytes during maternal ageing; proteins associated with oxidative stress/damage (e.g. superoxide dismutase) were infrequent. These quantitative alterations are either impoverishing or enriching. Using GO analysis, these alterations do not relate in any simple way to the classic signature of ageing known from somatic tissues. Given the lack of correlation, we conclude that proteome analysis of mouse oocytes may not be surrogated with transcriptome analysis. Furthermore, we conclude that the classic features of ageing may not be transposed from somatic tissues to oocytes in a one-to-one fashion. Overall, there is more to the maternal ageing of oocytes than mere cellular deterioration exemplified by the notorious increase of meiotic aneuploidy.

Systems Analysis of the Response of Photosynthesis, Metabolism, and Growth to an Increase in Irradiance in the Photosynthetic Model Organism Chlamydomonas reinhardtii

We investigated the systems response of metabolism and growth after an increase in irradiance in the nonsaturating range in the algal model Chlamydomonas reinhardtii. In a three-step process, photosynthesis and the levels of metabolites increased immediately, growth increased after 10 to 15 min, and transcript and protein abundance responded by 40 and 120 to 240 min, respectively. In the first phase, starch and metabolites provided a transient buffer for carbon until growth increased. This uncouples photosynthesis from growth in a fluctuating light environment. In the first and second phases, rising metabolite levels and increased polysome loading drove an increase in fluxes. Most Calvin-Benson cycle (CBC) enzymes were substrate-limited in vivo, and strikingly, many were present at higher concentrations than their substrates, explaining how rising metabolite levels stimulate CBC flux. Rubisco, fructose-1,6-biosphosphatase, and seduheptulose-1,7-bisphosphatase were close to substrate saturation in vivo, and flux was increased by posttranslational activation. In the third phase, changes in abundance of particular proteins, including increases in plastidial ATP synthase and some CBC enzymes, relieved potential bottlenecks and readjusted protein allocation between different processes. Despite reasonable overall agreement between changes in transcript and protein abundance (R2 = 0.24), many proteins, including those in photosynthesis, changed independently of transcript abundance.

Ionizing radiation-inducible miR-494 promotes glioma cell invasion through EGFR stabilization by targeting p190B rhoGAP

MicroRNAs (miRNAs) play an important role in various stages of tumor progression. miR-494, which we had previously identified as a miRNA induced by ionizing radiation (IR) in the glioma cell line U-251, was observed to enhance invasion of U-251 cells by activating MMP-2. The miR-494-induced invasive potential was accompanied by, and dependent on, epidermal growth factor receptor (EGFR) upregulation and the activation of its downstream signaling constituents, Akt and ERK. The upregulation of EGFR by miR-494 involved the suppression of lysosomal protein turnover. Among the putative target proteins tested, p190B RhoGAP (p190B) was downregulated by miR-494, and its reduced expression was responsible for the increase in EGFR expression. A reporter assay using a luciferase construct containing p190B 3'-untranslated region (3'UTR) confirmed that p190B is a direct target of miR-494. Downregulation of p190B by small interfering RNA (siRNA) transfection closely mimicked the outcomes of miR-494 transfection, and showed increased EGFR expression, MMP-2 secretion, and invasion. Ectopic expression of p190B suppressed the miR-494-induced EGFR upregulation and invasion promotion, thereby suggesting that p190B depletion is critical for the invasion-promoting action of miR-494. Collectively, our results suggest a novel function for miR-494 and its potential application as a target to control invasiveness in cancer therapy.

HiRIEF LC-MS enables deep proteome coverage and unbiased proteogenomics

We present a liquid chromatography-mass spectrometry (LC-MS)-based method permitting unbiased (gene prediction-independent) genome-wide discovery of protein-coding loci in higher eukaryotes. Using high-resolution isoelectric focusing (HiRIEF) at the peptide level in the 3.7-5.0 pH range and accurate peptide isoelectric point (pI) prediction, we probed the six-reading-frame translation of the human and mouse genomes and identified 98 and 52 previously undiscovered protein-coding loci, respectively. The method also enabled deep proteome coverage, identifying 13,078 human and 10,637 mouse proteins.

Resolving chromosome-centric human proteome with translating mRNA analysis: a strategic demonstration

Chromosome-centric human proteome project (C-HPP) aims at differentiating chromosome-based and tissue-specific protein compositions in terms of protein expression, quantification, and modification. We previously found that the analysis of translating mRNA (mRNA attached to ribosome-nascent chain complex, RNC-mRNA) can explain over 94% of mRNA-protein abundance. Therefore, we propose here to use full-length RNC-mRNA information to illustrate protein expression both qualitatively and quantitatively. We performed RNA-seq on RNC-mRNA (RNC-seq) and detected 12,758 and 14,113 translating genes in human normal bronchial epithelial (HBE) cells and human colorectal adenocarcinoma Caco-2 cells, respectively. We found that most of these genes were mapped with >80% of coding sequence coverage. In Caco-2 cells, we provided translating evidence on 4180 significant single-nucleotide variations. While using RNC-mRNA data as a standard for proteomic data integration, both translating and protein evidence of 7876 genes can be acquired from four interlaboratory data sets with different MS platforms. In addition, we detected 1397 noncoding mRNAs that were attached to ribosomes, suggesting a potential source of new protein explorations. By comparing the two cell lines, a total of 677 differentially translated genes were found to be nonevenly distributed across chromosomes. In addition, 2105 genes in Caco-2 and 750 genes in HBE cells are expressed in a cell-specific manner. These genes are significantly and specifically clustered on multiple chromosomes, such as chromosome 19. We conclude that HPP/C-HPP investigations can be considerably improved by integrating RNC-mRNA analysis with MS, bioinformatics, and antibody-based verifications.

Proteome changes induced by c-myb silencing in human chronic myeloid leukemia cells suggest molecular mechanisms and putative biomarkers of hematopoietic malignancies

To shed light on the molecular mechanisms associated with aberrant accumulation of c-Myb in chronic myeloid leukemia, comparative proteomic analysis was performed on c-myb RNAi-specifically silenced K562 cells, sampled on a time-course basis. 2D-DIGE technology highlighted 37 differentially-represented proteins that were further characterized by nLC-ESI-LIT-MS/MS and validated by western blotting and qRT-PCR analysis. Most of the deregulated proteins were related to protein folding, energy/primary metabolism, transcription/translation regulation and oxidative stress response. Protein network analysis suggested that glycolysis, gluconeogenesis and protein ubiquitination biosynthesis pathways were highly represented, confirming also the pivotal role of c-Myc. A specific reduced representation was observed for glyceraldehyde-3-phosphate-dehydrogenase and α-enolase, suggesting a possible role of c-Myb in the activation of aerobic glycolysis. A reduced amount was also observed for stress responsive heat shock 70kDa protein and 78kDa glucose-regulated protein, previously identified as direct targets of c-Myb. Among over-represented proteins, worth mentioning is the chromatin modifier chromobox protein homolog 3 that contributes to silencing of E2F- and Myc-responsive genes in quiescent G0 cells. Data here presented, while providing novel insights onto the molecular mechanisms underlying c-Myb activity, indicate potential protein biomarkers for monitoring the progression of chronic myeloid leukemia.

BIOLOGICAL SIGNIFICANCE

Myeloid leukemia is a malignant disease of the hematopoietic system in which cells of myeloid lineages accumulate to an undifferentiated state. In particular, it was shown that an aberrant accumulation of the c-Myb transcriptional factor is associated with the suppression of normal differentiation processes promoting the development of the hematopoietic malignancies. Many efforts have been recently made to identify novel genes directly targeted by c-Myb at a transcriptome level. In this work, we originally describe a differential proteomic approach to facilitate the comprehension of the regulation of the protein networks exerted by c-Myb. Our study reveals a complex network of proteins regulated by c-Myb. The functional heterogeneity of these proteins emphasizes the pleiotropic role of c-Myb as a regulator of genes that are crucial for energy production and stress response in leukemia. In fact, variations in glyceraldehyde-3-phosphate-dehydrogenase and α-enolase suggest a possible role of c-Myb in the activation of aerobic glycolysis. Moreover, significant differences were found for heat shock 70kDa protein and 78kDa glucose-regulated protein known as direct c-Myb targets. This work highlights potential protein biomarkers to look into disease progression and to develop translational medicine approaches in myeloid leukemia.

Translating mRNAs strongly correlate to proteins in a multivariate manner and their translation ratios are phenotype specific

As a well-known phenomenon, total mRNAs poorly correlate to proteins in their abundances as reported. Recent findings calculated with bivariate models suggested even poorer such correlation, whereas focusing on the translating mRNAs (ribosome nascent-chain complex-bound mRNAs, RNC-mRNAs) subset. In this study, we analysed the relative abundances of mRNAs, RNC-mRNAs and proteins on genome-wide scale, comparing human lung cancer A549 and H1299 cells with normal human bronchial epithelial (HBE) cells, respectively. As discovered, a strong correlation between RNC-mRNAs and proteins in their relative abundances could be established through a multivariate linear model by integrating the mRNA length as a key factor. The R² reached 0.94 and 0.97 in A549 versus HBE and H1299 versus HBE comparisons, respectively. This correlation highlighted that the mRNA length significantly contributes to the translational modulation, especially to the translational initiation, favoured by its correlation with the mRNA translation ratio (TR) as observed. We found TR is highly phenotype specific, which was substantiated by both pathway analysis and biased TRs of the splice variants of BDP1 gene, which is a key transcription factor of transfer RNAs. These findings revealed, for the first time, the intrinsic and genome-wide translation modulations at translatomic level in human cells at steady-state, which are tightly correlated to the protein abundance and functionally relevant to cellular phenotypes.

Interpretation, stratification and evidence for sequence variants affecting mRNA splicing in complete human genome sequences

Information theory-based methods have been shown to be sensitive and specific for predicting and quantifying the effects of non-coding mutations in Mendelian diseases. We present the Shannon pipeline software for genome-scale mutation analysis and provide evidence that the software predicts variants affecting mRNA splicing. Individual information contents (in bits) of reference and variant splice sites are compared and significant differences are annotated and prioritized. The software has been implemented for CLC-Bio Genomics platform. Annotation indicates the context of novel mutations as well as common and rare SNPs with splicing effects. Potential natural and cryptic mRNA splicing variants are identified, and null mutations are distinguished from leaky mutations. Mutations and rare SNPs were predicted in genomes of three cancer cell lines (U2OS, U251 and A431), which were supported by expression analyses. After filtering, tractable numbers of potentially deleterious variants are predicted by the software, suitable for further laboratory investigation. In these cell lines, novel functional variants comprised 6–17 inactivating mutations, 1–5 leaky mutations and 6–13 cryptic splicing mutations. Predicted effects were validated by RNA-seq analysis of the three aforementioned cancer cell lines, and expression microarray analysis of SNPs in HapMap cell lines.

Comparison of total and cytoplasmic mRNA reveals global regulation by nuclear retention and miRNAs

BACKGROUND

The majority of published gene-expression studies have used RNA isolated from whole cells, overlooking the potential impact of including nuclear transcriptome in the analyses. In this study, mRNA fractions from the cytoplasm and from whole cells (total RNA) were prepared from three human cell lines and sequenced using massive parallel sequencing.

RESULTS

For all three cell lines, of about 15000 detected genes approximately 400 to 1400 genes were detected in different amounts in the cytoplasmic and total RNA fractions. Transcripts detected at higher levels in the total RNA fraction had longer coding sequences and higher number of miRNA target sites. Transcripts detected at higher levels in the cytoplasmic fraction were shorter or contained shorter untranslated regions. Nuclear retention of transcripts and mRNA degradation via miRNA pathway might contribute to this differential detection of genes. The consequence of the differential detection was further investigated by comparison to proteomics data. Interestingly, the expression profiles of cytoplasmic and total RNA correlated equally well with protein abundance levels indicating regulation at a higher level.

CONCLUSIONS

We conclude that expression levels derived from the total RNA fraction be regarded as an appropriate estimate of the amount of mRNAs present in a given cell population, independent of the coding sequence length or UTRs.