Scalable noninvasive amplicon-based precision sequencing (SNAPseq) for genetic diagnosis and screening of β-thalassemia and sickle cell disease using a next-generation sequencing platform

Introduction

β-thalassemia (BT) and sickle cell disease (SCD) are mendelian monogenic red blood cell disorders and are among the most prevalent disorders worldwide. They are caused due to genetic mutations in the β-globin (HBB) gene and have an estimated incidence of 40-82 per 1,000 live births ^{1, 2}. β-hemoglobinopathies are a prevalent cause of health problems across multiple groups of ethnicities. SCD results from a single amino acid change at the 6^th position (GAG to GTG) in the β-globin gene to a replacement of glutamic acid to valine, whereas there are more than 200 different mutations reported for BT. These are majorly point mutations or deletions present across the HBB gene. An estimated 7% of the world’s population carries a defective β-globin gene, and around 300,000 babies are born each year with severe hemoglobin abnormalities ³. Recent estimates suggest that this number could rise to more than 400,000 by the year 2050 ⁴. Treatment options for β-hemoglobinopathies are primarily interventional rather than providing a permanent cure. Hydroxyurea, pain management and intravenous hydration are used in SCD whereas frequent blood transfusion coupled with iron chelation therapy are the treatment options for BT. Currently, allogeneic hematopoietic stem cell transplantation (HSCT) is the only curative option available for both the above diseases but success is primarily dependent on the availability of human leukocyte antigen (HLA)-matched donor ⁵.

Although advancements in disease management and care have been made, there is still a high incidence rate, which has a profound impact on the public health system demanding significant investment in patient care. Therefore, early detection of the disease remains the realistic approach for disease prevention, decreasing the burden on the healthcare system and aiding better disease management with the ultimate goal of reducing patient mortality and morbidity rates. At present, hemoglobin electrophoresis, high-performance liquid chromatography (HPLC) and isoelectric focusing are the predominantly used diagnostic methods. In general, HPLC tests require further confirmation through Sanger sequencing-based DNA analysis^{6, 7}. Currently, most commercially available diagnostic methods rely on antibody-based recognition of Hb variants, however, these tests may provide deceptive results in blood-transfused individuals. Moreover, commonly used antibody tests in newborn screening, may lead to inconclusive results as they often undergo blood transfusion due to multiple birth-related complications ⁸. Hence, subsequent validation after 90 days of transfusion is necessary for such cases ^{9, 10}. Additionally, the basal level of fetal hemoglobin expression is high in the first few months of the newborn thus limiting the applicability of antibody-based diagnostic tests. Despite advancement, the genetic diagnosis of BT continues to be challenging due to the prevalence of mutations and deletion throughout the HBB gene. Protein-based diagnostic tools cannot detect BT deletions, where 619 bp deletion is common in Southeast Asian population ¹¹. Hence, a reliable diagnostic method is required that can differentiate the genotype of β-hemoglobinopathies efficiently. Besides early detection of patients, carrier detection will unveil the scope of genetic counseling in reducing the frequency of the disease. World Health Organization (WHO) indexed hemoglobin testing as one of the most crucial in vitro diagnostics (IVD) tests for primary care use in developing and underdeveloped countries¹²

The classification of genetic diseases for accurate diagnosis can be provided by allele-specific detection of disease mutation. Molecular approaches to identify point mutations are mainly dependent on methods such as gold standard traditional Sanger sequencing ¹³ and real-time PCR ¹⁴ which are very low throughput and need multiple amplicons and several sequencing runs with different primers pairs in order to cover the entire HBB gene and to identify pathogenic variants. In addition, the read quality of Sanger sequencing is often not good in the first 50-80 bases where the primers bind to start amplification. Further, different combinations of compound heterozygous mutations in different ethnic groups demand the need for an integrated diagnosis system that can identify most of the mutations associated with the disease. Moreover, it is challenging to develop a molecular diagnostic test that does not necessitate nucleic acid isolation and purification because current genotyping methods based on DNA amplification use relatively high DNA purity. Thus, establishing a molecular test that analyzes the samples without DNA extraction and purification has remained elusive. The goal of the current study is to create a low-cost, low-complexity buffer system using common laboratory chemicals that enables direct PCR amplification using a straightforward process that effectively releases DNA from buccal swabs or finger prick blood.

In the last several years, next-generation sequencing (NGS) has become a cost-effective tool for the high-throughput identification of genetic variants, thereby unveiling new opportunities in the field of molecular diagnosis¹⁵. When compared to Sanger sequencing, NGS-based amplicon sequencing has several advantages such as higher sequencing depth with enhanced sensitivity, enhanced discovery power and higher mutation resolution, high throughput and maximum data that can be generated with the same quantity of DNA. The diagnosis of aneuploidy in embryo biopsy specimens is made possible by this technique, which has previously been neglected for use in assisted reproductive technology^16–18. Recently, single gene mutations in embryos have also been diagnosed using NGS; nevertheless, the clinically employed techniques still call for extensive customization in the design and optimization of multiplex-PCR for the amplification of mutant sites and polymorphisms ^{14, 19, 20}. Further, Hallem et al. and Haque et al. have validated that NGS in patient-based and population-based carrier screening of inherited recessive disorders and demonstrated acceptable rates of false-positive and false-negative as well as cost-effectiveness of the tool ^{21, 22}. Besides, NGS has higher sequencing depth, which enables enhanced sensitivity and high mutation resolution ²³, which is not possible in other methods.

In the current pilot study, we demonstrated a unique strategy based on the NGS approach, for the detection of virtually all HBB mutations accountable for SCD and BT. A simple non-invasive buccal swab specimen or finger prick blood has been adopted, which is economical and can be used as a direct sample matrix avoiding any DNA extraction process. For a proof-of-concept, the robustness of the SNAPseq assay in the detection of allele-specific BT and SCD genotypes from extraction-free non-invasive crude lysates of buccal swab samples with no additional PCR cleanup has been shown with patients, carriers and wild-type samples. A systematic SNAPseq data analysis pipeline was established and validated to prioritize and predict the pathogenicity of all HBB genetic mutations identified in each specimen. The present study shows that a simplified sampling procedure combined with NGS has enormous potential and clinical utility in the molecular diagnosis of genetically heterogeneous diseases such as BT. Our study demonstrated that the NGS-based SNAPseq assay, compared to conventional sequencing approaches, can identify virtually all β-globin mutations, including compound heterozygous and large deletions in a precise manner. These results suggest that SNAPseq is a precise and most efficient platform for large-scale genetic screening and clinical genotyping in subjects with SCD and BT.

2. Materials and Methods

3. Results

3.1 SNAPseq strategy

The design strategy adopted for the development of scalable noninvasive amplicon-based precision sequencing (SNAPseq) for genetic diagnosis and screening of β-hemoglobinopathies is shown in Fig.1. Clinical samples were collected in the forms of blood or buccal swab or finger prick. Genomic DNA was isolated from the whole blood and the other two sets of samples (buccal swab or finger prick) were treated with a lysis buffer to release the nucleic acids. PCR amplification was performed to amplify the entire HBB gene as a single amplicon. The PCR products were directly subjected to NGS library preparation as previously described. All the samples were processed in a 96-well plate. Subsequently, 96 samples were pooled together in a single tube. The pooled libraries were denatured and neutralized and a sequencing run was performed in either MiSeq or MiSeq or iSeq platform. The raw data generated in binary base call (BCL) format was analyzed using in-house developed SANPseq pipeline (Fig.1). In this study, the entire bioinformatics pipeline from data collection to clinical reporting is completely automated to reduce human intervention.

• Download figure
• Open in new tab

Figure 1:

Workflow of the SNAPseq methodology

3.2 Standardization of primers and polymerases

The main objective of the study was to develop a high-throughput integrative assay system that can detect virtually all the mutations associated with SCD and BT in a single amplicon amplified from the HBB gene. In order to enhance the versatility of the assay, designing a primer combination that could efficiently amplify using a wider range of proofreading enzymes was required. To accomplish this, Multiple primer combinations spanning the HBB gene were designed. All primer combinations were validated for amplification using multiple proofreading enzymes including NEBNext, Phusion Polymerase, Amplitaq GOLD, Q5 polymerase, PrimeStar max polymerase, LA Taq polymerase (Supplementary Fig. 1. Our results demonstrated a consistent amplification of the HBB gene using the HBB-F1 and HBB-R3 primer combination with all the above-mentioned polymerases. The primer combination HBB-F1 and HBB-R3 covers the entire HBB gene including promoter, 5’ and 3’ untranslated regions, all three exons and both the introns of the gene (Fig. 2A, 2B). Henceforth, PrimeStar polymerase was used in the study due to its fast extension rate and high fidelity. We further investigated the reliability of NGS outcomes obtained from direct amplicons or column purified amplicons. The elimination of this step would reduce human intervention, cost and a much wider scope of automation. Hence, PCR amplification was performed using genomic DNA obtained from 8 volunteers. Our data indicate a similar median coverage 12521.45 and 12226.15 in direct amplicons and purified amplicons. Hence, these findings indicate the feasibility of direct amplicons for NGS (Fig. 2C).

• Download figure
• Open in new tab

Figure 2:

A) Schematics of the HBB gene highlighting all types of BT mutation spread across the gene. B) Tracks of the amplicon in integrative genome viewer demonstrating coverage of complete HBB gene. C) Violin plot demonstrating the mean coverage of PCR purified and direct PCR samples, n=10. D) Violin plot demonstrating the mean coverage of direct PCR samples amplified from blood lysate, buccal swab lysate and purified genomic DNA, n=16.

3.4 Limit of detection

Next, we investigated the limit of detection for all three template types (purified genomic DNA, buccal swab lysate and blood lysate). In order to assess the limit of detention in purified DNA samples, a range of 10⁸-10¹ copies/ul dilution of pHBB plasmid was used. The results suggested that polymerase is efficient enough to amplify from as low as 10 target DNA copies (Fig.3A, Supplementary Fig 3A). Direct lysates contain genomic DNA in crude form; hence it becomes challenging to calculate the copy number accurately. To address this issue, we conducted additional tests to determine the maximum percentage of direct lysate that could potentially hinder PCR reaction due to the presence of salts. Additionally, we also tested the minimum amount of lysate required for amplification. To achieve this, we performed PCR using a range of 30%-0.01% of direct lysate as a template obtained from 3 volunteers. Our results indicate that as high as 30% of direct blood lysate was not inhibitory for PCR amplification (Fig. 3B, Supplementary 3B). On the other hand, 20% of buccal swab lysates in PCR reaction was inhibitory, possibly due to the presence of fibers from buccal swab that may have been left behind during processing (Fig. 3C, Supplementary Fig. 3C). Notably, a much higher variation in amplification was observed in buccal swab lysate as compared to blood lysate due to individual collection variability. Hence, the amount of genomic DNA may vary from individual to individual due to sampling variability. Nevertheless, our data demonstrate that as low as 0.1% of both the lysates was sufficient for amplification (Fig 3B and 3C, Supplementary Fig. 3D and 3E). Overall, 1%-15% of direct lysates are sufficient to amplify enough from both lysates.

• Download figure
• Open in new tab

Figure 3:

A) Minimum copy number required for the amplification using purified plasmid DNA. Densitometry analysis was done using ImageJ software and data from 10¹¹ dilution was considered as 100%. Remaining dilutions have been plotted as percentage relative to 100% (10¹¹ dilution), n=3. B) Box plot demonstrating minimum and maximum percentage of blood lysate required for amplification. Densitometry analysis was done using ImageJ software and the highest data value was considered as 100%. Remaining samples have been plotted as percentage relative to 100% (10¹¹ dilution), n=3. C) Box plot demonstrating minimum and maximum percentage of buccal swab lysate required for amplification. Densitometry analysis was done using ImageJ software and the highest data value was considered as 100%. Remaining samples were plotted as percentage relative to 100% (10¹¹ dilution), n=3.

3.5 Effect of time and temperature on crude genomic DNA lysates

As buccal swab lysates are non-invasive and can serve as a source of crude genomic DNA, we further evaluated the stability of unprocessed buccal swabs samples under variable temperatures. We have demonstrated previously that storing buccal swab samples at different temperatures (25°C, 37°C, and 42°C) prior to processing results in consistent amplification for short amplicons²⁴. We adopted the same methodology to assess its impact on the amplification of large amplicons. The buccal swab samples were collected and stored at all three temperatures for 5 days followed by further processing and PCR amplification. Our results indicated consistent amplification with buccal swab sample lysates at all different temperatures (Supplementary Fig. 4A) (Table 2). To further extend the scope of sample availability, we tested if dry spots of blood can also be used as a genomic source. In a similar manner to buccal swab samples, dry spots of blood samples were stored at 25°C, 37°C, and 42°C for 5 days. Afterwards, lysates were prepared and our amplification data demonstrated consistent amplification with all the samples (Supplementary Fig.4B) (Table 2). In summary, our findings indicate that both buccal swab and blood (dry spot) lysate results in consistent amplification for NGS. We further assessed the sensitivity of sequencing by employing three different Illumina sequencing platforms: MiSeq, iSeq and MiniSeq. Our data demonstrated 2848.78, 12010.85 and 2690.19 median coverage for MiSeq, iSeq and MiniSeq respectively (Fig. 4).

• Download figure
• Open in new tab

Figure 4:

Violin plot demonstrating the mean coverage from MiSeq and iSeq sequencing platform.

View this table:

• View inline
• View popup
• Download powerpoint

Table 2:

A comparison of PCR amplification using direct lysates stored at different temperature for 5 days

3.6 Comprehensive and precise identification of genetic mutations associated with β-hemoglobinopathies

Next, to validate the accuracy of our pipeline for precisely identifying the genetic mutation associated with β-hemoglobinopathies (Fig. 5), we collected samples in the form of blood (for genomic DNA isolation), or buccal swab or finger prick blood. We ensured coverage of BT mutation present in the different regions of the gene as well as the compound heterozygotes including Hb S/BT, BT/BT. We amplified the target region for over 100 samples with known mutations and found that the results generated by our pipeline were in concordance with the known genotype of the sample. This indicates 100% sensitivity and 100% specificity of the SNAPseq methodology (Table 3). We further sought to test the genotype of patients clinically declared as thalassemia based on their clinical profile. We were able to identify the mutation in the patients. Further, our data indicates out of all BT patients, 25% patients were compound heterozygous for the mutations. All the different mutations from known and unknown samples identified using SNAPseq methodology for this study are listed in Table 4.

• Download figure
• Open in new tab

Figure 5:

Flowchart demonstrating the pipeline for the identification of SCD and BT mutations. Sequencing data was analyzed using the inhouse pipeline. Fastq files were trimmed using the trimming algorithm and the resultant fastq files were aligned to the human reference genome. Subsequently, PCR duplicates were removed and the variants were called using Varscan and annotated to RefSeq using Annovar. Finally, the variants are matched against the clinvar likely pathogenic/ pathogenic list for SCD and BT as output from this pipeline.

View this table:

• View inline
• View popup
• Download powerpoint

Table 3:

Sensitivity of SNAP-seq metholod

View this table:

• View inline
• View popup

Table 4:

List of mutation identified in the study

In order to evaluate the robustness and reproducibility of the SNAPseq assay, we assessed the above protocol at a secondary health care center, where sample collection, sample processing, PCR amplification, NGS library preparation and sequencing run was performed by the laboratory technicians and good quality data was obtained. This clearly demonstrates the reproducibility and robustness of the assay and how easily it can be adopted in a healthcare set-up for high-throughput screening.

Discussion

Hemoglobinopathies are the most common monogenic disorder and impose a significant global health burden on families and the health care system across the globe. Hence the early diagnosis of hemoglobinopathies is imperative for both prevention and appropriate treatment for patients. It is also crucial to provide counseling to couples and families who may be at risk to suffer severe hematological symptoms. Despite extensive studies at biochemical, molecular and hematological levels, the differential diagnosis of hemoglobinopathies remains challenging.

At present, the diagnosis of β-hemoglobin disorders is an extensive process that involves multiple stages of screening: evaluation of the blood panel and specialized hemoglobin tests for characterization. The majority of newborn screening (NBS) programs utilize protein-based hemoglobin separation methods such as isoelectric focusing (IEF), HPLC, and gel or liquid-based electrophoresis ²⁵. Though protein-based techniques are the basis of hemoglobin diagnostics, they may not be able to identify BT mutations or to distinguish between HbSS and compound heterozygosity for HbS and hereditary persistence of fetal hemoglobin (HbS/HPFH) ²⁶. Hence, additional confirmation of mutation is required by DNA sequencing. A low MCH (Mean corpuscular hemoglobin) and MCV (Mean corpuscular volume) are common symptoms in hemoglobinopathies, but these symptoms may often be misdiagnosed for vitamin B12, iron and folic acid deficiency ²⁷. If β-hemoglobinopathy is suspected with the above symptoms, a series of specialized hemoglobin tests are recommended. Thus, multiple stages of diagnoses may lead to misdiagnosis, delayed treatment and missing out carriers of mutations who require counseling ²⁶. In addition, the diagnostic differentiation is now time-sensitive and crucial due to the rising use of early, pre-symptomatic disease-modifying medication (e.g., initiation of hydroxyurea therapy at 6 months of age)²⁸.

The rising prevalence of compound heterozygotes of BT and SCD possesses an additional challenge for point-of care or simple molecular based genetic tests. Similar to variable phenotypes associated with BT, a compound heterozygote of BT and SCD results in variable clinical presentation depending on the mutation. Around 10% to 15% of sickle cell disease is caused by a combination of the sickle cell and BT mutations, also known as compound heterozygosity²⁹, which is not possible to detect in point-of-care screening methods. The distinction between HbSS and HbS/β genotypes using current hemoglobin-based diagnostic methods can be difficult ^{30, 31}. As a result, there is a need for a scalable and comprehensive genetic diagnostic method that can identify virtually all HBB mutations associated with β-hemoglobinopathies.

NGS has transformed genomics research and created better sequencing methods. Although Sanger sequencing is competent for sequencing amplicons up to 1 kb ³², screening of longer templates needs primer walking or shotgun sequencing, which are more complicated and expensive ^{33, 34}. The ability of SNAPseq to evaluate long amplicons is a significant advantage in such cases. SNAPseq also has a notable advantage over Sanger sequencing as they frequently provided full coverage for the target amplicon. Generation of entire coverage through Sanger sequencing necessitates bidirectional analysis in order to avoid the interpretational complications brought on by “dye blobs” and ambiguity in base calling at the start of each sequence, which can vary from 60 to 120 bases of the sequence³⁵. Though Sanger sequencing can sequence up to 1 kb, the quality of the reads from 1 kb would be only around 700-800 bp due to the above reason.

NGS based molecular diagnosis has been successfully applied for various genomic applications ³⁶. Recently, Kubikova et. Al, have used an NGS based approach for the diagnosis of β-hemoglobinopathies in preimplantation embryos ³⁷. In this study, 8 pairs of primers were used to amplify 12 different PCR fragments spanning the coding regions of the HBB gene and splice acceptor and donor sites using purified genomic DNA extracted from participant samples. There are many BT mutations which are present in the intronic regions that are not completely covered in the study and amplification of multiple fragments with different primer sets using purified genomic DNA is an expensive, time consuming and cumbersome process, which are the major limitations of the above study. Another study by He et al reported a thalassemia carrier screening using NGS based approach but the study was designed to identify only the most common disease-causing point mutations or selected region ³⁸. SNAPseq, an NGS-based approach is designed to identify both common and rare, annotated and novel variants in carriers with and without BT and sickle cell trait phenotypes, significantly improving the detection of carrier status and subsequently improving the detection rate of at-risk couples. This is in contrast with conventional carrier screening assays, which are designed to search for only the most prevalent mutations within a gene and coverage of less comprehensive due to various factors²¹.

Advancements in the field of NGS and its high-depth sequencing to identify mutations and variants have fast-tracked the diagnosis of many human genetic disorders and have been more reliable in characterizing the disease genotype than other molecular diagnostic tests available ²². Significant progress in this field in recent years has found numerous applications in research and genomic medicine. However, an extraction-free, universal single amplicon sequencing strategy for identifying virtually all HBB mutations has not been established. In this pilot study, we designed, optimized and furnished a SNAPseq methodology by utilizing the advantage of inexpensive tools, to identify all the associated mutations with SCD and BT in the HBB gene, using crude lysates of buccal swabs and/or blood and analyzed using an established pipeline to prioritize pathogenic mutations with allele-specific sensitivity. In addition, in this study, we evaluated the ability of two major Illumina sequencing platforms, MiSeq, MiniSeq and iSeq, which covers both high throughput and low throughput samples and obtained similar results. Our optimized assay pipeline system provides clinically relevant sensitivity in mutation calling across an HBB gene.

The current study focuses on developing a single amplicon-based strategy for screening of β-hemoglobinopathies. During the development of this assay system, one of the main challenges was to address the method of sample collection and processing, ensuring that the samples were collected in the quickest and most efficient way without compromising the accuracy and sensitivity of the test. Establishing a genetic test that does not require genomic DNA isolation and purification is challenging since existing molecular tests based on DNA amplification require relatively high purity of genomic DNA ²⁴. Hence, developing a molecular test that analyzes the samples without genomic DNA isolation has remained evasive. Here, we sought to establish a cost-effective and low-complexity buffer system with common laboratory reagents that make direct genomic DNA amplification feasible by employing a simple protocol that enables the release of DNA most effectively. We have also compared our extraction-free protocol with PCR products amplified using purified genomic DNA of the same samples and obtained a similar outcome. In addition, the SNAPseq assay system has the additional advantage of no additional PCR cleanup steps required for the downstream applications (NGS) so that it reduces the overall assay cost. These cost reductions will have more benefits from the above steps in high throughput settings than from further advances in sequencing capacity.

Blood samples have been traditionally used for molecular studies and clinical diagnostics. Nevertheless, blood collection for genomic DNA isolation is an invasive procedure and demands skilled healthcare professionals³⁹. Hence, in this study, we have demonstrated a crude lysate prepared from a simple finger prick blood, which does not necessitate professional assistance, results in consistent amplification. Furthermore, saliva is a more suitable and minimally intrusive alternative to blood for molecular diagnostics. However, some study participants have trouble spitting into a tube, while others have trouble producing enough saliva for genetic testing ⁴⁰. In addition, saliva collection can be challenging due to the common medication side effect of dry mouth. However, buccal swab sample collection is rapid, non-invasive, and cost-effective, taking all this into account we indicated the feasibility of genotyping SCD and BT using buccal swab samples in this study. Additionally, it would make it simple for patients and families to collect their samples at their homes or clinics.

The maintenance of sample integrity is very crucial, especially in remote areas. Ideally, biological samples should be stored in temperature controlled and properly qualified storage and should be processed as quickly as possible after collection⁴¹. Hence, sample storage is another obstacle during sample collection since most of the collected samples need to be stored and transported before processing in a realistic scenario. Consequently, we have demonstrated that buccal swab samples in the collection buffer are stable at 25°C, 37°C, and 42°C for at least up to 5 days. In addition, our study also demonstrated that finger prick blood samples are also stable for at least up to 5 days as dry spots. Our study highlights the crucial significance of our findings for low-and middle-income nations with a high prevalence of β-hemoglobinopathies. These countries often face limited access to expensive refrigerated transport and adequate storage facilities for biological samples.

In this study, we provide a general workflow for the interpretation of the enormous amount of data produced through NGS-based molecular diagnosis. Using a SNAPseq approach, we analyzed over 250 samples with 100% sensitivity using a single amplification covering virtually all the HBB mutations with a simple, cost-effective and extraction-free platform. In addition to SCD and BT mutations, SNAseq can identify other genotypes such as hemoglobin type C genotypes AC, CC (HBC) and SC (compound heterozygous for HbC and HbS) and hemoglobin E (HBE) β-thalassemia mutation (point mutation in codon 26 GAG to AAG) on first allele and any other BT mutation in the second allele. Here, we have outlined a highly sensitive NGS-based carrier screening system that has a number of advantages over conventional genotyping techniques and exhibits excellent levels of precision, reproducibility, and resilience for clinical usage. The immense potential and clinical use of NGS-based SNAPseq is emphasized by a high diagnostic yield attained by the identification of genetic mutations of all forms of inheritance, particularly de novo mutations.

A similar diagnostic workflow can also be easily adopted for screening other genetic diseases with very less customization. It can offer a quick, simple, and potentially affordable high-throughput approach for screening genetic disease to detect a wide range of genotypes including point mutations, large deletions and indels. It provides a great opportunity for better assessment, disease management and genetic counseling.

In summary, we developed a highly sensitive comprehensive and precise molecular diagnostic tool SNAPseq with the combination of minimal work-up, high-throughput large-scale screening with systematic analysis and interpretation of all identified genetic mutations for β-hemoglobinopathies, a most common monogenic disease. It might be appropriate for large-scale screening where the prevalence of β-hemoglobinopathies is high.