Short-read Next-generation Sequencing Research Articles

Next-Generation Anchor Based Phylogeny (NexABP): Constructing phylogeny from Next-generation sequencing data

Whole genome sequences are ideally suited for deriving evolutionary relationship among organisms. With the availability of Next Generation sequencing (NGS) datasets in an unprecedented scale, it will be highly desirable if phylogenetic analysis can be carried out using short read NGS data. We described here an anchor based approach NexABP for phylogenetic construction of closely related strains/isolates from NGS data. This approach can be used even in the absence of a fully assembled reference genome and works by reducing the complexity of the datasets without compromising results. NexABP was used for constructing phylogeny of different strains of some of the common pathogens, such as Mycobacterium tuberculosis, Vibrio cholera and Escherichia coli. In addition to classification into distinct lineages, NexABP could resolve inner branches and also allow statistical testing using bootstrap analysis. We believe that there are some clear advantages of using NexABP based phylogenetic analysis as compared to other methods.

Frogs and estivation: transcriptional insights into metabolism and cell survival in a natural model of extended muscle disuse

Green-striped burrowing frogs (Cyclorana alboguttata) survive in arid environments by burrowing underground and entering into a deep, prolonged metabolic depression known as estivation. Throughout estivation, C. alboguttata is immobilized within a cast-like cocoon of shed skin and ceases feeding and moving. Remarkably, these frogs exhibit very little muscle atrophy despite extended disuse and fasting. Little is known about the transcriptional regulation of estivation or associated mechanisms that may minimize degradative pathways of atrophy. To investigate transcriptional pathways associated with metabolic depression and maintenance of muscle function in estivating burrowing frogs, we assembled a skeletal muscle transcriptome using next-generation short read sequencing and compared gene expression patterns between active and 4 mo estivating C. alboguttata. This identified a complex suite of gene expression changes that occur in muscle during estivation and provides evidence that estivation in burrowing frogs involves transcriptional regulation of genes associated with cytoskeletal remodeling, avoidance of oxidative stress, energy metabolism, the cell stress response, and apoptotic signaling. In particular, the expression levels of genes encoding cell cycle and prosurvival proteins, such as serine/threonine-protein kinase Chk1, cell division protein kinase 2, survivin, and vesicular overexpressed in cancer prosurvival protein 1, were upregulated during estivation. These data suggest that estivating C. alboguttata are able to regulate the expression of genes in several major cellular pathways critical to the survival and viability of cells, thus preserving muscle function while avoiding the deleterious consequences often seen in laboratory models of muscle disuse.

Viral population analysis and minority-variant detection using short read next-generation sequencing

RNA viruses within infected individuals exist as a population of evolutionary-related variants. Owing to evolutionary change affecting the constitution of this population, the frequency and/or occurrence of individual viral variants can show marked or subtle fluctuations. Since the development of massively parallel sequencing platforms, such viral populations can now be investigated to unprecedented resolution. A critical problem with such analyses is the presence of sequencing-related errors that obscure the identification of true biological variants present at low frequency. Here, we report the development and assessment of the Quality Assessment of Short Read (QUASR) Pipeline (http://sourceforge.net/projects/quasr) specific for virus genome short read analysis that minimizes sequencing errors from multiple deep-sequencing platforms, and enables post-mapping analysis of the minority variants within the viral population. QUASR significantly reduces the error-related noise in deep-sequencing datasets, resulting in increased mapping accuracy and reduction of erroneous mutations. Using QUASR, we have determined influenza virus genome dynamics in sequential samples from an in vitro evolution of 2009 pandemic H1N1 (A/H1N1/09) influenza from samples sequenced on both the Roche 454 GSFLX and Illumina GAIIx platforms. Importantly, concordance between the 454 and Illumina sequencing allowed unambiguous minority-variant detection and accurate determination of virus population turnover in vitro.

Telomeric gene deletion and intrachromosomal amplification in antimony‐resistant Leishmania

Antimonials are still the mainstay of treatment against leishmaniasis but drug resistance is increasing. We carried out short read next-generation sequencing (NGS) and comparative genomic hybridization (CGH) of three independent Leishmania major antimony-resistant mutants. Copy number variations were consistently detected with both NGS and CGH. A major attribute of antimony resistance was a novel terminal deletion of variable length (67 kb to 204 kb) of the polyploid chromosome 31 in the three mutants. Terminal deletions in two mutants occurred at the level of inverted repeated sequences. The AQP1 gene coding for an aquaglyceroporin was part of the deleted region and its transfection into resistant mutants reverted resistance to SbIII. We also highlighted an intrachromosomal amplification of a subtelomeric locus on chromosome 34 in one mutant. This region encoded for ascorbate-dependent peroxidase (APX) and glucose-6-phosphate dehydrogenase (G6PDH). Overexpression of these genes in revertant backgrounds demonstrated resistance to SbIII and protection from reactive oxygen species (ROS). Generation of a G6PDH null mutant in one revertant exhibited SbIII sensitivity and a decreased protection of ROS. Our genomic analyses and functional validation highlighted novel genomic rearrangements, functionally important resistant loci and the implication of new genes in antimony resistance in Leishmania.

Rapid Genome Mapping in Nanochannel Arrays for Highly Complete and Accurate De Novo Sequence Assembly of the Complex Aegilops tauschii Genome

Next-generation sequencing (NGS) technologies have enabled high-throughput and low-cost generation of sequence data; however, de novo genome assembly remains a great challenge, particularly for large genomes. NGS short reads are often insufficient to create large contigs that span repeat sequences and to facilitate unambiguous assembly. Plant genomes are notorious for containing high quantities of repetitive elements, which combined with huge genome sizes, makes accurate assembly of these large and complex genomes intractable thus far. Using two-color genome mapping of tiling bacterial artificial chromosomes (BAC) clones on nanochannel arrays, we completed high-confidence assembly of a 2.1-Mb, highly repetitive region in the large and complex genome of Aegilops tauschii, the D-genome donor of hexaploid wheat (Triticum aestivum). Genome mapping is based on direct visualization of sequence motifs on single DNA molecules hundreds of kilobases in length. With the genome map as a scaffold, we anchored unplaced sequence contigs, validated the initial draft assembly, and resolved instances of misassembly, some involving contigs <2 kb long, to dramatically improve the assembly from 75% to 95% complete.

MaSC: mappability-sensitive cross-correlation for estimating mean fragment length of single-end short-read sequencing data

Motivation: Reliable estimation of the mean fragment length for next-generation short-read sequencing data is an important step in next-generation sequencing analysis pipelines, most notably because of its impact on the accuracy of the enriched regions identified by peak-calling algorithms. Although many peak-calling algorithms include a fragment-length estimation subroutine, the problem has not been adequately solved, as demonstrated by the variability of the estimates returned by different algorithms.Results: In this article, we investigate the use of strand cross-correlation to estimate mean fragment length of single-end data and show that traditional estimation approaches have mixed reliability. We observe that the mappability of different parts of the genome can introduce an artificial bias into cross-correlation computations, resulting in incorrect fragment-length estimates. We propose a new approach, called mappability-sensitive cross-correlation (MaSC), which removes this bias and allows for accurate and reliable fragment-length estimation. We analyze the computational complexity of this approach, and evaluate its performance on a test suite of NGS datasets, demonstrating its superiority to traditional cross-correlation analysis.Availability: An open-source Perl implementation of our approach is available at http://www.perkinslab.ca/Software.html.Contact: tperkins@ohri.caSupplementary information: Supplementary data are available at Bioinformatics online.

Characterisation of the wheat (Triticum aestivum L.) transcriptome by de novo assembly for the discovery of phosphate starvation-responsive genes: gene expression in Pi-stressed wheat.

BackgroundPhosphorus (P) is an essential macronutrient for plant growth and development. To modulate their P homeostasis, plants must balance P uptake, mobilisation, and partitioning to various organs. Despite the worldwide importance of wheat as a cultivated food crop, molecular mechanisms associated with phosphate (Pi) starvation in wheat remain unclear. To elucidate these mechanisms, we used RNA-Seq methods to generate transcriptome profiles of the wheat variety ‘Chinese Spring’ responding to 10 days of Pi starvation.ResultsWe carried out de novo assembly on 73.8 million high-quality reads generated from RNA-Seq libraries. We then constructed a transcript dataset containing 29,617 non-redundant wheat transcripts, comprising 15,047 contigs and 14,570 non-redundant full-length cDNAs from the TriFLDB database. When compared with barley full-length cDNAs, 10,656 of the 15,047 contigs were unalignable, suggesting that many might be distinct from barley transcripts. The average expression level of the contigs was lower than that of the known cDNAs, implying that these contigs included transcripts that were rarely represented in the full-length cDNA library. Within the non-redundant transcript set, we identified 892–2,833 responsive transcripts in roots and shoots, corresponding on average to 23.4% of the contigs not covered by cDNAs in TriFLDB under Pi starvation. The relative expression level of the wheat IPS1 (Induced by Phosphate Starvation 1) homologue, TaIPS1, was 341-fold higher in roots and 13-fold higher in shoots; this finding was further confirmed by qRT-PCR analysis. A comparative analysis of the wheat- and rice-responsive transcripts for orthologous genes under Pi-starvation revealed commonly upregulated transcripts, most of which appeared to be involved in a general response to Pi starvation, namely, an IPS1-mediated signalling cascade and its downstream functions such as Pi remobilisation, Pi uptake, and changes in Pi metabolism.ConclusionsOur transcriptome profiles demonstrated the impact of Pi starvation on global gene expression in wheat. This study revealed that enhancement of the Pi-mediated signalling cascade using IPS1 is a potent adaptation mechanism to Pi starvation that is conserved in both wheat and rice and validated the effectiveness of using short-read next-generation sequencing data for wheat transcriptome analysis in the absence of reference genome information.

Use of bifidobacterial specific terminal restriction fragment length polymorphisms to complement next generation sequence profiling of infant gut communities

Bifidobacteria are intestinal anaerobes often associated with gut health. Specific bifidobacterial species are particularly common in the gastrointestinal tract of breast-fed infants. Current short read next-generation sequencing approaches to profile fecal microbial ecologies do not discriminate bifidobacteria to the species level. Here we describe a low-cost terminal restriction fragment length polymorphism (TRFLP) procedure to distinguish between the common infant-associated bifidobacterial species. An empirical database of TRF sizes was created from both common reference strains and well-identified isolates from infant feces. Species-specific quantitative PCR validated bifidobacterial-specific TRFLP profiles from infant feces. These results indicate that bifidobacterial-specific TRFLP is a useful method to monitor intestinal bifidobacterial populations from infant fecal samples. When used alongside next generation sequencing methods that detect broader population levels at lower resolution, this high-throughput, low-cost tool can help clarify the role of bifidobacteria in health and disease.

Paired-End Sequencing of Long-Range DNA Fragments for De Novo Assembly of Large, Complex Mammalian Genomes by Direct Intra-Molecule Ligation

BackgroundThe relatively short read lengths from next generation sequencing (NGS) technologies still pose a challenge for de novo assembly of complex mammal genomes. One important solution is to use paired-end (PE) sequence information experimentally obtained from long-range DNA fragments (>1 kb). Here, we characterize and extend a long-range PE library construction method based on direct intra-molecule ligation (or molecular linker-free circularization) for NGS.ResultsWe found that the method performs stably for PE sequencing of 2- to 5- kb DNA fragments, and can be extended to 10–20 kb (and even in extremes, up to ∼35 kb). We also characterized the impact of low quality input DNA on the method, and develop a whole-genome amplification (WGA) based protocol using limited input DNA (<1 µg). Using this PE dataset, we accurately assembled the YanHuang (YH) genome, the first sequenced Asian genome, into a scaffold N50 size of >2 Mb, which is over100-times greater than the initial size produced with only small insert PE reads(17 kb). In addition, we mapped two 7- to 8- kb insertions in the YH genome using the larger insert sizes of the long-range PE data.ConclusionsIn conclusion, we demonstrate here the effectiveness of this long-range PE sequencing method and its use for the de novo assembly of a large, complex genome using NGS short reads.

Next-generation sequencing and large genome assemblies.

The next-generation sequencing (NGS) revolution has drastically reduced time and cost requirements for sequencing of large genomes, and also qualitatively changed the problem of assembly. This article reviews the state of the art in de novo genome assembly, paying particular attention to mammalian-sized genomes. The strengths and weaknesses of the main sequencing platforms are highlighted, leading to a discussion of assembly and the new challenges associated with NGS data. Current approaches to assembly are outlined and the various software packages available are introduced and compared. The question of whether quality assemblies can be produced using short-read NGS data alone, or whether it must be combined with more expensive sequencing techniques, is considered. Prospects for future assemblers and tests of assembly performance are also discussed.

Scaffolder - software for manual genome scaffolding

BackgroundThe assembly of next-generation short-read sequencing data can result in a fragmented non-contiguous set of genomic sequences. Therefore a common step in a genome project is to join neighbouring sequence regions together and fill gaps. This scaffolding step is non-trivial and requires manually editing large blocks of nucleotide sequence. Joining these sequences together also hides the source of each region in the final genome sequence. Taken together these considerations may make reproducing or editing an existing genome scaffold difficult.MethodsThe software outlined here, “Scaffolder,” is implemented in the Ruby programming language and can be installed via the RubyGems software management system. Genome scaffolds are defined using YAML - a data format which is both human and machine-readable. Command line binaries and extensive documentation are available.ResultsThis software allows a genome build to be defined in terms of the constituent sequences using a relatively simple syntax. This syntax further allows unknown regions to be specified and additional sequence to be used to fill known gaps in the scaffold. Defining the genome construction in a file makes the scaffolding process reproducible and easier to edit compared with large FASTA nucleotide sequences.ConclusionsScaffolder is easy-to-use genome scaffolding software which promotes reproducibility and continuous development in a genome project. Scaffolder can be found at http://next.gs.

Improved variant discovery through local re-alignment of short-read next-generation sequencing data using SRMA

A primary component of next-generation sequencing analysis is to align short reads to a reference genome, with each read aligned independently. However, reads that observe the same non-reference DNA sequence are highly correlated and can be used to better model the true variation in the target genome. A novel short-read micro re-aligner, SRMA, that leverages this correlation to better resolve a consensus of the underlying DNA sequence of the targeted genome is described here.