Genomic Reads Research Articles

Prevalence and genetic variation of porcine circovirus type 2 in Hebei, China from 2004 to 2014

Porcine circovirus type 2 (PCV2), the primary causative agent of porcine circovirus-associated disease (PCVAD), causes severe economic losses to the pig industry in China since 2002. To investigate the molecular epidemic characteristics and genetic evolution of PCV2, 12 PCV2 isolates obtained from different pig farms with various clinical symptoms of PCVAD in Hebei, China from 2004 to 2014 were sequenced and analyzed. The phylogenetic analysis showed that the 12 isolates were divided into two distinct genotypes, PCV2b (7/12) and PCV2d (5/12), based on the sequences of either viral complete genome or open reading frame 2 (ORF2). Of the 7 PCV2b strains, 5 were isolated from 2004 to 2008 while all PCV2d were isolated from 2009 to 2014. This exhibited that PCV2b isolates were the most common before 2009 and then PCV2d isolates became predominant and widely distributed in pig farms. Sequence comparisons among total isolates indicated that the nucleotide identity ranged from 95.5% to 100% for complete genome and 93.1%–100% for ORF2. Compared with seven PCV2b isolates, there were thirteen amino-acid substitutions in the ORF2 region and one additional amino-acid K at this region terminal for five PCV2d isolates. The results suggest that a higher genetic variation and a distinct genotype shift occurred among the PCV2 isolates collected from 2004 to 2014 in Hebei.

The complete mitochondrial genome of the Anabas testudineus (Perciformes, Anabantidae) and its comparison with other related fish species

In the present study, the complete mitochondrial genome sequence of Anabas testudineusis reported using PGM sequencer (Ion Torrent, Life Technologies, La Jolla, CA). The complete mitogenome of climbing perch, A. testudineusis obtained by the de novo sequences assembly of genomic reads using the Torrent Mapping Alignment Program (TMAP), which is 16 603 bp in length. The mitogenome of A. testudineus composed of 13 protein- coding genes, two rRNA, and 22 tRNAs. Here, 20 tRNAs genes showed typical clover leaf model, and D-Loop as the control region along with gene order and organization, being closely similar to Osphronemidae and most of other Perciformes fish mitogenomes of NCBI databases. The mitogenome in the present study has 99% similarity to the complete mitogenome sequence of earlier reported A. testudineus. The phylogenetic analysis of Anabantidae depicted that their mitogenomes are closely related to each other. The complete mitogenome sequence of A. testudineus would be helpful in understanding the population genetics, phylogenetics, and evolution of Anabantidae.

Complete mitochondrial genome sequence of Indian medium carp, Labeo gonius (Hamilton, 1822) and its comparison with other related carp species

In the present study, the complete mitochondrial genome sequence of Labeo gonius is reported using PGM sequencer (Ion Torrent). The complete mitogenome of L. gonius is obtained by the de novo sequences assembly of genomic reads using the Torrent Mapping Alignment Program (TMAP) which is 16 614 bp in length. The mitogenome of L. gonius comprised of 13 protein-coding genes, 22 tRNAs, 2 rRNA genes, and D-loop as control region along with gene order and organization, being similar to most of other fish mitogenomes of NCBI databases. The mitogenome in the present study has 99% similarity to the complete mitogenome sequence of Labeo fimbriatus, as reported earlier. The phylogenetic analysis of Cypriniformes depicted that their mitogenomes are closely related to each other. The complete mitogenome sequence of L. gonius would be helpful in understanding the population genetics, phylogenetics, and evolution of Indian Carps.

Mitogenome assembly from genomic multiplex libraries: comparison of strategies and novel mitogenomes for five species of frogs.

Next-generation sequencing continues to revolutionize biodiversity studies by generating unprecedented amounts of DNA sequence data for comparative genomic analysis. However, these data are produced as millions or billions of short reads of variable quality that cannot be directly applied in comparative analyses, creating a demand for methods to facilitate assembly. We optimized an in silico strategy to efficiently reconstruct high-quality mitochondrial genomes directly from genomic reads. We tested this strategy using sequences from five species of frogs: Hylodes meridionalis (Hylodidae), Hyloxalus yasuni (Dendrobatidae), Pristimantis fenestratus (Craugastoridae), and Melanophryniscus simplex and Rhinella sp. (Bufonidae). These are the first mitogenomes published for these species, the genera Hylodes, Hyloxalus, Pristimantis, Melanophryniscus and Rhinella, and the families Craugastoridae and Hylodidae. Sequences were generated using only half of one lane of a standard Illumina HiqSeq 2000 flow cell, resulting in fewer than eight million reads. We analysed the reads of Hylodes meridionalis using three different assembly strategies: (1) reference-based (using bowtie2); (2) de novo (using abyss, soapdenovo2 and velvet); and (3) baiting and iterative mapping (using mira and mitobim). Mitogenomes were assembled exclusively with strategy 3, which we employed to assemble the remaining mitogenomes. Annotations were performed with mitos and confirmed by comparison with published amphibian mitochondria. In most cases, we recovered all 13 coding genes, 22 tRNAs, and two ribosomal subunit genes, with minor gene rearrangements. Our results show that few raw reads can be sufficient to generate high-quality scaffolds, making any Illumina machine run using genomic multiplex libraries a potential source of data for organelle assemblies as by-catch.

Selecting soluble/foldable protein domains through single-gene or genomic ORF filtering: structure of the head domain of Burkholderia pseudomallei antigen BPSL2063.

The 1.8 Å resolution crystal structure of a conserved domain of the potential Burkholderia pseudomallei antigen and trimeric autotransporter BPSL2063 is presented as a structural vaccinology target for melioidosis vaccine development. Since BPSL2063 (1090 amino acids) hosts only one conserved domain, and the expression/purification of the full-length protein proved to be problematic, a domain-filtering library was generated using β-lactamase as a reporter gene to select further BPSL2063 domains. As a result, two domains (D1 and D2) were identified and produced in soluble form in Escherichia coli. Furthermore, as a general tool, a genomic open reading frame-filtering library from the B. pseudomallei genome was also constructed to facilitate the selection of domain boundaries from the entire ORFeome. Such an approach allowed the selection of three potential protein antigens that were also produced in soluble form. The results imply the further development of ORF-filtering methods as a tool in protein-based research to improve the selection and production of soluble proteins or domains for downstream applications such as X-ray crystallography.

The rumen microbial metagenome associated with high methane production in cattle.

BackgroundMethane represents 16 % of total anthropogenic greenhouse gas emissions. It has been estimated that ruminant livestock produce ca. 29 % of this methane. As individual animals produce consistently different quantities of methane, understanding the basis for these differences may lead to new opportunities for mitigating ruminal methane emissions. Metagenomics is a powerful new tool for understanding the composition and function of complex microbial communities. Here we have applied metagenomics to the rumen microbial community to identify differences in the microbiota and metagenome that lead to high- and low-methane-emitting cattle phenotypes.MethodsFour pairs of beef cattle were selected for extreme high and low methane emissions from 72 animals, matched for breed (Aberdeen-Angus or Limousin cross) and diet (high or medium concentrate). Community analysis was carried out by qPCR of 16S and 18S rRNA genes and by alignment of Illumina HiSeq reads to the GREENGENES database. Total genomic reads were aligned to the KEGG genes databasefor functional analysis.ResultsDeep sequencing produced on average 11.3 Gb per sample. 16S rRNA gene abundances indicated that archaea, predominantly Methanobrevibacter, were 2.5× more numerous (P = 0.026) in high emitters, whereas among bacteria Proteobacteria, predominantly Succinivibrionaceae, were 4-fold less abundant (2.7 vs. 11.2 %; P = 0.002). KEGG analysis revealed that archaeal genes leading directly or indirectly to methane production were 2.7-fold more abundant in high emitters. Genes less abundant in high emitters included acetate kinase, electron transport complex proteins RnfC and RnfD and glucose-6-phosphate isomerase. Sequence data were assembled de novo and over 1.5 million proteins were annotated on the subsequent metagenome scaffolds. Less than half of the predicted genes matched matched a domain within Pfam. Amongst 2774 identified proteins of the 20 KEGG orthologues that correlated with methane emissions, only 16 showed 100 % identity with a publicly available protein sequence.ConclusionsThe abundance of archaeal genes in ruminal digesta correlated strongly with differing methane emissions from individual animals, a finding useful for genetic screening purposes. Lower emissions were accompanied by higher Succinovibrionaceae abundance and changes in acetate and hydrogen production leading to less methanogenesis, as similarly postulated for Australian macropods. Large numbers of predicted protein sequences differed between high- and low-methane-emitting cattle. Ninety-nine percent were unknown, indicating a fertile area for future exploitation.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2032-0) contains supplementary material, which is available to authorized users.

Draft genome of a commonly misdiagnosed multidrug resistant pathogen Candida auris.

BackgroundCandida auris is a multidrug resistant, emerging agent of fungemia in humans. Its actual global distribution remains obscure as the current commercial methods of clinical diagnosis misidentify it as C. haemulonii. Here we report the first draft genome of C. auris to explore the genomic basis of virulence and unique differences that could be employed for differential diagnosis.ResultsMore than 99.5 % of the C. auris genomic reads did not align to the current whole (or draft) genome sequences of Candida albicans, Candida lusitaniae, Candida glabrata and Saccharomyces cerevisiae; thereby indicating its divergence from the active Candida clade. The genome spans around 12.49 Mb with 8527 predicted genes. Functional annotation revealed that among the sequenced Candida species, it is closest to the hemiascomycete species Clavispora lusitaniae. Comparison with the well-studied species Candida albicans showed that it shares significant virulence attributes with other pathogenic Candida species such as oligopeptide transporters, mannosyl transfersases, secreted proteases and genes involved in biofilm formation. We also identified a plethora of transporters belonging to the ABC and major facilitator superfamily along with known MDR transcription factors which explained its high tolerance to antifungal drugs.ConclusionsOur study emphasizes an urgent need for accurate fungal screening methods such as PCR and electrophoretic karyotyping to ensure proper management of fungemia. Our work highlights the potential genetic mechanisms involved in virulence and pathogenicity of an important emerging human pathogen namely C. auris. Owing to its diversity at the genomic scale; we expect the genome sequence to be a useful resource to map species specific differences that will help develop accurate diagnostic markers and better drug targets.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1863-z) contains supplementary material, which is available to authorized users.

Cis-Regulatory Changes Associated with a Recent Mating System Shift and Floral Adaptation in Capsella.

The selfing syndrome constitutes a suite of floral and reproductive trait changes that have evolved repeatedly across many evolutionary lineages in response to the shift to selfing. Convergent evolution of the selfing syndrome suggests that these changes are adaptive, yet our understanding of the detailed molecular genetic basis of the selfing syndrome remains limited. Here, we investigate the role of cis-regulatory changes during the recent evolution of the selfing syndrome in Capsella rubella, which split from the outcrosser Capsella grandiflora less than 200 ka. We assess allele-specific expression (ASE) in leaves and flower buds at a total of 18,452 genes in three interspecific F1 C. grandiflora x C. rubella hybrids. Using a hierarchical Bayesian approach that accounts for technical variation using genomic reads, we find evidence for extensive cis-regulatory changes. On average, 44% of the assayed genes show evidence of ASE; however, only 6% show strong allelic expression biases. Flower buds, but not leaves, show an enrichment of cis-regulatory changes in genomic regions responsible for floral and reproductive trait divergence between C. rubella and C. grandiflora. We further detected an excess of heterozygous transposable element (TE) insertions near genes with ASE, and TE insertions targeted by uniquely mapping 24-nt small RNAs were associated with reduced expression of nearby genes. Our results suggest that cis-regulatory changes have been important during the recent adaptive floral evolution in Capsella and that differences in TE dynamics between selfing and outcrossing species could be important for rapid regulatory divergence in association with mating system shifts.

TEtranscripts: a package for including transposable elements in differential expression analysis of RNA-seq datasets.

Most RNA-seq data analysis software packages are not designed to handle the complexities involved in properly apportioning short sequencing reads to highly repetitive regions of the genome. These regions are often occupied by transposable elements (TEs), which make up between 20 and 80% of eukaryotic genomes. They can contribute a substantial portion of transcriptomic and genomic sequence reads, but are typically ignored in most analyses. Here, we present a method and software package for including both gene- and TE-associated ambiguously mapped reads in differential expression analysis. Our method shows improved recovery of TE transcripts over other published expression analysis methods, in both synthetic data and qPCR/NanoString-validated published datasets. The source code, associated GTF files for TE annotation, and testing data are freely available at http://hammelllab.labsites.cshl.edu/software. mhammell@cshl.edu. Supplementary data are available at Bioinformatics online.

Genome-Wide Search Identifies 1.9 Mb from the Polar Bear Y Chromosome for Evolutionary Analyses.

The male-inherited Y chromosome is the major haploid fraction of the mammalian genome, rendering Y-linked sequences an indispensable resource for evolutionary research. However, despite recent large-scale genome sequencing approaches, only a handful of Y chromosome sequences have been characterized to date, mainly in model organisms. Using polar bear (Ursus maritimus) genomes, we compare two different in silico approaches to identify Y-linked sequences: 1) Similarity to known Y-linked genes and 2) difference in the average read depth of autosomal versus sex chromosomal scaffolds. Specifically, we mapped available genomic sequencing short reads from a male and a female polar bear against the reference genome and identify 112 Y-chromosomal scaffolds with a combined length of 1.9 Mb. We verified the in silico findings for the longer polar bear scaffolds by male-specific in vitro amplification, demonstrating the reliability of the average read depth approach. The obtained Y chromosome sequences contain protein-coding sequences, single nucleotide polymorphisms, microsatellites, and transposable elements that are useful for evolutionary studies. A high-resolution phylogeny of the polar bear patriline shows two highly divergent Y chromosome lineages, obtained from analysis of the identified Y scaffolds in 12 previously published male polar bear genomes. Moreover, we find evidence of gene conversion among ZFX and ZFY sequences in the giant panda lineage and in the ancestor of ursine and tremarctine bears. Thus, the identification of Y-linked scaffold sequences from unordered genome sequences yields valuable data to infer phylogenomic and population-genomic patterns in bears.

Targeted single molecule sequencing methodology for ovarian hyperstimulation syndrome.

BackgroundOne of the most significant issues surrounding next generation sequencing is the cost and the difficulty assembling short read lengths. Targeted capture enrichment of longer fragments using single molecule sequencing (SMS) is expected to improve both sequence assembly and base-call accuracy but, at present, there are very few examples of successful application of these technologic advances in translational research and clinical testing. We developed a targeted single molecule sequencing (T-SMS) panel for genes implicated in ovarian response to controlled ovarian hyperstimulation (COH) for infertility.ResultsTarget enrichment was carried out using droplet-base multiplex polymerase chain reaction (PCR) technology (RainDance®) designed to yield amplicons averaging 1 kb fragment size from candidate 44 loci (99.8% unique base-pair coverage). The total targeted sequence was 3.18 Mb per sample. SMS was carried out using single molecule, real-time DNA sequencing (SMRT® Pacific Biosciences®), average raw read length = 1178 nucleotides, 5% of the amplicons >6000 nucleotides). After filtering with circular consensus (CCS) reads, the mean read length was 3200 nucleotides (97% CCS accuracy). Primary data analyses, alignment and filtering utilized the Pacific Biosciences® SMRT portal. Secondary analysis was conducted using the Genome Analysis Toolkit for SNP discovery l and wANNOVAR for functional analysis of variants. Filtered functional variants 18 of 19 (94.7%) were further confirmed using conventional Sanger sequencing. CCS reads were able to accurately detect zygosity. Coverage within GC rich regions (i.e.VEGFR; 72% GC rich) was achieved by capturing long genomic DNA (gDNA) fragments and reading into regions that flank the capture regions. As proof of concept, a non-synonymous LHCGR variant captured in two severe OHSS cases, and verified by conventional sequencing.ConclusionsCombining emulsion PCR-generated 1 kb amplicons and SMRT DNA sequencing permitted greater depth of coverage for T-SMS and facilitated easier sequence assembly. To the best of our knowledge, this is the first report combining emulsion PCR and T-SMS for long reads using human DNA samples, and NGS panel designed for biomarker discovery in OHSS.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1451-2) contains supplementary material, which is available to authorized users.

HyLiTE: accurate and flexible analysis of gene expression in hybrid and allopolyploid species.

BackgroundForming a new species through the merger of two or more divergent parent species is increasingly seen as a key phenomenon in the evolution of many biological systems. However, little is known about how expression of parental gene copies (homeologs) responds following genome merger. High throughput RNA sequencing now makes this analysis technically feasible, but tools to determine homeolog expression are still in their infancy.ResultsHere we present HyLiTE – a single-step analysis to obtain tables of homeolog expression in a hybrid or allopolyploid and its parent species directly from raw mRNA sequence files. By implementing on-the-fly detection of diagnostic parental polymorphisms, HyLiTE can perform SNP calling and read classification simultaneously, thus allowing HyLiTE to be run as parallelized code. HyLiTE accommodates any number of parent species, multiple data sources (including genomic DNA reads to improve SNP detection), and implements a statistical framework optimized for genes with low to moderate expression.ConclusionsHyLiTE is a flexible and easy-to-use program designed for bench biologists to explore patterns of gene expression following genome merger. HyLiTE offers practical advantages over manual methods and existing programs, has been designed to accommodate a wide range of genome merger systems, can identify SNPs that arose following genome merger, and offers accurate performance on non-model organisms.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-014-0433-8) contains supplementary material, which is available to authorized users.

Increased L-Ornithine Production in Corynebacterium glutamicum by Overexpression of a Gene Encoding a Putative Aminotransferase

Overexpression of the NCgl0462 open reading frame, encoding a class II aminotransferase, was studied in conjunction with other enzymes in <smlcap>L</smlcap>-ornithine biosynthesis in an <smlcap>L</smlcap>-ornithine-producing strain. Expression of the wild-type NCgl0462 open reading frame, which displayed aminotransferase activity, was amplified by placing it under the control of the glyceraldehyde 3-phosphate dehydrogenase gene promoter in the pEK0 plasmid and in the genome. <smlcap>L</smlcap>-Ornithine production in Corynebacteriumglutamicum SJC8260 harboring plasmid and the genomic NCgl0462 open reading frame was increased by 8.8 and 21.6%, respectively. In addition, the combined overexpression of the NCgl0462 open reading frame within the genome along with the mutated <smlcap>L</smlcap>-ornithine biosynthesis genes (argCJBD) placed in the pEK0 plasmid in C. glutamicum SJC8260 resulted in significant improvement in <smlcap>L</smlcap>-ornithine production (12.48 g/l for combined overexpression compared with 8.42 g/l for the control). These results suggest that overexpression of the aminotransferase-encoding NCgl0462 open reading frame plays an unequivocal role in the <smlcap>L</smlcap>-ornithine biosynthetic pathway, with overlapping substrate specificity in C. glutamicum.

Characterization of three active transposable elements recently inserted in three independent DFR-A alleles and one high-copy DNA transposon isolated from the Pink allele of the ANS gene in onion (Allium cepa L.).

Intact retrotransposon and DNA transposons inserted in a single gene were characterized in onions (Allium cepa) and their transcription and copy numbers were estimated in this study. While analyzing diverse onion germplasm, large insertions in the DFR-A gene encoding dihydroflavonol 4-reductase (DFR) involved in the anthocyanin biosynthesis pathway were found in two accessions. A 5,070-bp long terminal repeat (LTR) retrotransposon inserted in the active DFR-A (R4) allele was identified from one of the large insertions and designated AcCOPIA1. An intact ORF encoded typical domains of copia-like LTR retrotransposons. However, AcCOPIA1 contained atypical 'TG' and 'TA' dinucleotides at the ends of the LTRs. A 4,615-bp DNA transposon was identified in the other large insertion. This DNA transposon, designated AcCACTA1, contained an ORF coding for a transposase showing homology with the CACTA superfamily transposable elements (TEs). Another 5,073-bp DNA transposon was identified from the DFR-A (TRN) allele. This DNA transposon, designated AchAT1, belonged to the hAT superfamily with short 4-bp terminal inverted repeats (TIRs). Finally, a 6,258-bp non-autonomous DNA transposon, designated AcPINK, was identified in the ANS-p allele encoding anthocyanidin synthase, the next downstream enzyme to DFR in the anthocyanin biosynthesis pathway. AcPINK also possessed very short 3-bp TIRs. Active transcription of AcCOPIA1, AcCACTA1, and AchAT1 was observed through RNA-Seq analysis and RT-PCR. The copy numbers of AcPINK estimated by mapping the genomic DNA reads produced by NextSeq 500 were predominantly high compared with the other TEs. A series of evidence indicated that these TEs might have transposed in these onion genes very recently, providing a stepping stone for elucidation of enormously large-sized onion genome structure.

Fast and accurate mapping of Complete Genomics reads

Many recent advances in genomics and the expectations of personalized medicine are made possible thanks to power of high throughput sequencing (HTS) in sequencing large collections of human genomes. There are tens of different sequencing technologies currently available, and each HTS platform have different strengths and biases. This diversity both makes it possible to use different technologies to correct for shortcomings; but also requires to develop different algorithms for each platform due to the differences in data types and error models. The first problem to tackle in analyzing HTS data for resequencing applications is the read mapping stage, where many tools have been developed for the most popular HTS methods, but publicly available and open source aligners are still lacking for the Complete Genomics (CG) platform. Unfortunately, Burrows-Wheeler based methods are not practical for CG data due to the gapped nature of the reads generated by this method. Here we provide a sensitive read mapper (sirFAST) for the CG technology based on the seed-and-extend paradigm that can quickly map CG reads to a reference genome. We evaluate the performance and accuracy of sirFAST using both simulated and publicly available real data sets, showing high precision and recall rates.

Abstract 3560: Genome and RNA sequencing identifies fusion transcripts in bladder cancer

Abstract Urothelial cell carcinoma of the bladder (UCC) is a common malignancy with the highest prevalence rates found in Europe and United States where it rate as the fifth and fourth most common malignancy respectively. Bladder cancer is divided into two clinically divergent groups, non-muscle invasive (NMIBC) and muscle invasive bladder tumors (MIBC) accounting at diagnosis for 75% and 25% respectively and necessitates very different treatments. NMIBC have recurrence rate of 50-70% and often requires lifelong surveillance, however the prognosis is relatively good. In contrast MIBC often progresses and develop systemic disease. The aim of this study was to identify genomic rearrangements producing fusion transcripts in bladder cancer. Translocations in coding regions are known to be associated with malignancy. Especially in liquid tumors translocations are frequent; however the increased throughput in sequencing now identifies translocation associated with many solid tumors. In this study we used next generation sequencing to annotate genomic rearrangements in paired bladder cancer tumors (&amp;gt;90% tumor cells) presenting as Ta tumors and progressed to T1 or T2. Tumors were sequenced to 40-50x whole genome read depth (blood to &amp;gt;50x) and 65 million paired end reads of whole transcriptome RNAseq. Combining DNA and RNAseq, genomic rearrangements leading to expressed fusion transcript were identified (nFUSE, McPherson AW et. al. Genome Res. 2012 Jun 28). Approximately 25% of the predicted fusions could be validated using Sanger sequencing of cDNA. Importantly more than 80% of the fusions were only present in one of the paired tumors and primarily in the advanced tumor. Most of the fusions were out of frame and created truncated proteins. Two examples are truncation of metastasis suppressor 1 (MTSS1) and polymerase beta (POLB) both frequently mutated in cancer. As an interesting example of in frame fusions we found CTNNBL1 (Beta-Catenin-Like Protein 1) to be fused to FAM110A (family with sequence similarity 110, member A). Importantly the fusion places the promoter of CTNNBL1 in front of FAM110A and because the reading frame is preserved it is likely that both the start codon of CTNNBL1 and FAM110A can serve as translation start giving rise to nearly identical proteins. Furthermore, expression analysis showed significant increase of FAM110A expression in the affected tumor pair sample. Little is known about the FAM110 gene family, however they are shown to localizes to centrosomes and spindle poles and could potentially play an essential role in chromosome segregation and thereby cancer development. All sequencing was performed on the Illumina HiSeq platform. Citation Format: Iver K. Nordentoft, Karin Birkenkamp-Demtroder, Philippe Lamy, Thomas Reinert, Søren Vang, Jakob Hedegaard, Kasper Thorsen, Lars Dyrskjøt, Jakob Skou Pedersen, Torben Falck Ørntoft. Genome and RNA sequencing identifies fusion transcripts in bladder cancer. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 3560. doi:10.1158/1538-7445.AM2014-3560

Abstract 5321: BreaKmer: Detection of structural rearrangements in targeted next-generation sequencing data using kmers

Abstract Targeted next-generation sequencing (NGS) to capture genes or regions of interest has proven to be a cost-effective alternative to whole genome sequencing (WGS), particularly for cancer research and clinical cancer care. In this context, biologically or clinically relevant selected genes/regions are sequenced to several hundred-fold coverage. However, current algorithms and tools for detecting large structural variants (SVs), such as translocations, fail to achieve either high specificity or sensitivity, due to the fact that most available methods were designed for WGS data, and thus do not take advantage of the reduced size and higher coverage of targeted sequencing to improve SV calling. We developed a novel method, BreaKmer, to detect indels, rearrangements, and translocations from single sample targeted genomic reads. The algorithm extracts mis-mapped single or paired-end NGS reads. From these reads, hypothesized to contain SV breakpoints, contigs are built with a kmer strategy: reads are broken into k-length substrings or kmers, and those occurring in the reference are filtered. The remaining kmers represent sequences containing any sequence variant from the reference, ranging from single nucleotides to larger variants. Contigs are assembled from reads containing sample specific kmers. SVs are called based on alignment of the contigs to the reference sequence. With paired-end (PE) reads, discordantly mapped paired reads are extracted and coupled with SV calls that are made. To demonstrate BreaKmer, we analyzed NGS data from 166 samples enriched using 3 different capture panels (ranging from 305-504 genes). Our dataset contained 25 cancer specimens with known translocation events verified by orthogonal clinical methods and a negative control set of 141 ‘normal’ samples with no known SVs. The samples represented DNA extracted from FFPE, fresh frozen, blood and cell lines. Among 25 samples, 15 had additional replicate samples. Mean target coverage over all the samples was on average 150x. Specimens were barcoded at library preparation and pooled, followed by hybrid-capture targeting cancer genes and sequenced 2x100bp PE. All translocation events from the 25 test samples and their replicates (46) were detected by BreaKmer. An additional 19 translocations were detected among all the samples and their replicates, while no translocation calls were made for the 141 negative control samples. Our novel kmer strategy to detect SVs displayed high sensitivity and specificity. We reliably detected rearrangements of ALK, BCL2-IGH, BCR-ABL1 from lung adenocarcinoma, B-cell lymphoma, and chronic myeloid leukemia samples, respectively; indicating real clinical utility of this algorithm. In addition, our tool effectively detects other SV types - such as indels in FLT3 and KIT, among other genes. Our algorithm thus serves a pressing need for improved SV detection in targeted NGS data, particularly in precision cancer medicine. Citation Format: Ryan P. Abo, Elizabeth P. Garcia, Matthew Ducar, Ravali Adusumilli, Marc Breneiser, Vanessa Rojas-Rudilla, Lynette M. Sholl, Neal I. Lindeman, Matthew L. Meyerson, William C. Hahn, Paul Van Hummelen, Laura E. MacConaill. BreaKmer: Detection of structural rearrangements in targeted next-generation sequencing data using kmers. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 5321. doi:10.1158/1538-7445.AM2014-5321

Parallel Compression and Decompression of DNA Sequence Reads in FASTQ Format

Large volumes of short reads of genomic data are generated by high-throughput sequencing instruments. The FASTQ format is widely accepted as the input format of genomic reads and has presented challenges in data storage, management, and transfer. The performance of this type of serial algorithms such as G-SQZ and DSRC is limited by the single processor and the memory in a single computer. Utilizing data parallelism, the circular dual queues of buffers, memory mapping integrated with superblocks, pipeline parallelism with multi threads, and so on, we present the parallel compression and decompression methods for DNA sequence reads in FASTQ format based on the parallel computer architectures of the cluster and the SMP. Experimental results for the parallel DSRC algorithm clearly show the efficiency of using the powerful computing resources from multi computing nodes and multi cores of each node. The speedups vary from 46 to 62 for parallel compression and vary from 40 to 58 for parallel decompression by using 10 nodes of the cluster in Tianhe-1A super computer. Test results on the SMP machine are also pleasant. The methods could be applied to any serial compressing algorithms of DNA sequence reads in FASTQ format only if they have the traits of index and superblocks.

Software for pre-processing Illumina next-generation sequencing short read sequences.

BackgroundWhen compared to Sanger sequencing technology, next-generation sequencing (NGS) technologies are hindered by shorter sequence read length, higher base-call error rate, non-uniform coverage, and platform-specific sequencing artifacts. These characteristics lower the quality of their downstream analyses, e.g. de novo and reference-based assembly, by introducing sequencing artifacts and errors that may contribute to incorrect interpretation of data. Although many tools have been developed for quality control and pre-processing of NGS data, none of them provide flexible and comprehensive trimming options in conjunction with parallel processing to expedite pre-processing of large NGS datasets.MethodsWe developed ngsShoRT (next-generation sequencing Short Reads Trimmer), a flexible and comprehensive open-source software package written in Perl that provides a set of algorithms commonly used for pre-processing NGS short read sequences. We compared the features and performance of ngsShoRT with existing tools: CutAdapt, NGS QC Toolkit and Trimmomatic. We also compared the effects of using pre-processed short read sequences generated by different algorithms on de novo and reference-based assembly for three different genomes: Caenorhabditis elegans, Saccharomyces cerevisiae S288c, and Escherichia coli O157 H7.ResultsSeveral combinations of ngsShoRT algorithms were tested on publicly available Illumina GA II, HiSeq 2000, and MiSeq eukaryotic and bacteria genomic short read sequences with the focus on removing sequencing artifacts and low-quality reads and/or bases. Our results show that across three organisms and three sequencing platforms, trimming improved the mean quality scores of trimmed sequences. Using trimmed sequences for de novo and reference-based assembly improved assembly quality as well as assembler performance. In general, ngsShoRT outperformed comparable trimming tools in terms of trimming speed and improvement of de novo and reference-based assembly as measured by assembly contiguity and correctness.ConclusionsTrimming of short read sequences can improve the quality of de novo and reference-based assembly and assembler performance. The parallel processing capability of ngsShoRT reduces trimming time and improves the memory efficiency when dealing with large datasets. We recommend combining sequencing artifacts removal, and quality score based read filtering and base trimming as the most consistent method for improving sequence quality and downstream assemblies.ngsShoRT source code, user guide and tutorial are available at http://research.bioinformatics.udel.edu/genomics/ngsShoRT/. ngsShoRT can be incorporated as a pre-processing step in genome and transcriptome assembly projects.

Visualization of nucleotide substitutions in the (micro)transcriptome.

BackgroundRNA-related applications of the next-generation sequencing (NGS) technologies require context-specific interpretations: e.g., sequence mismatches may indicate sites of RNA editing, or uneven read coverage often points to mature form of microRNA. Existing visualization tools traditionally show RNA molecules in two dimensions, with their base pairing and the resulting secondary structure. However, it is not straightforward to combine a linear NGS data display with the 2-D RNA depictions.ResultsWe present a novel approach for interactive representation of nucleotide substitutions and modifications in the transcribed genome. With the focus on RNA secondary structure in the context of NGS data, it provides intuitive visualization of genomic environment, sequence reads, nucleotide polymorphisms and editing events integrated with the structural and functional elements of both coding and non-coding RNA molecules. Using our approach we present and discuss examples and general trends of polymorphisms and editing in the context of the secondary structure of microRNAs. As expected, most of the substitutions comprised A to G and C to T events, consistent with typical RNA editing patterns. However, we did not observe prevalence of editing in double-stranded regions of the microRNA stem-loop. We describe novel prominent editing event candidates, observed across several small RNA libraries of Drosophila melanogaster.ConclusionsIn contrast to the existing general tools for NGS data visualization, the power of our approach is not only in the display of read alignments and their counts, but the integration of RNA secondary structure, sequencing depth, and rates/patterns of editing or other modifications. It provides a comprehensive picture, important for large-scale studies and detailed analyses, helping to gain insight into the intricate relationships between different events in RNA biogenesis.