Long-read Sequencing Methods Research Articles

Abstract 5136: Adapting long read sequencing technologies for targeted and single cell applications

Abstract Structural variations (SV), a hallmark of genomic instability in cancer, include insertions, deletions, duplications, inversions, or translocations that can either activate oncogenes or inactivate tumor suppressor genes. These variants contribute to the complexity of the cancer genome but have not been successfully resolved yet. While short-read sequencing has aided cancer genomics, it has performed poorly in SV detection, with false positive and false negative rates of 50% or more. Long-read sequencing generates reads that are >10 kb long and has helped identify thousands of genomic features pertinent to cancer that were previously missed by short-read sequencing. Long reads can span SV with a single continuous read, giving a clearer idea of the variation, its position, and size. Currently, long-read sequencing methods are greatly limited by the DNA preparation step, with fragment lengths and yield often compromised by standard extraction methods. The potential to generate reads at the megabase level relies greatly on library preps using large amounts (several μg) of DNA input, which prevents its application to samples like patient tumors that are often limited in DNA. Our overarching goal is to develop a targeted approach that generates large contiguous reads from moderate to low input DNA masses. To this end, we use different extraction, amplification, fragmentation, and enrichment strategies to efficiently target and resolve complex regions of the genome that are often associated with disease conditions. Development of a targeted long-read sequencing strategy using low DNA input would give us the potential to explore the landscape of genetic variation across several individuals and heterogeneous cell populations. We are currently employing and adapting CRISPR-based targeting techniques to enrich for specific regions of the genome from intact cells and/or high-molecular-weight DNA. Citation Format: Shruti Iyer, Sara Goodwin, W. Richard McCombie. Adapting long read sequencing technologies for targeted and single cell applications [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 5136.

Draft Genome of the Rice Coral Montipora capitata Obtained from Linked-Read Sequencing.

The rice coral, Montipora capitata, is widely distributed throughout the Indo-Pacific and comprises one of the most important reef-building species in the Hawaiian Islands. Here, we describe a de novo assembly of its genome based on a linked-read sequencing approach developed by 10x Genomics. The final draft assembly consisted of 27,870 scaffolds with a N50 size of 186 kb and contained a fairly complete set (81%) of metazoan benchmarking (BUSCO) genes. Based on haploid assembly size (615 Mb) and read k-mer profiles, we estimated the genome size to fall between 600 and 700 Mb, although the high fraction of repetitive sequence introduced considerable uncertainty. Repeat analysis indicated that 42% of the assembly consisted of interspersed, mostly unclassified repeats, and almost 3% tandem repeats. We also identified 36,691 protein-coding genes with a median coding sequence length of 807 bp, together spanning 7% of the assembly. The high repeat content and heterozygosity of the genome proved a challenging scenario for assembly, requiring additional steps to merge haplotypes and resulting in a higher than expected fragmentation at the scaffold level. Despite these challenges, the assembly turned out to be comparable in most quality measures to that of other available coral genomes while being considerably more cost-effective, especially with respect to long-read sequencing methods. Provided high-molecular-weight DNA is available, linked-read technology may thus serve as a valuable alternative capable of providing quality genome assemblies of nonmodel organisms.

Presence of Newcastle disease viruses of sub-genotypes Vc and VIn in backyard chickens and in apparently healthy wild birds from Mexico in 2017.

Virulent Newcastle disease viruses (NDV) have been present in Mexico since 1946, and recently, multiple outbreaks have been reported in the country. Here, we characterized eleven NDV isolated from apparently healthy wild birds and backyard chickens in three different locations of Jalisco, Mexico in 2017. Total RNA from NDV was reverse-transcribed, and 1285 nucleotides, which includes 3/4 of the fusion gene, was amplified and sequenced using a long-read MinION sequencing method. The sequences were 99.99-100% identical to the corresponding region obtained using the Illumina MiSeq. Phylogenetic analysis using MinION sequences demonstrated that nine virulent NDV from wild birds belonged to sub-genotypes Vc and VIn, and two backyard chicken isolates were of sub-genotype Vc. The sub-genotype Vc viruses had nucleotide sequence identity that ranged from 97.7 to 98% to a virus of the same sub-genotype isolated from a chicken in Mexico in 2010. Three viruses from pigeons had 96.3-98.7% nucleotide identity to sub-genotype VIn pigeon viruses, commonly referred to as pigeon paramyxovirus, isolated in the USA during 2000-2016. This study demonstrates that viruses of sub-genotype Vc are still present in Mexico, and the detection of this sub-genotype in both chickens and wild birds suggests that transmission among these species may represent a biosecurity risk. This is the first detection and complete genome sequencing of genotype VI NDV from Mexico. In addition, the utilization of an optimized long-read sequencing method for rapid virulence and genotype identification using the Oxford nanopore MinION system is demonstrated.

Supervised Classification of CYP2D6 Genotype and Metabolizer Phenotype With Postmortem Tramadol-Exposed Finns.

Cytochrome p450 family 2, subfamily D, polypeptide 6 (CYP2D6) may be used to infer the metabolizer phenotype (MP) of an individual as poor, intermediate, extensive/normal, or ultrarapid. Metabolizer phenotypes may suggest idiosyncratic drug responses as contributing factors to cause and/or manner of death in postmortem investigations. Application of CYP2D6 has used long-range amplification of the locus and restriction enzyme digestion to detect single-nucleotide variants (SNVs) associated with MPs. This process can be cumbersome and requires knowledge of genotype phase. Phase may be achieved using long-read DNA sequencing and/or computational methods; however, both can be error prone, which may make it difficult or impractical for implementation into medicolegal practice. CYP2D6 was interrogated in postmortem autopsied Finns using supervised machine learning and feature selection to identify SNVs indicative of MP and/or rate of tramadol O-demethylation (T:M1). A subset of 18 CYP2D6 SNVs could predict MP/T:M1 with up to 96.3% accuracy given phased data. These data indicate that phase contributes to classification accuracy when using CYP2D6 data. Of these 18 SNVs, 3 are novel loci putatively associated with T:M1. These findings may enable design of small multiplexes for easy forensic application of MP prediction when cause and/or manner of death is unknown.

Long-Read Sequencing - A Powerful Tool in Viral Transcriptome Research.

Long-read sequencing (LRS) has become increasingly popular due to its strengths in de novo assembly and in resolving complex DNA regions as well as in determining full-length RNA molecules. Two important LRS technologies have been developed during the past few years, including single-molecule, real-time sequencing by Pacific Biosciences, and nanopore sequencing by Oxford Nanopore Technologies. Although current LRS methods produce lower coverage, and are more error prone than short-read sequencing, these methods continue to be superior in identifying transcript isoforms including multispliced RNAs and transcript-length variants as well as overlapping transcripts and alternative polycistronic RNA molecules. Viruses have small, compact genomes and therefore these organisms are ideal subjects for transcriptome analysis with the relatively low-throughput LRS techniques. Recent LRS studies have multiplied the number of previously known transcripts and have revealed complex networks of transcriptional overlaps in the examined viruses.

Improving nanopore read accuracy with the R2C2 method enables the sequencing of highly multiplexed full-length single-cell cDNA

High-throughput short-read sequencing has revolutionized how transcriptomes are quantified and annotated. However, while Illumina short-read sequencers can be used to analyze entire transcriptomes down to the level of individual splicing events with great accuracy, they fall short of analyzing how these individual events are combined into complete RNA transcript isoforms. Because of this shortfall, long-distance information is required to complement short-read sequencing to analyze transcriptomes on the level of full-length RNA transcript isoforms. While long-read sequencing technology can provide this long-distance information, there are issues with both Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) long-read sequencing technologies that prevent their widespread adoption. Briefly, PacBio sequencers produce low numbers of reads with high accuracy, while ONT sequencers produce higher numbers of reads with lower accuracy. Here, we introduce and validate a long-read ONT-based sequencing method. At the same cost, our Rolling Circle Amplification to Concatemeric Consensus (R2C2) method generates more accurate reads of full-length RNA transcript isoforms than any other available long-read sequencing method. These reads can then be used to generate isoform-level transcriptomes for both genome annotation and differential expression analysis in bulk or single-cell samples.

De novo assembly of the complex genome of Nippostrongylus brasiliensis using MinION long reads

BackgroundEukaryotic genome assembly remains a challenge in part due to the prevalence of complex DNA repeats. This is a particularly acute problem for holocentric nematodes because of the large number of satellite DNA sequences found throughout their genomes. These have been recalcitrant to most genome sequencing methods. At the same time, many nematodes are parasites and some represent a serious threat to human health. There is a pressing need for better molecular characterization of animal and plant parasitic nematodes. The advent of long-read DNA sequencing methods offers the promise of resolving complex genomes.ResultsUsing Nippostrongylus brasiliensis as a test case, applying improved base-calling algorithms and assembly methods, we demonstrate the feasibility of de novo genome assembly matching current community standards using only MinION long reads. In doing so, we uncovered an unexpected diversity of very long and complex DNA sequences repeated throughout the N. brasiliensis genome, including massive tandem repeats of tRNA genes.ConclusionBase-calling and assembly methods have improved sufficiently that de novo genome assembly of large complex genomes is possible using only long reads. The method has the added advantage of preserving haplotypic variants and so has the potential to be used in population analyses.

Pacbio sequencing of copper-tolerant Xanthomonas citri reveals presence of a chimeric plasmid structure and provides insights into reassortment and shuffling of transcription activator-like effectors among X. citri strains

BackgroundXanthomonas citri, a causal agent of citrus canker, has been a well-studied model system due to recent availability of whole genome sequences of multiple strains from different geographical regions. Major limitations in our understanding of the evolution of pathogenicity factors in X. citri strains sequenced by short-read sequencing methods have been tracking plasmid reshuffling among strains due to inability to accurately assign reads to plasmids, and analyzing repeat regions among strains. X. citri harbors major pathogenicity determinants, including variable DNA-binding repeat region containing Transcription Activator-like Effectors (TALEs) on plasmids. The long-read sequencing method, PacBio, has allowed the ability to obtain complete and accurate sequences of TALEs in xanthomonads. We recently sequenced Xanthomonas citri str. Xc-03-1638-1-1, a copper tolerant A group strain isolated from grapefruit in 2003 from Argentina using PacBio RS II chemistry. We analyzed plasmid profiles, copy number and location of TALEs in complete genome sequences of X. citri strains.ResultsWe utilized the power of long reads obtained by PacBio sequencing to enable assembly of a complete genome sequence of strain Xc-03-1638-1-1, including sequences of two plasmids, 249 kb (plasmid harboring copper resistance genes) and 99 kb (pathogenicity plasmid containing TALEs). The pathogenicity plasmid in this strain is a hybrid plasmid containing four TALEs. Due to the intriguing nature of this pathogenicity plasmid with Tn3-like transposon association, repetitive elements and multiple putative sites for origins of replication, we might expect alternative structures of this plasmid in nature, illustrating the strong adaptive potential of X. citri strains. Analysis of the pathogenicity plasmid among completely sequenced X. citri strains, coupled with Southern hybridization of the pathogenicity plasmids, revealed clues to rearrangements of plasmids and resulting reshuffling of TALEs among strains.ConclusionsWe demonstrate in this study the importance of long-read sequencing for obtaining intact sequences of TALEs and plasmids, as well as for identifying rearrangement events including plasmid reshuffling. Rearrangement events, such as the hybrid plasmid in this case, could be a frequent phenomenon in the evolution of X. citri strains, although so far it is undetected due to the inability to obtain complete plasmid sequences with short-read sequencing methods.

A survey of localized sequence rearrangements in human DNA.

Genomes mutate and evolve in ways simple (substitution or deletion of bases) and complex (e.g. chromosome shattering). We do not fully understand what types of complex mutation occur, and we cannot routinely characterize arbitrarily-complex mutations in a high-throughput, genome-wide manner. Long-read DNA sequencing methods (e.g. PacBio, nanopore) are promising for this task, because one read may encompass a whole complex mutation. We describe an analysis pipeline to characterize arbitrarily-complex ‘local’ mutations, i.e. intrachromosomal mutations encompassed by one DNA read. We apply it to nanopore and PacBio reads from one human cell line (NA12878), and survey sequence rearrangements, both real and artifactual. Almost all the real rearrangements belong to recurring patterns or motifs: the most common is tandem multiplication (e.g. heptuplication), but there are also complex patterns such as localized shattering, which resembles DNA damage by radiation. Gene conversions are identified, including one between hemoglobin gamma genes. This study demonstrates a way to find intricate rearrangements with any number of duplications, deletions, and repositionings. It demonstrates a probability-based method to resolve ambiguous rearrangements involving highly similar sequences, as occurs in gene conversion. We present a catalog of local rearrangements in one human cell line, and show which rearrangement patterns occur.

Analysis of Multiallelic CNVs by Emulsion Haplotype Fusion PCR.

Emulsion-fusion PCR recovers long-range sequence information by combining products in cis from individual genomic DNA molecules. Emulsion droplets act as very numerous small reaction chambers in which different PCR products from a single genomic DNA molecule are condensed into short joint products, to unite sequences in cis from widely separated genomic sites. These products can therefore provide information about the arrangement of sequences and variants at a larger scale than established long-read sequencing methods. The method has been useful in defining the phase of variants in haplotypes, the typing of inversions, and determining the configuration of sequence variants in multiallelic CNVs. In this description we outline the rationale for the application of emulsion-fusion PCR methods to the analysis of multiallelic CNVs, and give practical details for our own implementation of the method in that context.

NONCODE 2016: an informative and valuable data source of long non-coding RNAs.

NONCODE (http://www.bioinfo.org/noncode/) is an interactive database that aims to present the most complete collection and annotation of non-coding RNAs, especially long non-coding RNAs (lncRNAs). The recently reduced cost of RNA sequencing has produced an explosion of newly identified data. Revolutionary third-generation sequencing methods have also contributed to more accurate annotations. Accumulative experimental data also provides more comprehensive knowledge of lncRNA functions. In this update, NONCODE has added six new species, bringing the total to 16 species altogether. The lncRNAs in NONCODE have increased from 210 831 to 527,336. For human and mouse, the lncRNA numbers are 167,150 and 130,558, respectively. NONCODE 2016 has also introduced three important new features: (i) conservation annotation; (ii) the relationships between lncRNAs and diseases; and (iii) an interface to choose high-quality datasets through predicted scores, literature support and long-read sequencing method support. NONCODE is also accessible through http://www.noncode.org/.