Subset Of Reads Research Articles

SVNN: an efficient PacBio-specific pipeline for structural variations calling using neural networks

BackgroundOnce aligned, long-reads can be a useful source of information to identify the type and position of structural variations. However, due to the high sequencing error of long reads, long-read structural variation detection methods are far from precise in low-coverage cases. To be accurate, they need to use high-coverage data, which in turn, results in an extremely time-consuming pipeline, especially in the alignment phase. Therefore, it is of utmost importance to have a structural variation calling pipeline which is both fast and precise for low-coverage data.ResultsIn this paper, we present SVNN, a fast yet accurate, structural variation calling pipeline for PacBio long-reads that takes raw reads as the input and detects structural variants of size larger than 50 bp. Our pipeline utilizes state-of-the-art long-read aligners, namely NGMLR and Minimap2, and structural variation callers, videlicet Sniffle and SVIM. We found that by using a neural network, we can extract features from Minimap2 output to detect a subset of reads that provide useful information for structural variation detection. By only mapping this subset with NGMLR, which is far slower than Minimap2 but better serves downstream structural variation detection, we can increase the sensitivity in an efficient way. As a result of using multiple tools intelligently, SVNN achieves up to 20 percentage points of sensitivity improvement in comparison with state-of-the-art methods and is three times faster than a naive combination of state-of-the-art tools to achieve almost the same accuracy.ConclusionSince prohibitive costs of using high-coverage data have impeded long-read applications, with SVNN, we provide the users with a much faster structural variation detection platform for PacBio reads with high precision and sensitivity in low-coverage scenarios.

VARUS: sampling complementary RNA reads from the sequence read archive

BackgroundVast amounts of next generation sequencing RNA data has been deposited in archives, accompanying very diverse original studies. The data is readily available also for other purposes such as genome annotation or transcriptome assembly. However, selecting a subset of available experiments, sequencing runs and reads for this purpose is a nontrivial task and complicated by the inhomogeneity of the data.ResultsThis article presents the software VARUS that selects, downloads and aligns reads from NCBI’s Sequence Read Archive, given only the species’ binomial name and genome. VARUS automatically chooses runs from among all archived runs to randomly select subsets of reads. The objective of its online algorithm is to cover a large number of transcripts adequately when network bandwidth and computing resources are limited. For most tested species VARUS achieved both a higher sensitivity and specificity with a lower number of downloaded reads than when runs were manually selected. At the example of twelve eukaryotic genomes, we show that RNA-Seq that was sampled with VARUS is well-suited for fully-automatic genome annotation with BRAKER.ConclusionsWith VARUS, genome annotation can be automatized to the extent that not even the selection and quality control of RNA-Seq has to be done manually. This introduces the possibility to have fully automatized genome annotation loops over potentially many species without incurring a loss of accuracy over a manually supervised annotation process.

Experimental approaches for ancient DNA extraction and sample preparation for next generation sequencing in ultra-clean conditions

In this study, we aimed to develop an approach for genetic material extraction and sample preparation of ancient DNA for next generation sequencing that minimizes intra-laboratory contamination. A special module system consisting of four sterile gloveboxes connected by antechambers, which maintain a high purity atmosphere and allow conducting all stages of sample preparation in an isolated clean environment, starting with the processing of archaeological material, was developed. Furthermore, a fast and cost-effective double-stranded library preparation protocol for massive parallel sequencing was developed using non-standard modified sequencing adaptors, which enabled obtaining a library of fragments from extremely low amounts of starting degraded material. Bioinformatics data processing showed that the distribution of substitution frequencies in a subset of reads, mapped to the human reference genome, completely coincides with the pattern of postmortem modifications expected when the double-stranded library preparation method is employed. This is a strong evidence for the endogenous origin of the sequenced fragments of the human genome. Thus, the results of ancient sample and negative control sample sequencing, indicated the absence of visible contamination during sample preparation procedure.

Norgal: extraction and de novo assembly of mitochondrial DNA from whole-genome sequencing data

BackgroundWhole-genome sequencing (WGS) projects provide short read nucleotide sequences from nuclear and possibly organelle DNA depending on the source of origin. Mitochondrial DNA is present in animals and fungi, while plants contain DNA from both mitochondria and chloroplasts. Current techniques for separating organelle reads from nuclear reads in WGS data require full reference or partial seed sequences for assembling.ResultsNorgal (de Novo ORGAneLle extractor) avoids this requirement by identifying a high frequency subset of k-mers that are predominantly of mitochondrial origin and performing a de novo assembly on a subset of reads that contains these k-mers. The method was applied to WGS data from a panda, brown algae seaweed, butterfly and filamentous fungus. We were able to extract full circular mitochondrial genomes and obtained sequence identities to the reference sequences in the range from 98.5 to 99.5%. We also assembled the chloroplasts of grape vines and cucumbers using Norgal together with seed-based de novo assemblers.ConclusionNorgal is a pipeline that can extract and assemble full or partial mitochondrial and chloroplast genomes from WGS short reads without prior knowledge. The program is available at: https://bitbucket.org/kosaidtu/norgal.

BioSpotlight / Citations

BioSpotlight / Citations

Genotyping-by-Sequencing on Pooled Samples and its Use in Measuring Segregation Bias during the Course of Androgenesis in Barley.

Estimation of allelic frequencies is often required in breeding but genotyping many individuals at many loci can be expensive. We have developed a genotyping-by-sequencing (GBS) approach for estimating allelic frequencies on pooled samples (Pool-GBS) and used it to examine segregation distortion in doubled haploid (DH) populations of barley ( L.). In the first phase, we genotyped each line individually and exploited these data to explore a strategy to call single nucleotide polymorphisms (SNPs) on pooled reads. We measured both the number of SNPs called and the variance of the estimated allelic frequencies at various depths of coverage on a subset of reads containing 5 to 25 million reads. We show that allelic frequencies could be cost-effectively and accurately estimated at a depth of 50 reads per SNP using 15 million reads. This Pool-GBS approach yielded 1984 SNPs whose allelic frequency estimates were highly reproducible (CV = 10.4%) and correlated ( = 0.9167) with the "true" frequency derived from analysis of individual lines. In a second phase, we used Pool-GBS to investigate segregation bias throughout androgenesis from microspores to a population of regenerated plants. No strong bias was detected among the microspores resulting from the meiotic divisions, whereas significant biases could be shown to arise during embryo formation and plant regeneration. In summary, this methodology provides an approach to estimate allelic frequencies more efficiently and on materials that are unsuitable for individual analysis. In addition, it allowed us to shed light on the process of androgenesis in barley.