Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2022 Jul;40(7):1035-1041.
doi: 10.1038/s41587-022-01221-5. Epub 2022 Mar 28.

Accelerated identification of disease-causing variants with ultra-rapid nanopore genome sequencing

Sneha D Goenka #  1 , John E Gorzynski #  1 , Kishwar Shafin #  2 , Dianna G Fisk  3 , Trevor Pesout  2 , Tanner D Jensen  1 , Jean Monlong  2 , Pi-Chuan Chang  4 , Gunjan Baid  4 , Jonathan A Bernstein  1 , Jeffrey W Christle  1 , Karen P Dalton  1 , Daniel R Garalde  5 , Megan E Grove  3 , Joseph Guillory  5 , Alexey Kolesnikov  4 , Maria Nattestad  4 , Maura R Z Ruzhnikov  1 , Mehrzad Samadi  6 , Ankit Sethia  6 , Elizabeth Spiteri  1 , Christopher J Wright  5 , Katherine Xiong  1 , Tong Zhu  6 , Miten Jain  2 , Fritz J Sedlazeck  7 , Andrew Carroll  4 , Benedict Paten  2 , Euan A Ashley  8
Affiliations

Accelerated identification of disease-causing variants with ultra-rapid nanopore genome sequencing

Sneha D Goenka et al. Nat Biotechnol. 2022 Jul.

Abstract

Whole-genome sequencing (WGS) can identify variants that cause genetic disease, but the time required for sequencing and analysis has been a barrier to its use in acutely ill patients. In the present study, we develop an approach for ultra-rapid nanopore WGS that combines an optimized sample preparation protocol, distributing sequencing over 48 flow cells, near real-time base calling and alignment, accelerated variant calling and fast variant filtration for efficient manual review. Application to two example clinical cases identified a candidate variant in <8 h from sample preparation to variant identification. We show that this framework provides accurate variant calls and efficient prioritization, and accelerates diagnostic clinical genome sequencing twofold compared with previous approaches.

PubMed Disclaimer

Conflict of interest statement

J.E.G. owns stock in INVITAE, Illumina and PacBio. K.S. has performed paid internships at NVIDIA Corp. and Google LLC, and presented a talk at an ONT-sponsored event. P.C., G.B, A.K., M.N. and A.C. are employees of Google LLC and own Alphabet stock as part of the standard compensation package. D.R.G., J.G. and C.J.W. are employees of Oxford Nanopore Technologies and share/share option holders. M.S., A.S. and T.Z. are employees of NVIDIA Corp. and own NVIDIA stock as part of the standard compensation package. M.J. has received reimbursement for travel, accommodation and conference fees to speak at events organized by ONT. F.J.S. received travel compensation from Pacific Biotechnology and Oxford Nanopore Technologies. E.A.A. is cofounder of Personalis, Deepcell, and Svexa, Advisor to Apple, and a Non-Executive Director of AstraZeneca. The remaining authors declare no competing interests.

Figures

Fig. 1
Fig. 1. Overview of ultra-rapid computational pipeline.
a, After the start of sequencing on the PromethION48 device, raw signal files are periodically uploaded to cloud storage. Our cloud-based pipeline scales compute-intensive base calling and alignment across 16 instances with 4× Tesla V100 GPUs each and runs concurrently with sequencing. The instances aim for maximum resource utilization, where base calling using Guppy runs on GPU and alignment using Minimap2 (ref. ) runs on 42 virtual CPUs in parallel. b, Once the alignment file is ready, small-variant calling performed using GPU-accelerated PEPPER–Margin–DeepVariant on 14 instances and SV calling using Sniffles on 2 instances. Each instance processes a specific set of contigs. Specific details about the Google Cloud Platform-based instance configurations are provided in Supplementary Table 21. c, These variant calls are annotated to aid in the subsequent variant filtration and prioritization. Our score-based variant filtration method takes in millions of variants reported by the variant caller to surface any deleterious variant for review using Alissa. The filtration method is designed such that it reports a tractable number of variants for manual curation.
Fig. 2
Fig. 2. Comparison of barcoded and nonbarcoded variant calling performance, and standard and ultra-rapid variant filtration performance.
a, Stratified variant calling performance comparison between barcoded and nonbarcoded HG002 samples in all benchmarking and exonic regions. The HG002 nonbarcoded run was the seventh sample run on the flow cells. The similarity in variant calling performance shows that barcoding is not necessary to achieve high-quality variant calls. b, Comparison between standard and ultra-rapid variant filtration pipelines. The standard pipeline picks variants if any of the subcategories are true for any variant, whereas the ultra-rapid pipeline uses a score-based method to surface mostly relevant variants for fast manual review. We see that the standard pipeline flags 147 variants compared with 31 variants proposed by the ultra-rapid pipeline, which substantially reduces the time required for manual review. Source data
Fig. 3
Fig. 3. Ultra-rapid-sequencing pipeline performance.
The detailed schema of the ultra-rapid-sequencing pipeline and end-to-end performance of the ultra-rapid-sequencing pipeline on two different clinical samples. a, The ultra-rapid-sequencing pipeline starts from sample collection on the far left to final diagnosis on the far right. The details of each step are presented in Online Methods. All the steps that run in parallel are vertically stacked. b, The first patient was the fourth sample sequenced on the set of flow cells, resulting in 2:16 h of sequencing. Variant calling completed 6:55 h from the start of sample preparation. Subsequently, variant filtration and manual review identified a candidate variant in gene TNNT2 in 7:18 h. c, The second patient was the sixth sample on the same set of flow cells. The sequencing completed in 2:46 h with 200 Gb followed by another 2:44 h to identify a candidate variant resulting in an end-to-end time of 7:48 h.
See this image and copyright information in PMC

Similar articles

See all similar articles

Cited by

See all "Cited by" articles

References

    1. Ashley EA, et al. Clinical assessment incorporating a personal genome. Lancet. 2010;375:1525–1535. doi: 10.1016/S0140-6736(10)60452-7. - DOI - PMC - PubMed
    1. Dewey FE, et al. Clinical interpretation and implications of whole-genome sequencing. JAMA. 2014;311:1035–1045. doi: 10.1001/jama.2014.1717. - DOI - PMC - PubMed
    1. Buchan JG, White S, Joshi R, Ashley EA. Rapid genome sequencing in the critically ill. Clin. Chem. 2019;65:723–726. doi: 10.1373/clinchem.2018.293506. - DOI - PubMed
    1. Farnaes, L. et al. Rapid whole-genome sequencing decreases infant morbidity and cost of hospitalization. NPJ Genom. Med. 3, 10 (2018). - PMC - PubMed
    1. Saunders CarolJean, et al. Rapid whole-genome sequencing for genetic disease diagnosis in neonatal intensive care units. Sci. Transl. Med. 2012;4:154ra135–154ra135. doi: 10.1126/scitranslmed.3004041. - DOI - PMC - PubMed

Publication types