Thousands of short open reading frames (sORFs) are translated outside of annotated coding sequences. Recent studies have pioneered searching for sORF-encoded microproteins in mass spectrometry (MS)-based proteomics and peptidomics datasets. Here, we assessed literature-reported MS-based identifications of unannotated human proteins. We find that studies vary by three orders of magnitude in the number of unannotated proteins they report. Of nearly 10,000 reported sORF-encoded peptides, 96% were unique to a single study, and 12% mapped to annotated proteins or proteoforms. Manual curation of a benchmark dataset of 406 manually evaluated spectra from 204 sORF-encoded proteins revealed large variation in peptide-spectrum match (PSM) quality between studies, with immunopeptidomics studies generally reporting higher quality PSMs than conventional enzymatic digests of whole cell lysates. We estimate that 65% of predicted sORF-encoded protein detections in immunopeptidomics studies were supported by high-quality PSMs versus 7.8% in non-immunopeptidomics datasets. Our work stresses the need for standardized protocols and analysis workflows to guide future advancements in microprotein detection by MS towards uncovering how many human microproteins exist.

Thousands of short open reading frames (sORFs) are translated outside of annotated coding sequences. Recent studies have pioneered searching for sORF-encoded microproteins in mass spectrometry (MS)-based proteomics and peptidomics datasets. Here, we assessed literature-reported MS-based identifications of unannotated human proteins. We find that studies vary by three orders of magnitude in the number of unannotated proteins they report. Of nearly 10,000 reported sORF-encoded peptides, 96% were unique to a single study, and 12% mapped to annotated proteins or proteoforms. Manual curation of a benchmark dataset of 406 manually evaluated spectra from 204 sORF-encoded proteins revealed large variation in peptide-spectrum match (PSM) quality between studies, with immunopeptidomics studies generally reporting higher quality PSMs than conventional enzymatic digests of whole cell lysates. We estimate that 65% of predicted sORF-encoded protein detections in immunopeptidomics studies were supported by high-quality PSMs versus 7.8% in non-immunopeptidomics datasets. Our work stresses the need for standardized protocols and analysis workflows to guide future advancements in microprotein detection by MS towards uncovering how many human microproteins exist.

Integrating cross-linking mass spectrometry (XL-MS) into structural biology workflows provides valuable information about the spatial arrangement of amino acid stretches, which can guide elucidation of protein assembly architecture. Additionally, the combination of XL-MS with peptide quantitation techniques is a powerful approach to delineate protein interface dynamics across diverse conditions. While XL-MS is increasingly effective with isolated proteins or small complexes, its application to whole-cell samples poses technical challenges related to analysis depth and throughput. The use of enrichable cross-linkers has greatly improved the detectability of protein interfaces in a proteome-wide scale, facilitating global protein-protein interaction mapping. Therefore, bringing together enrichable cross-linking and multiplexed peptide quantification is an appealing approach to enable comparative characterization of structural attributes of proteins and protein interactions. Here, we combined phospho-enrichable cross-linking with TMT labeling to develop a streamline workflow (PhoXplex) for the detection of differential structural features across a panel of cell lines in a global scale. We achieved deep coverage with quantification of over 9000 cross-links and long loop-links in total including potentially novel interactions. Overlaying AlphaFold predictions and disorder protein annotations enables exploration of the quantitative cross-linking data set, to reveal possible associations between mutations and protein structures. Lastly, we discuss current shortcomings and perspectives for deep whole-cell profiling of protein interfaces at large-scale.

Understanding how genetic variants impact molecular phenotypes is a key goal of functional genomics, currently hindered by reliance on a single haploid reference genome. Here, we present the EN-TEx resource of 1,635 open-access datasets from four donors (∼30 tissues × ∼15 assays). The datasets are mapped to matched, diploid genomes with long-read phasing and structural variants, instantiating a catalog of >1 million allele-specific loci. These loci exhibit coordinated activity along haplotypes and are less conserved than corresponding, non-allele-specific ones. Surprisingly, a deep-learning transformer model can predict the allele-specific activity based only on local nucleotide-sequence context, highlighting the importance of transcription-factor-binding motifs particularly sensitive to variants. Furthermore, combining EN-TEx with existing genome annotations reveals strong associations between allele-specific and GWAS loci. It also enables models for transferring known eQTLs to difficult-to-profile tissues (e.g., from skin to heart). Overall, EN-TEx provides rich data and generalizable models for more accurate personal functional genomics.

Abstract SF3B1 hotspot mutations are highly prevalent in several tumor types, where they lead to a global disruption of canonical splicing and are associated with a poor clinical outcome. Whilst evidence suggests that spliceosomal inhibitors have a therapeutic benefit for specific SF3B1 mutant cancers, recent trials have failed to report any clinical responses. Here we used a myriad of in vitro and in vivo techniques to identify and explore a synthetic lethal vulnerability in SF3B1 mutant cancers. These included a high throughput drug screen with an in-house curated library of 80 small-molecule inhibitors in an isogenic cell line model with the most prevalent SF3B1 mutation, a resistance CRISPR knockout screen and a range of DNA fibre analyses to gain mechanistic insights into the observed sensitivity. We identified that SF3B1 mutant cells are selectively sensitive to multiple PARP inhibitors (PARPi), which validated in numerous cell lines and patient derived models in vitro and in vivo, whereby they prevent metastasis. SF3B1 selective PARPi sensitivity occurs via a previously undescribed mechanism, as a result of excessive dormant origin firing and altered fork dynamics, due to the loss of the canonical replication stress response (pCHK1, pRPA, 53BP1), when SF3B1 mutant cells are treated with PARPi. This was further substantiated by the lack of MUS81 recruitment to under-replicated DNA intermediates, leading to the accumulation of DNA damage and cell cycle stalling, via the induction of the G2/M checkpoint. These findings demonstrate that SF3B1 mutant cells have wild-type-like replication under normal growth conditions, but fail to safeguard replication when exposed to PARPi. High content proteomic profiling identified that the defective ATR response is a direct consequence of the downregulation of CINP in SF3B1 mutant cells, which was confirmed in multiple isogenic models and in primary tumors. The defective ATR response was rescued with CINP re-expression and phenocopied through CINP silencing. In SF3B1 mutant cells G2/M checkpoint activation under PARP inhibition can be further abrogated by the addition of ATM inhibitors, but not by the inhibition of ATR or Chk1. Our discovery defines a unique mechanism underpinning PARPi sensitivity in SF3B1 mutant cancers. This synthetic-lethal relationship is further validated in a Phase 1 clinical trial, which shows that those patients with SF3B1 mutant relapsed leukemia have the longest progression free survival when treated with single agent PARPi. This novel synthetic lethality extends the utility of PARPi beyond BRCA1/2 mutant cancers and introduces a new therapeutic approach ready for proof-of-concept clinical trials. Citation Format: Philip Bland, Harry Saville, Patty T Wai, Lucinda Curnow, Jadwiga Nieminuszczy, Nivedita Ravindran, Holly E Barker, James Wright, Lu Yu, Ioanna Mavrommati, Barrie Peck, Patrycja Gazinska, Helen N Pemberton, Aditi Gulati, Naomi Guppy, Ioannis Roxanis, Guy Pratt, Ceri Oldreive, Tatjana Stankovic, Samantha Barlow, Helen Kalirai, Sarah E Coupland, Ronan Broderick, Samar Alsafadi, Alexandre Houy, Marc-Henri Stern, Stephen Pettit, Jyoti S Choudhary, Syed Haider, Wojciech Niedzwiedz, Christopher J Lord, Rachael Natrajan. SF3B1 hotspot mutations confer sensitivity to PARP inhibition through a defective replication stress response [abstract]. In: Proceedings of the AACR-NCI-EORTC Virtual International Conference on Molecular Targets and Cancer Therapeutics; 2023 Oct 11-15; Boston, MA. Philadelphia (PA): AACR; Mol Cancer Ther 2023;22(12 Suppl):Abstract nr B020.

Background: Heterozygous hotspot mutations in the RNA splicing factor SF3B1, occur in 3% of unselected breast cancers and are associated with oestrogen receptor (ER+) breast cancer (BC) where they are enriched in metastatic disease and are associated with a poor clinical outcome. SF3B1 mutations drive distinct signatures of alternative splicing through cryptic 3’ splice site selection leading to global transcriptomic and proteomic changes. The functional consequences of the mis-splicing events and resultant genetic vulnerabilities are poorly understood and precision medicine approaches that exploit these characteristics are not clinically available (Table 1). Methods: To understand the role of SF3B1 mutations in ER+ BC, we generated a series of SF3B1 mutant (SF3B1MUT) isogenic cell lines which were characterised using RNA-sequencing and high content mass-spectrometry proteomic profiling. SF3B1 interactome analysis was also performed using immunoprecipitation of SF3B1 followed by mass-spectrometry. The molecular consequences of aberrant splicing were investigated using a targeted screening approach of 280 genes predicted to be alternatively spliced in SF3B1MUT BC, while high-throughput drug screens were used to identify novel therapeutic options for patients with SF3B1MUT breast cancer using isogenic cells. Hits were validated in vitro and in vivo using cell line and patient derived xenografts. Results: Transcriptomic and proteomic profiling of SF3B1MUT cells identified global alternative 3’ splice site selection and subsequent proteomic changes induced by the mutations. Investigation of the SF3B1K700E interactome identified an enrichment of SF3B1K700E binding with ER, aberrant splicing of ER target genes, global rewiring of ER chromatin binding and resistance to endocrine therapy. Silencing of the aberrantly spliced candidate genes PPIH, TRIM37, HIGD1A, BRD9, and PHKG2 significantly enhanced the growth of the SF3B1 mutant cells, suggestive of a dose dependent tumour suppressive effect. Through synthetic-lethal drug screens we found that SF3B1MUT cells are selectively sensitive to PARP inhibitors. SF3B1MUT cells display a defective response to PARPi induced replication stress. Mechanistically, this occurs via defective ATR signalling in SF3B1MUT cells, which upon PARPi exposure leads to increased replication origin firing and loss of pChk1 (S317) induction. The resultant replication stress leads to failure to resolve DNA replication intermediates via the endonuclease MUS81 and cell cycle stalling at the G2/M checkpoint. These defects can be further targeted by ATM, CDK7 or FACT inhibition, when used in combination with PARPi treatment. This SF3B1MUT selective PARPi sensitivity is preserved across multiple cell lines and patient derived tumour models. In vivo, PARPi produce profound anti-tumour effects in multiple SF3B1MUT cancer models and eliminate distant metastases. Conclusions: Our integrative analysis reveals mechanistic insight into the role of SF3B1 mutations in endocrine therapy response in ER+ breast cancers, where altered SF3B1 induces ER-transcriptional re-programming. We further identified a robust synthetic-lethal relationship of mutant SF3B1 with PARP inhibition that is caused by a defective response to PARPi induced replication stress. Furthermore, we identified several potential selective combination strategies together with PARPi that are selective for SF3B1MUT cells. Together, these data provide the pre-clinical and mechanistic rationale for assessing already-approved PARPi in a biomarker-defined subset of advanced ER+ BC. Table 1. Identified potential therapies for SF3B1 mutant cancers from this study and the literature Citation Format: Phil Bland, Harry Saville, Abigail Read, Patty Wai, Gareth Muirhead, Lucinda Curnow, Jadwiga Nieminuszczy, Nivedita Ravindran, Marie John, Somaieh Hedayat, Holly Barker, James Wright, Lu Yu, Ioanna Mavrommati, Barrie Peck, Mark Allen, Patrycja Gazinska, Helen Pemberton, Aditi Gulati, Sarah Nash, Farzana Noor, Naomi Guppy, Ioannis Roxanis, Samantha Barlow, Helen Kalirai, Sarah Coupland, Ronan Broderick, Samar Alsafadi, Alexandre Houy, Marc-Henri Stern, Stephen Pettit, Jyoti Choudhary, Syed Haider, Wojciech Niedzwiedz, Christopher Lord, Rachael Natrajan. Mutations in the RNA Splicing Factor SF3B1 drive endocrine therapy resistance and confer a targetable replication stress response defect through PARP inhibition. [abstract]. In: Proceedings of the 2022 San Antonio Breast Cancer Symposium; 2022 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2023;83(5 Suppl):Abstract nr P6-10-05.

Oncolytic Reovirus type 3 Dearing (RT3D), is a naturally occurring double-stranded (ds) RNA virus that is under development as an oncolytic immunotherapy We used an unbiased high-throughput cytotoxicity screen of different targeted therapeutic agents with the aim of identifying potential drug-viral sensitizers to enhance RT3D tumour killing. Talazoparib, a clinical poly(ADP)-ribose polymerase 1 (PARP-1) inhibitor, was identified as a top hit and found to cause profound sensitisation to RT3D. This effect was not seen with other classes of oncolytic virus and was not mediated by enhanced viral replication or PARP inhibitor-related effects on the DNA damage response.RT3D interacts with retinoic acid-induced gene-1 (RIG-I) and activates PARP-1, with consequent PARylation of components of the extrinsic apoptosis pathway. Pharmacological and genetic inhibition of PARP-1 abrogates this PARylation and increases levels of extrinsic apoptosis, NF-kB signalling and pro-inflammatory cell death. Direct interaction between PARP-1 and RIG-I following RT3D/talazoparib treatment is a key factor in activating downstream signaling pathways that lead to IFN-β and TNF-α/TRAIL production which, in turn, amplify the therapeutic effect through positive feedback. Critically, it was possible to phenocopy the effect of RT3D through the use of non-viral ds-RNA therapy and RIG-I agonism. In in vivo studies, we demonstrated profound combinatorial efficacy of RT3D and talazoparib in human A375 melanoma in immunodeficient mice. More impressively, in immunocompetent mouse models of 4434 murine melanoma, we achieved 100% tumour control and protection from subsequent tumour rechallenge with the combination regimen. Correlative immunophenotyping confirmed significant innate and adaptive immune activation with the combination of RT3D and PARP inhibition. Taken together, these data provide a clear line of sight to clinical translation of combined regimens of PARP inhibition or ds-RNA agonism, with either viral or non-viral agents, in tumour types beyond the relatively narrow confines of current licensed indications for PARP inhibition.### Competing Interest StatementThe authors received RT3D from Oncolytics Biotech. Inc., but no financial support for the conduct of these studies

SF3B1 hotspot mutations are associated with a poor prognosis in several tumor types and lead to global disruption of canonical splicing. Through synthetic lethal drug screens, we identify that SF3B1 mutant ( SF3B1 MUT ) cells are selectively sensitive to poly (ADP-ribose) polymerase inhibitors (PARPi), independent of hotspot mutation and tumor site. SF3B1 MUT cells display a defective response to PARPi-induced replication stress that occurs via downregulation of the cyclin-dependent kinase 2 interacting protein (CINP), leading to increased replication fork origin firing and loss of phosphorylated CHK1 (pCHK1; S317) induction. This results in subsequent failure to resolve DNA replication intermediates and G 2 /M cell cycle arrest. These defects are rescued through CINP overexpression, or further targeted by a combination of ataxia-telangiectasia mutated and PARP inhibition. In vivo, PARPi produce profound antitumor effects in multiple SF3B1 MUT cancer models and eliminate distant metastases. These data provide the rationale for testing the clinical efficacy of PARPi in a biomarker-driven, homologous recombination proficient, patient population.

Cytoplasmic pattern recognition receptors (PRRs) for double-stranded RNA (RIG-I/MDA5) are key mediators of anti-viral responses. PRR agonists, such as dsRNA oncolytic Reovirus type 3 Dearing (Rt3D), potently activate RNA sensors. We used an unbiased cytotoxicity screen to reveal synergistic drug-virotherapy combinations and found potent effects of Rt3D combined with the CDK4/6 inhibitor, palbociclib. The combination augmented oncolytic virus-induced endoplasmic reticulum (ER) stress/unfolded protein response (UPR) and the expression and activation/signaling of RNA sensors. Combined Rt3D-palbociclib treatment potently increased interferon production and signaling, and knockdown studies implicated key UPR proteins and the RNA sensor, RIG-I, as essential to the phenotype observed. Further experiments, using canonical RIG-I agonists and an ER stress inducer, thapsigargin, confirmed cross-talk between RNA sensing and ER stress pathways that augmented cancer cell death and interferon production. Combined Rt3D-palbociclib also increased innate immune activation within tumour cells and IFN-induced HLA expression. Analysis of the immunopeptidome revealed changes to HLA-captured peptides with Rt3D-palbociclib, including altered expression of peptides from cancer/testis antigens (CTA) and endogenous retroviral elements (ERVs). Our findings highlight cross-talk between UPR signaling and RNA-mediated PRR activation as a means of enhancing anti-cancer efficacy with potential pro-immunogenic consequences. This has implications for future clinical development of PRR agonists and oncolytic viruses, and broadens the therapeutic remit of CDK4/6 inhibitors to include roles as both ER stress and dsRNA PRR sensitizers.

Multi-omics approaches including proteomics analyses are becoming an integral component of precision medicine. As clinical proteomics studies gain momentum and their sensitivity increases, research on identifying individuals based on their proteomics data is here examined for risks and ethics-related issues. A great deal of work has already been done on this topic for DNA/RNA sequencing data, but it has yet to be widely studied in other omics fields. The current state-of-the-art for the identification of individuals based solely on proteomics data is explained. Protein sequence variation analysis approaches are covered in more detail, including the available analysis workflows and their limitations. We also outline some previous forensic and omics proteomics studies that are relevant for the identification of individuals. Following that, we discuss the risks of patient reidentification using other proteomics data types such as protein expression abundance and post-translational modification (PTM) profiles. In light of the potential identification of individuals through proteomics data, possible legal and ethical implications are becoming increasingly important in the field.

Abstract Cytoplasmic nucleic acid sensors for double-stranded (ds) RNA (RIG-I/MDA5) and DNA (cGAS-STING) are pattern recognition receptors (PRRs) key to intracellular anti-viral responses. Recent research has highlighted roles for PRR agonists, including oncolytic virotherapy agents, in anti-tumor immunotherapy. Reovirus type 3 Dearing (Rt3D) is an oncolytic dsRNA virus with limited single-agent activity in clinical studies, but potential for use in combination regimens. We sought synergistic drug-virotherapy combinations using an unbiased screening approach that highlighted the CDK4/6 inhibitor, palbociclib, as a leading hit. We found that, when combined with Rt3D, palbociclib augmented oncolytic virus-induced endoplasmic reticulum (ER) stress/unfolded protein response (UPR) signaling. Combined Rt3D-palbociclib treatment potently increased interferon signaling and endogenous retroviral transcripts. Knockdown (siRNA) studies indicated key UPR proteins and the RNA sensor, RIG-I, were essential to the phenotype observed. Mechanistically independent experiments, using canonical RIG-I agonists and the ER stress inducer (thapsigargin), confirmed cross-talk between RNA sensing and ER stress pathways that augment cancer cell death and interferon production. Combined Rt3D-palbociclib increased innate immune activation and effector function. Our findings demonstrate that UPR signaling and innate immune RNA sensor crosstalk can be exploited to enhance anti-cancer efficacy with pro-immunogenic consequences. This has implications for future clinical development of PRR agonists and oncolytic viruses, as well as broadening the therapeutic remit of CDK4/6 inhibitors to include their role as ER stress sensitizers. Citation Format: Victoria Roulstone, Joan Kyula, Richard Elliot, Christopher J. Lord, Nik Matthews, Vicki Jennings, Harriet Whittock, David Mansfield, Jyoti Choudhary, James Wright, Lu Yu, Alan Melcher, Richard Vile, Matt Coffey, Martin McLaughlin, Kevin Harrington. Mechanisms of therapeutic synergy between pattern recognition response agonists and cdk4 inhibitors [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 1960.

Abstract The GENCODE project annotates human and mouse genes and transcripts supported by experimental data with high accuracy, providing a foundational resource that supports genome biology and clinical genomics. GENCODE annotation processes make use of primary data and bioinformatic tools and analysis generated both within the consortium and externally to support the creation of transcript structures and the determination of their function. Here, we present improvements to our annotation infrastructure, bioinformatics tools, and analysis, and the advances they support in the annotation of the human and mouse genomes including: the completion of first pass manual annotation for the mouse reference genome; targeted improvements to the annotation of genes associated with SARS-CoV-2 infection; collaborative projects to achieve convergence across reference annotation databases for the annotation of human and mouse protein-coding genes; and the first GENCODE manually supervised automated annotation of lncRNAs. Our annotation is accessible via Ensembl, the UCSC Genome Browser and https://www.gencodegenes.org.

Abstract Reovirus (RT3D) is a naturally occurring double-stranded RNA oncolytic virus that has shown preclinical efficacy in a wide range of tumor types. Early phase clinical studies have shown that this agent has modest monotherapy efficacy and can safely be combined with cytotoxic chemotherapy regimens. In the current studies, we used a high-throughput drug screen approach of different targeted therapeutic agents with the aim of looking for potential viral sensitizers that could enhance RT3D tumor killing. BMN-673 (talazoparib), a clinically approved poly(ADP)-ribose polymerase 1 (PARP-1) inhibitor was identified as a top hit and found to sensitize profoundly to RT3D both in vitro and in vivo in human xenograft tumors in a nude mouse model. We found that RT3D activated cellular PARP1 and was associated with PARylation of cellular proteins, including components of the DISC-associated cell death machinery. Combined treatment with RT3D and talazoparib enhanced extrinsic apoptosis (amplified by autocrine/paracrine TNF-α and TRAIL signaling), NF-κB pathway activity and pro-inflammatory cytokine production (CCL5/RANTES, CXCL8/IL8, CXCL1/GRO and CXCL10/IP10). Signaling was shown to be dependent on nucleic acid sensing mechanisms mediated by RIG-I and TLR3. We also found anti-tumour efficacy in an immunocompetent mouse model and this correlated with an increase in an immune response following combination treatment of RT3D and talazoparib. Our data provide a strong rationale for the combination of oncolytic RT3D with PARP1 inhibitors to exploit immunogenic response in cancer treatment. Citation Format: Joan N. Kyula, Victoria Roulstone, Richard Elliott, Harriet Whittock, Galabina Bozhanova, Martin McLaughlin, Malin Pedersen, Dragomir Krastev, Stephen Pettitt, Arnaud Legrand, Tencho Tenev, James Wright, Lu Yu, Jyoti Choudhary, Pascal Meier, Christopher J. Lord, Alan Melcher, Grey Wilkinson, Matt Coffey, Kevin J. Harrington. Talazoparib interacts with oncolytic reovirus to enhance death-inducing signaling complex (DISC)-mediated apoptosis and immune response [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 1932.

Evaluating the impact of genetic variants on transcriptional regulation is a central goal in biological science that has been constrained by reliance on a single reference genome. To address this, we constructed phased, diploid genomes for four cadaveric donors (using long-read sequencing) and systematically charted noncoding regulatory elements and transcriptional activity across more than 25 tissues from these donors. Integrative analysis revealed over a million variants with allele-specific activity, coordinated, locus-scale allelic imbalances, and structural variants impacting proximal chromatin structure. We relate the personal genome analysis to the ENCODE encyclopedia, annotating allele- and tissue-specific elements that are strongly enriched for variants impacting expression and disease phenotypes. These experimental and statistical approaches, and the corresponding EN-TEx resource, provide a framework for personalized functional genomics. decorations. The heatmap (center) shows the GWAS tag SNP enrichment of distal active CTCF+ cCREs. The tissues (colored according to the GTEx convention from Fig. 1) run along the bottom, and phenotypes (usually diseases, not explicitly labeled) run along the vertical axis. Around the heatmap, we show zoomed-in views highlighting various comparisons. At the bottom, we show that the tissue specificity of the enrichment for the trait Granulocyte % of Myeloid WBC is much stronger for AS versus non-AS cCREs. On the left, we show higher LDSC enrichment for AS elements compared to the corresponding non-AS ones for one tissue (coronary artery) across many associated traits. Left-bottom extends this comparison (for a single tissue-trait pair) to include roadmap annotations.

BACKGROUND:POLG, located on nuclear chromosome 15, encodes the DNA polymerase γ(Pol γ). Pol γ is responsible for the replication and repair of mitochondrial DNA (mtDNA). Pol γ is the only DNA polymerase found in mitochondria for most animal cells. Mutations in POLG are the most common single-gene cause of diseases of mitochondria and have been mapped over the coding region of the POLG ORF.RESULTS:Using PhyloCSF to survey alternative reading frames, we found a conserved coding signature in an alternative frame in exons 2 and 3 of POLG, herein referred to as ORF-Y that arose de novo in placental mammals. Using the synplot2 program, synonymous site conservation was found among mammals in the region of the POLG ORF that is overlapped by ORF-Y. Ribosome profiling data revealed that ORF-Y is translated and that initiation likely occurs at a CUG codon. Inspection of an alignment of mammalian sequences containing ORF-Y revealed that the CUG codon has a strong initiation context and that a well-conserved predicted RNA stem-loop begins 14 nucleotides downstream. Such features are associated with enhanced initiation at near-cognate non-AUG codons. Reanalysis of the Kim et al. (2014) draft human proteome dataset yielded two unique peptides that map unambiguously to ORF-Y. An additional conserved uORF, herein referred to as ORF-Z, was also found in exon 2 of POLG. Lastly, we surveyed Clinvar variants that are synonymous with respect to the POLG ORF and found that most of these variants cause amino acid changes in ORF-Y or ORF-Z.CONCLUSIONS:We provide evidence for a novel coding sequence, ORF-Y, that overlaps the POLG ORF. Ribosome profiling and mass spectrometry data show that ORF-Y is expressed. PhyloCSF and synplot2 analysis show that ORF-Y is subject to strong purifying selection. An abundance of disease-correlated mutations that map to exons 2 and 3 of POLG but also affect ORF-Y provides potential clinical significance to this finding.

Mass spectrometry (MS)-based quantitative proteomics experiments typically assay a subset of up to 60% of the ≈20 000 human protein coding genes. Computational methods for imputing the missing values using RNA expression data usually allow only for imputations of proteins measured in at least some of the samples. In silico methods for comprehensively estimating abundances across all proteins are still missing. Here, a novel method is proposed using deep learning to extrapolate the observed protein expression values in label-free MS experiments to all proteins, leveraging gene functional annotations and RNA measurements as key predictive attributes. This method is tested on four datasets, including human cell lines and human and mouse tissues. This method predicts the protein expression values with average R2 scores between 0.46 and 0.54, which is significantly better than predictions based on correlations using the RNA expression data alone. Moreover, it is demonstrated that the derived models can be "transferred" across experiments and species. For instance, the model derived from human tissues gave a R2=0.51 when applied to mouse tissue data. It is concluded that protein abundances generated in label-free MS experiments can be computationally predicted using functional annotated attributes and can be used to highlight aberrant protein abundance values.

The accurate identification and description of the genes in the human and mouse genomes is a fundamental requirement for high quality analysis of data informing both genome biology and clinical genomics. Over the last 15 years, the GENCODE consortium has been producing reference quality gene annotations to provide this foundational resource. The GENCODE consortium includes both experimental and computational biology groups who work together to improve and extend the GENCODE gene annotation. Specifically, we generate primary data, create bioinformatics tools and provide analysis to support the work of expert manual gene annotators and automated gene annotation pipelines. In addition, manual and computational annotation workflows use any and all publicly available data and analysis, along with the research literature to identify and characterise gene loci to the highest standard. GENCODE gene annotations are accessible via the Ensembl and UCSC Genome Browsers, the Ensembl FTP site, Ensembl Biomart, Ensembl Perl and REST APIs as well as https://www.gencodegenes.org.

The most widely appreciated role of DNA is to encode protein, yet the exact portion of the human genome that is translated remains to be ascertained. We previously developed PhyloCSF, a widely used tool to identify evolutionary signatures of protein-coding regions using multispecies genome alignments. Here, we present the first whole-genome PhyloCSF prediction tracks for human, mouse, chicken, fly, worm, and mosquito. We develop a workflow that uses machine learning to predict novel conserved protein-coding regions and efficiently guide their manual curation. We analyze more than 1000 high-scoring human PhyloCSF regions and confidently add 144 conserved protein-coding genes to the GENCODE gene set, as well as additional coding regions within 236 previously annotated protein-coding genes, and 169 pseudogenes, most of them disabled after primates diverged. The majority of these represent new discoveries, including 70 previously undetected protein-coding genes. The novel coding genes are additionally supported by single-nucleotide variant evidence indicative of continued purifying selection in the human lineage, coding-exon splicing evidence from new GENCODE transcripts using next-generation transcriptomic data sets, and mass spectrometry evidence of translation for several new genes. Our discoveries required simultaneous comparative annotation of other vertebrate genomes, which we show is essential to remove spurious ORFs and to distinguish coding from pseudogene regions. Our new coding regions help elucidate disease-associated regions by revealing that 118 GWAS variants previously thought to be noncoding are in fact protein altering. Altogether, our PhyloCSF data sets and algorithms will help researchers seeking to interpret these genomes, while our new annotations present exciting loci for further experimental characterization.

The early infantile epileptic encephalopathies (EIEE) are a group of rare, severe neurodevelopmental disorders, where even the most thorough sequencing studies leave 60-65% of patients without a molecular diagnosis. Here, we explore the incompleteness of transcript models used for exome and genome analysis as one potential explanation for lack of current diagnoses. Therefore, we have updated the GENCODE gene annotation for 191 epilepsy-associated genes, using human brain-derived transcriptomic libraries and other data to build 3,550 novel putative transcript models. The extended transcriptional footprint of these genes allowed for 294 intronic or intergenic variants, found in human mutation databases, to be reclassified as exonic, while a further 70 intronic variants were reclassified as splice-site proximal. Using SCN1A as a case study due to its close phenotype/genotype correlation with Dravet syndrome, we screened 122 people with Dravet syndrome, or a similar phenotype, with a panel of novel exon sequences representing eight established genes and identified two de novo SCN1A variants that now, through improved gene annotation can be ascribed to residing among novel exons. These two (from 122 screened patients, 1.6%) new molecular diagnoses carry significant clinical implications. Furthermore, we identified a previously-classified SCN1A intronic Dravet-associated variant that now lies within a deeply conserved novel exon. Our findings illustrate the potential gains of thorough gene annotation in improving diagnostic yields for genetic disorders. We would expect to find new molecular diagnoses in our 191 genes that were originally suspected by clinicians for patients, with a negative diagnosis.

Infection with Citrobacter rodentium triggers robust tissue damage repair responses, manifested by secretion of IL-22, in the absence of which mice succumbed to the infection. Of the main hallmarks of C. rodentium infection are colonic crypt hyperplasia (CCH) and dysbiosis. In order to colonize the host and compete with the gut microbiota, C. rodentium employs a type III secretion system (T3SS) that injects effectors into colonic intestinal epithelial cells (IECs). Once injected, the effectors subvert processes involved in innate immune responses, cellular metabolism and oxygenation of the mucosa. Importantly, the identity of the effector/s triggering the tissue repair response is/are unknown. Here we report that the effector EspO ,an orthologue of OspE found in Shigella spp, affects proliferation of IECs 8 and 14 days post C. rodentium infection as well as secretion of IL-22 from colonic explants. While we observed no differences in the recruitment of group 3 innate lymphoid cells (ILC3s) and T cells, which are the main sources of IL-22 at the early and late stages of C. rodentium infection respectively, infection with ΔespO was characterized by diminished recruitment of sub-mucosal neutrophils, which coincided with lower abundance of Mmp9 and chemokines (e.g. S100a8/9) in IECs. Moreover, mice infected with ΔespO triggered significantly lesser nutritional immunity (e.g. calprotectin, Lcn2) and expression of antimicrobial peptides (Reg3β, Reg3γ) compared to mice infected with WT C. rodentium. This overlapped with a decrease in STAT3 phosphorylation in IECs. Importantly, while the reduced CCH and abundance of antimicrobial proteins during ΔespO infection did not affect C. rodentium colonization or the composition of commensal Proteobacteria, they had a subtle consequence on Firmicutes subpopulations. EspO is the first bacterial virulence factor that affects neutrophil recruitment and secretion of IL-22, as well as expression of antimicrobial and nutritional immunity proteins in IECs.