Matthew Loose

Establishment and cryptic transmission of Zika virus in Brazil and the Americas

Transmission of Zika virus (ZIKV) in the Americas was first confirmed in May 2015 in northeast Brazil. Brazil has had the highest number of reported ZIKV cases worldwide (more than 200,000 by 24 December 2016) and the most cases associated with microcephaly and other birth defects (2,366 confirmed by 31 December 2016). Since the initial detection of ZIKV in Brazil, more than 45 countries in the Americas have reported local ZIKV transmission, with 24 of these reporting severe ZIKV-associated disease. However, the origin and epidemic history of ZIKV in Brazil and the Americas remain poorly understood, despite the value of this information for interpreting observed trends in reported microcephaly. Here we address this issue by generating 54 complete or partial ZIKV genomes, mostly from Brazil, and reporting data generated by a mobile genomics laboratory that travelled across northeast Brazil in 2016. One sequence represents the earliest confirmed ZIKV infection in Brazil. Analyses of viral genomes with ecological and epidemiological data yield an estimate that ZIKV was present in northeast Brazil by February 2014 and is likely to have disseminated from there, nationally and internationally, before the first detection of ZIKV in the Americas. Estimated dates for the international spread of ZIKV from Brazil indicate the duration of pre-detection cryptic transmission in recipient regions. The role of northeast Brazil in the establishment of ZIKV in the Americas is further supported by geographic analysis of ZIKV transmission potential and by estimates of the basic reproduction number of the virus.

MinION Analysis and Reference Consortium: Phase 1 data release and analysis

The advent of a miniaturized DNA sequencing device with a high-throughput contextual sequencing capability embodies the next generation of large scale sequencing tools. The MinION™ Access Programme (MAP) was initiated by Oxford Nanopore Technologies™ in April 2014, giving public access to their USB-attached miniature sequencing device. The MinION Analysis and Reference Consortium (MARC) was formed by a subset of MAP participants, with the aim of evaluating and providing standard protocols and reference data to the community. Envisaged as a multi-phased project, this study provides the global community with the Phase 1 data from MARC, where the reproducibility of the performance of the MinION was evaluated at multiple sites. Five laboratories on two continents generated data using a control strain of Escherichia coli K-12, preparing and sequencing samples according to a revised ONT protocol. Here, we provide the details of the protocol used, along with a preliminary analysis of the characteristics of typical runs including the consistency, rate, volume and quality of data produced. Further analysis of the Phase 1 data presented here, and additional experiments in Phase 2 of E. coli from MARC are already underway to identify ways to improve and enhance MinION performance.

Nanopore sequencing and assembly of a human genome with ultra-long reads

We report the sequencing and assembly of a reference genome for the human GM12878 Utah/Ceph cell line using the MinION (Oxford Nanopore Technologies) nanopore sequencer. 91.2 Gb of sequence data, representing ∼30× theoretical coverage, were produced. Reference-based alignment enabled detection of large structural variants and epigenetic modifications. De novo assembly of nanopore reads alone yielded a contiguous assembly (NG50 ∼3 Mb). We developed a protocol to generate ultra-long reads (N50 > 100 kb, read lengths up to 882 kb). Incorporating an additional 5× coverage of these ultra-long reads more than doubled the assembly contiguity (NG50 ∼6.4 Mb). The final assembled genome was 2,867 million bases in size, covering 85.8% of the reference. Assembly accuracy, after incorporating complementary short-read sequencing data, exceeded 99.8%. Ultra-long reads enabled assembly and phasing of the 4-Mb major histocompatibility complex (MHC) locus in its entirety, measurement of telomere repeat length, and closure of gaps in the reference human genome assembly GRCh38.

Same-day diagnostic and surveillance data for tuberculosis via whole genome sequencing of direct respiratory samples.

ABSTRACT Routine full characterization of Mycobacterium tuberculosis is culture based, taking many weeks. Whole-genome sequencing (WGS) can generate antibiotic susceptibility profiles to inform treatment, augmented with strain information for global surveillance; such data could be transformative if provided at or near the point of care. We demonstrate a low-cost method of DNA extraction directly from patient samples for M. tuberculosis WGS. We initially evaluated the method by using the Illumina MiSeq sequencer (40 smear-positive respiratory samples obtained after routine clinical testing and 27 matched liquid cultures). M. tuberculosis was identified in all 39 samples from which DNA was successfully extracted. Sufficient data for antibiotic susceptibility prediction were obtained from 24 (62%) samples; all results were concordant with reference laboratory phenotypes. Phylogenetic placement was concordant between direct and cultured samples. With Illumina MiSeq/MiniSeq, the workflow from patient sample to results can be completed in 44/16 h at a reagent cost of £96/£198 per sample. We then employed a nonspecific PCR-based library preparation method for sequencing on an Oxford Nanopore Technologies MinION sequencer. We applied this to cultured Mycobacterium bovis strain BCG DNA and to combined culture-negative sputum DNA and BCG DNA. For flow cell version R9.4, the estimated turnaround time from patient to identification of BCG, detection of pyrazinamide resistance, and phylogenetic placement was 7.5 h, with full susceptibility results 5 h later. Antibiotic susceptibility predictions were fully concordant. A critical advantage of MinION is the ability to continue sequencing until sufficient coverage is obtained, providing a potential solution to the problem of variable amounts of M. tuberculosis DNA in direct samples.

Multiplex PCR method for MinION and Illumina sequencing of Zika and other virus genomes directly from clinical samples

Genome sequencing has become a powerful tool for studying emerging infectious diseases; however, genome sequencing directly from clinical samples (i.e., without isolation and culture) remains challenging for viruses such as Zika, for which metagenomic sequencing methods may generate insufficient numbers of viral reads. Here we present a protocol for generating coding-sequence-complete genomes, comprising an online primer design tool, a novel multiplex PCR enrichment protocol, optimized library preparation methods for the portable MinION sequencer (Oxford Nanopore Technologies) and the Illumina range of instruments, and a bioinformatics pipeline for generating consensus sequences. The MinION protocol does not require an Internet connection for analysis, making it suitable for field applications with limited connectivity. Our method relies on multiplex PCR for targeted enrichment of viral genomes from samples containing as few as 50 genome copies per reaction. Viral consensus sequences can be achieved in 1–2 d by starting with clinical samples and following a simple laboratory workflow. This method has been successfully used by several groups studying Zika virus evolution and is facilitating an understanding of the spread of the virus in the Americas. The protocol can be used to sequence other viral genomes using the online Primal Scheme primer designer software. It is suitable for sequencing either RNA or DNA viruses in the field during outbreaks or as an inexpensive, convenient method for use in the lab.

Transcriptional networks regulating hematopoietic cell fate decisions.

Purpose of reviewWe provide a summary of the temporal cascade of transcriptional networks giving rise to the hematopoietic stem cell (HSC) and controlling differentiation of the erythroid lineage from it. We focus on the mechanisms by which cell fate decisions are made and comment on recent developments and additions to the networks. Recent findingsA role for an SCL/LMO2 complex in HSC emergence, as well as in subsequent erythroid differentiation, has received support. Connections between the transcriptional networks and signaling molecules are being made but more work is needed in this area. Evidence that transcriptional cross-antagonistic switches underlie the choice between lineage pathways is increasing, and we highlight how the dynamics of earlier lineage decisions can influence later ones. Mathematical models are being built and reveal a surprising degree of power in these simple motifs to explain lineage choices. SummaryNew links in the transcriptional networks underlying cell-fate decisions are constantly emerging, and their incorporation into the evolving networks will make mathematical modeling more precise in its predictions of cell behavior, which can be tested experimentally.

Real-time selective sequencing using nanopore technology

The Oxford Nanopore Technologies MinION sequencer enables the selection of specific DNA molecules for sequencing by reversing the driving voltage across individual nanopores. To directly select molecules for sequencing, we used dynamic time warping to match reads to reference sequences. We demonstrate our open-source Read Until software in real-time selective sequencing of regions within small genomes, individual amplicon enrichment and normalization of an amplicon set.

Acquisition of germ plasm accelerates vertebrate evolution.

Tangling Evolutionary Trees Evolutionary rates tend to vary among taxa and may result in phylogenetic trees that do not reflect the true relationships among taxa, depending on the sequences input into the analysis. Examining vertebrate trees, Evans et al. (p. 200) demonstrate that differences in evolutionary rates, leading to phylogenetic distortions, are correlated with the mechanisms underlying germ cell formation. Evolutionary rate is faster in cases where germ cells are established by maternal molecules (“preformed”) relative to those that are induced during embryogenesis (“epigenesis”) in slowly evolving and, presumably, ancestral lineages. For example, frogs evolve more rapidly than salamanders, and teleosts more rapidly than ascipenseriform fishes. Thus, epigenesis constrains the ability of gene regulatory networks to change, with the repeated and convergent evolution of preformation eliminating this constraint. Evolutionary rates are phylogenetically correlated with how germ cells are specified. Primordial germ cell (PGC) specification occurs either by induction from pluripotent cells (epigenesis) or by a cell-autonomous mechanism mediated by germ plasm (preformation). Among vertebrates, epigenesis is basal, whereas germ plasm has evolved convergently across lineages and is associated with greater speciation. We compared protein-coding sequences of vertebrate species that employ preformation with their sister taxa that use epigenesis and demonstrate that genes evolve more rapidly in species containing germ plasm. Furthermore, differences in rates of evolution appear to cause phylogenetic incongruence in protein-coding sequence comparisons between vertebrate taxa. Our results support the hypothesis that germ plasm liberates constraints on somatic development and that enhanced evolvability drives the evolution of germ plasm.

A conserved mechanism for vertebrate mesoderm specification in urodele amphibians and mammals

Understanding how mesoderm is specified during development is a fundamental issue in biology, and it has been studied intensively in embryos from Xenopus. The gene regulatory network (GRN) for Xenopus is surprisingly complex and is not conserved in vertebrates, including mammals, which have single copies of the key genes Nodal and Mix. Why the Xenopus GRN should express multiple copies of Nodal and Mix genes is not known. To understand how these expanded gene families evolved, we investigated mesoderm specification in embryos from axolotls, representing urodele amphibians, since urodele embryology is basal to amphibians and was conserved during the evolution of amniotes, including mammals. We show that single copies of Nodal and Mix are required for mesoderm specification in axolotl embryos, suggesting the ancestral vertebrate state. Furthermore, we uncovered a novel genetic interaction in which Mix induces Brachyury expression, standing in contrast to the relationship of these molecules in Xenopus. However, we demonstrate that this functional relationship is conserved in mammals by showing that it is involved in the production of mesoderm from mouse embryonic stem cells. From our results, we produced an ancestral mesoderm (m)GRN, which we suggest is conserved in vertebrates. The results are discussed within the context of a theory in which the evolution of mechanisms governing early somatic development is constrained by the ancestral germ line-soma relationship, in which germ cells are produced by epigenesis.

5-hydroxymethyl-cytosine enrichment of non-committed cells is not a universal feature of vertebrate development

5-hydroxymethyl-cytosine (5-hmC) is a cytosine modification that is relatively abundant in mammalian pre-implantation embryos and embryonic stem cells (ESC) derived from mammalian blastocysts. Recent observations imply that both 5-hmC and Tet1/2/3 proteins, catalyzing the conversion of 5-methyl-cytosine to 5-hmC, may play an important role in self renewal and differentiation of ESCs. Here we assessed the distribution of 5-hmC in zebrafish and chick embryos and found that, unlike in mammals, 5-hmC is immunochemically undetectable in these systems before the onset of organogenesis. In addition, Tet1/2/3 transcripts are either low or undetectable at corresponding stages of zebrafish development. However, 5-hmC is enriched in later zebrafish and chick embryos and exhibits tissue-specific distribution in adult zebrafish. Our findings show that 5-hmC enrichment of non-committed cells is not a universal feature of vertebrate development and give insights both into evolution of embryonic pluripotency and the potential role of 5-hmC in its regulation.