A Comprehensive Biomarker Analysis of Microsatellite Unstable/Mismatch Repair Deficient Colorectal Cancer Cohort Treated with Immunotherapy.

The search for immunotherapy biomarkers in Microsatellite Instability High/Deficient Mismatch Repair system (MSI-H/dMMR) metastatic colorectal cancer (mCRC) is an unmet need. Sixteen patients with mCRC and MSI-H/dMMR (determined by either immunohistochemistry or polymerase chain reaction) treated with PD-1/PD-L1 inhibitors at our institution were included. According to whether the progression-free survival with PD-1/PD-L1 inhibitors was longer than 6 months or shorter, patients were clustered into the IT-responder group (n: 9 patients) or IT-resistant group (n: 7 patients), respectively. In order to evaluate determinants of benefit with PD-1/PD-L1 inhibitors, we performed multimodal analysis including genomics (through NGS panel tumour-only with 431 genes) and the immune microenvironment (using CD3, CD8, FOXP3 and PD-L1 antibodies). The following mutations were more frequent in IT-resistant compared with IT-responder groups: B2M (4/7 versus 2/9), CTNNB1 (2/7 versus 0/9), and biallelic PTEN (3/7 versus 1/9). Biallelic ARID1A mutations were found exclusively in the IT-responder group (4/9 patients). Tumour mutational burden did not correlate with immunotherapy benefit, neither the rate of indels in homopolymeric regions. Of note, biallelic ARID1A mutated tumours had the highest immune infiltration and PD-L1 scores, contrary to tumours with CTNNB1 mutation. Immune microenvironment analysis showed higher densities of different T cell subpopulations and PD-L1 expression in IT-responders. Misdiagnosis of MSI-H/dMMR inferred by discordances between immunohistochemistry and polymerase chain reaction was only found in the IT-resistant population (3/7 patients). Biallelic ARID1A mutations and Wnt signalling activation through CTNNB1 mutation were associated with high and low T cell immune infiltrates, respectively, and deserve special attention as determinants of response to PD-1/PD-L1 inhibitors. The non-MSI-H phenotype in dMMR is associated with poor benefit to immunotherapy. Our results suggest that mechanisms of resistance to immunotherapy are multi-factorial.

Molecular profiling of long-term responders to immune checkpoint inhibitors in advanced non-small cell lung cancer.

Immunotherapy has transformed advanced non-small cell lung cancer (NSCLC) treatment strategies and has led to unprecedented long-lasting responses in some patients. However, the molecular determinants driving these long-term responses remain elusive. To address this issue, we performed an integrative analysis of genomic and transcriptomic features of long-term immune checkpoint inhibitors (ICIs)-associated responders. We assembled a cohort of 47 patients with NSCLC receiving ICIs that was enriched in long-term responders [>18 months of progression-free survival (PFS)]. We performed whole-exome sequencing from tumor samples, estimated the tumor mutational burden (TMB), and inferred the somatic copy number alterations (SCNAs). We also obtained gene transcription data for a subset of patients using Nanostring, which we used to assess the tumor immune infiltration status and PD-L1 expression. Our results indicate that there is an association between TMB and benefit to ICIs, which is driven by those patients with long-term response. Additionally, high SCNAs burden is associated with poor response and negatively correlates with the presence of several immune cell types (B cells, natural killers, regulatory T cells or effector CD8 T cells). Also, CD274 (PD-L1) expression is increased in patients with benefit, mainly in those with long-term response. In our cohort, combined assessment of TMB and SCNAs burden enabled identification of long-term responders (considering PFS and overall survival). Notably, the association between TMB, SCNAs burden, and PD-L1 expression with the outcomes of ICIs treatment was validated in two public datasets of ICI-treated patients with NSCLC. Thus, our data indicate that TMB is associated with long-term benefit following ICIs treatment in NSCLC and that TMB, SCNAs burden, and PD-L1 are complementary determinants of response to ICIs.

Clinical value of next generation sequencing of plasma cell-free DNA in gastrointestinal stromal tumors.

BACKGROUND: Gastrointestinal stromal tumor (GIST) initiation and evolution is commonly framed by KIT/PDGFRA oncogenic activation, and in later stages by the polyclonal expansion of resistant subpopulations harboring KIT secondary mutations after the onset of imatinib resistance. Thus, circulating tumor (ct)DNA determination is expected to be an informative non-invasive dynamic biomarker in GIST patients. METHODS: We performed amplicon-based next-generation sequencing (NGS) across 60 clinically relevant genes in 37 plasma samples from 18 GIST patients collected prospectively. ctDNA alterations were compared with NGS of matched tumor tissue samples (obtained either simultaneously or at the time of diagnosis) and cross-validated with droplet digital PCR (ddPCR). RESULTS: We were able to identify cfDNA mutations in five out of 18 patients had detectable in at least one timepoint. Overall, NGS sensitivity for detection of cell-free (cf)DNA mutations in plasma was 28.6%, showing high concordance with ddPCR confirmation. We found that GIST had relatively low ctDNA shedding, and mutations were at low allele frequencies. ctDNA was detected only in GIST patients with advanced disease after imatinib failure, predicting tumor dynamics in serial monitoring. KIT secondary mutations were the only mechanism of resistance found across 10 imatinib-resistant GIST patients progressing to sunitinib or regorafenib. CONCLUSIONS: ctDNA evaluation with amplicon-based NGS detects KIT primary and secondary mutations in metastatic GIST patients, particularly after imatinib progression. GIST exhibits low ctDNA shedding, but ctDNA monitoring, when positive, reflects tumor dynamics.

Identification of Expression Profiles Defining Distinct Prognostic Subsets of Radioactive-Iodine Refractory Differentiated Thyroid Cancer from the DECISION Trial.

Several biomarkers have been suggested to have prognostic value in differentiated thyroid carcinomas (DTC) with no validation in the refractory setting, including all tumor subtypes. We aim to correlate RNA expression profiles with survival based on patients included in the DECISION trial. We obtained 247 samples from the 417 patients included in the DECISION study and performed RNAseq analysis (77 million paired-end reads for each sample on HiSeq2000). After quality control, 125 samples were included in the secondary analysis and mapped against the human reference genome (GRCh38) with STAR (v2.5.1b) using ENCODE parameter. Survival analysis was calculated using the Kaplan-Meier method and log-rank test was used for statistical comparison. In this post hoc analysis, we identified three groups of tumors based on their gene expression profile: BRAF-like, RAS-like, and non-BRAF-non-RAS-like (NoBRaL). No significant correlation with sorafenib responders was observed. However, we identified a statistically significant correlation between the RNA-expression profiles and progression-free survival. The BRAF-like profile had a significantly better outcome compared with RAS-like and NoBRaL (11.8, 6.2, and 5.5 months, respectively) [HR: 0.31, 95% confidence interval (CI), 0.17-0.60; P < 0.001 and HR: 0.36 (95% CI, 0.21-0.63); P < 0.001] and HR: 0.36 (95% CI, 0.21-0.63; P < 0.001) and maintained significance as an independent prognostic factor for overall survival in the multivariate analysis for papillary thyroid cancers. To our knowledge, this is the first comprehensive RNA-seq analysis of all histologic subtypes of DTC. The RNA expression profiles identified may suggest a new prognostic parameter to be considered before recommendation of systemic therapies or the design of stratification factors for future clinical trials.

Epigenetic EGFR Gene Repression Confers Sensitivity to Therapeutic BRAFV600E Blockade in Colon Neuroendocrine Carcinomas.

PURPOSE: The limited knowledge of the molecular alterations that characterize poorly differentiated neuroendocrine carcinomas has limited the clinical development of targeted agents directed to driver mutations. Here we aim to identify new molecular targets in colon neuroendocrine carcinomas (co-NEC) and proof the efficacy of matching drugs. EXPERIMENTAL DESIGN: We performed a multi-omic analysis of co-NEC to identify genetic or epigenetic alterations that could be exploited as effective drug targets. We compared co-NEC samples with colorectal carcinomas (CRC) to identify neuroendocrine-specific traits. Patients with co-NEC and patient-derived xenografts were treated with a BRAFV600E-blocking drug to demonstrate sensitivity. RESULTS: co-NEC and CRC are similar in their mutational repertoire, although co-NECs are particularly enriched in BRAFV600E mutations. We report for the first time that V600EBRAF-mutant co-NECs may benefit from BRAF inhibition in monotherapy and how EGFR status is essential to predict innate sensitivity and acquired resistance by a differential methylation of its gene regulatory regions. CONCLUSIONS: The identification of V600E BRAF mutations in high-grade co-NECs has allowed the description of radiological responses to combination therapy of BRAF and MEK inhibitors in basket clinical trials. However, the molecular rationale for this treatment combination was based on the presence of the BRAF mutation and the efficacy observed in other cancer types such as melanoma. Future drug development in this setting should test BRAF inhibitors upfront and the addition of anti-EGFR antibodies instead of MEK inhibitors for an efficient blockade of acquired resistance.

Signatures of Evolutionary Adaptation in Quantitative Trait Loci Influencing Trace Element Homeostasis in Liver.

Essential trace elements possess vital functions at molecular, cellular, and physiological levels in health and disease, and they are tightly regulated in the human body. In order to assess variability and potential adaptive evolution of trace element homeostasis, we quantified 18 trace elements in 150 liver samples, together with the expression levels of 90 genes and abundances of 40 proteins involved in their homeostasis. Additionally, we genotyped 169 single nucleotide polymorphism (SNPs) in the same sample set. We detected significant associations for 8 protein quantitative trait loci (pQTL), 10 expression quantitative trait loci (eQTLs), and 15 micronutrient quantitative trait loci (nutriQTL). Six of these exceeded the false discovery rate cutoff and were related to essential trace elements: 1) one pQTL for GPX2 (rs10133290); 2) two previously described eQTLs for HFE (rs12346) and SELO (rs4838862) expression; and 3) three nutriQTLs: The pathogenic C282Y mutation at HFE affecting iron (rs1800562), and two SNPs within several clustered metallothionein genes determining selenium concentration (rs1811322 and rs904773). Within the complete set of significant QTLs (which involved 30 SNPs and 20 gene regions), we identified 12 SNPs with extreme patterns of population differentiation (FST values in the top 5% percentile in at least one HapMap population pair) and significant evidence for selective sweeps involving QTLs at GPX1, SELENBP1, GPX3, SLC30A9, and SLC39A8. Overall, this detailed study of various molecular phenotypes illustrates the role of regulatory variants in explaining differences in trace element homeostasis among populations and in the human adaptive response to environmental pressures related to micronutrients.

Rescuing discarded spectra: Full comprehensive analysis of a minimal proteome.

A common problem encountered when performing large-scale MS proteome analysis is the loss of information due to the high percentage of unassigned spectra. To determine the causes behind this loss we have analyzed the proteome of one of the smallest living bacteria that can be grown axenically, Mycoplasma pneumoniae (729 ORFs). The proteome of M. pneumoniae cells, grown in defined media, was analyzed by MS. An initial search with both Mascot and a species-specific NCBInr database with common contaminants (NCBImpn), resulted in around 79% of the acquired spectra not having an assignment. The percentage of non-assigned spectra was reduced to 27% after re-analysis of the data with the PEAKS software, thereby increasing the proteome coverage of M. pneumoniae from the initial 60% to over 76%. Nonetheless, 33,413 spectra with assigned amino acid sequences could not be mapped to any NCBInr database protein sequence. Approximately, 1% of these unassigned peptides corresponded to PTMs and 4% to M. pneumoniae protein variants (deamidation and translation inaccuracies). The most abundant peptide sequence variants (Phe-Tyr and Ala-Ser) could be explained by alterations in the editing capacity of the corresponding tRNA synthases. About another 1% of the peptides not associated to any protein had repetitions of the same aromatic/hydrophobic amino acid at the N-terminus, or had Arg/Lys at the C-terminus. Thus, in a model system, we have maximized the number of assigned spectra to 73% (51,453 out of the 70,040 initial acquired spectra). All MS data have been deposited in the ProteomeXchange with identifier PXD002779 (http://proteomecentral.proteomexchange.org/dataset/PXD002779).

The propagation of perturbations in rewired bacterial gene networks.

What happens to gene expression when you add new links to a gene regulatory network? To answer this question, we profile 85 network rewirings in E. coli. Here we report that concerted patterns of differential expression propagate from reconnected hub genes. The rewirings link promoter regions to different transcription factor and σ-factor genes, resulting in perturbations that span four orders of magnitude, changing up to ∼ 70% of the transcriptome. Importantly, factor connectivity and promoter activity both associate with perturbation size. Perturbations from related rewirings have more similar transcription profiles and a statistical analysis reveals ∼ 20 underlying states of the system, associating particular gene groups with rewiring constructs. We examine two large clusters (ribosomal and flagellar genes) in detail. These represent alternative global outcomes from different rewirings because of antagonism between these major cell states. This data set of systematically related perturbations enables reverse engineering and discovery of underlying network interactions.

miR-17-92 fine-tunes MYC expression and function to ensure optimal B cell lymphoma growth.

The synergism between c-MYC and miR-17-19b, a truncated version of the miR-17-92 cluster, is well-documented during tumor initiation. However, little is known about miR-17-19b function in established cancers. Here we investigate the role of miR-17-19b in c-MYC-driven lymphomas by integrating SILAC-based quantitative proteomics, transcriptomics and 3' untranslated region (UTR) analysis upon miR-17-19b overexpression. We identify over one hundred miR-17-19b targets, of which 40% are co-regulated by c-MYC. Downregulation of a new miR-17/20 target, checkpoint kinase 2 (Chek2), increases the recruitment of HuR to c-MYC transcripts, resulting in the inhibition of c-MYC translation and thus interfering with in vivo tumor growth. Hence, in established lymphomas, miR-17-19b fine-tunes c-MYC activity through a tight control of its function and expression, ultimately ensuring cancer cell homeostasis. Our data highlight the plasticity of miRNA function, reflecting changes in the mRNA landscape and 3' UTR shortening at different stages of tumorigenesis.

Network-based proteomic approaches reveal the neurodegenerative, neuroprotective and pain-related mechanisms involved after retrograde axonal damage.

Neurodegenerative processes are preceded by neuronal dysfunction and synaptic disconnection. Disconnection between spinal motoneuron (MN) soma and synaptic target leads either to a retrograde degenerative process or to a regenerative reaction, depending injury proximity among other factors. Distinguished key events associated with one or other processes may give some clues towards new therapeutical approaches based on boosting endogenous neuroprotective mechanisms. Root mechanical traction leads to retrograde MN degeneration, but share common initial molecular mechanisms with a regenerative process triggered by distal axotomy and suture. By 7 days post-injury, key molecular events starts to diverge and sign apart each destiny. We used comparative unbiased proteomics to define these signatures, coupled to a novel network-based analysis to get biological meaning. The procedure implicated the previous generation of combined topological information from manual curated 19 associated biological processes to be contrasted with the proteomic list using gene enrichment analysis tools. The novel and unexpected results suggested that motoneurodegeneration is better explained mainly by the concomitant triggering of anoikis, anti-apoptotic and neuropathic-pain related programs. In contrast, the endogenous neuroprotective mechanisms engaged after distal axotomy included specifically rather anti-anoikis and selective autophagy. Validated protein-nodes and processes are highlighted across discussion.

Profiling the secretome and extracellular proteome of the potato late blight pathogen Phytophthora infestans.

Oomycetes are filamentous organisms that cause notorious diseases, several of which have a high economic impact. Well known is Phytophthora infestans, the causal agent of potato late blight. Previously, in silico analyses of the genome and transcriptome of P. infestans resulted in the annotation of a large number of genes encoding proteins with an N-terminal signal peptide. This set is collectively referred to as the secretome and comprises proteins involved in, for example, cell wall growth and modification, proteolytic processes, and the promotion of successful invasion of plant cells. So far, proteomic profiling in oomycetes was primarily focused on subcellular, intracellular or cell wall fractions; the extracellular proteome has not been studied systematically. Here we present the first comprehensive characterization of the in vivo secretome and extracellular proteome of P. infestans. We have used mass spectrometry to analyze P. infestans proteins present in seven different growth media with mycelial cultures and this resulted in the consistent identification of over two hundred proteins. Gene ontology classification pinpointed proteins involved in cell wall modifications, pathogenesis, defense responses, and proteolytic processes. Moreover, we found members of the RXLR and CRN effector families as well as several proteins lacking an obvious signal peptide. The latter were confirmed to be bona fide extracellular proteins and this suggests that, similar to other organisms, oomycetes exploit non-conventional secretion mechanisms to transfer certain proteins to the extracellular environment.

Lifelong exercise training modulates cardiac mitochondrial phosphoproteome in rats.

Moderate physical activity has traditionally been associated with the improvement of cardiac function and, consequently, with the extension of life span. Mitochondria play a key role in the adaptation of heart muscle to exercise-related metabolic demands. In order to disclose the molecular mechanisms underlying the beneficial effect of lifelong physical activity in cardiac function, we performed label-free quantitative mass spectrometry-based proteomics of Sprague-Dawley rat heart mitochondrial proteome and phosphoproteome. Our data revealed that 54 weeks of moderate treadmill exercise modulates the abundance of proteins involved in the generation of precursor metabolites and cellular respiration, suggesting an increase in carbohydrate oxidation-based metabolism. Moreover, from the 1335 phosphopeptides identified in this study, 6 phosphosites were exclusively assigned to heart mitochondria from sedentary rats and 17 to exercised animals, corresponding to 6 and 16 proteins, respectively. Most proteins exhibiting significant alterations in specific phosphorylation sites were involved in metabolism. Analysis of the acquired data led to the identification of several kinases potentially modulated by exercise training, which were selected for further validation. Indeed, higher protein abundance levels of RAF and p38 in mitochondria were confirmed to be modulated by sustained exercise. Our work describes the plasticity of heart mitochondria in response to long exercise programs manifested by the reprogramming of phosphoproteome and provides evidence for the kinases involved in the regulation of metabolic pathways and mitochondrial maintenance.

Integrative quantitation enables a comprehensive proteome comparison of two Mycoplasma pneumoniae genetic perturbations.

Quantitative proteomics is an essential tool in proteome research since it enables measuring changes in protein abundance in response to biological perturbations. During the last few years, different quantitative strategies have been developed in proteomics to compare different experimental conditions, including label-free and isobaric chemical labeling approaches. Here we show that different quantitation techniques have an important influence on detected sample variability, and we use the combination of six different quantitation strategies to perform a proteome comparison of three different Mycoplasma pneumoniae strains (ldh knockdown, Δprkc, and wild-type). The integration of the different datasets indicates that the ldh knockdown strongly affects the abundance of ribosomal proteins and enzymes involved in the regulation of the cellular redox state, whereas the prkc deletion affects key cellular physiological processes such as protein and DNA synthesis, and cytoadherence.

Widespread generation of alternative UTRs contributes to sex-specific RNA binding by UNR.

Upstream of N-ras (UNR) is a conserved RNA-binding protein that regulates mRNA translation and stability by binding to sites generally located in untranslated regions (UTRs). In Drosophila, sex-specific binding of UNR to msl2 mRNA and the noncoding RNA roX is believed to play key roles in the control of X-chromosome dosage compensation in both sexes. To investigate broader sex-specific functions of UNR, we have identified its RNA targets in adult male and female flies by high-throughput RNA binding and transcriptome analysis. Here we show that UNR binds to a large set of protein-coding transcripts and to a smaller set of noncoding RNAs in a sex-specific fashion. The analyses also reveal a strong correlation between sex-specific binding of UNR and sex-specific differential expression of UTRs in target genes. Validation experiments indicate that UNR indeed recognizes sex-specifically processed transcripts. These results suggest that UNR exploits the transcript diversity generated by alternative processing and alternative promoter usage to bind and regulate target genes in a sex-specific manner.

HumMeth27QCReport: an R package for quality control and primary analysis of Illumina Infinium methylation data.

BACKGROUND: The study of the human DNA methylome has gained particular interest in the last few years. Researchers can nowadays investigate the potential role of DNA methylation in common disorders by taking advantage of new high-throughput technologies. Among these, Illumina Infinium assays can interrogate the methylation levels of hundreds of thousands of CpG sites, offering an ideal solution for genome-wide methylation profiling. However, like for other high-throughput technologies, the main bottleneck remains at the stage of data analysis rather than data production. FINDINGS: We have developed HumMeth27QCReport, an R package devoted to researchers wanting to quickly analyse their Illumina Infinium methylation arrays. This package automates quality control steps by generating a report including sample-independent and sample-dependent quality plots, and performs primary analysis of raw methylation calls by computing data normalization, statistics, and sample similarities. This package is available at CRAN repository, and can be integrated in any Galaxy instance through the implementation of ad-hoc scripts accessible at Galaxy Tool Shed. CONCLUSIONS: Our package provides users of the Illumina Infinium Methylation assays with a simplified, automated, open-source quality control and primary analysis of their methylation data. Moreover, to enhance its use by experimental researchers, the tool is being distributed along with the scripts necessary for its implementation in the Galaxy workbench. Finally, although it was originally developed for HumanMethylation27, we proved its compatibility with data generated with the HumanMethylation450 Bead Chip.

Substantial histone reduction modulates genomewide nucleosomal occupancy and global transcriptional output.

The basic unit of genome packaging is the nucleosome, and nucleosomes have long been proposed to restrict DNA accessibility both to damage and to transcription. Nucleosome number in cells was considered fixed, but recently aging yeast and mammalian cells were shown to contain fewer nucleosomes. We show here that mammalian cells lacking High Mobility Group Box 1 protein (HMGB1) contain a reduced amount of core, linker, and variant histones, and a correspondingly reduced number of nucleosomes, possibly because HMGB1 facilitates nucleosome assembly. Yeast nhp6 mutants lacking Nhp6a and -b proteins, which are related to HMGB1, also have a reduced amount of histones and fewer nucleosomes. Nucleosome limitation in both mammalian and yeast cells increases the sensitivity of DNA to damage, increases transcription globally, and affects the relative expression of about 10% of genes. In yeast nhp6 cells the loss of more than one nucleosome in four does not affect the location of nucleosomes and their spacing, but nucleosomal occupancy. The decrease in nucleosomal occupancy is non-uniform and can be modelled assuming that different nucleosomal sites compete for available histones. Sites with a high propensity to occupation are almost always packaged into nucleosomes both in wild type and nucleosome-depleted cells; nucleosomes on sites with low propensity to occupation are disproportionately lost in nucleosome-depleted cells. We suggest that variation in nucleosome number, by affecting nucleosomal occupancy both genomewide and gene-specifically, constitutes a novel layer of epigenetic regulation.