Next Article in Journal
Serenoa repens and Urtica dioica Fixed Combination: In-Vitro Validation of a Therapy for Benign Prostatic Hyperplasia (BPH)
Next Article in Special Issue
Functional Genomics in Health and Disease
Previous Article in Journal
RB Regulates DNA Double Strand Break Repair Pathway Choice by Mediating CtIP Dependent End Resection
Previous Article in Special Issue
Multi-Omics Data Integration in Extracellular Vesicle Biology—Utopia or Future Reality?
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 23
10.3390/ijms21239177
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Long-Read Sequencing Approach for Direct Haplotype Phasing in Clinical Settings

by
Simone Maestri
,
Maria Giovanna Maturo
,
Emanuela Cosentino
,
Luca Marcolungo
,
Barbara Iadarola
,
Elisabetta Fortunati
,
Marzia Rossato
* and
Massimo Delledonne
Department of Biotechnology, University of Verona, 37134 Verona, Italy
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(23), 9177; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21239177
Submission received: 13 October 2020 / Revised: 20 November 2020 / Accepted: 27 November 2020 / Published: 1 December 2020
(This article belongs to the Special Issue Functional Genomics in Health and Disease)
Browse Figures
Versions Notes

Abstract

:
The reconstruction of individual haplotypes can facilitate the interpretation of disease risks; however, high costs and technical challenges still hinder their assessment in clinical settings. Second-generation sequencing is the gold standard for variant discovery but, due to the production of short reads covering small genomic regions, allows only indirect haplotyping based on statistical methods. In contrast, third-generation methods such as the nanopore sequencing platform developed by Oxford Nanopore Technologies (ONT) generate long reads that can be used for direct haplotyping, with fewer drawbacks. However, robust standards for variant phasing in ONT-based target resequencing efforts are not yet available. In this study, we presented a streamlined proof-of-concept workflow for variant calling and phasing based on ONT data in a clinically relevant 12-kb region of the APOE locus, a hotspot for variants and haplotypes associated with aging-related diseases and longevity. Starting with sequencing data from simple amplicons of the target locus, we demonstrated that ONT data allow for reliable single-nucleotide variant (SNV) calling and phasing from as little as 60 reads, although the recognition of indels is less efficient. Even so, we identified the best combination of ONT read sets (600) and software (BWA/Minimap2 and HapCUT2) that enables full haplotype reconstruction when both SNVs and indels have been identified previously using a highly-accurate sequencing platform. In conclusion, we established a rapid and inexpensive workflow for variant phasing based on ONT long reads. This allowed for the analysis of multiple samples in parallel and can easily be implemented in routine clinical practice, including diagnostic testing.

1. Introduction

When a group of genetic variants is inherited together on a chromosome from one parent because of their genetic linkage, it defines a haplotype [1]. Reconstructing individual haplotypes, a process known as haplotype phasing or haplotyping, has important implications for the interpretation of disease risks. Multiple studies have identified haplotypes (rather than individual variants) as protective or risk factors for various conditions, such as Alzheimer’s disease [2,3,4], thrombosis [5], and breast cancer [6].
The region on chromosome 19, including the genes encoding apolipoprotein E (APOE), translocase of outer mitochondrial membrane 40 homolog (TOMM40), and apolipoprotein C1 (APOC1), is a hotspot for single nucleotide variants (SNVs) and haplotypes associated with Alzheimer’s disease and other age-related diseases, or with longevity [3,4,7,8,9,10]. The APOE ε4 haplotype (rs429358C–rs7412C) is an established risk factor for late-onset Alzheimer’s disease (LOAD) and other age-related diseases [3,4,11,12]. The risk of Alzheimer’s disease is also influenced by different variants within the APOE promoter when they are present in phase (in cis on the same chromosome) with the ε4 haplotype [13]. In particular, the variant –219T rs405509T is located approximately 3000 bp upstream of ε4 and significantly reduces the average age of LOAD diagnosis [4]. Many genome-wide association studies for human longevity have identified significant variants in the TOMM40-APOE-APOC1 region [14,15,16], including a specific haplotype comprising 10 SNVs spanning 26 kb in a Chinese population [17]. However, it is currently impossible to identify carriers of such haplotypes directly by Sanger sequencing or the most common next-generation sequencing (NGS) platforms because they cannot generate reads longer than ~1500 bp from a single DNA molecule. Methods that can directly report haplotype phasing in this genomic region, and in others associated with disease risks, may allow earlier diagnosis and more effective therapy to limit the disease burden.
Technical issues and high costs currently hinder the assessment of haplotypes in clinical settings. Haplotype phasing traditionally involves the analysis of multiple individuals from the same family, to indicate the chromosomal origin of each variant [18]. However, pedigree-based haplotype inference depends on the availability of family DNA samples and cannot phase de novo variations that have arisen in the last generation. Long-range PCR and fluorescence in situ hybridization allow the direct haplotyping of particular homozygous loci, but are not useful for genome-wide analysis [19]. Second-generation NGS technologies, such as Illumina sequencing by synthesis, produce short reads (up to 300 bp) that usually span no more than a single variant, and thus allow only indirect haplotype reconstruction based on statistical inferences from population genotyping data. [18]. To overcome these issues, several library preparation methods that preserve haplotype information in short-read sequencing data have been proposed. These include mate-pair libraries, the single-cell DNA template strand sequencing (Strand-seq) [20], the proximity-ligation libraries of the Chicago [21], and Chromosome Conformation Capture (Hi-C) [22]. More recently, the company 10× Genomics launched a novel microfluidics-based library preparation approach generating linked reads that can be assembled into long haplotypes [23]. Although these approaches allow for accurate haplotyping, ad hoc library preparation protocols require additional sample preparation time and attract higher costs.
Third-generation sequencing technologies developed by companies such as Oxford Nanopore Technologies (ONT) and Pacific Biosciences (PacBio) can facilitate rapid and direct haplotyping due to the production of long sequencing reads (>10 kb) [18,24,25,26]. Differently to short-read sequencing platforms, which rely on amplification to enlarge clusters of a given DNA template, third-generation sequencing platforms directly analyze single DNA molecules in real-time. In particular, the ONT platform identifies DNA bases by measuring the changes in electrical conductivity generated while DNA strands pass through a biological pore [27]. However, ONT technology has a higher error rate than other NGS methods (~10%) and specific adjustments, namely bioinformatic pipelines that deal with this drawback, are therefore required to fully exploit the opportunities provided by ONT sequencing [28]. Some phasing tools that enable haplotype phasing based on ONT sequencing have been developed [29,30,31], but performance benchmarking analysis has yet to be published and there is no standard pipeline for this type of analysis. Moreover, the impact of sequencing coverage on the accuracy of reconstructed haplotypes is unclear.
To address the lack of guidelines for users wishing to perform haplotype phasing based on ONT sequencing reads, we presented a workflow for SNV calling and haplotype phasing of variants localized in a region of interest. The method exploits available and validated software and requires only a modest number of ONT reads, thus enabling efficient haplotype phasing at low costs in clinical settings.

2. Results

2.1. Optimization of APOE Locus Amplification

The genomic region chr19:44,902,730–44,915,006, comprising the whole APOE gene as well as its promoter and small parts of TOMM40 and APOC1, was selected as a target region because it includes several SNVs and haplotypes previously linked to Alzheimer’s disease and other age-related diseases (Figure 1 and Supplementary File S1). According to Illumina whole-genome sequencing data and the GIAB dataset, five available samples showed numerous heterozygous SNVs within this genomic locus (Table 1) and were thus used as case studies for the analysis of variant calling and haplotype phasing in the target region.
To amplify this ~12-kb genomic region, a long-range PCR protocol was tested using either commercially available or custom-designed primer pairs (Figure 2A). Both primer pairs produced amplicons >10 kb that included the target locus, as demonstrated by the positive amplification of APOE exon 3 by nested PCR (Figure 2B). However, the commercial primers produced multiple amplicons including nonspecific ones with a molecular weight <12 kb, whereas the custom primers generated a cleaner but consistently weaker product (Figure 2A). To increase the amplification efficiency of the custom primers, allowing the production of enough DNA for downstream analysis, the reactions were supplemented with different amounts of DMSO and tested on three distinct biological samples. At 4% DMSO, the custom primers achieved the best compromise between efficiency and specificity (Figure 2C), despite yielding minor nonspecific products with lower molecular weights than expected (Figure 2C,E). Therefore, we used two different approaches to purify the specific 12-kb PCR products: Manual gel excision followed by on-column purification, and automatic selection using the BluePippin device. Both methods purified similar amounts of DNA (36% and 39%, respectively) but the BluePippin method yielded a significantly cleaner product, verified by manual and automatic electrophoretic analysis (Figure 2D,E). The TOMM40-APOE-APOC1 target region was therefore amplified in all samples using the custom primers in the presence of 4% DMSO with subsequent BluePippin selection.

2.2. Variant Calling Performance Using ONT Data

The five samples with available genotyping data (Supplementary File S1 and Table 2) were amplified in the target locus, multiplexed, and sequenced on a R9.4 flow-cell, yielding an average of 47,835 reads per sample (Table 3). An average of 99.99% MinION PASS reads mapped properly to the target region and 98.07% of these reads spanned the entire 12-kb target.
Variant calling was carried out using sets composed of 10–10,000 ONT reads and we compared the results with the set of ground-truth variants derived from Illumina data. ONT-based variant calling showed poor results for indels, regardless of the sequence coverage, generating an average number of true positive (TP), false negative (FN), and false positive (FP) calls of 0.2, 1.4, and 2.2, respectively, at 60× coverage, and no improvement at greater sequencing depth. In contrast, SNV calling was already highly accurate when using 60 reads, generating an average number of TP, FN, and FP calls of 9.8, 0.2, and 0.0, respectively (Table 4). Only one homozygous SNV from sample V1 (found within a homopolymer run) was not properly called (Figure 3). This demonstrates that variant calling based on ONT reads is reliable only for SNVs, as long as they are not embedded in long homopolymer runs. Interestingly, the SNV calling performance reached a plateau from 60× to 300× coverage, and started to decline at 600× coverage, as shown by the decrease in the F1 score (Figure 4)

2.3. Optimization of a Haplotype Phasing Pipeline Based on ONT Data

The reference haplotypes for three of the five samples (V1, V2, and V3) were identified using the 10× Genomics Chromium platform by whole-genome sequencing at an average coverage of 33× (Supplementary File S2 and S3). The 10× Genomics analysis provided phasing information for 39 of the 40 heterozygous variants identified across these samples (Table 2). For the two remaining samples, haplotypes were determined based on trio analysis, starting with the publicly available VCF file of individuals from the 1000 Genomes Project (Table 2 and Supplementary File S3).
Compared to reference haplotypes, SNVs called at 60× coverage based on ONT data were all correctly phased by the variant calling pipeline, which used Minimap2 for alignment and WhatsHap for phasing (Table 2). In contrast, we could not evaluate indel phasing because this type of variant was not properly called by ONT. In order to assess the capability of ONT data to phase both SNVs and indels, the set of SNVs properly called by ONT (all true positives) was integrated with the truth set of indels (derived from reference VCF data) and phased solely using ONT data. At this aim, we assessed the performance of different software combinations with different numbers of ONT reads. In particular, we compared the haplotypes obtained when using Minimap2 [32] or BWA mem [33] for alignment, HapCUT2 [29] or WhatsHap [30] for phasing, and read sets ranging from 10 to 10,000 (Figure 5). In all cases, accuracy was estimated over 100 iterations. When starting from 600 read groups, all software combinations using HapCUT2 for phasing identified the correct haplotype with 100% accuracy for all samples, and the alignment software had a minimal impact on phasing accuracy (Figure 6). These results demonstrated that phasing both SNVs and indels is possible with no residual errors when exploiting at least 600 ONT reads using the combination of Minimap2 or BWA with HapCUT2.

3. Discussion

Despite the importance of haplotype reconstruction for the interpretation of disease risks and therapeutic responses, haplotype assessment in clinical settings is hindered by technical difficulties and high costs. Short-read sequencing (e.g., the Illumina platform) is regarded as the state-of-the-art approach for SNV discovery in the clinic. However, direct haplotyping is restricted to ~600 bp due to the maximum read length; thus, haplotypes in longer regions can only be inferred by statistical analysis. Emerging approaches such as the Hi-C method or 10× Genomics linked-read sequencing provide additional opportunities for haplotype analysis, but they require laborious and costly sample preparation methods and/or expensive equipment. Moreover, 10× Genomics announced the discontinuation of its chromium genome and exome product lines at the end of June 2020 [34].
Third-generation sequencing methods such as the ONT platform offer a promising alternative for direct haplotyping thanks to the production of long reads. Although its per-base accuracy and cost per base have improved significantly, ONT is not yet on par with short-read platforms [35,36,37]. However, in contrast to other long-read technologies such as PacBio SMRT sequencing, ONT provides flexibility in terms of data acquisition because the user can stop the sequencing run having reached the desired throughput (real-time data analysis). Moreover, the ONT platform is based on affordable instruments such as the MinION sequencer ($US 1000) and it can exploit sequencing flow cells with different output capacities, allowing a wide range of applications. ONT sequencing may therefore represent a cost-effective solution for direct SNV calling and phasing or may complement short-read technologies for the phasing of variants previously identified using these platforms.
One drawback of ONT sequencing for variant phasing is the absence of a well-defined and easy-to-use workflow from sample preparation to ONT data analysis. We therefore established a simple laboratory workflow and bioinformatics pipeline to call and phase heterozygous SNVs in any region of interest based on sole ONT sequencing. In addition, we showed that ONT reads can also be used to reconstruct correct haplotypes for both SNVs and indels, when the latter have been previously identified by high-accuracy methods such as Illumina or Sanger. The work identified the optimal coverage and the best combination of software required for variant phasing based on ONT reads, aspects not dissected in the available literature yet. Compared to currently available approaches, our workflow produces reliable results in a short time and with minimal additional costs, due to the straightforward multiplexing of samples in a single flow cell.
Attempts of using the ONT platform for direct haplotype phasing have been reported since its release to the community through the MinION Access Program (MAP). For example, Ammar et al. [38] generated nanopore reads to detect variants and haplotypes in genes involved in drug responses, whereas Stancu et al. [39] used whole-genome nanopore sequencing reads to phase heterozygous variants previously identified using the Illumina platform in a trio-based analysis. Leija-Salazar [40] evaluated the ability of nanopore sequencing to phase variants in the GBA gene in order to identify patients suffering from Gaucher’s disease or Parkinson’s disease. However, these feasibility studies did not determine the accuracy of haplotyping. Recently, the performance of four SNV-calling tools was evaluated on nanopore data, but the authors did not consider the impact of phasing tools [35].
We addressed the lack of benchmarking data by assessing the impact of alignment and phasing tools on variant calling and phasing accuracy with reference to sequence coverage. By calculating the percentage of properly reconstructed haplotypes, we showed that the phasing software has a much greater impact than the alignment software. Specifically, the medaka_variant workflow provided by ONT (integrating Medaka, Minimap2, and WhatsHap) called SNVs with high accuracy in all samples, with the exception of one SNV embedded in a homopolymer run. Optimal performance was achieved using 60–600 ONT reads. Interestingly, the SNV calling accuracy declined significantly when we used 1000 and 10,000 reads, which is important because high sequence coverage is very common in the context of target resequencing studies [41,42,43]. Users of the Medaka tool should pay particular attention to this drawback. In contrast, the variant calling pipeline based on sole ONT data performed less well for the identification of indels. This is likely to reflect the known bias of ONT reads, which are especially prone to errors in homopolymer runs [26,41,44]. Improvements in base-calling accuracy are constantly contributing to reduce this type of errors, with Bonito now representing a very promising base-caller [45,46]. Moreover, the R10.3 chemistry recently released by ONT is expected to improve the accuracy of indel calling because it uses pores with a longer barrel and a dual reader head, which achieves better raw read accuracy even in homopolymer regions [44,47]. Although indel calling from R9.4 data is not yet accurate, we showed that ONT reads can be used to accurately phase both SNVs and indels when these have already been identified using a more accurate sequencing approach (such as Illumina or Sanger). For this purpose, we have established a second-step pipeline where ONT reads are aligned again to the reference genome using BWA or Minimap2 and phasing of variants from a given VCF file is defined with HapCUT2. For this phasing step, 600 ONT reads allowed the recovery of fully correct haplotypes.
Our proof-of-concept study focused on the APOE locus, a clinically relevant genomic region, allowing for the analysis of numerous SNVs and haplotypes that have been linked to Alzheimer’s disease and other age-related diseases. For example, the proposed workflow may be used to assess the phasing of SNVs >3000 bp upstream of the APOE gene that may affect the disease outcome when present in the cis configuration [17]. Moreover, because SNV frequencies within the APOE locus are highly skewed to population-specific haplotypes [12], our workflow will allow wider population studies to assess how the ancestral background within different sites of the APOE gene shapes the disease phenotype. Our workflow can also be modified easily for the analysis of other genomic regions containing clinically relevant haplotypes. The enrichment of the APOE target locus was achieved by direct amplification of the selected region using long-range PCR, and this approach could be applied to any target of interest. Simple experimental parameters (i.e., PCR primer position and buffer composition) were tuned to maximize amplicon production and both tested methods for amplicon purification yielded enough DNA for ONT sequencing. Even so, the generation of long amplicons can be challenging in repetitive/highly homologous or GC-rich genomic targets. Innovative target-enrichment approaches such as the CRISPR/Cas9-based approach [35,48] and X-Drop technology [49], recently coupled to nanopore sequencing, may help to overcome the limitations of a PCR-based workflow and thus facilitate the broader implementation of our bioinformatics pipeline.
Given that the established pipeline requires only 600 reads per sample, multiplexing to 96 samples using the Native Barcoding Kit can be achieved in a single MinION flow cell, reducing costs to ~€20 per sample. Costs could be reduced even further by using the recent ONT Flongle flow cells, which are specifically developed for applications requiring rapid, low-cost analysis such as diagnostic applications. The simplicity and versatility of our proposed workflow therefore makes it ideal for clinical settings, where haplotype information could facilitate decision-making processes related to disease management [50].

4. Materials and Methods

4.1. DNA Extraction

DNA of donors V1, V2, and V3 was extracted from 2 mL of whole blood using the Invisorb Mini Kit (Stratec Molecular, Berlin, Germany) following the manufacturer’s instructions. DNA purity was assessed using a Nanodrop 1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and DNA integrity was assessed by capillary electrophoresis using a TapeStation 2200 system (Agilent Technologies, Santa Clara, CA, USA). DNA was quantified using the Qubit dsDNA HS Assay kit (Thermo Fisher Scientific).

4.2. PCR Amplification of The Target Region

Long-range PCR was carried out in a 25-µl reaction containing 100 ng of genomic template DNA, 0%, 2.5%, 4%, or 5% DMSO, LongAmp Taq Polymerase 2× master mix (New England Biolabs, Ipswich, MA, USA), and 0.4 µM of one of the following primer pairs: Qiagen forward primer (Hs_APOE_01_LR 5′-CCG TCC ACT TTC CCA TCT CCT CGG TAT AAA TC-3′) and Qiagen reverse primer (St000357847 5′-AGA GTC AGG AAG ATT GAG AGG TGA GAG TGC-3′), or custom forward primer (5′-CTT CAC GGG ACT CAG TGA GAG GAA CAG ATT C-3′) and custom reverse primer (5′-TAG CCA CCT TCA CTC TTC CAA ATC TCC AA-3′). The reaction was heated to 94 °C for 3 min, followed by 35 cycles of 94 °C for 30 s, 64 °C for 60 s, and 65 °C for 9 min 20 s, and a final extension at 65 °C for 10 min. PCR products generated using the custom primers were separated by agarose gel electrophoresis and purified using the Zymoclean Large Fragment DNA Recovery Kit after manually excising the gel band (Zymo Research, Irvine, CA, USA). Alternatively, we applied the high pass protocol for fragments larger than 10 kb on BluePippin (Sage Science, Beverly, MA, USA) using a 0.75% agarose gel cassette. The appropriate size of purified amplicons was verified by capillary electrophoresis on the TapeStation 2200 system. An average of four reactions per sample was necessary to yield enough amplified DNA for nanopore sequencing.

4.3. 10× Genomics Library Preparation, Sequencing, and Data Analysis

Whole-genome sequencing libraries were constructed using the Chromium Gel Bead and Library Kit (10× Genomics, Pleasanton, CA, USA) with the Chromium instrument (10× Genomics) at the Functional Genomics Center of Zurich. Sequencing was performed at Macrogen (Seoul, Korea) on an Illumina (San Diego, CA, USA) Hiseq X producing 400 million paired-end reads of 150 p. LongRanger v2.2.2 was used to map Illumina reads to the hg38 human reference genome, and to call and phase the variants. Only variants with the “PASS” flag were retained for subsequent analysis. Unphased and phased heterozygous variants present in the target region were stored in two separate VCF files. The VCF files of the reference samples NA12878 and NA12892 were downloaded from ftp://ftp-trace.ncbi.nih.gov/1000genomes/ftp/technical/working/20140625_high_coverage_trios_broad/CEU.wgs.consensus.20131118.snps_indels.high_coverage_pcr_free_v2.genotypes.vcf.gz. Variants in the region of interest were extracted from the VCF files and genomic coordinates were lifted from the hg19 to the hg38 reference genome using CrossMap v0.3.0 [51]. The phasing of selected variants was determined using GATK PhaseByTransmission v3.8.1 [52].

4.4. ONT Sequencing and Data Analysis

ONT MinION 1D libraries were prepared from purified amplicons (~1 μg for each sample) using the Native Barcoding Genomic DNA Kit (SQK-LSK108) in accordance with the manufacturer’s instructions (Oxford Nanopore Technologies, Oxford, UK). The five samples were sequenced on a single R9.4 flow-cell in a 48-h sequencing run using the NC_48Hr_sequencing_Run_FLO_MIN106_SQK-LSK108 protocol and MinKNOW v1.2.8. Base-calling of raw fast5 files was carried out using Guppy v3.6.0 with high-accuracy mode. Reads were then demultiplexed using Guppy v3.6.0 by requiring the presence of barcodes at both ends of the read (--require_barcode_both_ends). Reads with quality <7 and length <11,276 nt or >13,276 nt were discarded using NanoFilt [53]. Reads were randomly subsampled using seqtk sample v1.3-r106 (https://github.com/lh3/seqtk). Variants were called using the medaka_variant subroutine of Medaka v0.13.0 (https://github.com/nanoporetech/medaka) after mapping ONT reads with Minimap2 v2.17-r941 [32] to hg38 region chr19:44,899,031–chr19:44,916,122, which we describe hereafter as the hg38_chr19_APOE_locus. Only variants overlapping hg38_chr19_APOE_locus:3730–15948 and with a “PASS” flag were retained using bedtools intersect [54]. Variant calling performance was evaluated using hap.py script (https://github.com/Illumina/hap.py).
Haplotyping was performed by mapping ONT reads to the hg38_chr19_APOE_locus using Minimap2 v2.17-r941 [32] or BWA mem v0.7.17-r1188 [33] for alignment to produce BAM files. HapCUT2 v1.3.2 [29] or WhatsHap v1.0 [30] were used to phase the heterozygous variants stored in the Illumina VCF files, relying on the information included in the BAM file. The accuracy of haplotype phasing was determined by counting the number of times that the entirely correct haplotype was reconstructed for each sample, over 100 iterations on different sizes of read sets, ranging from 10 to 10,000 reads. All the scripts are available at https://github.com/MaestSi/ONT_preprocessing and https://github.com/MaestSi/ONToHap.

5. Conclusions

The work presented provides broadly applicable guidelines, namely number of reads and bioinformatic software, to exploit ONT data for the phasing of genetic variants in a selected region of interest. The rapid and inexpensive workflow proposed in this work allows for the analysis of multiple samples in parallel and can easily be implemented in routine clinical practice, including diagnostic testing.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/21/23/9177/s1.

Author Contributions

Conceptualization, M.D. and M.R.; methodology, S.M., M.G.M., M.R., and M.D.; software, S.M.; validation, S.M. and M.G.M.; formal analysis, S.M. and M.R.; investigation, S.M., M.G.M., L.M., B.I., E.F., and E.C.; writing–original draft preparation, S.M. and M.R.; writing–review and editing, M.R. and M.D.; visualization, S.M.; supervision, M.R. and M.D.; project administration, M.R. and M.D.; funding acquisition, M.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Scitable by Nature Education. Available online: https://0-www-nature-com.brum.beds.ac.uk/scitable/definition/haplotype-haplotypes-142/ (accessed on 30 November 2020).
  2. Allen, M.; Kachadoorian, M.; Quicksall, Z.; Zou, F.; Chai, H.S.; Younkin, C.; E Crook, J.; Pankratz, V.S.; Carrasquillo, M.M.; Krishnan, S.; et al. Association of MAPT haplotypes with Alzheimer’s disease risk and MAPT brain gene expression levels. Alzheimer’s Res. Ther. 2014, 6, 39. [Google Scholar] [CrossRef] [Green Version]
  3. Williams, M.A.; McKay, G.J.; Carson, R.; Craig, D.; Silvestri, G.; Passmore, P. Age-Related Macular Degeneration-Associated Genes in Alzheimer Disease. Am. J. Geriatr. Psychiatry 2015, 23, 1290–1296. [Google Scholar] [CrossRef]
  4. Lescai, F.; Chiamenti, A.M.; Codemo, A.; Pirazzini, C.; D’Agostino, G.; Ruaro, C.; Ghidoni, R.; Benussi, L.; Galimberti, D.; Esposito, F.; et al. An APOE Haplotype Associated with Decreased epsilon4 Expression Increases the Risk of Late Onset Alzheimer’s Disease. J. Alzheimer’s Dis. 2011, 24, 235–245. [Google Scholar] [CrossRef]
  5. Navarro, S.; Medina, P.; Mira, Y.; Estelles, A.; Villa, P.; Ferrando, F.; Vaya, A.; Bertina, R.M.; España, F. Haplotypes of the EPCR gene, prothrombin levels, and the risk of venous thrombosis in carriers of the prothrombin G20210A mutation. Haematol. 2008, 93, 885–891. [Google Scholar] [CrossRef] [Green Version]
  6. Vymetalkova, V.P.; Soucek, P.; Kunická, T.; Jiraskova, K.; Brynychová, V.; Pardini, B.; Novosadova, V.; Polivkova, Z.; Kubáčková, K.; Kozevnikovova, R.; et al. Genotype and Haplotype Analyses of TP53 Gene in Breast Cancer Patients: Association with Risk and Clinical Outcomes. PLOS ONE 2015, 10, e0134463. [Google Scholar] [CrossRef]
  7. Schächter, F.; Faure-Delanef, L.; Guénot, F.; Rouger, H.; Froguel, P.; Lesueur-Ginot, L.; Cohen, D. Genetic associations with human longevity at the APOE and ACE loci. Nat. Genet. 1994, 6, 29–32. [Google Scholar] [CrossRef]
  8. Soerensen, M.; Dato, S.; Tan, Q.; Thinggaard, M.; Kleindorp, R.; Beekman, M.; Suchiman, H.E.D.; Jacobsen, R.; McGue, M.; Stevnsner, T.; et al. Evidence from case-control and longitudinal studies supports associations of genetic variation in APOE, CETP, and IL6 with human longevity. AGE 2012, 35, 487–500. [Google Scholar] [CrossRef]
  9. Ferrari, R.; Wang, Y.; Vandrovcova, J.; Guelfi, S.; Witeolar, A.; Karch, C.M.; Schork, A.J.; Fan, C.C.; Brewer, J.B.; Momeni, P.; et al. Genetic architecture of sporadic frontotemporal dementia and overlap with Alzheimer’s and Parkinson’s diseases. J. Neurol. Neurosurg. Psychiatry 2017, 88, 152–164. [Google Scholar] [CrossRef]
  10. Coon, K.D.; Myers, A.J.; Craig, D.W.; Webster, J.A.; Pearson, J.V.; Lince, D.H.; Zismann, V.L.; Beach, T.G.; Leung, D.; Bryden, L.; et al. A High-Density Whole-Genome Association Study Reveals That APOE Is the Major Susceptibility Gene for Sporadic Late-Onset Alzheimer’s Disease. J. Clin. Psychiatry 2007, 68, 613–618. [Google Scholar] [CrossRef]
  11. Corder, E.H.; Saunders, A.M.; Strittmatter, W.J.; E Schmechel, D.; Gaskell, P.C.; Small, G.W.; Roses, A.D.; Haines, J.L.; A Pericak-Vance, M. Gene dose of apolipoprotein E type 4 allele and the risk of Alzheimer’s disease in late onset families. Science 1993, 261, 921–923. [Google Scholar] [CrossRef]
  12. Babenko, V.N.; Afonnikov, D.A.; Ignatieva, E.V.; Klimov, A.V.; Gusev, F.E.; Rogaev, E.I. Haplotype analysis of APOE intragenic SNPs. BMC Neurosci. 2018, 19, 29–40. [Google Scholar] [CrossRef]
  13. Xiao, H.; Gao, Y.; Liu, L.; Li, Y. Association between polymorphisms in the promoter region of the apolipoprotein E (APOE) gene and Alzheimer’s disease: A meta-analysis. EXCLI J. 2017, 16, 921–938. [Google Scholar]
  14. Deelen, J.; Beekman, M.; Uh, H.-W.; Helmer, Q.; Kuningas, M.; Christiansen, L.; Kremer, D.; Van Der Breggen, R.; Suchiman, H.E.D.; Lakenberg, N.; et al. Genome-wide association study identifies a single major locus contributing to survival into old age; the APOE locus revisited. Aging Cell 2011, 10, 686–698. [Google Scholar] [CrossRef]
  15. Deelen, J.; Beekman, M.; Uh, H.-W.; Broer, L.; Ayers, K.L.; Tan, Q.; Kamatani, Y.; Bennet, A.M.; Tamm, R.; Trompet, S.; et al. Genome-wide association meta-analysis of human longevity identifies a novel locus conferring survival beyond 90 years of age. Hum. Mol. Genet. 2014, 23, 4420–4432. [Google Scholar] [CrossRef]
  16. Nebel, A.; Kleindorp, R.; Caliebe, A.; Nothnagel, M.; Blanché, H.; Junge, O.; Wittig, M.; Ellinghaus, D.; Flachsbart, F.; Wichmann, H.-E.; et al. A genome-wide association study confirms APOE as the major gene influencing survival in long-lived individuals. Mech. Ageing Dev. 2011, 132, 324–330. [Google Scholar] [CrossRef]
  17. Lin, R.; Zhang, Y.; Yan, D.; Liao, X.; Gong, G.; Hu, J.; Fu, Y.; Cai, W. Association of common variants in TOMM40/APOE/APOC1 region with human longevity in a Chinese population. J. Hum. Genet. 2016, 61, 323–328. [Google Scholar] [CrossRef]
  18. Snyder, M.W.; Adey, A.; Kitzman, J.O.; Shendure, J. Haplotype-resolved genome sequencing: Experimental methods and applications. Nat. Rev. Genet. 2015, 16, 344–358. [Google Scholar] [CrossRef]
  19. Huang, M.; Tu, J.; Lu, Z. Recent Advances in Experimental Whole Genome Haplotyping Methods. Int. J. Mol. Sci. 2017, 18, 1944. [Google Scholar] [CrossRef] [Green Version]
  20. Porubský, D.; Sanders, A.D.; Van Wietmarschen, N.; Falconer, E.; Hills, M.; Spierings, D.C.; Bevova, M.R.; Guryev, V.; Lansdorp, P.M. Direct chromosome-length haplotyping by single-cell sequencing. Genome Res. 2016, 26, 1565–1574. [Google Scholar] [CrossRef] [Green Version]
  21. Putnam, N.H.; O’Connell, B.L.; Stites, J.C.; Rice, B.J.; Blanchette, M.; Calef, R.; Troll, C.J.; Fields, A.; Hartley, P.D.; Sugnet, C.W.; et al. Chromosome-scale shotgun assembly using an in vitro method for long-range linkage. Genome Res. 2016, 26, 342–350. [Google Scholar] [CrossRef] [Green Version]
  22. Lieberman-Aiden, E.; Van Berkum, N.L.; Williams, L.; Imakaev, M.; Ragoczy, T.; Telling, A.; Amit, I.; Lajoie, B.R.; Sabo, P.J.; Dorschner, M.O.; et al. Comprehensive Mapping of Long-Range Interactions Reveals Folding Principles of the Human Genome. Science 2009, 326, 289–293. [Google Scholar] [CrossRef] [Green Version]
  23. Zheng, G.X.; Lau, B.T.; Schnall-Levin, M.; Jarosz, M.; Bell, J.M.; Hindson, C.M.; Kyriazopoulou-Panagiotopoulou, S.; Masquelier, D.A.; Merrill, L.; Terry, J.M.; et al. Haplotyping germline and cancer genomes with high-throughput linked-read sequencing. Nat. Biotechnol. 2016, 34, 303–311. [Google Scholar] [CrossRef] [PubMed]
  24. Jain, M.; Koren, S.; Miga, K.H.; Quick, J.; Rand, A.C.; A Sasani, T.; Tyson, J.R.; Beggs, A.D.; Dilthey, A.T.; Fiddes, I.T.; et al. Nanopore sequencing and assembly of a human genome with ultra-long reads. Nat. Biotechnol. 2018, 36, 338–345. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Porubsky, D.; Garg, S.; Sanders, A.D.; Korbel, J.O.; Guryev, V.; Lansdorp, P.M.; Marschall, T. Dense and accurate whole-chromosome haplotyping of individual genomes. Nat. Commun. 2017, 8, 1293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Rang, F.J.; Kloosterman, W.P.; De Ridder, J. From squiggle to basepair: Computational approaches for improving nanopore sequencing read accuracy. Genome Biol. 2018, 19, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Lu, H.; Giordano, F.; Ning, Z. Oxford Nanopore MinION Sequencing and Genome Assembly. Genom. Proteom. Bioinform. 2016, 14, 265–279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Laver, T.W.; Caswell, R.C.; Moore, K.A.; Poschmann, J.; Johnson, M.B.; Owens, M.M.; Ellard, S.; Paszkiewicz, K.H.; Weedon, M.N. Pitfalls of haplotype phasing from amplicon-based long-read sequencing. Sci. Rep. 2016, 6, 21746. [Google Scholar] [CrossRef] [Green Version]
  29. Edge, P.; Bafna, V.; Bansal, V. HapCUT2: Robust and accurate haplotype assembly for diverse sequencing technologies. Genome Res. 2017, 27, 801–812. [Google Scholar] [CrossRef]
  30. Patterson, M.; Marschall, T.; Pisanti, N.; Van Iersel, L.; Stougie, L.; Klau, G.W.; Hu, Y.-J. WhatsHap: Weighted Haplotype Assembly for Future-Generation Sequencing Reads. J. Comput. Biol. 2015, 22, 498–509. [Google Scholar] [CrossRef]
  31. Kuleshov, V. Probabilistic single-individual haplotyping. Bioinform. 2014, 30, i379–i385. [Google Scholar] [CrossRef] [Green Version]
  32. Li, H. Minimap2: Pairwise alignment for nucleotide sequences. Bioinform. 2018, 34, 3094–3100. [Google Scholar] [CrossRef]
  33. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. 10x Genomics. Discontinuation of Linked-Reads. Available online: https://www.10xgenomics.com/products/linked-reads (accessed on 30 November 2020).
  35. Gilpatrick, T.; Lee, I.; Graham, J.E.; Raimondeau, E.; Bowen, R.; Heron, A.; Downs, B.; Sukumar, S.; Sedlazeck, F.J.; Timp, W. Targeted nanopore sequencing with Cas9-guided adapter ligation. Nat. Biotechnol. 2020, 38, 433–438. [Google Scholar] [CrossRef]
  36. Edge, P.; Bansal, V. Longshot enables accurate variant calling in diploid genomes from single-molecule long read sequencing. Nat. Commun. 2019, 10, 1–10. [Google Scholar] [CrossRef] [Green Version]
  37. Luo, R.; Wong, C.-L.; Wong, Y.-S.; Tang, C.-I.; Liu, C.-M.; Leung, C.-M.; Lam, T.-W. Exploring the limit of using a deep neural network on pileup data for germline variant calling. Nat. Mach. Intell. 2020, 2, 220–227. [Google Scholar] [CrossRef]
  38. Ammar, R.; Paton, T.A.; Torti, D.; Shlien, A.; Bader, G.D. Long read nanopore sequencing for detection of HLA and CYP2D6 variants and haplotypes. F1000Research 2015, 4, 17. [Google Scholar] [CrossRef] [Green Version]
  39. Stancu, M.C.; Van Roosmalen, M.J.; Renkens, I.; Nieboer, M.M.; Middelkamp, S.; De Ligt, J.; Pregno, G.; Giachino, D.; Mandrile, G.; Valle-Inclan, J.E.; et al. Mapping and phasing of structural variation in patient genomes using nanopore sequencing. Nat. Commun. 2017, 8, 1–13. [Google Scholar] [CrossRef] [Green Version]
  40. Leija-Salazar, M.; Sedlazeck, F.; Toffoli, M.; Mullin, S.; Mokretar, K.; Athanasopoulou, M.; Donald, A.; Sharma, R.; Hughes, D.; Schapira, A.H.V.; et al. Detection of GBA missense mutations and other variants using the Oxford Nanopore MinION. Mol. Genet. Genomic. Med. 2019, 7, 3. [Google Scholar] [CrossRef] [Green Version]
  41. Maestri, S.; Cosentino, E.; Paterno, M.; Freitag, H.; Garces, J.M.; Marcolungo, L.; Alfano, M.; Njunjić, I.; Schilthuizen, M.; Slik, F.J.; et al. A Rapid and Accurate MinION-Based Workflow for Tracking Species Biodiversity in the Field. Genes 2019, 10, 468. [Google Scholar] [CrossRef] [Green Version]
  42. Menegon, M.; Cantaloni, C.; Rodriguez-Prieto, A.; Centomo, C.; Abdelfattah, A.; Rossato, M.; Bernardi, M.; Xumerle, L.; Loader, S.; Delledonne, M. On site DNA barcoding by nanopore sequencing. PLoS ONE 2017, 12, e0184741. [Google Scholar] [CrossRef] [Green Version]
  43. Knot, I.E.; Zouganelis, G.D.; Weedall, G.D.; Wich, S.A.; Rae, R. DNA Barcoding of Nematodes Using the MinION. Front. Ecol. Evol. 2020, 8, 8. [Google Scholar] [CrossRef]
  44. Karst, S.; Ziels, R.; Kirkegaard, R.; Sørensen, E.; McDonald, D.; Zhu, Q.; Knight, R.; Albertsen, M. Enabling high-accuracy long-read amplicon sequences using unique molecular identifiers with Nanopore or PacBio sequencing. bioRxiv 2020. Available online: https://www.biorxiv.org/content/10.1101/645903v3.full (accessed on 30 November 2020).
  45. Wick, R.R.; Judd, L.M.; E Holt, K. Performance of neural network basecalling tools for Oxford Nanopore sequencing. Genome Biol. 2019, 20, 1–10. [Google Scholar] [CrossRef] [Green Version]
  46. Vereecke, N.; Bokma, J.; Haesebrouck, F.; Nauwynck, H.; Boyen, F.; Pardon, B.; Theuns, S. High quality genome assemblies of Mycoplasma bovis using a taxon-specific Bonito basecaller for MinION and Flongle long-read nanopore sequencing. BMC Bioinform. 2020, 21, 1–16. [Google Scholar] [CrossRef]
  47. Tytgat, O.; Gansemans, Y.; Weymaere, J.; Rubben, K.; Deforce, D.; Van Nieuwerburgh, F. Nanopore Sequencing of a Forensic STR Multiplex Reveals Loci Suitable for Single-Contributor STR Profiling. Genes 2020, 11, 381. [Google Scholar] [CrossRef] [Green Version]
  48. Gabrieli, T.; Sharim, H.; Michaeli, Y.; Ebenstein, Y. Cas9-Assisted Targeting of CHromosome segments (CATCH) for targeted nanopore sequencing and optical genome mapping. BioRxiv 2017. Available online: https://www.biorxiv.org/content/10.1101/110163v3 (accessed on 30 November 2020).
  49. Madsen, E.B.; Höijer, I.; Kvist, T.; Ameur, A.; Mikkelsen, M.J. Xdrop: Targeted sequencing of long DNA molecules from low input samples using droplet sorting. Hum. Mutat. 2020, 41, 1671–1679. [Google Scholar] [CrossRef]
  50. Mantere, T.; Kersten, S.; Hoischen, A. Long-Read Sequencing Emerging in Medical Genetics. Front. Genet. 2019, 10, 426. [Google Scholar] [CrossRef] [Green Version]
  51. Zhao, H.; Sun, Z.; Wang, J.; Huang, H.; Kocher, J.-P.; Wang, L. CrossMap: A versatile tool for coordinate conversion between genome assemblies. Bioinform. 2014, 30, 1006–1007. [Google Scholar] [CrossRef]
  52. Van Der Auwera, G.A.; Carneiro, M.O.; Hartl, C.; Poplin, R.; Del Angel, G.; Levy-Moonshine, A.; Jordan, T.; Shakir, K.; Roazen, D.; Thibault, J.; et al. From FastQ Data to High-Confidence Variant Calls: The Genome Analysis Toolkit Best Practices Pipeline. Curr. Protoc. Bioinform. 2013, 43, 11.10.1–11.10.33. [Google Scholar] [CrossRef]
  53. De Coster, W.; D’Hert, S.; Schultz, D.T.; Cruts, M.; Van Broeckhoven, C. NanoPack: Visualizing and processing long-read sequencing data. Bioinform. 2018, 34, 2666–2669. [Google Scholar] [CrossRef]
  54. Quinlan, A.R.; Hall, I.M. BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics 2010, 26, 841–842. [Google Scholar] [CrossRef] [Green Version]
Ijms 21 09177 g001 550
Figure 1. Target region selected for this the study. The genomic region chr19:44,902,730–44,915,006 (hg38) visualized in the UCSC Genome Browser, identifying the single-nucleotide polymorphisms (SNVs) included in dbSNP. SNVs indicated with gray ticks represent synonymous variants in coding regions. SNVs indicated with red ticks represent non-synonymous variants in coding regions. SNVs indicated with blue ticks represent variants in untranslated regions. SNVs indicated with black ticks represent either variants in intronic regions or upstream/downstream gene variants. The positions of the PCR primers are shown in green.
Figure 1. Target region selected for this the study. The genomic region chr19:44,902,730–44,915,006 (hg38) visualized in the UCSC Genome Browser, identifying the single-nucleotide polymorphisms (SNVs) included in dbSNP. SNVs indicated with gray ticks represent synonymous variants in coding regions. SNVs indicated with red ticks represent non-synonymous variants in coding regions. SNVs indicated with blue ticks represent variants in untranslated regions. SNVs indicated with black ticks represent either variants in intronic regions or upstream/downstream gene variants. The positions of the PCR primers are shown in green.
Ijms 21 09177 g001
Ijms 21 09177 g002 550
Figure 2. Amplification of the APOE locus. (A) Agarose gel electrophoresis of amplicons generated from sample V1 genomic DNA using LongAmp Taq Polymerase (NEB) and commercial primers (Qiagen) or custom primers. The expected amplicon size is ~12 kb. (B) Agarose gel electrophoresis of amplicons generated from sample V1 genomic DNA or from PCR products (V1 amplicons) generated in (A), using nested primers designed to amplify a 412-bp product from APOE exon 3. (C) Agarose gel electrophoresis of long-range PCR amplicons generated using custom primers in the presence of 2.5%, 4%, or 5% DMSO. (D) Agarose-gel electrophoresis of amplicons generated using custom primers and 4% DMSO and further selected by gel excision and on-column purification (manual) or the BluePippin device. (E) Capillary electrophoresis of the same products from (D). Starting and recovered μg refer to the amount of amplicon before and after size selection and the relative yield. NTC = no-template control.
Figure 2. Amplification of the APOE locus. (A) Agarose gel electrophoresis of amplicons generated from sample V1 genomic DNA using LongAmp Taq Polymerase (NEB) and commercial primers (Qiagen) or custom primers. The expected amplicon size is ~12 kb. (B) Agarose gel electrophoresis of amplicons generated from sample V1 genomic DNA or from PCR products (V1 amplicons) generated in (A), using nested primers designed to amplify a 412-bp product from APOE exon 3. (C) Agarose gel electrophoresis of long-range PCR amplicons generated using custom primers in the presence of 2.5%, 4%, or 5% DMSO. (D) Agarose-gel electrophoresis of amplicons generated using custom primers and 4% DMSO and further selected by gel excision and on-column purification (manual) or the BluePippin device. (E) Capillary electrophoresis of the same products from (D). Starting and recovered μg refer to the amount of amplicon before and after size selection and the relative yield. NTC = no-template control.
Ijms 21 09177 g002
Ijms 21 09177 g003 550
Figure 3. Integrative Genomics Viewer (IGV) visualization of ONT mapped reads from sample V1 at position 11,658. IGV visualization showing a region where a homozygous SNV at position 11,658 was not identified based on ONT data. The blue box shows that the SNV occurs within a long homopolymer run.
Figure 3. Integrative Genomics Viewer (IGV) visualization of ONT mapped reads from sample V1 at position 11,658. IGV visualization showing a region where a homozygous SNV at position 11,658 was not identified based on ONT data. The blue box shows that the SNV occurs within a long homopolymer run.
Ijms 21 09177 g003
Ijms 21 09177 g004 550
Figure 4. Accuracy of SNV calling from ONT reads. The graph shows the harmonic mean of precision and recall for SNV calling as an F1 score averaged across samples, when different numbers of ONT reads (10, 30, 60, 100, 300, 600, 1000, or 10,000) are used. Error bars represent the standard error.
Figure 4. Accuracy of SNV calling from ONT reads. The graph shows the harmonic mean of precision and recall for SNV calling as an F1 score averaged across samples, when different numbers of ONT reads (10, 30, 60, 100, 300, 600, 1000, or 10,000) are used. Error bars represent the standard error.
Ijms 21 09177 g004
Ijms 21 09177 g005 550
Figure 5. Experimental pipeline to assess the performance of alignment and phasing software. (A) Sets ranging from 10 to 10,000 ONT reads were considered. (B) Variants in the region of interest were identified using the highly accurate Illumina platform and stored in a VCF file (ground-truth unphased variants). ONT reads were mapped to the reference sequence using either BWA or Minimap2, producing a BAM file. (C) VCF and BAM files constitute the input files for the phasing software (WhatsHap or HapCUT2), which generated a phased VCF file allowing haplotype reconstruction. The ONT-phased VCF file was then compared to a reference-phased VCF file (ground-truth phased variants). Analysis was carried out 100 times to obtain more robust estimates of accuracy.
Figure 5. Experimental pipeline to assess the performance of alignment and phasing software. (A) Sets ranging from 10 to 10,000 ONT reads were considered. (B) Variants in the region of interest were identified using the highly accurate Illumina platform and stored in a VCF file (ground-truth unphased variants). ONT reads were mapped to the reference sequence using either BWA or Minimap2, producing a BAM file. (C) VCF and BAM files constitute the input files for the phasing software (WhatsHap or HapCUT2), which generated a phased VCF file allowing haplotype reconstruction. The ONT-phased VCF file was then compared to a reference-phased VCF file (ground-truth phased variants). Analysis was carried out 100 times to obtain more robust estimates of accuracy.
Ijms 21 09177 g005
Ijms 21 09177 g006 550
Figure 6. Accuracy of full haplotype reconstruction. Haplotypes were called 100 times using all possible combinations of two alignment and phasing software and different read numbers (10, 30, 60, 100, 300, 600, 1000, and 10,000). The graph shows the frequency at which the fully correct haplotype was called (averaged across samples) with error bars representing the standard error.
Figure 6. Accuracy of full haplotype reconstruction. Haplotypes were called 100 times using all possible combinations of two alignment and phasing software and different read numbers (10, 30, 60, 100, 300, 600, 1000, and 10,000). The graph shows the frequency at which the fully correct haplotype was called (averaged across samples) with error bars representing the standard error.
Ijms 21 09177 g006
Table
Table 1. Genetic variants in the region of interest. The columns show the identification code in dbSNP (dbSNP ID), the hg38-based coordinate positions (Pos.), and the reference (Ref.) and alternative (Alt.) nucleotides of each variant.
Table 1. Genetic variants in the region of interest. The columns show the identification code in dbSNP (dbSNP ID), the hg38-based coordinate positions (Pos.), and the reference (Ref.) and alternative (Alt.) nucleotides of each variant.
dbSNP IDPos.Ref.Alt.
rs893292251chr19:44902764CT
rs34215622chr19:44903281CCG
-chr19:44903416GA
-chr19:44904398CT
rs7259620chr19:44904531GA,C
rs449647chr19:44905307AT
rs769446chr19:44905371TC
rs405509chr19:44905579TG
-chr19:44905832CT
rs440446chr19:44905910CG
-chr19:44906337GA
rs769450chr19:44907187GA
rs429358chr19:44908684TC
rs7412chr19:44908822CT
rs747519137chr19:44909521CTC
-chr19:44909967TGT
rs1065853chr19:44909976GC,T
rs75627662chr19:44910319CT
-chr19:44910393AC
-chr19:44910397TC
-chr19:44910405AC
-chr19:44910410AC
rs72654469chr19:44910678TC
rs72654473chr19:44911142CA
rs439401chr19:44911194TC
rs5828224chr19:44911609ATA
rs445925chr19:44912383GA,C
rs483082chr19:44912921GT
rs59325138chr19:44913034CT
rs584007chr19:44913221AG
rs438811chr19:44913484CT
rs390082chr19:44913574TA,C,G
rs72654445chr19:44913943GA
-chr19:44914318AT
rs34954997chr19:44914381CCTTCG
Table
Table 2. Number of total and phased ground-truth variants. For each sample, the table shows the total number of homozygous and heterozygous variants identified by Illumina data and the number of heterozygous variants phased by 10× Genomics or trio analysis. SNVs = single nucleotide variants, Indels = insertion/deletion variants.
Table 2. Number of total and phased ground-truth variants. For each sample, the table shows the total number of homozygous and heterozygous variants identified by Illumina data and the number of heterozygous variants phased by 10× Genomics or trio analysis. SNVs = single nucleotide variants, Indels = insertion/deletion variants.
Sample NameHomozygous VariantsHeterozygous VariantsHeterozygous Phased Variants
SNVsIndelsSNVsIndelsSNVsIndels
NA1287800101101
NA12892004040
V163111111
V2406261
V3414141
Table
Table 3. ONT sequencing statistics. The total number of ONT reads and the number and percentage of PASS reads, along with the read length of each category.
Table 3. ONT sequencing statistics. The total number of ONT reads and the number and percentage of PASS reads, along with the read length of each category.
Sample NameNumber of ReadsReads Mean LengthNumber of PASS Reads (%)PASS Reads Mean Length (bp)Number of PASS Reads AlignedNumber of PASS Reads Covering the Whole Region
NA1287828,25111,87124,335 (86%)12,04624,330 (99.99%)23,799 (97.82%)
NA1289227,41011,71023,541 (86%)12,04623,541 (100%)23,022 (97.80%)
V139,99311,55233,489 (84%)12,04233,486 (99.99%)32,941 (98.37%)
V238,50111,69033,373 (87%)12,04433,371 (99.99%)32,755 (98.15%)
V332,08411,58827,165 (85%)12,04627,162 (99.99%)26,673 (98.20%)
Table
Table 4. Number of variants identified by the ONT platform. The total number of correct and spurious variants identified in the samples based solely on ONT data. SNVs = single nucleotide variants, indels = insertion/deletion variants, TP = true positives (identified by ONT and Illumina), FN = false negatives (only Illumina), FP = false positives (only ONT).
Table 4. Number of variants identified by the ONT platform. The total number of correct and spurious variants identified in the samples based solely on ONT data. SNVs = single nucleotide variants, indels = insertion/deletion variants, TP = true positives (identified by ONT and Illumina), FN = false negatives (only Illumina), FP = false positives (only ONT).
NA12878NA12892V1V2V3
SNVsIndelsSNVsIndelsSNVsIndelsSNVsIndelsSNVsIndels
Number of ReadsTPFNFPTPFNFPTPFNFPTPFNFPTPFNFPTPFNFPTPFNFPTPFNFPTPFNFPTPFNFP
10100901140300016200339110022623023
3010000114000031620122910024800022
60100001240000317100321000022800112
100100001240000217100321000022800022
300100001140000017100321000020800020
60010000104000001710032910020800020
100010000104000001620031820020800021
100008200114000011080031730023710022
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Maestri, S.; Maturo, M.G.; Cosentino, E.; Marcolungo, L.; Iadarola, B.; Fortunati, E.; Rossato, M.; Delledonne, M. A Long-Read Sequencing Approach for Direct Haplotype Phasing in Clinical Settings. Int. J. Mol. Sci. 2020, 21, 9177. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21239177

AMA Style

Maestri S, Maturo MG, Cosentino E, Marcolungo L, Iadarola B, Fortunati E, Rossato M, Delledonne M. A Long-Read Sequencing Approach for Direct Haplotype Phasing in Clinical Settings. International Journal of Molecular Sciences. 2020; 21(23):9177. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21239177

Chicago/Turabian Style

Maestri, Simone, Maria Giovanna Maturo, Emanuela Cosentino, Luca Marcolungo, Barbara Iadarola, Elisabetta Fortunati, Marzia Rossato, and Massimo Delledonne. 2020. "A Long-Read Sequencing Approach for Direct Haplotype Phasing in Clinical Settings" International Journal of Molecular Sciences 21, no. 23: 9177. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21239177

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop