Nanopore Long-read Sequencing Research Articles

Cas9-based enrichment and single-molecule sequencing for precise characterization of genomic duplications

The widespread use of genome-wide diagnostic screening methods has greatly increased the frequency with which incidental (but possibly pathogenic) copy number changes affecting single genes are detected. These findings require validation to allow appropriate clinical management. Deletion variants can usually be readily validated using a range of short-read next-generation sequencing (NGS) strategies, but the characterization of duplication variants at nucleotide resolution remains challenging. This presents diagnostic problems, since pathogenicity cannot generally be assessed without knowing the structure of the variant. We have used a novel Cas9 enrichment strategy, in combination with long-read single-molecule nanopore sequencing, to address this need. We describe the nucleotide-level resolution of two problematic cases, both of whom presented with neurodevelopmental problems and were initially investigated by array CGH. In the first case, an incidental 1.7-kb imbalance involving a partial duplication of VHL exon 3 was detected. This variant was inherited from the patient's father, who had a history of renal cancer at 38 years. In the second case, an incidental ~200-kb de novo duplication that included DMD exons 30–44 was resolved. In both cases, the long-read data yielded sufficient information to enable Sanger sequencing to define the rearrangement breakpoints, and creation of breakpoint-spanning PCR assays suitable for testing of relatives. Our Cas9 enrichment and nanopore sequencing approach can be readily adopted by molecular diagnostic laboratories for cost-effective and rapid characterization of challenging duplication-containing alleles. We also anticipate that in future this method may prove useful for characterizing acquired translocations in tumor cells, and for precisely identifying transgene integration sites in mouse models.

Inside the USCAP Journals

Inside the USCAP Journals

Generation of a Transcriptional Radiation Exposure Signature in Human Blood Using Long-Read Nanopore Sequencing.

In the event of a large-scale event leading to acute ionizing radiation exposure, high-throughput methods would be required to assess individual dose estimates for triage purposes. Blood-based gene expression is a broad source of biomarkers of radiation exposure which have great potential for providing rapid dose estimates for a large population. Time is a crucial component in radiological emergencies and the shipment of blood samples to relevant laboratories presents a concern. In this study, we performed nanopore sequencing analysis to determine if the technology can be used to detect radiation-inducible genes in human peripheral blood mononuclear cells (PBMCs). The technology offers not only long-read sequencing but also a portable device which can overcome issues involving sample shipment, and provide faster results. For this goal, blood from nine healthy volunteers was 2 Gy ex vivo X irradiated. After PBMC isolation, irradiated samples were incubated along with the controls for 24 h at 37°C. RNA was extracted, poly(A)+ enriched and reverse-transcribed before sequencing. The data generated was analyzed using a Snakemake pipeline modified to handle paired samples. The sequencing analysis identified a radiation signature consisting of 46 differentially expressed genes (DEGs) which included 41 protein-coding genes, a long non-coding RNA and four pseudogenes, five of which have been identified as radiation-responsive transcripts for the first time. The genes in which transcriptional expression is most significantly modified after radiation exposure were APOBEC3H and FDXR, presenting a 25- and 28-fold change on average, respectively. These levels of transcriptional response were comparable to results we obtained by quantitative polymerase chain reaction (qPCR) analysis. In vivo exposure analyses showed a transcriptional radioresponse at 24 h postirradiation for both genes together with a strong dose-dependent response in blood irradiated ex vivo. Finally, extrapolating from the data we obtained, the minimum sequencing time required to detect an irradiated sample using APOBEC3H transcripts would be less than 3 min for a total of 50,000 reads. Future improvements, in sample processing and bioinformatic pipeline for specific radiation-responsive transcript identification, will allow the provision of a portable, rapid, real-time biodosimetry platform based on this new sequencing technology. In summary, our data show that nanopore sequencing can identify radiation-responsive genes and can also be used for identification of new transcripts.

Rapid, multiplexed, whole genome and plasmid sequencing of foodborne pathogens using long-read nanopore technology

U.S. public health agencies have employed next-generation sequencing (NGS) as a tool to quickly identify foodborne pathogens during outbreaks. Although established short-read NGS technologies are known to provide highly accurate data, long-read sequencing is still needed to resolve highly-repetitive genomic regions and genomic arrangement, and to close the sequences of bacterial chromosomes and plasmids. Here, we report the use of long-read nanopore sequencing to simultaneously sequence the entire chromosome and plasmid of Salmonella enterica subsp. enterica serovar Bareilly and Escherichia coli O157:H7. We developed a rapid and random sequencing approach coupled with de novo genome assembly within a customized data analysis workflow that uses publicly-available tools. In sequencing runs as short as four hours, using the MinION instrument, we obtained full-length genomes with an average identity of 99.87% for Salmonella Bareilly and 99.89% for E. coli in comparison to the respective MiSeq references. These nanopore-only assemblies provided readily available information on serotype, virulence factors, and antimicrobial resistance genes. We also demonstrate the potential of nanopore sequencing assemblies for rapid preliminary phylogenetic inference. Nanopore sequencing provides additional advantages as very low capital investment and footprint, and shorter (10 hours library preparation and sequencing) turnaround time compared to other NGS technologies.

Complete Genome Sequence of Staphylococcus arlettae Strain P2, Isolated from a Laboratory Environment.

Staphylococcus arlettae is one coagulase-negative species in the bacterial genus Staphylococcus Here, we describe the closed complete genome sequence of S. arlettae strain P2, which was obtained using a hybrid approach combining Oxford Nanopore long-read and Illumina MiSeq short-read sequencing data.

Molecular characterization of carbapenem-resistant Klebsiella pneumoniae isolates with focus on antimicrobial resistance

BackgroundThe enhancing incidence of carbapenem-resistant Klebsiella pneumoniae (CRKP)-mediated infections in Mengchao Hepatobiliary Hospital of Fujian Medical University in 2017 is the motivation behind this investigation to study gene phenotypes and resistance-associated genes of emergence regarding the CRKP strains. In current study, seven inpatients are enrolled in the hospital with complete treatments. The carbapenem-resistant K. pneumoniae whole genome is sequenced using MiSeq short-read and Oxford Nanopore long-read sequencing technology. Prophages are identified to assess genetic diversity within CRKP genomes.ResultsThe investigation encompassed eight CRKP strains that collected from the patients enrolled as well as the environment, which illustrate that blaKPC-2 is responsible for phenotypic resistance in six CRKP strains that K. pneumoniae sequence type (ST11) is informed. The plasmid with IncR, ColRNAI and pMLST type with IncF[F33:A-:B-] co-exist in all ST11 with KPC-2-producing CRKP strains. Along with carbapenemases, all K. pneumoniae strains harbor two or three extended spectrum β-lactamase (ESBL)-producing genes. fosA gene is detected amongst all the CRKP strains. The single nucleotide polymorphisms (SNP) markers are indicated and validated among all CRKP strains, providing valuable clues for distinguishing carbapenem-resistant strains from conventional K. pneumoniae.ConclusionsST11 is the main CRKP type, and blaKPC-2 is the dominant carbapenemase gene harbored by clinical CRKP isolates from current investigations. The SNP markers detected would be helpful for characterizing CRKP strain from general K. pneumoniae. The data provides insights into effective strategy developments for controlling CRKP and nosocomial infection reductions.

Complete Genome Sequence of a Virulent African Swine Fever Virus from a Domestic Pig in Ukraine.

Here, we report the complete genome sequence of an African swine fever (ASF) virus (ASFV/Kyiv/2016/131) isolated from the spleen of a domestic pig in Ukraine with a lethal case of African swine fever. Using only long-read Nanopore sequences, we assembled a full-length genome of 191,911 base pairs in a single contig.

Detection of Complex Structural Variations of X-Linked Deafness-2 (DFNX2) by Single-Molecule Sequencing

Purpose: To investigate genetic causes of X-linked deafness-2 (DFNX2) and analyse the efficiencies and advantages/disadvantages of different sequencing methods. Methods: From 2016 to 2018, 12 unrelated patients with characteristic high-resolution computed tomographic images showing incomplete partition type III (IP-III) were enrolled in this study. Targeted next-generation sequencing (NGS) and Nanopore long-read sequencing were used to analyse the genetic etiology. Results: Six variants of POU3F4 were identified in this cohort by NGS. The six variants were all located upstream of the POU-specific domain and were considered pathogenic due to truncated protein. Among six patients with negative NGS results, three patients were volunteered to choose further Nanopore sequencing and were all found to have structural variations (SVs). We detected an 870 kb deletion in J0012, located at ChrX:81,140,313-82,011,100, approximately 750 kb upstream of the POU3F4 gene. Non-homologous end joining may be a deletion mechanism. An 8 Mb inversion was found in J0011 at ChrX:82,146,071-90,188,492. The sequences of the breakpoint analysis suggested non-allelic homologous recombination occurred within a single chromatid. A 6 kb deletion at ChrX:81,096,652-81,102,705 was found in J0007, which is a region with a high density of retrotransposons. Conclusions: The most common mutations in patients with IP-III were point mutations and small INDELs or exons that can be detected by Sanger sequencing and NGS. Single-molecule long-fragment sequencing technology provides a new technique to detect the pathogenic SVs in patients with negative NGS results. Thus, the DNA diagnostic strategy for IP-III and its efficiency have been greatly improved. At the same time, the sequence characteristics around the breakpoint and the possible mechanism of SV formation should be explored further. Funding Statement: This work was supported by National Key Research and Development Program of China (2016YFC1000700, 2016YFC1000704 and 2017YFC1001800), Shanghai Key Laboratory of Translational Medicine on Ear and Nose diseases (14DZ2260300). Declaration of Interests: The authors have declared that no competing interests exist. Ethics Approval Statement: This investigation was performed with the approval of the ethics committees of the Chinese PLA General Hospital. All of the subjects or their guardians signed the informed consent forms prior to blood sampling.

Hybrid Resistance and Virulence Plasmids in "High-Risk" Clones of Klebsiella pneumoniae, Including Those Carrying blaNDM-5.

Virulence plasmids are associated with hypervirulent types of Klebsiella pneumoniae, which generally do not carry antibiotic resistance genes. In contrast, nosocomial isolates are often associated with resistance, but rarely with virulence plasmids. Here, we describe virulence plasmids in nosocomial isolates of “high-risk” clones of sequence types (STs) 15, 48, 101, 147 and 383 carrying carbapenemase genes. The whole genome sequences were determined by long-read nanopore sequencing. The 12 isolates all contained hybrid plasmids containing both resistance and virulence genes. All carried rmpA/rmpA2 and the aerobactin cluster, with the virulence plasmids of two of three representatives of ST383 carrying blaNDM-5 and seventeen other resistance genes. Representatives of ST48 and ST15 had virulence plasmid-associated genes distributed between two plasmids, both containing antibiotic resistance genes. Representatives of ST101 were remarkable in all sharing virulence plasmids in which iucC and terAWXYZ were missing and iucB and iucD truncated. The combination of resistance and virulence in plasmids of high-risk clones is extremely worrying. Virulence plasmids were often notably consistent within a lineage, even in the absence of epidemiological links, suggesting they are not moving between types. However, there was a common segment containing multiple resistance genes in virulence plasmids of representatives of both STs 48 and 383.

Rapid detection of chromosomal translocation and precise breakpoint characterization in acute myeloid leukemia by nanopore long-read sequencing

Detection of chromosomal translocation is a key component in diagnosis and management of acute myeloid leukemia (AML). Targeted RNA next-generation sequencing (NGS) is emerging as a powerful and clinically practical tool, but it depends on expression of RNA transcript from the underlying DNA translocation. Here, we show the clinical utility of nanopore long-read sequencing in rapidly detecting DNA translocation with exact breakpoints. In a newly diagnosed patient with AML, conventional karyotyping showed translocation t(10;12)(q22;p13) but RNA NGS detected NUP98-NSD1 fusion transcripts from a known cryptic translocation t(5;11)(q35;p15). Rapid PCR-free nanopore whole-genome sequencing yielded a 26,194 bp sequencing read and revealed the t(10;12) breakpoint to be DUSP13 and GRIN2B in head-to-head configuration. This translocation was then classified as a passenger structural variant. The sequencing also yielded a 20,709 bp sequencing read and revealed the t(5;11) breakpoint of the driver NUP98-NSD1 fusion. The identified DNA breakpoints also served as markers for molecular monitoring, in addition to fusion transcript expression by digital PCR and sequence mutations by NGS. We illustrate that third-generation nanopore sequencing is a simple and low-cost workflow for DNA translocation detection.

Complete Genome Sequence of Rubrobacter xylanophilus Strain AA3-22, Isolated from Arima Onsen in Japan.

Rubrobacter xylanophilus strain AA3-22, belonging to the phylum Actinobacteria, was isolated from nonvolcanic Arima Onsen (hot spring) in Japan. Here, we report the complete genome sequence of this organism, which was obtained by combining Oxford Nanopore long-read and Illumina short-read sequencing data.

Identification of intermediate-sized deletions and inference of their impact on gene expression in a human population

BackgroundNext-generation sequencing has allowed for the identification of different genetic variations, which are known to contribute to diseases. Of these, insertions and deletions are the second most abundant type of variations in the genome, but their biological importance or disease association is not well-studied, especially for deletions of intermediate sizes.MethodsWe identified intermediate-sized deletions from whole-genome sequencing (WGS) data of Japanese samples (n = 174) with a novel deletion calling method which considered multiple samples. These deletions were used to construct a reference panel for use in imputation. Imputation was then conducted using the reference panel and data from 82 publically available Japanese samples with gene expression data. The accuracy of the deletion calling and imputation was examined with Nanopore long-read sequencing technology. We also conducted an expression quantitative trait loci (eQTL) association analysis using the deletions to infer their functional impacts on genes, before characterizing the deletions causal for gene expression level changes.ResultsWe obtained a set of polymorphic 4378 high-confidence deletions and constructed a reference panel. The deletions were successfully imputed into the Japanese samples with high accuracy (97.3%). The eQTL analysis identified 181 deletions (4.1%) suggested as causal for gene expression level changes. The causal deletion candidates were significantly enriched in promoters, super-enhancers, and transcription elongation chromatin states. Generation of deletions in a cell line with the CRISPR-Cas9 system confirmed that they were indeed causative variants for gene expression change. Furthermore, one of the deletions was observed to affect the gene expression levels of a gene it was not located in.ConclusionsThis paper reports an accurate deletion calling method for genotype imputation at the whole genome level and shows the importance of intermediate-sized deletions in the human population.

Complete Genomes of Symbiotic Cyanobacteria Clarify the Evolution of Vanadium-Nitrogenase.

Plant endosymbiosis with nitrogen-fixing cyanobacteria has independently evolved in diverse plant lineages, offering a unique window to study the evolution and genetics of plant–microbe interaction. However, very few complete genomes exist for plant cyanobionts, and therefore little is known about their genomic and functional diversity. Here, we present four complete genomes of cyanobacteria isolated from bryophytes. Nanopore long-read sequencing allowed us to obtain circular contigs for all the main chromosomes and most of the plasmids. We found that despite having a low 16S rRNA sequence divergence, the four isolates exhibit considerable genome reorganizations and variation in gene content. Furthermore, three of the four isolates possess genes encoding vanadium (V)-nitrogenase (vnf), which is uncommon among diazotrophs and has not been previously reported in plant cyanobionts. In two cases, the vnf genes were found on plasmids, implying possible plasmid-mediated horizontal gene transfers. Comparative genomic analysis of vnf-containing cyanobacteria further identified a conserved gene cluster. Many genes in this cluster have not been functionally characterized and would be promising candidates for future studies to elucidate V-nitrogenase function and regulation.

Complete Genome Sequence of Bacillus mycoides TH26, Isolated from a Permafrost Thaw Gradient.

Here, we describe the complete genome assembly of a Bacillus mycoides isolate collected from a boreal forest soil associated with permafrost thaw. Using only long-read Nanopore sequences and VolTRAX library preparation, we assembled two circular contigs totaling 5,789,722 bp (N 50, 5,306,036), a complete chromosome with one plasmid.

Nanopype: a modular and scalable nanopore data processing pipeline.

Long-read third-generation nanopore sequencing enables researchers to now address a range of questions that are difficult to tackle with short read approaches. The rapidly expanding user base and continuously increasing throughput have sparked the development of a growing number of specialized analysis tools. However, streamlined processing of nanopore datasets using reproducible and transparent workflows is still lacking. Here we present Nanopype, a nanopore data processing pipeline that integrates a diverse set of established bioinformatics software while maintaining consistent and standardized output formats. Seamless integration into compute cluster environments makes the framework suitable for high-throughput applications. As a result, Nanopype facilitates comparability of nanopore data analysis workflows and thereby should enhance the reproducibility of biological insights. https://github.com/giesselmann/nanopype, https://nanopype.readthedocs.io. Supplementary data are available at Bioinformatics online.

Detection of DNA base modifications by deep recurrent neural network on Oxford Nanopore sequencing data

DNA base modifications, such as C5-methylcytosine (5mC) and N6-methyldeoxyadenosine (6mA), are important types of epigenetic regulations. Short-read bisulfite sequencing and long-read PacBio sequencing have inherent limitations to detect DNA modifications. Here, using raw electric signals of Oxford Nanopore long-read sequencing data, we design DeepMod, a bidirectional recurrent neural network (RNN) with long short-term memory (LSTM) to detect DNA modifications. We sequence a human genome HX1 and a Chlamydomonas reinhardtii genome using Nanopore sequencing, and then evaluate DeepMod on three types of genomes (Escherichia coli, Chlamydomonas reinhardtii and human genomes). For 5mC detection, DeepMod achieves average precision up to 0.99 for both synthetically introduced and naturally occurring modifications. For 6mA detection, DeepMod achieves ~0.9 average precision on Escherichia coli data, and have improved performance than existing methods on Chlamydomonas reinhardtii data. In conclusion, DeepMod performs well for genome-scale detection of DNA modifications and will facilitate epigenetic analysis on diverse species.

Accessing a Russian Wheat Aphid Resistance Gene in Bread Wheat by Long-Read Technologies.

Russian wheat aphid (RWA) ( Kurdjumov) is a serious invasive pest of small-grain cereals and many grass species. An efficient strategy to defy aphid attacks is to identify sources of natural resistance and transfer resistance genes into susceptible crop cultivars. Revealing the genes helps understand plant defense mechanisms and engineer plants with durable resistance to the pest. To date, more than 15 RWA resistance genes have been identified in wheat ( L.) but none of them has been cloned. Previously, we genetically mapped the RWA resistance gene into an interval of 0.83 cM on the short arm of chromosome 7D and spanned it with five bacterial artificial chromosome (BAC) clones. Here, we used a targeted strategy combining traditional approaches toward gene cloning (genetic mapping and sequencing of BAC clones) with novel technologies, including optical mapping and long-read nanopore sequencing. The latter, with reads spanning the entire length of a BAC insert, enabled us to assemble the whole region, a task that was not achievable with short reads. Long-read optical mapping validated the DNA sequence in the interval and revealed a difference in the locus organization between resistant and susceptible genotypes. The complete and accurate sequence of the region facilitated the identification of new markers and precise annotation of the interval, revealing six high-confidence genes. Identification of as the most likely candidate opens an avenue for its validation through functional genomics approaches.

Ultra-deep, long-read nanopore sequencing of mock microbial community standards.

BackgroundLong sequencing reads are information-rich: aiding de novo assembly and reference mapping, and consequently have great potential for the study of microbial communities. However, the best approaches for analysis of long-read metagenomic data are unknown. Additionally, rigorous evaluation of bioinformatics tools is hindered by a lack of long-read data from validated samples with known composition.FindingsWe sequenced 2 commercially available mock communities containing 10 microbial species (ZymoBIOMICS Microbial Community Standards) with Oxford Nanopore GridION and PromethION. Both communities and the 10 individual species isolates were also sequenced with Illumina technology. We generated 14 and 16 gigabase pairs from 2 GridION flowcells and 150 and 153 gigabase pairs from 2 PromethION flowcells for the evenly distributed and log-distributed communities, respectively. Read length N50 ranged between 5.3 and 5.4 kilobase pairs over the 4 sequencing runs. Basecalls and corresponding signal data are made available (4.2 TB in total). Alignment to Illumina-sequenced isolates demonstrated the expected microbial species at anticipated abundances, with the limit of detection for the lowest abundance species below 50 cells (GridION). De novo assembly of metagenomes recovered long contiguous sequences without the need for pre-processing techniques such as binning.ConclusionsWe present ultra-deep, long-read nanopore datasets from a well-defined mock community. These datasets will be useful for those developing bioinformatics methods for long-read metagenomics and for the validation and comparison of current laboratory and software pipelines.

Interchromosomal template-switching as a novel molecular mechanism for imprinting perturbations associated with Temple syndrome

BackgroundIntrachromosomal triplications (TRP) can contribute to disease etiology via gene dosage effects, gene disruption, position effects, or fusion gene formation. Recently, post-zygotic de novo triplications adjacent to copy-number neutral genomic intervals with runs of homozygosity (ROH) have been shown to result in uniparental isodisomy (UPD). The genomic structure of these complex genomic rearrangements (CGRs) shows a consistent pattern of an inverted triplication flanked by duplications (DUP-TRP/INV-DUP) formed by an iterative DNA replisome template-switching mechanism during replicative repair of a single-ended, double-stranded DNA (seDNA), the ROH results from an interhomolog or nonsister chromatid template switch. It has been postulated that these CGRs may lead to genetic abnormalities in carriers due to dosage-sensitive genes mapping within the copy-number variant regions, homozygosity for alleles at a locus causing an autosomal recessive (AR) disease trait within the ROH region, or imprinting-associated diseases.MethodsHere, we report a family wherein the affected subject carries a de novo 2.2-Mb TRP followed by 42.2 Mb of ROH and manifests clinical features overlapping with those observed in association with chromosome 14 maternal UPD (UPD(14)mat). UPD(14)mat can cause clinical phenotypic features enabling a diagnosis of Temple syndrome. This CGR was then molecularly characterized by high-density custom aCGH, genome-wide single-nucleotide polymorphism (SNP) and methylation arrays, exome sequencing (ES), and the Oxford Nanopore long-read sequencing technology.ResultsWe confirmed the postulated DUP-TRP/INV-DUP structure by multiple orthogonal genomic technologies in the proband. The methylation status of known differentially methylated regions (DMRs) on chromosome 14 revealed that the subject shows the typical methylation pattern of UPD(14)mat. Consistent with these molecular findings, the clinical features overlap with those observed in Temple syndrome, including speech delay.ConclusionsThese data provide experimental evidence that, in humans, triplication can lead to segmental UPD and imprinting disease. Importantly, genotype/phenotype analyses further reveal how a post-zygotically generated complex structural variant, resulting from a replication-based mutational mechanism, contributes to expanding the clinical phenotype of known genetic syndromes. Mechanistically, such events can distort transmission genetics resulting in homozygosity at a locus for which only one parent is a carrier as well as cause imprinting diseases.

Long-read direct RNA sequencing by 5’-Cap capturing reveals the impact of Piwi on the widespread exonization of transposable elements in locusts

ABSTRACTThe large genome of the migratory locust (Locusta migratoria) genome accumulates massive amount of accumulated transposable elements (TEs), which show intrinsic transcriptional activities. Hampering the ability to precisely determine full-length RNA transcript sequences are exonized TEs, which produce numerous highly similar fragments that are difficult to resolve using short-read sequencing technology. Here, we applied a 5ʹ-Cap capturing method using Nanopore long-read direct RNA sequencing to characterize full-length transcripts in their native RNA form and to analyze the TE exonization pattern in the locust transcriptome. Our results revealed the widespread establishment of TE exonization and a substantial contribution of TEs to RNA splicing in the locust transcriptome. The results of the transcriptomic spectrum influenced by Piwi expression indicated that TE-derived sequences were the main targets of Piwi-mediated repression. Furthermore, our study showed that Piwi expression regulates the length of RNA transcripts containing TE-derived sequences, creating an alternative UTR usage. Overall, our results reveal the transcriptomic characteristics of TE exonization in the species characterized by large and repetitive genomes.