- Technical advance
- Open Access
Rapid detection of pathological mutations and deletions of the haemoglobin beta gene (HBB) by High Resolution Melting (HRM) analysis and Gene Ratio Analysis Copy Enumeration PCR (GRACE-PCR)
BMC Medical Genetics volume 17, Article number: 75 (2016)
Inherited disorders of haemoglobin are the world’s most common genetic diseases, resulting in significant morbidity and mortality. The large number of mutations associated with the haemoglobin beta gene (HBB) makes gene scanning by High Resolution Melting (HRM) PCR an attractive diagnostic approach. However, existing HRM-PCR assays are not able to detect all common point mutations and have only a very limited ability to detect larger gene rearrangements. The aim of the current study was to develop a HBB assay, which can be used as a screening test in highly heterogeneous populations, for detection of both point mutations and larger gene rearrangements.
The assay is based on a combination of conventional HRM-PCR and a novel Gene Ratio Analysis Copy Enumeration (GRACE) PCR method. HRM-PCR was extensively optimised, which included the use of an unlabelled probe and incorporation of universal bases into primers to prevent interference from common non-pathological polymorphisms. GRACE-PCR was employed to determine HBB gene copy numbers relative to a reference gene using melt curve analysis to detect rearrangements in the HBB gene. The performance of the assay was evaluated by analysing 410 samples.
A total of 44 distinct pathological genotypes were detected. In comparison with reference methods, the assay has a sensitivity of 100 % and a specificity of 98 %.
We have developed an assay that detects both point mutations and larger rearrangements of the HBB gene. This assay is quick, sensitive, specific and cost effective making it suitable as an initial screening test that can be used for highly heterogeneous cohorts.
According to the World Health Organization (WHO), inherited disorders of haemoglobin are the most common monogenic disorders in the world. The WHO estimates that 300,000 to 400,000 severely affected infants are born each year and that greater than 95 % of these severe cases arise from genetic abnormalities of the β-globin gene (HBB) . Globally, these conditions result in significant childhood morbidity and mortality.
The HBB gene is 1606 base pairs long, contains three exons and is located on the short arm of chromosome 11 towards the 3’ end of a 70 kb region known as the β-globin locus. The HBB gene codes for β-globin, a protein that is an essential component of normal adult haemoglobin . To date, 868 mutations of the HBB gene have been described that have either a qualitative or quantitative effect on β-globin synthesis. These include point mutations, small insertions / deletions (indels) and larger gene rearrangements . Qualitative defects are normally the result of point mutations and lead to the formation of abnormal β-globin chains, which in turn give rise to variant haemoglobins. Quantitative defects of the HBB gene either reduce, or totally stop, the synthesis of β-globin resulting in β-thalassaemia. Unlike α-thalassaemia, which is mostly caused by large deletions, the majority of β-thalassaemia alleles are the result of point mutations or small indels. However, in recent years several novel deletions of the HBB gene have been identified and it has become apparent that deletional forms of β-thalassaemia are more prevalent than was previously thought. Indeed, it has been suggested that large deletions may account for as much as 10 % of β-thalassaemia mutations .
Premarital screening programs specifically targeting β-thalassaemia were introduced in Greece  and Cyprus [6, 7] in the 1970s. These programs significantly reduced the number of children being born with β-thalassaemia, and as a result similar programs were subsequently introduced in a number of countries where the disease is prevalent . The success of such premarital screening programs is dependent on reliable laboratory tests for the detection of carriers.
In general, carriers for the common variant haemoglobins are readily detected by routine biochemical techniques, such as cation exchange High Performance Liquid Chromatography (HPLC) or capillary electrophoresis (CE). However, reliable detection of β-thalassaemia carriers can be more challenging . Haemoglobin A2 (Hb A2) quantification is widely used for screening, since Hb A2 is easily quantified by HPLC or CE and shows a modest increase in the majority of individuals with β-thalassaemia trait . However, some β-thalassaemia mutations result in only minimal increases in Hb A2 levels and are easily missed during screening . Interpretation of β-thalassaemia screening results can be further complicated by external factors that may influence the Hb A2 level. For example, Hb A2 is lowered by both iron deficiency  and α-thalassaemia , but is elevated by mutations in the Krüppel-like Factor 1 (KLF1) gene . In addition, co-inheritance a of β-thalassaemia trait with a δ-globin gene (HBD) mutation can lead to Hb A2 levels that are normal or even reduced . Consequently, although Hb A2 quantification remains a useful screening tool, confirmation by molecular analysis is frequently required .
A variety of PCR techniques, such as Amplification Refractory Mutation System (ARMS) PCR, Allele Specific Oligonucleotide (ASO) PCR and real time PCR have been used to detect β-thalassaemia resulting from point mutations or indels [16–18], while Gap-PCR is normally used to detect larger gene rearrangements . These assays are designed to detect either a single specific mutation or a small set of mutations. This approach works well in populations where a limited number of β-thalassaemia mutations are prevalent, but is less useful when testing highly heterogeneous populations where the range of mutations is less predictable. In such situations gene scanning methods may be useful and assays using various technologies including Denaturing HPLC (DHPLC), Single Strand Conformational Polymorphism (SSCP) and Denaturing Gradient Gel Electrophoresis (DGGE) have been described . These assays offer an alternative means of detecting point mutations or indels of the HBB gene, but are not appropriate for the detection of larger deletions and have the disadvantage of being open tube techniques that require significant post amplification processing.
An alternative gene scanning technology called High Resolution Melting PCR (HRM-PCR) offers potential advantages over the other techniques since it does not require complicated post amplification processing by chromatography or polyacrylamide gel electrophoresis . HRM-PCR is a closed tube technique that allows both amplification and detection to be performed on the same instrument without requiring further hands-on time after the initial PCR setup . Thus, HRM-PCR is well suited to clinical applications, since it permits a simple and efficient workflow, with minimal risk of amplicon contamination. In addition to being quicker and easier to perform, HRM-PCR has specificity and sensitivity equal to, or better than, other gene scanning technologies such as DHPLC [21, 22]. These advantages have resulted in the application of HRM-PCR in a wide range of clinical settings [23–28].
HRM-PCR assays for scanning the HBB gene mutations have been described, targeting the most common mutations found in Chinese  and South East Asian [30, 31] populations. These assays are well suited for screening the specific cohorts targeted. However, the approaches used to overcome interference from high frequency single nucleotide polymorphisms (SNPs) made these assays insensitive to some important pathological mutations that are common in other populations including Hb S and Hb C. Additionally, with the exception of the 3.4 kb deletion, these assays cannot detect larger gene rearrangements .
We have recently developed a Gene Ratio Copy Enumeration PCR (GRACE-PCR) to determine copy numbers for the α-globin genes  and we describe a similar approach for the detection of deletions of the HBB gene. We also developed a HRM-PCR method that was optimized with the use of an unlabelled probe and primers containing universal bases to address the issue of SNP interference. The combined use of GRACE-PCR and HRM-PCR allowed for the detection of point mutations, indels and larger deletions of the HBB gene, thus resulting in an assay that is universally applicable.
Materials and methods
Samples and cohort selection
This study was conducted using anonymous, archived material from blood specimens submitted to the Sheikh Khalifa Medical City (SKMC) laboratory, Abu Dhabi, United Arab Emirates (UAE) for thalassaemia/haemoglobinopathy screening. Ethical clearance was obtained from the SKMC Institutional Research and Ethics Committee to use this material for the current study. Since the study used anonymous archived material it was not possible to identify the subjects, therefore specific informed consent was not possible. The material used contained no patient identifiers, but the results of the haemoglobin analysis by HPLC had been assigned prior to the anonymization of samples. A total of 410 samples were used in the study, these included 305 with elevated Hb A2 (Hb A2 > 3.5 %), 27 with low or normal Hb A2, 76 with abnormal haemoglobins and two samples known to be positive for the 619 base pair deletion (NM_000518.4:c.316-149_*342delinsAAGTAGA) . The initial assessment of the assay was conducted using 342 of these samples and an additional 68 were used for a blinded validation.
Genomic DNA was extracted from whole blood using the QIAamp DNA Blood Mini kit (Qiagen, Germany) in accordance with the manufacturer’s instructions. The concentration of the extracted DNA was measured using a NanoDrop 2000 spectrophotometer (Thermo Scientific, USA) and adjusted to 10 ng/μl with 10 mM Tris, 0.5 mM EDTA (pH 9.0).
Primer design and synthesis
Analysis of the HbVar online database (http://globin.bx.psu.edu/hbvar/menu.html) was performed to identify the regions of the HBB gene of potential clinical relevance and to identify the breakpoints for known gene deletions. The final version of the assay used a total of 13 PCR reactions, 11 HRM-PCRs to scan for point mutations and indels of the HBB gene and two GRACE-PCRs for the detection of larger gene rearrangements (Table 1 and Fig. 1).
The previously described primer set H7  was included to test for the HBB:c.316-197C > T mutation found in intervening sequence two (IVS2). The reverse primers in primer sets H6 and H7 incorporated the universal base inosine to eliminate the risk of potential interference from common non-pathological SNPs (Table 1). In addition, an unlabelled probe was designed to overlap the region of the first exon containing the HBB:c.19G > A and HBB:c.20A > T mutations (Fig. 1 and Table 1). This probe was designed to have a melting temperature approximately 8 °C lower than the corresponding amplicon and was synthesized with a phosphate block at the 3’ prime end to prevent it from acting as a primer.
Primer sets G1 and G2 were used to amplify targets in the HBB gene promoter and third exon, respectively. These primers were designed to be used in GRACE-PCR with primer set G3, which amplifies a target in the CLCN7 reference gene .
All primers were checked for specificity using the NCBI Primer BLAST program (http://www.ncbi.nlm.nih.gov/tools/primer-blast). HPLC purified primers and probes were commercially synthesized (Metabion, Germany).
PCR reaction conditions
For scanning the HBB gene by HRM with primers sets H1 to H11 each 10 μl HRM-PCR reaction included 5 μL of LightCycler® 480 High Resolution Melting Master (Roche Diagnostics, Germany), 3.75 mM MgCl2, 25 ng of genomic DNA and 0.2 μM of each primer. When the unlabelled probe was included, it was used at a concentration of 0.2 μM and the forward primer concentration was reduced to 0.02 μM. PCR was performed using a Rotor-Gene Q-5plex HRM thermocycler (Qiagen). This system can accommodate up to 72 PCR reactions per run. After a 10 min hold at 95 °C, 45 cycles of PCR were performed as follows: 15 s at 95 °C, 15 s at the temperature given in Table 1, 15 s at 72 °C. Cycling was followed by a 2 min hold at 72 °C, a 1 min hold at 95 °C and then a 1 min hold at 40 °C. Amplicon melting was performed at a rate of 0.1 °C/s over the temperature ranges indicated in Table 1.
For the detection of deletions in the promoter of the HBB gene promoter, each 12.5 μL GRACE-PCR reaction included 0.25 units of Kappa HiFi HotStart polymerase (Kappa Biosystems, USA), 0.25 units of Platinum Taq polymerase (Invitrogen, USA), 1x Kappa GC buffer (Kappa Biosystems), 0.3 mM of each dNTP (Kappa Biosystems), 0.25 μL Resolite dye (Roche Diagnostics), 25 ng of genomic DNA, 0.16 μM of each G1 primer and 0.1 μM of each G3 primer (Table 1). After a 2 min hold at 98 °C, 26 cycles of PCR were performed as follows: 15 s at 95 °C, 15 s at 55 °C, 15 s at 72 °C. Melting was performed from 75 to 90 °C at a rate of 0.2 °C/s.
For the detection of deletions in the third exon of the HBB gene, each 10 μL GRACE-PCR reaction included 5 μL of LightCycler® 480 High Resolution Melting Master (Roche Diagnostics), 3.0 mM MgCl2, 25 ng of genomic DNA, 0.2 μM of each G2 primer and 0.08 μM of each G3 primer (Table 1). After a 10 min hold at 95 °C, 26 cycles of PCR were performed as follows: 15 s at 95 °C, 15 s at 55 °C, 15 s at 72 °C. Melting was performed from 80 to 90 °C at a rate of 0.2 °C/s.
Acquisition and analysis of data
Data acquisition and analysis were performed using the Rotor-Gene Q series software version 1.7 (Qiagen). For gene scanning by HRM, the raw melting curves were normalized and then converted to difference plots (Figs. 2 and 3). For the determination of gene copy number by GRACE-PCR, the melt curves were normalized before the melting of the CLCN7 control gene and after the HBB gene (Fig. 4). This resulted in a plot which allowed the CLCN7:HBB ratio to be determined by simple visual examination (Fig. 4b). Since the copy number of the CLCN7 gene is assumed to be two, the copy number of the HBB gene is easily derived.
Confirmation of deletions in the HBB gene
The Filipino β-thalassaemia deletion (NG_000007.3:g.66258_184734del118477) [3, 35] was confirmed by Gap-PCR using previously described primers . Each 12.5 μL PCR reaction included 0.25 units of Platinum Taq, 1x PCR buffer (Invitrogen), 0.3 mM of each dNTP, 2.0 mM MgCl2, 0.1 μM of each primer and 25 ng of genomic DNA (Table 2). After an initial 4 min hold at 94 °C, 38 cycles of PCR was performed as follows: 30 s at 94 °C, 30 s at 60 °C, 60 s at 72 °C. PCR products were evaluated by electrophoresis using 1 % agarose gels, with detection of a 920 bp product indicating the presence of the deletion.
The presence of the 619 bp deletion of the HBB gene was confirmed using the commercial β-Globin StripAssay (ViennaLab Diagnostics, Austria) in accordance with the manufacturer’s instructions.
Sequencing of the HBB gene
Amplicons for sequencing were obtained using the primers and reaction conditions given in Table 2. The PCR products were diluted 1:150 with distilled water prior to cycle sequencing. Cycle sequencing was performed with the BigDye Terminator v3.1 Cycle Sequencing kit (Applied Biosystems, USA) in accordance with the manufacturer’s instructions using M13 primers. Cycle sequencing products were cleaned by an ethanol and ethylenediaminetetraacetic acid (EDTA) precipitation and sequenced on a 3130 Genetic Analyser (Applied Biosystems).
Amplicon and primer design
When a common non-pathological SNP was present in a region that was optimal for primer placement, we substituted the nucleotide at the SNP position with the universal base inosine. Incorporating inosine into the primer allows equal amplification of both the wild-type and variant alleles . This proved to be a good strategy for obtaining clear HRM-PCR profiles without interference by the non-pathological SNP. The SNPs concerned were rs10768683 (HBB:c.315 + 16G > C) under the reverse H6 primer and rs1609812 (HBB:c.316-185C > T) under the reverse H7 primer (Table 1). An example showing the detection of the HBB:c.315 + 1G > A mutant and wild type alleles where inosine in the primer was used to mask the non-pathological SNP rs10768683 is shown in Fig. 2, plots 1A and 1B, and Table 3 (genotype 37).
The CLCN7 gene was selected as a reference gene for GRACE-PCR . As previously described, the CLCN7 gene is well suited as a reference gene since it has no pseudo genes and defects in the gene result in an obvious phenotype [34, 38]. The G3 primer set was designed to target the CLCN7 gene and generate a product with a melting temperature approximately 8 °C lower that the products of the G1 and G2 primer sets. This resulted in GRACE-PCR melting curves with two distinct steps corresponding to the melting of the reference gene product and the HBB product (Fig. 4). The drop in fluorescence as each product melts is proportional to the amount of the product present. Limiting the number of cycles to 26 ensured that the amount of products remained proportional to the gene copy numbers. Thus, normalized melt curve data allowed the number of gene copies to be determined by simple visual analysis (Fig. 4b).
Initial assessment of the assay
The assay was initially assessed using 342 samples for which the HBB gene had been previously sequenced and included two samples positive for the 619 bp deletion. For 26 of these samples the assay did not yield the expected result, with 23 due to interference from the common synonymous SNP HBB:c.9C > T, two due to deletions and one was an abnormal melt curve obtained with primer set H1.
The initial version of the assay was not able to distinguish heterozygous Hb C (Table 3; Genotype 7, HBB:c.19G > A; n = 5) from the heterozygous form of the common synonymous SNP HBB:c.9C > T (Fig. 3a and b). Likewise, homozygous Hb S (Table 3; Genotype 10, HBB:c.20A > T; n = 18) could not be distinguished reliably from the homozygous form of the HBB:c.9C > T SNP. This issue was resolved by the use of the unlabelled probe [37, 39] designed to cover the HBB:c.19G > A and HBB:c.20A > T mutations, but not the nearby HBB:c.9C > T SNP (Fig. 1 and Table 1). This approach resulted in clear separation of the pathological mutations in the probe melt region (Fig. 3c). Repeating the HRM-PCR analysis in the presence of the unlabelled probe for the 23 samples in which a mutation was initially missed, yielded results consistent with sequencing (Table 3; columns Sequencing versus Sample set 1).
GRACE-PCR correctly detected the two samples known to be positive for the 619 bp deletion as having deletions of the third exon of the HBB gene (Table 3, genotypes 42 and 43). GRACE-PCR also indicated that two additional samples had large deletions of the HBB gene involving both the promoter and the third exon. Both of these unexpected cases were subsequently confirmed by Gap-PCR to be positive for the Filipino β-thalassaemia deletion (Table 3, genotype 44).
After modifying the assay to include the probe and confirming the presence of the deletions, only one of the 342 samples in the initial assessment failed to give the expected result. This sample produced an abnormal melting curve with the H1 primer set the first time it was tested. On retesting the sample was found to be wild type, this sample is listed as a false positive in Table 3.
Validation of the assay
An additional 68 samples without prior sequence information were analysed by the assay. Sixty three samples had HRM melting curves that indicated the presence of sequence variations, other than the common HBB:c.9C > T SNP, and for one sample GRACE-PCR indicated a gene deletion. For the remaining four samples there was no evidence of a sequence change. It was noted that the HRM-PCR difference plots were highly reproducible between runs and between samples (Fig. 2; panels a versus b). This raised the possibility that HRM-PCR could be used not only to detect sequence variations but also reveal the actual nucleotide change. For 53 of the 63 samples in which a sequence variation was detected, it was possible to assign a provisional genotype based only on the HRM-PCR data. Subsequent sequencing of the HBB gene confirmed that the four samples without variation by HRM-PCR had not contained any sequence changes, and that the 63 samples identified with variations all carried pathological mutations. The sample, which tested positive by GRACE-PCR, was confirmed by Gap-PCR to be another example of the Filipino β-thalassaemia deletion. Furthermore, all 53 of the provisional genotypes assigned by analysis of the HRM-PCR data were proved to be correct. The 10 samples, for which no provisional genotype could be assigned, were found to have genotypes that were either absent, or only present in a very few cases in the initial 342 samples (Table 3; columns Sample set 2 versus Sample set 1).
A combination of sequencing, Gap-PCR and the commercial β-globin StripAssay showed that 361 of the 410 samples included in the study contained one or more mutations in the HBB gene. The assay was able to detect all 361 abnormal samples and generated just one false positive result (Table 3; columns Sequencing versus Sample sets 1&2). This equates to a sensitivity of 100.0 % (95 % confidence interval 98.7 - 100.0 %) and a specificity of 98.0 % (95 % confidence interval 87.8 – 99.9 %). In total, the 361 positive samples in the study accounted for 44 distinct pathological genotypes (Table 3).
Quantification of Hb A2 is routinely used for the detection of β-thalassaemia carriers. The normal range for Hb A2 is 2.0–3.3 % with levels above 3.5 % being considered indicative of β-thalassaemia trait . The current study included 22 samples with Hb A2 levels greater than 3.5 %, which did not have detectable mutations or deletions of the HBB gene. This is in line with previous observations that false positive results can be a significant problem when screening for β-thalassaemia carriers based on Hb A2 levels . Indeed, elevated Hb A2 in the absence of β-thalassaemia has recently been linked to mutations in the KLF1 gene , suggesting that molecular confirmation may be required. In our laboratory, HBB gene mutations are routinely detected with a commercial reverse hybridization assay (ViennaLab Diagnostics). This commercial kit was designed to detect the 22 most prevalent mutations found in Indian and Middle Eastern populations. The current study identified 63 samples (15.4 %) with mutations not targeted by this kit (Table 3). This further illustrates the limitations of conventional PCR techniques when screening for β-thalassaemia trait in highly heterogeneous populations, and highlights the advantages of using a more universal assay, such as the one described here.
The first application of HRM-PCR to the HBB gene was described by Wittwer et al. . This was a small-scale study designed to detect the mutations responsible for the abnormal haemoglobins Hb S and Hb C . Despite the fact that they amplified approximately the same section of the gene covered by the H2 primer set in the current study, the authors did not report any problems with the HBB:c.9C > T SNP. This is presumably due to the relatively small number of samples tested (n = 12). In contrast, in the current study the HBB:c.9C > T SNP was encountered frequently and it was not possible to resolve all genotypes without the addition of the unlabelled probe.
Three previous studies have described HRM-PCR assays to scan the HBB gene in different populations. Shih et al. , Saetung et al.  and Lin et al.  used the technique to detect the most prevalent mutations in Taiwanese, Thai and Chinese populations respectively [29–31]. These assays scanned selected areas of the HBB gene and the assay described by Saetung et al. also included a tube for the detection of the 3.4 kb deletion. Both Shih et al.  and Lin et al.  reported problems with the HBB:c.9C > T SNP, which were resolved by placing the primer over the SNP. A possible concern with this approach is that it might result in allele dropout, particularly if a SNP is located close to the 3’ end of a primer. The primers described by Shih et al.  put the SNP close to the 3’ end, while the primer used by Lin et al.  place the SNP in the middle of the primer where it is less likely to have a significant effect. However, the HBB:c.19G > A and HBB:c.20A > T mutations are located at the 3’ end of the primer used by Lin et al. , which would probably affect the amplification of these common pathologically significant alleles. Since these mutations are rare in China  and were not included in the specific mutations being targeted this would not have been a significant concern when screening that population. However, these mutations are very common throughout Africa and the Middle East , thus the assay described by Lin et al.  is not well suited for screening in these regions. In contrast to overlaying the HBB:c.9C > T SNP with a primer, the approach described here of using the unlabelled probe, allowed detection of all the mutations we encountered in this area of the HBB gene.
Saetung et al.  avoided complications from the HBB:c9.C > T SNP by simply not scanning this region of the HBB gene. This approach limits the usefulness of their assay in populations where mutations in this area of the gene are common.
Other high frequency SNPs also complicated HRM analysis of the HBB gene. These were HBB:c.315 + 16G > C (rs10768683) and HBB:c.316-185C > T (rs1609812). We addressed this issue by substituting the respective SNP with inosine in each primer designed to cover these bases. This approach, as previously demonstrated for the BRCA genes , eliminated the detection of this high frequency SNP and provided clean HRM profiles (Fig. 2; 1a and 1b).
Until recently, deletions of the HBB gene were considered to be rare with only the 619 base deletion being common, accounting for around 20 % of β-thalassaemia alleles in some parts of India and Pakistan . However, it has now become apparent that HBB gene deletions are more frequent than was previously believed. The 3.4 kb deletion has a prevalence of up to 4 % in some parts of Thailand  and the Filipino β-thalassaemia deletion has been described as being a common cause of β-thalassaemia in Filipinos [36, 45]. Indeed, it has been suggested that deletions involving the β-globin gene cluster account for as many as 10 % of all β-thalassaemia mutations .
Conventional HRM-PCR works well for the detection of point mutations and small indels, but is not a useful technique for the detection of larger gene rearrangements. Melting curve assays have been described to detect the 3.4 kb deletion of the HBB gene [31, 44]. However, for other deletions Gap-PCR is still commonly used. Although Gap-PCR is a simple and robust technique, it is not ideal for high volume screening, since the use of agarose gel electrophoresis for detection makes it a time consuming open tube technique. An important limitation of both Gap-PCR and the existing melt curve assays is that they can only be used to detect specific deletions with well-defined breakpoints and consequently are not suitable for the detection of novel or rare deletions. In contrast, the GRACE-PCRs included in the current study determine the copy number ratio between the HBB gene and the reference gene. Thus GRACE-PCR can be used to detect large rearrangements of the HBB gene without prior knowledge of the breakpoints. In order to maximise the number of deletions that could be detected our assay used two GRACE-PCR reactions, one targeting the promoter and the other targeting the third exon.
Allele dropout due to a SNP within a priming site is a potential risk with any PCR based assay. The use of the probe and primers containing inosine in the current assay addressed this concern for the most common high frequency SNPs associated with the HBB gene. However, with 868 known pathological mutations of the HBB gene, as well as a number of non-pathological ones, it is inevitable that there are SNPs associated with most priming sites. The risk of misinterpretation due to allele dropout is significantly mitigated by the overlapping nature of the HRM amplicons and by the use of two GRACE-PCR reactions. Furthermore, the fact that no problems due to allele dropout were encountered in the 410 samples evaluated in the study indicated that the risk is indeed small. Nevertheless, it is recommended that the possibility of allele dropout be considered when results of the assay do not correlate with other laboratory findings or the clinical picture.
The general purpose of gene scanning is to identify amplicons carrying variant sequences. The HRM-PCR profiles of the amplicons not only indicated the presence of a variation in the sequence, but also often allowed the identification of the exact nucleotide change. In the 68 blindly analysed samples all variant sequences were correctly detected. Furthermore, in 84 % of these cases the melting curves were sufficiently distinctive for the correct genotype to be assigned prior to sequencing (Table 3; columns Sample set 2 versus Sample set 1). We believe, as more experience is gained with this assay, provisional direct genotyping of almost all samples should be possible. The observed genotyping ability of HRM-PCR for the HBB gene contrasts favourably with denaturing HPLC (DHPLC) assays, as DHPLC may require either a two-step assay , or may require to be repeated after spiking with control DNA from a known variant .
When a mutation occurred in overlapping amplicon regions (Fig. 1 and Table 1) it was detected in both amplicons (Table 3). This provided additional information assisting in the assignment of the correct genotype. For example, a sample found to be positive with both primer sets H4 and H5 was likely to be either HBB:c.126_129delCTTT, HBB:c.135delC or HBB:c.140G > A, since they were the only mutations identified in the current study that occurred in this region of overlap.
Interestingly, SNPs in linkage disequilibrium with pathological mutations also provide additional information which assists in assigning the correct genotype. For example, an association between Hb Monroe (HBB:c.92G > C) and the SNP (HBB:c.-92C > G) in the HBB gene promoter has been reported . In the current study three cases of Hb Monroe in association with HBB:c.-92C > G were found (Table 3; Genotype 22). The Hb Monroe mutation and the SNP resulted in distinctive melting curves in amplicons obtained with primer sets H3 and H1, respectively.
No mutations were found in the region of the HBB gene scanned by primer sets H7, H10 and H11, indicating that mutations in this part of the gene are comparatively rare in the population studied (Table 3). It is interesting to note that even in the highly heterogeneous population of the UAE, 91 % of positive cases could be detected using just 5 of the 13 primer sets. This suggests that a stepwise approach could be used starting with high mutation frequency amplicons, thus reducing the number of PCR reactions required to obtain results.
The combination of the HRM-PCR and GRACE-PCR assays theoretically cover 97 % of the mutations listed in the HbVar database . This compares favourably with previously described assays that range from 35 to 90 % theoretical coverage [29, 30]. An example of a mutation that would not be detected by the current assays is HBB:c.203_204delTG  due to lack of overlap between primer sets H5 and H6. Similarly, HBB:c.316-146 T > G and HBB:c.316-125A > G  would not be detected due to lack of overlap between primer sets H7 and H8.
Only negative (wild type) controls were included in the current study, since all samples screened were also sequenced. However, when using the assays for clinical testing, the inclusion of appropriate positive and negative controls with each run is highly recommended. Such controls could be genomic DNA from confirmed cases or plasmids .
In our laboratory, the costs for scanning the HBB gene with the combined HRM-PCR/GRACE-PCR assay were approximately one eighth the cost of sequencing. Given that in at least 84 % of cases HRM-PCR allowed the assignment of a provisional genotype, this represents an opportunity for significant cost savings. Should sequencing be required to evaluate an unexpected melt curve, this can be performed by directly sequencing the positive HRM amplicon , which is still cheaper and faster than sequencing the entire HBB gene. In addition, the HRM-PCR / GRACE-PCR approach allows the detection of large deletions which are not readily detected by sequencing. In the current study, the samples were manually pipetted and processed in batches of 72 PCR reactions. Since both HRM-PCR and GRACE-PCR lend themselves to automation, robotic reaction setup is possible, and if a micro-titre plate based PCR system was used, up to 384 reactions could be run simultaneously on a single plate. Although we used different annealing temperatures, all the primers used for HRM-PCR will perform adequately at 55 °C, which would simplify batch processing if micro-titre plates were adopted.
In conclusion, the current study has resulted in the development of a robust assay for the detection of pathological HBB gene mutations. The approaches used for handling the high frequency SNPs mean that the assay is universally applicable, rather than being population specific. The inclusion of GRACE-PCR targeting both the gene promoter and third exon means that the assay has the ability to detect a wide range of large deletions of the HBB gene. The assay has sensitivity and specificity comparable to sequencing and was able to detect all mutations in the samples tested. Furthermore, this assay can be run at a fraction of the cost of a full sequencing approach as 84 % of positive samples could be provisionally genotyped. The assay has provided useful data on the spectrum of HBB gene mutations in the highly heterogeneous population of the United Arab Emirates. In future it should be able to serve as a screening test for suspected β-thalassaemia carriers and variant haemoglobins.
Amplification Refractory Mutation System
Allele specific oligonucleotide
- CLCN7 :
Chloride channel voltage sensitive 7 gene
Denaturing gradient gel electrophoresis
Denaturing high performance liquid chromatography
Gene ratio assay copy enumeration
- HBB :
Haemoglobin beta gene
- HBD :
Haemoglobin delta gene
High performance liquid chromatography
High resolution melting
- KLF1 :
Krüppel-like Factor 1 gene
Polymerase chain reaction
Single nucleotide polymorphism
Single Strand Conformational Polymorphism
World Health Organization
Modell B. Global epidemiology of haemoglobin disorders and derived service indicators. Bull World Health Organ. 2008;86:480–7.
Thein S. Genetic modifiers of β-thalassemia. Haematologica. 2005;90:649–60.
Hbvar. Database of Human Hemoglobin Variants and Thalassemia Mutations. http://globin.bx.psu.edu/hbvar/menu.html. Accessed 1 Sept 2015.
Harteveld CL. Nine unknown rearrangements in 16p13.3 and 11p15.4 causing α- and β-thalassaemia characterised by high resolution multiplex ligation-dependent probe amplification. J Med Genet. 2005;42:922–31.
Theodoridou S, Alemayehou M, Prappas N, Karakasidou O, Aletra V, Plata E, Tsaftaridis P, Karababa P, Boussiou M, Sinopoulou K, Hatzi A, Voskaridou E, Loutradi A, Manitsa A. Carrier screening and prenatal diagnosis of hemoglobinopathies. A study of indigenous and immigrant couples in northern Greece, Over the last 5 years. Hemoglobin. 2008;32:434–9.
Angastiniotis M, Kyriakidou S, Hadjiminas M. How thalassaemia was controlled in Cyprus. World Health Forum. 1986;7:291–7.
Bozkurt G. Results from the North Cyprus Thalassemia prevention program. Hemoglobin. 2007;31:257–64.
Cousens NE, Gaff CL, Metcalfe SA, Delatycki MB. Carrier screening for Beta-thalassaemia: a review of international practice. Eur J Hum Genet. 2010;18:1077–83.
Bain BJ. Haemoglobinopathy diagnosis: Algorithms, lessons and pitfalls. Blood Rev. 2011;25:205–13.
Ryan K, Bain BJ, Worthington D, James J, Plews D, Mason A, Roper D, Rees DC, la Salle de B, Streetly A. Significant haemoglobinopathies: guidelines for screening and diagnosis. British J Haematol. 2010;149:35–49.
Stephens AD, Angastiniotis M, Baysal E, Chan V, Fucharoen S, Giordano PC, Hoyer JD, Mosca A, Wild B. ICSH recommendations for the measurement of Haemoglobin A2. Int J Lab Hematol. 2011;34:1–13.
Galanello R, Ruggeri R, Addis M, Paglietti E, Cao A. Hemoglobin A2 in iron deficient β-thalassemia heterozygotes. Hemoglobin. 1981;5:613–8.
Harteveld CL, Higgs DR. α-thalassaemia. Orphanet J Rare Dis. 2010;5:13.
Perseu L, Satta S, Moi P, Demartis FR, Manunza L, Sollaino MC, Barella S, Cao A, Galanello R. KLF1 gene mutations cause borderline HbA2. Blood. 2011;118:4454–8.
Cao A, Galanello R. Beta thalassemia. Genet Med. 2010;12:61–76.
Old J, Petrou M, Varnavides L, Layton M, Modell B. Accuracy of prenatal diagnosis for haemoglobin disorders in the UK: 25 years’ experience. Prenat Diagn. 2000;20:986–91.
Old JM. Screening and genetic diagnosis of haemoglobin disorders. Blood Rev. 2003;17:43–53.
Vrettou C, Traeger-Synodinos J, Tzetis M, Malamis G, Kanavakis E. Rapid Screening of multiple β-globin gene mutations by real-time PCR on the LightCycler: application to carrier screening and prenatal diagnosis of thalassemia syndromes. Clin Chem. 2003;49:769–76.
Patrinos GP, Kollia P, Papadakis MN. Molecular diagnosis of inherited disorders: lessons from hemoglobinopathies. Hum Mutat. 2005;26:399–412.
Erali M, Wittwer CT. High resolution melting analysis for gene scanning. Methods. 2010;50:250–61.
Chou L-S, Lyon E, Wittwer CT. A comparison of high-resolution melting analysis with denaturing high-performance liquid chromatography for mutation scanning: cystic fibrosis transmembrane conductance regulator gene as a model. Am J Clin Pathol. 2005;124:330–8.
Sadr-Nabavi A, Hoffmann M, Weilke C, Sun Y-L, Nevinny-Stickel-Hinzpeter C. Mutation scanning using high resolution melting or dHPLC: a performance comparison study. Biochemica. 2009;4:30–3.
Payne MS, Tabone T, Kemp MW, Keelan JA, Spiller OB, Newnham JP, Carroll KC. High-Resolution Melt PCR analysis for genotyping of ureaplasma parvum Isolates directly from clinical samples. J Clin Microbiol. 2014;52:599–606.
Pickering J, Binks MJ, Beissbarth J, Hare KM, Kirkham LAS, Smith-Vaughan H, Forbes BA. A PCR-high-resolution melt assay for rapid differentiation of Nontypeable Haemophilus influenzae and Haemophilus haemolyticus. J Clin Microbiol. 2014;52:663–7.
Ihle MA, Fassunke J, Konig K, Grunewald I, Schlaak M, Kreuzberg N, Tietze L, Schildhaus H-U, Buttner R, Merkelbauch-Brause S. Comparison of high resolution melting analysis, pyrosequencing, next generation sequencing and immunohistochemistry to conventional Sanger sequencing for the detection of p.V600E and non-p. 600E BRAF mutations. BMC Cancer. 2014;14:1–13.
Chen Y, Wu Q, Shao Y, Zhang J, Guan M, Wan J, Yu B, Zhang W. Identify the association between polymorphisms of BLK and systemic lupus erythematosus through unlabelled probe-based high-resolution melting analysis. Int J Immunogenetics. 2012;39:321–7.
Blombery PA, Wong SQ, Hewitt CA, Dobrovic A, Maxwell EL, Juneja S, Grigoriadis G, Westerman DA. Detection of BRAF mutations in patients with hairy cell leukemia and related lymphoproliferative disorders. Haematologica. 2012;97:780–3.
Santos S, Marques V, Pires M, Silveira L, Oliveira H, Lança V, Brito D, Madeira H, Esteves JF, Freitas A, Carreira IM, Gaspar IM, Monteiro C, Fernandes AR. High resolution melting: improvements in the genetic diagnosis of hypertrophic cardiomyopathy in a Portuguese cohort. BMC Med Genet. 2012;13:17.
Lin M, Jiao J-W, Zhan X-H, Zhan X-F, Pan M-C, Wang J-L, Wang C-F, Zhong T-Y, Zhang Q, Yu X, Wu J-R, Yang H-T, Lin F, Tong X, Yang H, Zha G-C, Wang Q, Zheng L, Wen Y-F, Yang L-Y. High resolution melting analysis: a rapid screening and typing tool for common β-Thalassemia mutation in Chinese population. PLoS One. 2014;9:e102243.
Shih H-C, Er T-K, Chang T-J, Chang Y-S, Liu T-C, Chang J-G. Rapid identification of HBB gene mutations by high-resolution melting analysis. Clin Biochem. 2009;42:1667–76.
Saetung R, Ongchai S, Charoenkwan P, Sanguansermsri T. Genotyping of beta thalassemia trait by high-reolution melting analysis. Southeast Asian J Trop Med Public Health. 2013;44:1–10.
Turner A, Sasse J, Varadi A. Development and validation of a high throughput, closed tube method for the determination of haemoglobin alpha gene (HBA1 and HBA2) numbers by gene ratio assay copy enumeration-PCR (GRACE-PCR). BMC Med Genet. 2015;16:115.
Pritchard CC, Tait JF, Buller-Burckle AM, Mikula M. Annotation error of a common β0-thalassemia has implications for molecular diagnosis. Am J Hematol. 2010;85:987.
Fallah M-S, Mahdian R, Aleyasin S-A, Jamali S, Hayat-Nosaeid M, Karimipour M, Raeisi M, Zeinali S. Development of a quantitative real-time PCR assay for detection of unknown α-globin gene deletions. Blood Cells Mol Dis. 2010;45:58–64.
Eng B, Chui DH, Saunderson J, Olivieri N. Identification of two novel beta zero-thalassemia mutations in a Filipino family: framshift codon 67 (-TG) and a beta-globin deletion. Hum Mutat. 1993;2:375–9.
Yamashiro Y, Hattori Y, Nitta T, Adhiyanto C, Matar M, Ferania M, Takagi F. Filipino-type βo-thalassemia has 116 kb deletion: its correct breakpoints and five cases found in Japan. Bull Yamaguchi Med Sch. 2012;59:53–9.
Krypuy M, Ahmed A, Etemadmoghadam D, Hyland SJ, deFazio A, Fox SB, Brenton JD, Bowtell DD, Dobrovic A. A high-throughput protocol for mutation scanning of the BRCA1 and BRCA2 genes. BMC Cancer. 2011;11:265.
Babashah S, Jamali S, Mahdian R, Nosaeid MH, Karimipoor M, Alimohammadi R, Raeisi M, Maryami F, Masoudifar M, Zeinali S. Detection of unknown deletions in β-globin gene cluster using relative quantitative PCR methods. Eur J Haematol. 2009;83:261–9.
Zhou L, Wang L, Palais R, Pryor R, Wittwer CT. High-resolution DNA melting analysis for simultaneous mutation scanning and genotyping in solution. Clin Chem. 2005;51:1770–7.
Giambona A, Passarello C, Vinciguerra M, Li Muli R, Teresi P, Anza M, Ruggeri G, Renda D, Maggio A. Significance of borderline hemoglobin A2 values in an Italian population with a high prevalence of β-thalassemia. Haematologica. 2008;93:1380–4.
Wittwer CT, Reed GH, Gundry CN, Vandersteen JG, Pryor RJ. High-resolution genotyping by amplicon melting analysis using LCGreen. Clin Chem. 2003;49:853–60.
Zeng Y-T, Huang S-Z. The studies of hemoglobinopathies and thalassemia in China—the experiences in Shanghai Institute of Medical Genetics. Clin Chim Acta. 2001;3:107–11.
Weatherall D, Clegg J. Inherited haemoglobin disorders: an increasing global health problem. Bull World Health Organ. 2001;79:704–12.
Prathomtanapong P, Pornprasert S, Phusua A, Suanta S, Saetung R, Sanguansermsri T. Detection and identification of β thalassemia 3.5kb deletion by SYBR Green1 and high resolution melting analysis. Eur J Haematol. 2009;82:159–60.
Waye JS, Eng B, Hunt JA, Chui DHK. Filipino β-thalassemia due to a large deletion: identification of the deletion endpoints and polymerase chain reaction (PCR)-based diagnosis. Hum Genet. 1994;94:530–2.
Su Y-N, Lee C-N, Hung C-C, Chen C-A, Cheng W-F, Tsao P-N, Yu C-L, Hsieh F-J. Rapid detection of β-globin gene (HBB) mutations coupling heteroduplex and primer-extension analysis by DHPLC. Hum Mutat. 2003;22:326–36.
Li Q, Li L-Y, Huang S-W, Li L, Chen X-W, Zhou W-J, Xu X-M. Rapid genotyping of known mutations and polymorphisms in β-globin gene based on the DHPLC profile patterns of homoduplexes and heteroduplexes. Clin Biochem. 2008;41:681–7.
Fedorov AN, Nasyrov FY, Smirnova EA, Bocharova TN, Limborska SA. IVS-1-1 (G → C) in Combination with -42 (C → G) in the promoter region of the β-globin gene in patients from Tajikistan. Hemoglobin. 1993;17:275–8.
AT conducted the experimental work and wrote the manuscript. JS and AV were both significantly involved in the study design, writing and critical revision of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
About this article
Cite this article
Turner, A., Sasse, J. & Varadi, A. Rapid detection of pathological mutations and deletions of the haemoglobin beta gene (HBB) by High Resolution Melting (HRM) analysis and Gene Ratio Analysis Copy Enumeration PCR (GRACE-PCR). BMC Med Genet 17, 75 (2016). https://doi.org/10.1186/s12881-016-0334-y
- Beta thalassaemia
- Copy number determination
- Gene quantification