Transcriptional variability accelerates preleukemia by cell diversification and perturbation of protein synthesis.

Transcriptional variability facilitates stochastic cell diversification and can in turn underpin adaptation to stress or injury. We hypothesize that it may analogously facilitate progression of premalignancy to cancer. To investigate this, we initiated preleukemia in mouse cells with enhanced transcriptional variability due to conditional disruption of the histone lysine acetyltransferase gene Kat2a. By combining single-cell RNA sequencing of preleukemia with functional analysis of transformation, we show that Kat2a loss results in global variegation of cell identity and accumulation of preleukemic cells. Leukemia progression is subsequently facilitated by destabilization of ribosome biogenesis and protein synthesis, which confer a transient transformation advantage. The contribution of transcriptional variability to early cancer evolution reflects a generic role in promoting cell fate transitions, which, in the case of well-adapted malignancies, contrastingly differentiates and depletes cancer stem cells. That is, transcriptional variability confers forward momentum to cell fate systems, with differential multistage impact throughout cancer evolution.

Novel open reading frames in human accelerated regions and transposable elements reveal new leads to understand schizophrenia and bipolar disorder.

Schizophrenia (SCZ) and bipolar disorder are debilitating neuropsychiatric disorders arising from a combination of environmental and genetic factors. Novel open reading frames (nORFs) are genomic loci that give rise to previously uncharacterized transcripts and protein products. In our previous work, we have shown that nORFs can be biologically regulated and that they may play a role in cancer and rare diseases. More importantly, we have shown that nORFs may emerge in accelerated regions of the genome giving rise to species-specific functions. We hypothesize that nORFs represent a potentially important group of biological factors that may contribute to SCZ and bipolar disorder pathophysiology. Human accelerated regions (HARs) are genomic features showing human-lineage-specific rapid evolution that may be involved in biological regulation and have additionally been found to associate with SCZ genes. Transposable elements (TEs) are another set of genomic features that have been shown to regulate gene expression. As with HARs, their relevance to SCZ has also been suggested. Here, nORFs are investigated in the context of HARs and TEs. This work shows that nORFs whose expression is disrupted in SCZ and bipolar disorder are in close proximity to HARs and TEs and that some of them are significantly associated with SCZ and bipolar disorder genomic hotspots. We also show that nORF encoded proteins can form structures and potentially constitute novel drug targets.

In silico identification of novel open reading frames in Plasmodium falciparum oocyte and salivary gland sporozoites using proteogenomics framework.

BACKGROUND: Plasmodium falciparum causes the deadliest form of malaria, which remains one of the most prevalent infectious diseases. Unfortunately, the only licensed vaccine showed limited protection and resistance to anti-malarial drug is increasing, which can be largely attributed to the biological complexity of the parasite's life cycle. The progression from one developmental stage to another in P. falciparum involves drastic changes in gene expressions, where its infectivity to human hosts varies greatly depending on the stage. Approaches to identify candidate genes that are responsible for the development of infectivity to human hosts typically involve differential gene expression analysis between stages. However, the detection may be limited to annotated proteins and open reading frames (ORFs) predicted using restrictive criteria. METHODS: The above problem is particularly relevant for P. falciparum; whose genome annotation is relatively incomplete given its clinical significance. In this work, systems proteogenomics approach was used to address this challenge, as it allows computational detection of unannotated, novel Open Reading Frames (nORFs), which are neglected by conventional analyses. Two pairs of transcriptome/proteome were obtained from a previous study where one was collected in the mosquito-infectious oocyst sporozoite stage, and the other in the salivary gland sporozoite stage with human infectivity. They were then re-analysed using the proteogenomics framework to identify nORFs in each stage. RESULTS: Translational products of nORFs that map to antisense, intergenic, intronic, 3' UTR and 5' UTR regions, as well as alternative reading frames of canonical proteins were detected. Some of these nORFs also showed differential expression between the two life cycle stages studied. Their regulatory roles were explored through further bioinformatics analyses including the expression regulation on the parent reference genes, in silico structure prediction, and gene ontology term enrichment analysis. CONCLUSION: The identification of nORFs in P. falciparum sporozoites highlights the biological complexity of the parasite. Although the analyses are solely computational, these results provide a starting point for further experimental validation of the existence and functional roles of these nORFs.

Pan-cancer analysis of transcripts encoding novel open-reading frames (nORFs) and their potential biological functions.

Uncharacterized and unannotated open-reading frames, which we refer to as novel open reading frames (nORFs), may sometimes encode peptides that remain unexplored for novel therapeutic opportunities. To our knowledge, no systematic identification and characterization of transcripts encoding nORFs or their translation products in cancer, or in any other physiological process has been performed. We use our curated nORFs database (nORFs.org), together with RNA-Seq data from The Cancer Genome Atlas (TCGA) and Genotype-Expression (GTEx) consortiums, to identify transcripts containing nORFs that are expressed frequently in cancer or matched normal tissue across 22 cancer types. We show nORFs are subject to extensive dysregulation at the transcript level in cancer tissue and that a small subset of nORFs are associated with overall patient survival, suggesting that nORFs may have prognostic value. We also show that nORF products can form protein-like structures with post-translational modifications. Finally, we perform in silico screening for inhibitors against nORF-encoded proteins that are disrupted in stomach and esophageal cancer, showing that they can potentially be targeted by inhibitors. We hope this work will guide and motivate future studies that perform in-depth characterization of nORF functions in cancer and other diseases.

A platform for curated products from novel open reading frames prompts reinterpretation of disease variants.

Recent evidence from proteomics and deep massively parallel sequencing studies have revealed that eukaryotic genomes contain substantial numbers of as-yet-uncharacterized open reading frames (ORFs). We define these uncharacterized ORFs as novel ORFs (nORFs). nORFs in humans are mostly under 100 codons and are found in diverse regions of the genome, including in long noncoding RNAs, pseudogenes, 3' UTRs, 5' UTRs, and alternative reading frames of canonical protein coding exons. There is therefore a pressing need to evaluate the potential functional importance of these unannotated transcripts and proteins in biological pathways and human disease on a larger scale, rather than one at a time. In this study, we outline the creation of a valuable nORFs data set with experimental evidence of translation for the community, use measures of heritability and selection that reveal signals for functional importance, and show the potential implications for functional interpretation of genetic variants in nORFs. Our results indicate that some variants that were previously classified as being benign or of uncertain significance may have to be reinterpreted.

Evolutionary divergence of novel open reading frames in cichlids speciation.

Novel open reading frames (nORFs) with coding potential may arise from noncoding DNA. Not much is known about their emergence, functional role, fixation in a population or contribution to adaptive radiation. Cichlids fishes exhibit extensive phenotypic diversification and speciation. Encounters with new environments alone are not sufficient to explain this striking diversity of cichlid radiation because other taxa coexistent with the Cichlidae demonstrate lower species richness. Wagner et al. analyzed cichlid diversification in 46 African lakes and reported that both extrinsic environmental factors and intrinsic lineage-specific traits related to sexual selection have strongly influenced the cichlid radiation, which indicates the existence of unknown molecular mechanisms responsible for rapid phenotypic diversification, such as emergence of novel open reading frames (nORFs). In this study, we integrated transcriptomic and proteomic signatures from two tissues of two cichlids species, identified nORFs and performed evolutionary analysis on these nORF regions. Our results suggest that the time scale of speciation of the two species and evolutionary divergence of these nORF genomic regions are similar and indicate a potential role for these nORFs in speciation of the cichlid fishes.

PhosphoEffect: Prioritizing Variants On or Adjacent to Phosphorylation Sites through Their Effect on Kinase Recognition Motifs.

Phosphorylation sites often have key regulatory functions and are central to many cellular signaling pathways, so mutations that modify them have the potential to contribute to pathological states such as cancer. Although many classifiers exist for prioritization of coding genomic variants, to our knowledge none of them explicitly account for the alteration or creation of kinase recognition motifs that alter protein structure, function, regulation of activity, and interaction networks through modifying the pattern of phosphorylation. We present a novel computational pipeline that uses a random forest classifier to predict the pathogenicity of a variant, according to its direct or indirect effect on local phosphorylation sites and the predicted functional impact of perturbing a phosphorylation event. We call this classifier PhosphoEffect and find that it compares favorably and with increased accuracy to the existing classifier PolyPhen 2.2.2 when tested on a dataset of known variants enriched for phosphorylation sites and their neighbors.

Big data in digital healthcare: lessons learnt and recommendations for general practice.

Big Data will be an integral part of the next generation of technological developments-allowing us to gain new insights from the vast quantities of data being produced by modern life. There is significant potential for the application of Big Data to healthcare, but there are still some impediments to overcome, such as fragmentation, high costs, and questions around data ownership. Envisioning a future role for Big Data within the digital healthcare context means balancing the benefits of improving patient outcomes with the potential pitfalls of increasing physician burnout due to poor implementation leading to added complexity. Oncology, the field where Big Data collection and utilization got a heard start with programs like TCGA and the Cancer Moon Shot, provides an instructive example as we see different perspectives provided by the United States (US), the United Kingdom (UK) and other nations in the implementation of Big Data in patient care with regards to their centralization and regulatory approach to data. By drawing upon global approaches, we propose recommendations for guidelines and regulations of data use in healthcare centering on the creation of a unique global patient ID that can integrate data from a variety of healthcare providers. In addition, we expand upon the topic by discussing potential pitfalls to Big Data such as the lack of diversity in Big Data research, and the security and transparency risks posed by machine learning algorithms.

Loss of Kat2a enhances transcriptional noise and depletes acute myeloid leukemia stem-like cells.

Acute Myeloid Leukemia (AML) is an aggressive hematological malignancy with abnormal progenitor self-renewal and defective white blood cell differentiation. Its pathogenesis comprises subversion of transcriptional regulation, through mutation and by hijacking normal chromatin regulation. Kat2a is a histone acetyltransferase central to promoter activity, that we recently associated with stability of pluripotency networks, and identified as a genetic vulnerability in AML. Through combined chromatin profiling and single-cell transcriptomics of a conditional knockout mouse, we demonstrate that Kat2a contributes to leukemia propagation through preservation of leukemia stem-like cells. Kat2a loss impacts transcription factor binding and reduces transcriptional burst frequency in a subset of gene promoters, generating enhanced variability of transcript levels. Destabilization of target programs shifts leukemia cell fate out of self-renewal into differentiation. We propose that control of transcriptional variability is central to leukemia stem-like cell propagation, and establish a paradigm exploitable in different tumors and distinct stages of cancer evolution.

Less than 30% of patients with acute myeloid leukaemia ­ an aggressive cancer of the white blood cells ­ survive five years post-diagnosis. This disease disrupts the maturation of white blood cells, resulting in the accumulation of immature cells that multiply and survive but are incapable of completing their maturation process. Amongst these, a group of cancer cells known as leukemic stem cells is responsible for continually replenishing the leukaemia, thus perpetuating its growth. Cancers develop when cells in the body acquire changes or mutations to their genetic makeup. The mutations that lead to acute myeloid leukaemia often affect the activity of genes known as epigenetic regulators. These genes regulate which proteins and other molecules cells make by controlling the way in which cells 'read' their genetic instructions. The epigenetic regulator Kat2a is thought to 'tune' the frequency at which cells read their genetic instructions. This tuning mechanism decreases random fluctuations in the execution of the instructions cells receive to make proteins and other molecules. In turn, this helps to ensure that individual cells of the same type behave in a similar way, for example by keeping leukaemia cells in an immature state. Here, Domingues, Kulkarni et al. investigated whether interfering with Kat2a can make acute myeloid leukaemia less aggressive by allowing the immature white blood cells to mature. Domingues, Kulkarni et al. genetically engineered mice to remove Kat2a from blood cells on demand and then inserted a mutation that causes acute myeloid leukaemia. The experiments showed that the loss of Kat2a delayed the development of leukaemia in the mice and progressively depleted leukaemia stem cells, causing the disease to become less aggressive. The results also showed that loss of Kat2a caused more fluctuations in how the white blood cells read their genetic code, which resulted in more variability in the molecules they produced and increased the tendency of the cells to mature. These findings establish that loss of Kat2a causes leukaemia stem cells to mature and stop multiplying by untuning the frequency at which the cells read their genetic instructions. In the future, it may be possible to develop drugs that target human KAT2A to treat acute myeloid leukaemia.

A Complex Hierarchy of Avoidance Behaviors in a Single-Cell Eukaryote.

Complex behavior is associated with animals with nervous systems, but decision-making and learning also occur in non-neural organisms [1], including singly nucleated cells [2-5] and multi-nucleate synctia [6-8]. Ciliates are single-cell eukaryotes, widely dispersed in aquatic habitats [9], with an extensive behavioral repertoire [10-13]. In 1906, Herbert Spencer Jennings [14, 15] described in the sessile ciliate Stentor roeseli a hierarchy of responses to repeated stimulation, which are among the most complex behaviors reported for a singly nucleated cell [16, 17]. These results attracted widespread interest [18, 19] and exert continuing fascination [7, 20-22] but were discredited during the behaviorist orthodoxy by claims of non-reproducibility [23]. These claims were based on experiments with the motile ciliate Stentor coeruleus. We acquired and maintained the correct organism in laboratory culture and used micromanipulation and video microscopy to confirm Jennings' observations. Despite significant individual variation, not addressed by Jennings, S. roeseli exhibits avoidance behaviors in a characteristic hierarchy of bending, ciliary alteration, contractions, and detachment, which is distinct from habituation or conditioning. Remarkably, the choice of contraction versus detachment is consistent with a fair coin toss. Such behavioral complexity may have had an evolutionary advantage in protist ecosystems, and the ciliate cortex may have provided mechanisms for implementing such behavior prior to the emergence of multicellularity. Our work resurrects Jennings' pioneering insights and adds to the list of exceptional features, including regeneration [24], genome rearrangement [25], codon reassignment [26], and cortical inheritance [27], for which the ciliate clade is renowned.

Behavioural analysis of single-cell aneural ciliate, Stentor roeseli, using machine learning approaches.

There is still a significant gap between our understanding of neural circuits and the behaviours they compute-i.e. the computations performed by these neural networks (Carandini 2012 Nat. Neurosci.15, 507-509. (doi:10.1038/nn.3043)). Cellular decision-making processes, learning, behaviour and memory formation-all that have been only associated with animals with neural systems-have also been observed in many unicellular aneural organisms, namely Physarum, Paramecium and Stentor (Tang & Marshall2018 Curr. Biol.28, R1180-R1184. (doi:10.1016/j.cub.2018.09.015)). As these are fully functioning organisms, yet being unicellular, there is a much better chance to elucidate the detailed mechanisms underlying these learning processes in these organisms without the complications of highly interconnected neural circuits. An intriguing learning behaviour observed in Stentor roeseli (Jennings 1902 Am. J. Physiol. Legacy Content8, 23-60. (doi:10.1152/ajplegacy.1902.8.1.23)) when stimulated with carmine has left scientists puzzled for more than a century. So far, none of the existing learning paradigm can fully encapsulate this particular series of five characteristic avoidance reactions. Although we were able to observe all responses described in the literature and in a previous study (Dexter et al. 2019), they do not conform to any particular learning model. We then investigated whether models inferred from machine learning approaches, including decision tree, random forest and feed-forward artificial neural networks could infer and predict the behaviour of S. roeseli. Our results showed that an artificial neural network with multiple 'computational' neurons is inefficient at modelling the single-celled ciliate's avoidance reactions. This has highlighted the complexity of behaviours in aneural organisms. Additionally, this report will also discuss the significance of elucidating molecular details underlying learning and decision-making processes in these unicellular organisms, which could offer valuable insights that are applicable to higher animals.

The Macronuclear Genome of Stentor coeruleus Reveals Tiny Introns in a Giant Cell.

The giant, single-celled organism Stentor coeruleus has a long history as a model system for studying pattern formation and regeneration in single cells. Stentor [1, 2] is a heterotrichous ciliate distantly related to familiar ciliate models, such as Tetrahymena or Paramecium. The primary distinguishing feature of Stentor is its incredible size: a single cell is 1 mm long. Early developmental biologists, including T.H. Morgan [3], were attracted to the system because of its regenerative abilities-if large portions of a cell are surgically removed, the remnant reorganizes into a normal-looking but smaller cell with correct proportionality [2, 3]. These biologists were also drawn to Stentor because it exhibits a rich repertoire of behaviors, including light avoidance, mechanosensitive contraction, food selection, and even the ability to habituate to touch, a simple form of learning usually seen in higher organisms [4]. While early microsurgical approaches demonstrated a startling array of regenerative and morphogenetic processes in this single-celled organism, Stentor was never developed as a molecular model system. We report the sequencing of the Stentor coeruleus macronuclear genome and reveal key features of the genome. First, we find that Stentor uses the standard genetic code, suggesting that ciliate-specific genetic codes arose after Stentor branched from other ciliates. We also discover that ploidy correlates with Stentor's cell size. Finally, in the Stentor genome, we discover the smallest spliceosomal introns reported for any species. The sequenced genome opens the door to molecular analysis of single-cell regeneration in Stentor.

The Study of Posttranslational Modifications of Tau Protein by Nuclear Magnetic Resonance Spectroscopy: Phosphorylation of Tau Protein by ERK2 Recombinant Kinase and Rat Brain Extract, and Acetylation by Recombinant Creb-Binding Protein.

Nuclear magnetic resonance (NMR) spectroscopy can be used as an analytical tool to investigate posttranslational modifications of protein. NMR is a valuable tool to map the interaction regions of protein partners. Here, we present protocols that have been developed in the course of our studies of the neuronal Tau protein. Tau is found aggregated in the neurons of Alzheimer's disease patients. Development of the disease is accompanied by increased, abnormal phosphorylation and acetylation of Tau. We have used NMR to investigate how these posttranslational modifications of Tau affect the interactions with its partners. We present here detailed protocols of in vitro phosphorylation of Tau by recombinant kinase, ERK2, or kinase activity of rat brain extracts, and acetylation by recombinant Creb-binding protein (CBP) acetyltransferase. The analytical characterization of the modified Tau by NMR spectroscopy is additionally described.

Characterization of Neuronal Tau Protein as a Target of Extracellular Signal-regulated Kinase.

Tau neuronal protein has a central role in neurodegeneration and is implicated in Alzheimer disease development. Abnormal phosphorylation of Tau impairs its interaction with other proteins and is associated with its dysregulation in pathological conditions. Molecular mechanisms leading to hyperphosphorylation of Tau in pathological conditions are unknown. Here, we characterize phosphorylation of Tau by extracellular-regulated kinase (ERK2), a mitogen-activated kinase (MAPK) that responds to extracellular signals. Analysis ofin vitrophosphorylated Tau by activated recombinant ERK2 with nuclear magnetic resonance spectroscopy (NMR) reveals phosphorylation of 15 Ser/Thr sites.In vitrophosphorylation of Tau using rat brain extract and subsequent NMR analysis identifies the same sites. Phosphorylation with rat brain extract is known to transform Tau into an Alzheimer disease-like state. Our results indicate that phosphorylation of Tau by ERK2 alone is sufficient to produce the same characteristics. We further investigate the mechanism of ERK2 phosphorylation of Tau. Kinases are known to recognize their protein substrates not only by their specificity for a targeted Ser or Thr phosphorylation site but also by binding to linear-peptide motifs called docking sites. We identify two main ERK2 docking sites in Tau sequence using NMR. Our results suggest that ERK2 dysregulation in Alzheimer disease could lead to abnormal phosphorylation of Tau resulting in the pathology of the disease.

Quantitative profiling of peptides from RNAs classified as noncoding.

Only a small fraction of the mammalian genome codes for messenger RNAs destined to be translated into proteins, and it is generally assumed that a large portion of transcribed sequences--including introns and several classes of noncoding RNAs (ncRNAs)--do not give rise to peptide products. A systematic examination of translation and physiological regulation of ncRNAs has not been conducted. Here we use computational methods to identify the products of non-canonical translation in mouse neurons by analysing unannotated transcripts in combination with proteomic data. This study supports the existence of non-canonical translation products from both intragenic and extragenic genomic regions, including peptides derived from antisense transcripts and introns. Moreover, the studied novel translation products exhibit temporal regulation similar to that of proteins known to be involved in neuronal activity processes. These observations highlight a potentially large and complex set of biologically regulated translational events from transcripts formerly thought to lack coding potential.

Paradoxical results in perturbation-based signaling network reconstruction.

Mathematical models are extensively employed to understand physicochemical processes in biological systems. In the absence of detailed mechanistic knowledge, models are often based on network inference methods, which in turn rely upon perturbations to nodes by biochemical means. We have discovered a potential pitfall of the approach underpinning such methods when applied to signaling networks. We first show experimentally, and then explain mathematically, how even in the simplest signaling systems, perturbation methods may lead to paradoxical conclusions: for any given pair of two components X and Y, and depending upon the specific intervention on Y, either an activation or a repression of X could be inferred. This effect is of a different nature from incomplete network identification due to underdetermined data and is a phenomenon intrinsic to perturbations. Our experiments are performed in an in vitro minimal system, thus isolating the effect and showing that it cannot be explained by feedbacks due to unknown intermediates. Moreover, our in vitro system utilizes proteins from a pathway in mammalian (and other eukaryotic) cells that play a central role in proliferation, gene expression, differentiation, mitosis, cell survival, and apoptosis. This pathway is the perturbation target of contemporary therapies for various types of cancers. The results presented here show that the simplistic view of intracellular signaling networks being made up of activation and repression links is seriously misleading, and call for a fundamental rethinking of signaling network analysis and inference methods.

Post-translational modification: nature's escape from genetic imprisonment and the basis for dynamic information encoding.

We discuss protein post-translational modification (PTM) from an information processing perspective. PTM at multiple sites on a protein creates a combinatorial explosion in the number of potential 'mod-forms', or global patterns of modification. Distinct mod-forms can elicit distinct downstream responses, so that the overall response depends partly on the effectiveness of a particular mod-form to elicit a response and partly on the stoichiometry of that mod-form in the molecular population. We introduce the 'mod-form distribution'-the relative stoichiometries of each mod-form-as the most informative measure of a protein's state. Distinct mod-form distributions may summarize information about distinct cellular and physiological conditions and allow downstream processes to interpret this information accordingly. Such information 'encoding' by PTMs may facilitate evolution by weakening the need to directly link upstream conditions to downstream responses. Mod-form distributions provide a quantitative framework in which to interpret ideas of 'PTM codes' that are emerging in several areas of biology, as we show by reviewing examples of ion channels, GPCRs, microtubules, and transcriptional co-regulators. We focus particularly on examples other than the well-known 'histone code', to emphasize the pervasive use of information encoding in molecular biology. Finally, we touch briefly on new methods for measuring mod-form distributions.

Comparative analysis of Erk phosphorylation suggests a mixed strategy for measuring phospho-form distributions.

The functional impact of multisite protein phosphorylation can depend on both the numbers and the positions of phosphorylated sites-the global pattern of phosphorylation or 'phospho-form'-giving biological systems profound capabilities for dynamic information processing. A central problem in quantitative systems biology, therefore, is to measure the 'phospho-form distribution': the relative amount of each of the 2(n) phospho-forms of a protein with n-phosphorylation sites. We compared four potential methods-western blots with phospho-specific antibodies, peptide-based liquid chromatography (LC) and mass spectrometry (MS; pepMS), protein-based LC/MS (proMS) and nuclear magnetic resonance spectroscopy (NMR)-on differentially phosphorylated samples of the well-studied mitogen-activated protein kinase Erk2, with two phosphorylation sites. The MS methods were quantitatively consistent with each other and with NMR to within 10%, but western blots, while highly sensitive, showed significant discrepancies with MS. NMR also uncovered two additional phosphorylations, for which a combination of pepMS and proMS yielded an estimate of the 16-member phospho-form distribution. This combined MS strategy provides an optimal mixture of accuracy and coverage for quantifying distributions, but positional isomers remain a challenging problem.