TY - JOUR AB - Background: Many cancer genomes are extensively rearranged with highly aberrant chromosomal karyotypes. Structural and copy number variations in cancer genomes can be determined via abnormal mapping of sequenced reads to the reference genome. Recently it became possible to reconcile both of these types of large-scale variations into a karyotype graph representation of the rearranged cancer genomes. Such a representation, however, does not directly describe the linear and/or circular structure of the underlying rearranged cancer chromosomes, thus limiting possible analysis of cancer genomes somatic evolutionary process as well as functional genomic changes brought by the large-scale genome rearrangements. Results: Here we address the aforementioned limitation by introducing a novel methodological framework for recovering rearranged cancer chromosomes from karyotype graphs. For a cancer karyotype graph we formulate an Eulerian Decomposition Problem (EDP) of finding a collection of linear and/or circular rearranged cancer chromosomes that are determined by the graph. We derive and prove computational complexities for several variations of the EDP. We then demonstrate that Eulerian decomposition of the cancer karyotype graphs is not always unique and present the Consistent Contig Covering Problem (CCCP) of recovering unambiguous cancer contigs from the cancer karyotype graph, and describe a novel algorithm CCR capable of solving CCCP in polynomial time. We apply CCR on a prostate cancer dataset and demonstrate that it is capable of consistently recovering large cancer contigs even when underlying cancer genomes are highly rearranged. Conclusions: CCR can recover rearranged cancer contigs from karyotype graphs thereby addressing existing limitation in inferring chromosomal structures of rearranged cancer genomes and advancing our understanding of both patient/cancer-specific as well as the overall genetic instability in cancer. AU - Aganezov, Sergey AU - Zban, Ilya AU - Aksenov, Vitalii AU - Alexeev, Nikita AU - Schatz, Michael C. ID - 7214 JF - BMC Bioinformatics TI - Recovering rearranged cancer chromosomes from karyotype graphs VL - 20 ER - TY - JOUR AB - This is a literature teaching resource review for biologically inspired microfluidics courses or exploring the diverse applications of microfluidics. The structure is around key papers and model organisms. While courses gradually change over time, a focus remains on understanding how microfluidics has developed as well as what it can and cannot do for researchers. As a primary starting point, we cover micro-fluid mechanics principles and microfabrication of devices. A variety of applications are discussed using model prokaryotic and eukaryotic organisms from the set of bacteria (Escherichia coli), trypanosomes (Trypanosoma brucei), yeast (Saccharomyces cerevisiae), slime molds (Physarum polycephalum), worms (Caenorhabditis elegans), flies (Drosophila melangoster), plants (Arabidopsis thaliana), and mouse immune cells (Mus musculus). Other engineering and biochemical methods discussed include biomimetics, organ on a chip, inkjet, droplet microfluidics, biotic games, and diagnostics. While we have not yet reached the end-all lab on a chip, microfluidics can still be used effectively for specific applications. AU - Merrin, Jack ID - 7225 IS - 4 JF - Bioengineering TI - Frontiers in microfluidics, a teaching resource review VL - 6 ER - TY - CONF AB - Traditional concurrent programming involves manipulating shared mutable state. Alternatives to this programming style are communicating sequential processes (CSP) and actor models, which share data via explicit communication. These models have been known for almost half a century, and have recently had started to gain significant traction among modern programming languages. The common abstraction for communication between several processes is the channel. Although channels are similar to producer-consumer data structures, they have different semantics and support additional operations, such as the select expression. Despite their growing popularity, most known implementations of channels use lock-based data structures and can be rather inefficient. In this paper, we present the first efficient lock-free algorithm for implementing a communication channel for CSP programming. We provide implementations and experimental results in the Kotlin and Go programming languages. Our new algorithm outperforms existing implementations on many workloads, while providing non-blocking progress guarantee. Our design can serve as an example of how to construct general communication data structures for CSP and actor models. AU - Koval, Nikita AU - Alistarh, Dan-Adrian AU - Elizarov, Roman ID - 7228 SN - 0302-9743 T2 - 25th Anniversary of Euro-Par TI - Scalable FIFO channels for programming via communicating sequential processes VL - 11725 ER - TY - CONF AB - We present LiveTraVeL (Live Transit Vehicle Labeling), a real-time system to label a stream of noisy observations of transit vehicle trajectories with the transit routes they are serving (e.g., northbound bus #5). In order to scale efficiently to large transit networks, our system first retrieves a small set of candidate routes from a geometrically indexed data structure, then applies a fine-grained scoring step to choose the best match. Given that real-time data remains unavailable for the majority of the world’s transit agencies, these inferences can help feed a real-time map of a transit system’s trips, infer transit trip delays in real time, or measure and correct noisy transit tracking data. This system can run on vehicle observations from a variety of sources that don’t attach route information to vehicle observations, such as public imagery streams or user-contributed transit vehicle sightings.We abstract away the specifics of the sensing system and demonstrate the effectiveness of our system on a "semisynthetic" dataset of all New York City buses, where we simulate sensed trajectories by starting with fully labeled vehicle trajectories reported via the GTFS-Realtime protocol, removing the transit route IDs, and perturbing locations with synthetic noise. Using just the geometric shapes of the trajectories, we demonstrate that our system converges on the correct route ID within a few minutes, even after a vehicle switches from serving one trip to the next. AU - Osang, Georg F AU - Cook, James AU - Fabrikant, Alex AU - Gruteser, Marco ID - 7216 SN - 9781538670248 T2 - 2019 IEEE Intelligent Transportation Systems Conference TI - LiveTraVeL: Real-time matching of transit vehicle trajectories to transit routes at scale ER - TY - CONF AB - Piecewise Barrier Tubes (PBT) is a new technique for flowpipe overapproximation for nonlinear systems with polynomial dynamics, which leverages a combination of barrier certificates. PBT has advantages over traditional time-step based methods in dealing with those nonlinear dynamical systems in which there is a large difference in speed between trajectories, producing an overapproximation that is time independent. However, the existing approach for PBT is not efficient due to the application of interval methods for enclosure-box computation, and it can only deal with continuous dynamical systems without uncertainty. In this paper, we extend the approach with the ability to handle both continuous and hybrid dynamical systems with uncertainty that can reside in parameters and/or noise. We also improve the efficiency of the method significantly, by avoiding the use of interval-based methods for the enclosure-box computation without loosing soundness. We have developed a C++ prototype implementing the proposed approach and we evaluate it on several benchmarks. The experiments show that our approach is more efficient and precise than other methods in the literature. AU - Kong, Hui AU - Bartocci, Ezio AU - Jiang, Yu AU - Henzinger, Thomas A ID - 7231 SN - 0302-9743 T2 - 17th International Conference on Formal Modeling and Analysis of Timed Systems TI - Piecewise robust barrier tubes for nonlinear hybrid systems with uncertainty VL - 11750 ER - TY - JOUR AB - Coupling of endoplasmic reticulum stress to dimerisation‑dependent activation of the UPR transducer IRE1 is incompletely understood. Whilst the luminal co-chaperone ERdj4 promotes a complex between the Hsp70 BiP and IRE1's stress-sensing luminal domain (IRE1LD) that favours the latter's monomeric inactive state and loss of ERdj4 de-represses IRE1, evidence linking these cellular and in vitro observations is presently lacking. We report that enforced loading of endogenous BiP onto endogenous IRE1α repressed UPR signalling in CHO cells and deletions in the IRE1α locus that de-repressed the UPR in cells, encode flexible regions of IRE1LD that mediated BiP‑induced monomerisation in vitro. Changes in the hydrogen exchange mass spectrometry profile of IRE1LD induced by ERdj4 and BiP confirmed monomerisation and were consistent with active destabilisation of the IRE1LD dimer. Together, these observations support a competition model whereby waning ER stress passively partitions ERdj4 and BiP to IRE1LD to initiate active repression of UPR signalling. AU - Amin-Wetzel, Niko Paresh AU - Neidhardt, Lisa AU - Yan, Yahui AU - Mayer, Matthias P. AU - Ron, David ID - 7340 JF - eLife TI - Unstructured regions in IRE1α specify BiP-mediated destabilisation of the luminal domain dimer and repression of the UPR VL - 8 ER - TY - JOUR AB - Biochemical reactions often occur at low copy numbers but at once in crowded and diverse environments. Space and stochasticity therefore play an essential role in biochemical networks. Spatial-stochastic simulations have become a prominent tool for understanding how stochasticity at the microscopic level influences the macroscopic behavior of such systems. While particle-based models guarantee the level of detail necessary to accurately describe the microscopic dynamics at very low copy numbers, the algorithms used to simulate them typically imply trade-offs between computational efficiency and biochemical accuracy. eGFRD (enhanced Green’s Function Reaction Dynamics) is an exact algorithm that evades such trade-offs by partitioning the N-particle system into M ≤ N analytically tractable one- and two-particle systems; the analytical solutions (Green’s functions) then are used to implement an event-driven particle-based scheme that allows particles to make large jumps in time and space while retaining access to their state variables at arbitrary simulation times. Here we present “eGFRD2,” a new eGFRD version that implements the principle of eGFRD in all dimensions, thus enabling efficient particle-based simulation of biochemical reaction-diffusion processes in the 3D cytoplasm, on 2D planes representing membranes, and on 1D elongated cylinders representative of, e.g., cytoskeletal tracks or DNA; in 1D, it also incorporates convective motion used to model active transport. We find that, for low particle densities, eGFRD2 is up to 6 orders of magnitude faster than conventional Brownian dynamics. We exemplify the capabilities of eGFRD2 by simulating an idealized model of Pom1 gradient formation, which involves 3D diffusion, active transport on microtubules, and autophosphorylation on the membrane, confirming recent experimental and theoretical results on this system to hold under genuinely stochastic conditions. AU - Sokolowski, Thomas R AU - Paijmans, Joris AU - Bossen, Laurens AU - Miedema, Thomas AU - Wehrens, Martijn AU - Becker, Nils B. AU - Kaizu, Kazunari AU - Takahashi, Koichi AU - Dogterom, Marileen AU - ten Wolde, Pieter Rein ID - 7422 IS - 5 JF - The Journal of Chemical Physics SN - 0021-9606 TI - eGFRD in all dimensions VL - 150 ER - TY - CONF AB - Simple drawings of graphs are those in which each pair of edges share at most one point, either a common endpoint or a proper crossing. In this paper we study the problem of extending a simple drawing D(G) of a graph G by inserting a set of edges from the complement of G into D(G) such that the result is a simple drawing. In the context of rectilinear drawings, the problem is trivial. For pseudolinear drawings, the existence of such an extension follows from Levi’s enlargement lemma. In contrast, we prove that deciding if a given set of edges can be inserted into a simple drawing is NP-complete. Moreover, we show that the maximization version of the problem is APX-hard. We also present a polynomial-time algorithm for deciding whether one edge uv can be inserted into D(G) when {u,v} is a dominating set for the graph G. AU - Arroyo Guevara, Alan M AU - Derka, Martin AU - Parada, Irene ID - 7230 SN - 0302-9743 T2 - 27th International Symposium on Graph Drawing and Network Visualization TI - Extending simple drawings VL - 11904 ER - TY - CONF AB - We present Mixed-time Signal Temporal Logic (STL−MX), a specification formalism which extends STL by capturing the discrete/ continuous time duality found in many cyber-physical systems (CPS), as well as mixed-signal electronic designs. In STL−MX, properties of components with continuous dynamics are expressed in STL, while specifications of components with discrete dynamics are written in LTL. To combine the two layers, we evaluate formulas on two traces, discrete- and continuous-time, and introduce two interface operators that map signals, properties and their satisfaction signals across the two time domains. We show that STL-mx has the expressive power of STL supplemented with an implicit T-periodic clock signal. We develop and implement an algorithm for monitoring STL-mx formulas and illustrate the approach using a mixed-signal example. AU - Ferrere, Thomas AU - Maler, Oded AU - Nickovic, Dejan ID - 7232 SN - 0302-9743 T2 - 17th International Conference on Formal Modeling and Analysis of Timed Systems TI - Mixed-time signal temporal logic VL - 11750 ER - TY - JOUR AB - β1-integrins mediate cell–matrix interactions and their trafficking is important in the dynamic regulation of cell adhesion, migration and malignant processes, including cancer cell invasion. Here, we employ an RNAi screen to characterize regulators of integrin traffic and identify the association of Golgi-localized gamma ear-containing Arf-binding protein 2 (GGA2) with β1-integrin, and its role in recycling of active but not inactive β1-integrin receptors. Silencing of GGA2 limits active β1-integrin levels in focal adhesions and decreases cancer cell migration and invasion, which is in agreement with its ability to regulate the dynamics of active integrins. By using the proximity-dependent biotin identification (BioID) method, we identified two RAB family small GTPases, i.e. RAB13 and RAB10, as novel interactors of GGA2. Functionally, RAB13 silencing triggers the intracellular accumulation of active β1-integrin, and reduces integrin activity in focal adhesions and cell migration similarly to GGA2 depletion, indicating that both facilitate active β1-integrin recycling to the plasma membrane. Thus, GGA2 and RAB13 are important specificity determinants for integrin activity-dependent traffic. AU - Sahgal, Pranshu AU - Alanko, Jonna H AU - Icha, Jaroslav AU - Paatero, Ilkka AU - Hamidi, Hellyeh AU - Arjonen, Antti AU - Pietilä, Mika AU - Rokka, Anne AU - Ivaska, Johanna ID - 7420 IS - 11 JF - Journal of Cell Science SN - 0021-9533 TI - GGA2 and RAB13 promote activity-dependent β1-integrin recycling VL - 132 ER - TY - JOUR AB - We compare finite rank perturbations of the following three ensembles of complex rectangular random matrices: First, a generalised Wishart ensemble with one random and two fixed correlation matrices introduced by Borodin and Péché, second, the product of two independent random matrices where one has correlated entries, and third, the case when the two random matrices become also coupled through a fixed matrix. The singular value statistics of all three ensembles is shown to be determinantal and we derive double contour integral representations for their respective kernels. Three different kernels are found in the limit of infinite matrix dimension at the origin of the spectrum. They depend on finite rank perturbations of the correlation and coupling matrices and are shown to be integrable. The first kernel (I) is found for two independent matrices from the second, and two weakly coupled matrices from the third ensemble. It generalises the Meijer G-kernel for two independent and uncorrelated matrices. The third kernel (III) is obtained for the generalised Wishart ensemble and for two strongly coupled matrices. It further generalises the perturbed Bessel kernel of Desrosiers and Forrester. Finally, kernel (II), found for the ensemble of two coupled matrices, provides an interpolation between the kernels (I) and (III), generalising previous findings of part of the authors. AU - Akemann, Gernot AU - Checinski, Tomasz AU - Liu, Dangzheng AU - Strahov, Eugene ID - 7423 IS - 1 JF - Annales de l'Institut Henri Poincaré, Probabilités et Statistiques SN - 0246-0203 TI - Finite rank perturbations in products of coupled random matrices: From one correlated to two Wishart ensembles VL - 55 ER - TY - JOUR AB - X and Y chromosomes can diverge when rearrangements block recombination between them. Here we present the first genomic view of a reciprocal translocation that causes two physically unconnected pairs of chromosomes to be coinherited as sex chromosomes. In a population of the common frog (Rana temporaria), both pairs of X and Y chromosomes show extensive sequence differentiation, but not degeneration of the Y chromosomes. A new method based on gene trees shows both chromosomes are sex‐linked. Furthermore, the gene trees from the two Y chromosomes have identical topologies, showing they have been coinherited since the reciprocal translocation occurred. Reciprocal translocations can thus reshape sex linkage on a much greater scale compared with inversions, the type of rearrangement that is much better known in sex chromosome evolution, and they can greatly amplify the power of sexually antagonistic selection to drive genomic rearrangement. Two more populations show evidence of other rearrangements, suggesting that this species has unprecedented structural polymorphism in its sex chromosomes. AU - Toups, Melissa A AU - Rodrigues, Nicolas AU - Perrin, Nicolas AU - Kirkpatrick, Mark ID - 7421 IS - 8 JF - Molecular Ecology SN - 0962-1083 TI - A reciprocal translocation radically reshapes sex‐linked inheritance in the common frog VL - 28 ER - TY - CONF AB - Proofs of sequential work (PoSW) are proof systems where a prover, upon receiving a statement χ and a time parameter T computes a proof ϕ(χ,T) which is efficiently and publicly verifiable. The proof can be computed in T sequential steps, but not much less, even by a malicious party having large parallelism. A PoSW thus serves as a proof that T units of time have passed since χ was received. PoSW were introduced by Mahmoody, Moran and Vadhan [MMV11], a simple and practical construction was only recently proposed by Cohen and Pietrzak [CP18]. In this work we construct a new simple PoSW in the random permutation model which is almost as simple and efficient as [CP18] but conceptually very different. Whereas the structure underlying [CP18] is a hash tree, our construction is based on skip lists and has the interesting property that computing the PoSW is a reversible computation. The fact that the construction is reversible can potentially be used for new applications like constructing proofs of replication. We also show how to “embed” the sloth function of Lenstra and Weselowski [LW17] into our PoSW to get a PoSW where one additionally can verify correctness of the output much more efficiently than recomputing it (though recent constructions of “verifiable delay functions” subsume most of the applications this construction was aiming at). AU - Abusalah, Hamza M AU - Kamath Hosdurg, Chethan AU - Klein, Karen AU - Pietrzak, Krzysztof Z AU - Walter, Michael ID - 7411 SN - 0302-9743 T2 - Advances in Cryptology – EUROCRYPT 2019 TI - Reversible proofs of sequential work VL - 11477 ER - TY - JOUR AB - Background Synaptic vesicles (SVs) are an integral part of the neurotransmission machinery, and isolation of SVs from their host neuron is necessary to reveal their most fundamental biochemical and functional properties in in vitro assays. Isolated SVs from neurons that have been genetically engineered, e.g. to introduce genetically encoded indicators, are not readily available but would permit new insights into SV structure and function. Furthermore, it is unclear if cultured neurons can provide sufficient starting material for SV isolation procedures. New method Here, we demonstrate an efficient ex vivo procedure to obtain functional SVs from cultured rat cortical neurons after genetic engineering with a lentivirus. Results We show that ∼108 plated cortical neurons allow isolation of suitable SV amounts for functional analysis and imaging. We found that SVs isolated from cultured neurons have neurotransmitter uptake comparable to that of SVs isolated from intact cortex. Using total internal reflection fluorescence (TIRF) microscopy, we visualized an exogenous SV-targeted marker protein and demonstrated the high efficiency of SV modification. Comparison with existing methods Obtaining SVs from genetically engineered neurons currently generally requires the availability of transgenic animals, which is constrained by technical (e.g. cost and time) and biological (e.g. developmental defects and lethality) limitations. Conclusions These results demonstrate the modification and isolation of functional SVs using cultured neurons and viral transduction. The ability to readily obtain SVs from genetically engineered neurons will permit linking in situ studies to in vitro experiments in a variety of genetic contexts. AU - Mckenzie, Catherine AU - Spanova, Miroslava AU - Johnson, Alexander J AU - Kainrath, Stephanie AU - Zheden, Vanessa AU - Sitte, Harald H. AU - Janovjak, Harald L ID - 7406 JF - Journal of Neuroscience Methods SN - 0165-0270 TI - Isolation of synaptic vesicles from genetically engineered cultured neurons VL - 312 ER - TY - CONF AB - Most of today's distributed machine learning systems assume reliable networks: whenever two machines exchange information (e.g., gradients or models), the network should guarantee the delivery of the message. At the same time, recent work exhibits the impressive tolerance of machine learning algorithms to errors or noise arising from relaxed communication or synchronization. In this paper, we connect these two trends, and consider the following question: Can we design machine learning systems that are tolerant to network unreliability during training? With this motivation, we focus on a theoretical problem of independent interest-given a standard distributed parameter server architecture, if every communication between the worker and the server has a non-zero probability p of being dropped, does there exist an algorithm that still converges, and at what speed? The technical contribution of this paper is a novel theoretical analysis proving that distributed learning over unreliable network can achieve comparable convergence rate to centralized or distributed learning over reliable networks. Further, we prove that the influence of the packet drop rate diminishes with the growth of the number of parameter servers. We map this theoretical result onto a real-world scenario, training deep neural networks over an unreliable network layer, and conduct network simulation to validate the system improvement by allowing the networks to be unreliable. AU - Yu, Chen AU - Tang, Hanlin AU - Renggli, Cedric AU - Kassing, Simon AU - Singla, Ankit AU - Alistarh, Dan-Adrian AU - Zhang, Ce AU - Liu, Ji ID - 7437 SN - 9781510886988 T2 - 36th International Conference on Machine Learning, ICML 2019 TI - Distributed learning over unreliable networks VL - 2019-June ER - TY - JOUR AB - We develop a framework for the rigorous analysis of focused stochastic local search algorithms. These algorithms search a state space by repeatedly selecting some constraint that is violated in the current state and moving to a random nearby state that addresses the violation, while (we hope) not introducing many new violations. An important class of focused local search algorithms with provable performance guarantees has recently arisen from algorithmizations of the Lovász local lemma (LLL), a nonconstructive tool for proving the existence of satisfying states by introducing a background measure on the state space. While powerful, the state transitions of algorithms in this class must be, in a precise sense, perfectly compatible with the background measure. In many applications this is a very restrictive requirement, and one needs to step outside the class. Here we introduce the notion of measure distortion and develop a framework for analyzing arbitrary focused stochastic local search algorithms, recovering LLL algorithmizations as the special case of no distortion. Our framework takes as input an arbitrary algorithm of such type and an arbitrary probability measure and shows how to use the measure as a yardstick of algorithmic progress, even for algorithms designed independently of the measure. AU - Achlioptas, Dimitris AU - Iliopoulos, Fotis AU - Kolmogorov, Vladimir ID - 7412 IS - 5 JF - SIAM Journal on Computing SN - 0097-5397 TI - A local lemma for focused stochastical algorithms VL - 48 ER - TY - JOUR AB - Multiple importance sampling (MIS) has become an indispensable tool in Monte Carlo rendering, widely accepted as a near-optimal solution for combining different sampling techniques. But an MIS combination, using the common balance or power heuristics, often results in an overly defensive estimator, leading to high variance. We show that by generalizing the MIS framework, variance can be substantially reduced. Specifically, we optimize one of the combined sampling techniques so as to decrease the overall variance of the resulting MIS estimator. We apply the approach to the computation of direct illumination due to an HDR environment map and to the computation of global illumination using a path guiding algorithm. The implementation can be as simple as subtracting a constant value from the tabulated sampling density done entirely in a preprocessing step. This produces a consistent noise reduction in all our tests with no negative influence on run time, no artifacts or bias, and no failure cases. AU - Karlík, Ondřej AU - Šik, Martin AU - Vévoda, Petr AU - Skrivan, Tomas AU - Křivánek, Jaroslav ID - 7418 IS - 6 JF - ACM Transactions on Graphics SN - 0730-0301 TI - MIS compensation: Optimizing sampling techniques in multiple importance sampling VL - 38 ER - TY - JOUR AB - We consider Bose gases consisting of N particles trapped in a box with volume one and interacting through a repulsive potential with scattering length of order N−1 (Gross–Pitaevskii regime). We determine the ground state energy and the low-energy excitation spectrum, up to errors vanishing as N→∞. Our results confirm Bogoliubov’s predictions. AU - Boccato, Chiara AU - Brennecke, Christian AU - Cenatiempo, Serena AU - Schlein, Benjamin ID - 7413 IS - 2 JF - Acta Mathematica SN - 0001-5962 TI - Bogoliubov theory in the Gross–Pitaevskii limit VL - 222 ER - TY - JOUR AB - The study of parallel ecological divergence provides important clues to the operation of natural selection. Parallel divergence often occurs in heterogeneous environments with different kinds of environmental gradients in different locations, but the genomic basis underlying this process is unknown. We investigated the genomics of rapid parallel adaptation in the marine snail Littorina saxatilis in response to two independent environmental axes (crab-predation versus wave-action and low-shore versus high-shore). Using pooled whole-genome resequencing, we show that sharing of genomic regions of high differentiation between environments is generally low but increases at smaller spatial scales. We identify different shared genomic regions of divergence for each environmental axis and show that most of these regions overlap with candidate chromosomal inversions. Several inversion regions are divergent and polymorphic across many localities. We argue that chromosomal inversions could store shared variation that fuels rapid parallel adaptation to heterogeneous environments, possibly as balanced polymorphism shared by adaptive gene flow. AU - Morales, Hernán E. AU - Faria, Rui AU - Johannesson, Kerstin AU - Larsson, Tomas AU - Panova, Marina AU - Westram, Anja M AU - Butlin, Roger K. ID - 7393 IS - 12 JF - Science Advances SN - 2375-2548 TI - Genomic architecture of parallel ecological divergence: Beyond a single environmental contrast VL - 5 ER - TY - JOUR AB - Polymer additives can substantially reduce the drag of turbulent flows and the upperlimit, the so called “maximum drag reduction” (MDR) asymptote is universal, i.e. inde-pendent of the type of polymer and solvent used. Until recently, the consensus was that,in this limit, flows are in a marginal state where only a minimal level of turbulence activ-ity persists. Observations in direct numerical simulations using minimal sized channelsappeared to support this view and reported long “hibernation” periods where turbu-lence is marginalized. In simulations of pipe flow we find that, indeed, with increasingWeissenberg number (Wi), turbulence expresses long periods of hibernation if the domainsize is small. However, with increasing pipe length, the temporal hibernation continuouslyalters to spatio-temporal intermittency and here the flow consists of turbulent puffs sur-rounded by laminar flow. Moreover, upon an increase in Wi, the flow fully relaminarises,in agreement with recent experiments. At even larger Wi, a different instability is en-countered causing a drag increase towards MDR. Our findings hence link earlier minimalflow unit simulations with recent experiments and confirm that the addition of polymersinitially suppresses Newtonian turbulence and leads to a reverse transition. The MDRstate on the other hand results from a separate instability and the underlying dynamicscorresponds to the recently proposed state of elasto-inertial-turbulence (EIT). AU - Lopez Alonso, Jose M AU - Choueiri, George H AU - Hof, Björn ID - 7397 JF - Journal of Fluid Mechanics SN - 0022-1120 TI - Dynamics of viscoelastic pipe flow at low Reynolds numbers in the maximum drag reduction limit VL - 874 ER -