Genome sequencing and analysis of the Tasmanian devil and its transmissible cancer.

The Tasmanian devil (Sarcophilus harrisii), the largest marsupial carnivore, is endangered due to a transmissible facial cancer spread by direct transfer of living cancer cells through biting. Here we describe the sequencing, assembly, and annotation of the Tasmanian devil genome and whole-genome sequences for two geographically distant subclones of the cancer. Genomic analysis suggests that the cancer first arose from a female Tasmanian devil and that the clone has subsequently genetically diverged during its spread across Tasmania. The devil cancer genome contains more than 17,000 somatic base substitution mutations and bears the imprint of a distinct mutational process. Genotyping of somatic mutations in 104 geographically and temporally distributed Tasmanian devil tumors reveals the pattern of evolution and spread of this parasitic clonal lineage, with evidence of a selective sweep in one geographical area and persistence of parallel lineages in other populations.

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.

Improved genome sequencing using an engineered transposase.

BACKGROUND: Next-generation sequencing (NGS) has transformed genomic research by reducing turnaround time and cost. However, no major breakthrough has been made in the upstream library preparation methods until the transposase-based Nextera method was invented. Nextera combines DNA fragmentation and barcoding in a single tube reaction and therefore enables a very fast workflow to sequencing-ready DNA libraries within a couple of hours. When compared to the traditional ligation-based methods, transposed-based Nextera has a slight insertion bias. RESULTS: Here we present the discovery of a mutant transposase (Tn5-059) with a lowered GC insertion bias through protein engineering. We demonstrate Tn5-059 reduces AT dropout and increases uniformity of genome coverage in both bacterial genomes and human genome. We also observe higher library diversity generated by Tn5-059 when compared to Nextera v2 for human exomes, which leads to less sequencing and lower cost per genome. In addition, when used for human exomes, Tn5-059 delivers consistent library insert size over a range of input DNA, allowing up to a tenfold variance from the 50 ng input recommendation. CONCLUSIONS: Enhanced DNA input tolerance of Tn5-059 can translate to flexibility and robustness of workflow. DNA input tolerance together with superior uniformity of coverage and lower AT dropouts extend the applications of transposase based library preps. We discuss possible mechanisms of improvements in Tn5-059, and potential advantages of using the new mutant in varieties of applications including microbiome sequencing and chromatin profiling.

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

Clonality of HTLV-2 in natural infection.

Human T-lymphotropic virus type 1 (HTLV-1) and type 2 (HTLV-2) both cause lifelong persistent infections, but differ in their clinical outcomes. HTLV-1 infection causes a chronic or acute T-lymphocytic malignancy in up to 5% of infected individuals whereas HTLV-2 has not been unequivocally linked to a T-cell malignancy. Virus-driven clonal proliferation of infected cells both in vitro and in vivo has been demonstrated in HTLV-1 infection. However, T-cell clonality in HTLV-2 infection has not been rigorously characterized. In this study we used a high-throughput approach in conjunction with flow cytometric sorting to identify and quantify HTLV-2-infected T-cell clones in 28 individuals with natural infection. We show that while genome-wide integration site preferences in vivo were similar to those found in HTLV-1 infection, expansion of HTLV-2-infected clones did not demonstrate the same significant association with the genomic environment of the integrated provirus. The proviral load in HTLV-2 is almost confined to CD8+ T-cells and is composed of a small number of often highly expanded clones. The HTLV-2 load correlated significantly with the degree of dispersion of the clone frequency distribution, which was highly stable over ∼8 years. These results suggest that there are significant differences in the selection forces that control the clonal expansion of virus-infected cells in HTLV-1 and HTLV-2 infection. In addition, our data demonstrate that strong virus-driven proliferation per se does not predispose to malignant transformation in oncoretroviral infections.

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.

Estimating abundances of retroviral insertion sites from DNA fragment length data.

MOTIVATION: The relative abundance of retroviral insertions in a host genome is important in understanding the persistence and pathogenesis of both natural retroviral infections and retroviral gene therapy vectors. It could be estimated from a sample of cells if only the host genomic sites of retroviral insertions could be directly counted. When host genomic DNA is randomly broken via sonication and then amplified, amplicons of varying lengths are produced. The number of unique lengths of amplicons of an insertion site tends to increase according to its abundance, providing a basis for estimating relative abundance. However, as abundance increases amplicons of the same length arise by chance leading to a non-linear relation between the number of unique lengths and relative abundance. The difficulty in calibrating this relation is compounded by sample-specific variations in the relative frequencies of clones of each length. RESULTS: A likelihood function is proposed for the discrete lengths observed in each of a collection of insertion sites and is maximized with a hybrid expectation-maximization algorithm. Patient data illustrate the method and simulations show that relative abundance can be estimated with little bias, but that variation in highly abundant sites can be large. In replicated patient samples, variation exceeds what the model implies-requiring adjustment as in Efron (2004) or using jackknife standard errors. Consequently, it is advantageous to collect replicate samples to strengthen inferences about relative abundance.

The host genomic environment of the provirus determines the abundance of HTLV-1-infected T-cell clones.

Human T-lymphotropic virus type 1 (HTLV-1) persists by driving clonal proliferation of infected T lymphocytes. A high proviral load predisposes to HTLV-1-associated diseases. Yet the reasons for the variation within and between persons in the abundance of HTLV-1-infected clones remain unknown. We devised a high-throughput protocol to map the genomic location and quantify the abundance of > 91,000 unique insertion sites of the provirus from 61 HTLV-1(+) persons and > 2100 sites from in vitro infection. We show that a typical HTLV-1-infected host carries between 500 and 5000 unique insertion sites. We demonstrate that negative selection dominates during chronic infection, favoring establishment of proviruses integrated in transcriptionally silenced DNA: this selection is significantly stronger in asymptomatic carriers. We define a parameter, the oligoclonality index, to quantify clonality. The high proviral load characteristic of HTLV-1-associated inflammatory disease results from a larger number of unique insertion sites than in asymptomatic carriers and not, as previously thought, from a difference in clonality. The abundance of established HTLV-1 clones is determined by genomic features of the host DNA flanking the provirus. HTLV-1 clonal expansion in vivo is favored by orientation of the provirus in the same sense as the nearest host gene.

Haplotype phasing of whole human genomes using bead-based barcode partitioning in a single tube.

Haplotype-resolved genome sequencing promises to unlock a wealth of information in population and medical genetics. However, for the vast majority of genomes sequenced to date, haplotypes have not been determined because of cumbersome haplotyping workflows that require fractions of the genome to be sequenced in a large number of compartments. Here we demonstrate barcode partitioning of long DNA molecules in a single compartment using "on-bead" barcoded tagmentation. The key to the method that we call "contiguity preserving transposition" sequencing on beads (CPTv2-seq) is transposon-mediated transfer of homogenous populations of barcodes from beads to individual long DNA molecules that get fragmented at the same time (tagmentation). These are then processed to sequencing libraries wherein all sequencing reads originating from each long DNA molecule share a common barcode. Single-tube, bulk processing of long DNA molecules with ∼150,000 different barcoded bead types provides a barcode-linked read structure that reveals long-range molecular contiguity. This technology provides a simple, rapid, plate-scalable and automatable route to accurate, haplotype-resolved sequencing, and phasing of structural variants of the genome.

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.

Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms.

DNA sequencing continues to decrease in cost with the Illumina HiSeq2000 generating up to 600 Gb of paired-end 100 base reads in a ten-day run. Here we present a protocol for community amplicon sequencing on the HiSeq2000 and MiSeq Illumina platforms, and apply that protocol to sequence 24 microbial communities from host-associated and free-living environments. A critical question as more sequencing platforms become available is whether biological conclusions derived on one platform are consistent with what would be derived on a different platform. We show that the protocol developed for these instruments successfully recaptures known biological results, and additionally that biological conclusions are consistent across sequencing platforms (the HiSeq2000 versus the MiSeq) and across the sequenced regions of amplicons.

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.

Predictors of institutionalization in patients with dementia in Korea.

BACKGROUND: Many studies have sought to determine the predictors of institutionalization of patients with dementia. Such studies, performed in developed western societies, have come to various conclusions which may not be supported in an East Asian culture such as that found in Korea. OBJECTIVES: This study aimed to determine the factors that predict institutionalization of patients in Korea diagnosed with dementia. METHODS: Seventy-nine cases (37 institutionalized, 42 community-dwelling) in the Kwangju area were evaluated for patient characteristics, severity of dementia symptoms, caregiver characteristics, burden and distress. Logistic regression was performed to determine predictors of actual institutionalization. RESULTS: Six predictors of institutionalization were identified. Of these, three were patient-related factors: higher score on the Clinical Dementia Rating, higher score on the Brief Psychiatric Rating Scale, and shorter duration of dementia. The other three were caregiver-related factors: younger age, higher education (formal schooling), and higher cost of home care. CONCLUSIONS: As seen in previous western studies, institutionalization of dementia sufferers was influenced by both patient and caregiver factors. But, the specific predictors and their relative influences might be explained best by the particular social, cultural and economic situation in Korea. This study was the first of its kind in Korea and, as such, could serve as a reference for future intra-cultural and cross-cultural comparisons.