Bead-linked transposomes enable a normalization-free workflow for NGS library preparation.

BACKGROUND: Transposome-based technologies have enabled the streamlined production of sequencer-ready DNA libraries; however, current methods are highly sensitive to the amount and quality of input nucleic acid. RESULTS: We describe a new library preparation technology (Nextera DNA Flex) that utilizes a known concentration of transposomes conjugated directly to beads to bind a fixed amount of DNA, and enables direct input of blood and saliva using an integrated extraction protocol. We further report results from libraries generated outside the standard parameters of the workflow, highlighting novel applications for Nextera DNA Flex, including human genome builds and variant calling from below 1 ng DNA input, customization of insert size, and preparation of libraries from short fragments and severely degraded FFPE samples. Using this bead-linked library preparation method, library yield saturation was observed at an input amount of 100 ng. Preparation of libraries from a range of species with varying GC levels demonstrated uniform coverage of small genomes. For large and complex genomes, coverage across the genome, including difficult regions, was improved compared with other library preparation methods. Libraries were successfully generated from amplicons of varying sizes (from 50 bp to 11 kb), however, a decrease in efficiency was observed for amplicons smaller than 250 bp. This library preparation method was also compatible with poor-quality DNA samples, with sequenceable libraries prepared from formalin-fixed paraffin-embedded samples with varying levels of degradation. CONCLUSIONS: In contrast to solution-based library preparation, this bead-based technology produces a normalized, sequencing-ready library for a wide range of DNA input types and amounts, largely obviating the need for DNA quantitation. The robustness of this bead-based library preparation kit and flexibility of input DNA facilitates application across a wide range of fields.

NxTrim: optimized trimming of Illumina mate pair reads.

MOTIVATION: Mate pair protocols add to the utility of paired-end sequencing by boosting the genomic distance spanned by each pair of reads, potentially allowing larger repeats to be bridged and resolved. The Illumina Nextera Mate Pair (NMP) protocol uses a circularization-based strategy that leaves behind 38-bp adapter sequences, which must be computationally removed from the data. While 'adapter trimming' is a well-studied area of bioinformatics, existing tools do not fully exploit the particular properties of NMP data and discard more data than is necessary. RESULTS: We present NxTrim, a tool that strives to discard as little sequence as possible from NMP reads. NxTrim makes full use of the sequence on both sides of the adapter site to build 'virtual libraries' of mate pairs, paired-end reads and single-ended reads. For bacterial data, we show that aggregating these datasets allows a single NMP library to yield an assembly whose quality compares favourably to that obtained from regular paired-end reads. AVAILABILITY AND IMPLEMENTATION: The source code is available at https://github.com/sequencing/NxTrim

A comprehensive catalogue of somatic mutations from a human cancer genome.

All cancers carry somatic mutations. A subset of these somatic alterations, termed driver mutations, confer selective growth advantage and are implicated in cancer development, whereas the remainder are passengers. Here we have sequenced the genomes of a malignant melanoma and a lymphoblastoid cell line from the same person, providing the first comprehensive catalogue of somatic mutations from an individual cancer. The catalogue provides remarkable insights into the forces that have shaped this cancer genome. The dominant mutational signature reflects DNA damage due to ultraviolet light exposure, a known risk factor for malignant melanoma, whereas the uneven distribution of mutations across the genome, with a lower prevalence in gene footprints, indicates that DNA repair has been preferentially deployed towards transcribed regions. The results illustrate the power of a cancer genome sequence to reveal traces of the DNA damage, repair, mutation and selection processes that were operative years before the cancer became symptomatic.

Rapid single-colony whole-genome sequencing of bacterial pathogens.

OBJECTIVES: As a result of the introduction of rapid benchtop sequencers, the time required to subculture a bacterial pathogen to extract sufficient DNA for library preparation can now exceed the time to sequence said DNA. We have eliminated this rate-limiting step by developing a protocol to generate DNA libraries for whole-genome sequencing directly from single bacterial colonies grown on primary culture plates. METHODS: We developed our protocol using single colonies of 17 bacterial pathogens responsible for severe human infection that were grown using standard diagnostic media and incubation conditions. We then applied this method to four clinical scenarios that currently require time-consuming reference laboratory tests: full identification and genotyping of salmonellae; identification of blaNDM-1, a highly transmissible carbapenemase resistance gene, in Klebsiella pneumoniae; detection of genes encoding staphylococcal toxins associated with specific disease syndromes; and monitoring of vaccine targets to detect vaccine escape in Neisseria meningitidis. RESULTS: We validated our single-colony whole-genome sequencing protocol for all 40 combinations of pathogen and selective, non-selective or indicator media tested in this study. Moreover, we demonstrated the clinical value of this method compared with current reference laboratory tests. CONCLUSIONS: This advance will facilitate the implementation of whole-genome sequencing into diagnostic and public health microbiology.

Accurate whole human genome sequencing using reversible terminator chemistry.

DNA sequence information underpins genetic research, enabling discoveries of important biological or medical benefit. Sequencing projects have traditionally used long (400-800 base pair) reads, but the existence of reference sequences for the human and many other genomes makes it possible to develop new, fast approaches to re-sequencing, whereby shorter reads are compared to a reference to identify intraspecies genetic variation. Here we report an approach that generates several billion bases of accurate nucleotide sequence per experiment at low cost. Single molecules of DNA are attached to a flat surface, amplified in situ and used as templates for synthetic sequencing with fluorescent reversible terminator deoxyribonucleotides. Images of the surface are analysed to generate high-quality sequence. We demonstrate application of this approach to human genome sequencing on flow-sorted X chromosomes and then scale the approach to determine the genome sequence of a male Yoruba from Ibadan, Nigeria. We build an accurate consensus sequence from >30x average depth of paired 35-base reads. We characterize four million single-nucleotide polymorphisms and four hundred thousand structural variants, many of which were previously unknown. Our approach is effective for accurate, rapid and economical whole-genome re-sequencing and many other biomedical applications.

DNA supercoiling enables the type IIS restriction enzyme BspMI to recognise the relative orientation of two DNA sequences.

Many proteins can sense the relative orientations of two sequences at distant locations in DNA: some require sites in inverted (head-to-head) orientation, others in repeat (head-to-tail) orientation. Like many restriction enzymes, the BspMI endonuclease binds two copies of its target site before cleaving DNA. Its target is an asymmetric sequence so two sites in repeat orientation differ from sites in inverted orientation. When tested against supercoiled plasmids with two sites 700 bp apart in either repeated or inverted orientations, BspMI had a higher affinity for the plasmid with repeated sites than the plasmid with inverted sites. In contrast, on linear DNA or on supercoiled DNA with sites 1605 bp apart, BspMI interacted equally with repeated or inverted sites. The ability of BspMI to detect the relative orientation of two DNA sequences thus depends on both the topology and the length of the intervening DNA. Supercoiling may restrain the juxtaposition of sites 700 bp apart to a particular alignment across the superhelical axis, but the juxtaposition of sites in linear DNA or far apart in supercoiled DNA may occur without restraint. BspMI can therefore act as a sensor of the conformational dynamics of supercoiled DNA.

Efficient de novo assembly of single-cell bacterial genomes from short-read data sets.

Whole genome amplification by the multiple displacement amplification (MDA) method allows sequencing of DNA from single cells of bacteria that cannot be cultured. Assembling a genome is challenging, however, because MDA generates highly nonuniform coverage of the genome. Here we describe an algorithm tailored for short-read data from single cells that improves assembly through the use of a progressively increasing coverage cutoff. Assembly of reads from single Escherichia coli and Staphylococcus aureus cells captures >91% of genes within contigs, approaching the 95% captured from an assembly based on many E. coli cells. We apply this method to assemble a genome from a single cell of an uncultivated SAR324 clade of Deltaproteobacteria, a cosmopolitan bacterial lineage in the global ocean. Metabolic reconstruction suggests that SAR324 is aerobic, motile and chemotaxic. Our approach enables acquisition of genome assemblies for individual uncultivated bacteria using only short reads, providing cell-specific genetic information absent from metagenomic studies.

Genomic diversity among drug sensitive and multidrug resistant isolates of Mycobacterium tuberculosis with identical DNA fingerprints.

BACKGROUND: Mycobacterium tuberculosis complex (MTBC), the causative agent of tuberculosis (TB), is characterized by low sequence diversity making this bacterium one of the classical examples of a genetically monomorphic pathogen. Because of this limited DNA sequence variation, routine genotyping of clinical MTBC isolates for epidemiological purposes relies on highly discriminatory DNA fingerprinting methods based on mobile and repetitive genetic elements. According to the standard view, isolates exhibiting the same fingerprinting pattern are considered direct progeny of the same bacterial clone, and most likely reflect ongoing transmission or disease relapse within individual patients. METHODOLOGY/PRINCIPAL FINDINGS: Here we further investigated this assumption and used massively parallel whole-genome sequencing to compare one drug-susceptible (K-1) and one multidrug resistant (MDR) isolate (K-2) of a rapidly spreading M. tuberculosis Beijing genotype clone from a high incidence region (Karakalpakstan, Uzbekistan). Both isolates shared the same IS6110 RFLP pattern and the same allele at 23 out of 24 MIRU-VNTR loci. We generated 23.9 million (K-1) and 33.0 million (K-2) paired 50 bp purity filtered reads corresponding to a mean coverage of 483.5 fold and 656.1 fold respectively. Compared with the laboratory strain H37Rv both Beijing isolates shared 1,209 SNPs. The two Beijing isolates differed by 130 SNPs and one large deletion. The susceptible isolate had 55 specific SNPs, while the MDR variant had 75 specific SNPs, including the five known resistance-conferring mutations. CONCLUSIONS: Our results suggest that M. tuberculosis isolates exhibiting identical DNA fingerprinting patterns can harbour substantial genomic diversity. Because this heterogeneity is not captured by traditional genotyping of MTBC, some aspects of the transmission dynamics of tuberculosis could be missed or misinterpreted. Furthermore, a valid differentiation between disease relapse and exogenous reinfection might be impossible using standard genotyping tools if the overall diversity of circulating clones is limited. These findings have important implications for clinical trials of new anti-tuberculosis drugs.

The type IIs restriction endonuclease BspMI is a tetramer that acts concertedly at two copies of an asymmetric DNA sequence.

Type IIs endonucleases recognize asymmetric DNA sequences and cleave both strands at fixed positions downstream of the sequence. Many type IIs enzymes, including BspMI, cleave substrates with two sites more rapidly than those with one site. They usually act sequentially on DNA with two sites, but BspMI converted such a substrate directly to the final products cut at both sites. The BspMI endonuclease was found to be a tetramer, in contrast to the monomeric structures for many type IIs enzymes. No change in subunit association occurred during the BspMI reaction. Plasmids with two BspMI sites were cleaved in cis, in reactions spanning sites in the same DNA, even when the sites were separated by just 38 bp. Plasmids with one BspMI site were cleaved in trans, with the enzyme bridging sites in separate DNA molecules: these slow reactions could be accelerated by adding a second DNA with the recognition sequence. Thus, whereas many type IIs enzymes dimerize before cleaving DNA, a process facilitated by two recognition sites in cis, the BspMI tetramer binds two copies of its recognition sequence before cleaving the DNA in both strands at both sites.

Many type IIs restriction endonucleases interact with two recognition sites before cleaving DNA.

Type IIs restriction endonucleases recognize asymmetric DNA sequences and cleave both DNA strands at fixed positions, typically several base pairs away from the recognition site. These enzymes are generally monomers that transiently associate to form dimers to cleave both strands. Their reactions could involve bridging interactions between two copies of their recognition sequence. To examine this possibility, several type IIs enzymes were tested against substrates with either one or two target sites. Some of the enzymes cleaved the DNA with two target sites at the same rate as that with one site, but most cut their two-site substrate more rapidly than the one-site DNA. In some cases, the two sites were cut sequentially, at rates that were equal to each other but that exceeded the rate on the one-site DNA. In another case, the DNA with two sites was cleaved rapidly at one site, but the residual site was cleaved at a much slower rate. In a further example, the two sites were cleaved concertedly to give directly the final products cut at both sites. Many type IIs enzymes thus interact with two copies of their recognition sequence before cleaving DNA, although via several different mechanisms.