Skip to main content
  • Research article
  • Open access
  • Published:

Lower reference limits of quantitative cord glucose-6-phosphate dehydrogenase estimated from healthy term neonates according to the clinical and laboratory standards institute guidelines: a cross sectional retrospective study

Abstract

Background

Previous studies have reported the lower reference limit (LRL) of quantitative cord glucose-6-phosphate dehydrogenase (G6PD), but they have not used approved international statistical methodology. Using common standards is expecting to yield more true findings. Therefore, we aimed to estimate LRL of quantitative G6PD detection in healthy term neonates by using statistical analyses endorsed by the International Federation of Clinical Chemistry (IFCC) and the Clinical and Laboratory Standards Institute (CLSI) for reference interval estimation.

Methods

This cross sectional retrospective study was performed at King Abdulaziz Hospital, Saudi Arabia, between March 2010 and June 2012. The study monitored consecutive neonates born to mothers from one Arab Muslim tribe that was assumed to have a low prevalence of G6PD-deficiency. Neonates that satisfied the following criteria were included: full-term birth (37 weeks); no admission to the special care nursery; no phototherapy treatment; negative direct antiglobulin test; and fathers of female neonates were from the same mothers’ tribe. The G6PD activity (Units/gram Hemoglobin) was measured spectrophotometrically by an automated kit. This study used statistical analyses endorsed by IFCC and CLSI for reference interval estimation. The 2.5th percentiles and the corresponding 95% confidence intervals (CI) were estimated as LRLs, both in presence and absence of outliers.

Results

207 males and 188 females term neonates who had cord blood quantitative G6PD testing met the inclusion criteria. Method of Horn detected 20 G6PD values as outliers (8 males and 12 females). Distributions of quantitative cord G6PD values exhibited a normal distribution in absence of the outliers only. The Harris-Boyd method and proportion criteria revealed that combined gender LRLs were reliable. The combined bootstrap LRL in presence of the outliers was 10.0 (95% CI: 7.5-10.7) and the combined parametric LRL in absence of the outliers was 11.0 (95% CI: 10.5-11.3).

Conclusion

These results contribute to the LRL of quantitative cord G6PD detection in full-term neonates. They are transferable to another laboratory when pre-analytical factors and testing methods are comparable and the IFCC-CLSI requirements of transference are satisfied. We are suggesting using estimated LRL in absence of the outliers as mislabeling G6PD-deficient neonates as normal is intolerable whereas mislabeling G6PD-normal neonates as deficient is tolerable.

Peer Review reports

Background

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is a common X-linked recessive enzymopathy [1]. Genotypically, males are either hemizygous for the G6PD gene with normal gene expression or G6PD deficient, while females can be normal, heterozygous, or homozygous. While about 186 variants have been recognized, not all of them are clinically significant [2]. The World Health Organization (WHO) has classified G6PD variants into five classes: Class I (severe enzyme deficiency with chronic non-spherocytic hemolytic anemia); Class II (<10% of normal); Class III (10–60% of normal); Class IV (60–100% of normal); and Class V (>200% of normal) [3]. The G6PD deficiency is well-known to cause hyperbilirubinemia that may be severe enough to cause kernicterus or neonatal death [1, 4]. Early recognition of G6PD deficiency can help prevent these serious complications. Thus, the WHO has endorsed screening cord blood samples from all neonates in populations with a prevalence of G6PD deficiency of 3 to 5% or more in males [5]. The WHO has endorsed 5 methods as screening or diagnostic tests for G6PD deficiency including quantitative measurement of G6PD activity in red blood cells (RBCs) [5, 6]. Reference limits are the cornerstones of interpretation of any laboratory result, including quantitative cord G6PD values. Previous studies have reported on the lower reference limit (LRL) or lower decision limit (LDL) of quantitative cord G6PD values, but these have methodological and/or statistical flaws [718]. These studies did not implement standard statistics used to estimate LRLs that have long been endorsed by the International Federation of Clinical Chemistry (IFCC) and the Clinical and Laboratory Standards Institute (CLSI) [1925]. None of these studies have addressed the detection and handling of outliers. Only two of these studies addressed the type of distribution of G6PD values and were in accordance with known methods of defining reference intervals (RIs), including the 95% central RIs [7, 8]. Instead of reporting non-parametric RIs that do not make any assumptions about the data distribution, [26] some researchers have based the LRL on the mean without addressing normal distribution testing [13]. Some reported LRLs/LDLs have been estimated from small sample sizes, [7, 11, 12] from populations with high G6PD deficiency prevalence, [810, 13, 15, 16] and from mixtures of G6PD-normal and G6PD-deficient preterm and full-term neonates [8, 15]. Therefore, we aimed to estimate the LRLs of quantitative cord G6PD activity from a large population of healthy term neonates in accordance with the standard method of IFCC-CLSI [1925]. In our opinion, using of these common standards would overcome limitations of the previous studies and yield true findings [27]. We elected not to estimate the upper reference limit as it has no clinical implication [28].

Methods

Setting

We conducted this cross sectional retrospective study at King Abdulaziz Hospital (KAH) in the Al-Ahsa area of the Eastern province of Saudi Arabia. Since 2009, KAH has been accredited by the Joint Commission of International Accreditation and its laboratory has been accredited by the College of American Pathologists, which implements the CLSI standards in its accreditation checklists.

Analytic method for G6PD

In early 2008, KAH began universal cord blood screening for G6PD deficiency coupled with direct antiglobulin testing and blood grouping [29]. A semi-qualitative fluorescence spot test (FST) with a cut-off point of 2.1 Units/gram Hemoglobin (U/g Hb) was used until 01 March 2010, at which point it was replaced with an automated commercial kit. This kit is Udilipse Auto Analyzer from United Diagnostics Industry, Dammam, Saudi Arabia that offers quantitative measurement of G6PD activity.

Just after delivery of the placenta, whole cord blood was collected in ethylenediaminetetraacetic acid (EDTA) Vacutainer tubes (Becton-Dickinson, Rutherford, NJ, USA). Blood specimens were transported to the laboratory by pneumatic tubes. Quantitative G6PD activity measurements were performed in batches every morning, 7 days a week, with blood samples stored in the vertical position at 2–8 degree Celsius (°C) until analysis. The principal method of the Udilipse kit is in accordance with the standardized WHO method for G6PD assay of the hemolysate at 25°C outlined as follows [6]:

Glucose 6 Phosphate + NADP G 6 PD 6 - Phosphogluconate + NADPH + H +

The reagents of this kit are outlined as follows:

Reagent 1 (G6PD Buffer): Ready to use 50 micromolar (mM)/Liter triethanolamine buffer, 5 mM/Liter EDTA, pH 7.6 ± 0.05 (25°C).

Reagent 2 (G6PD NADP): Reconstituted 30 mM nicotinamide adenine dinucleotide phosphate (NADP).

Reagent 3 (G6PD Substrate): Reconstituted 17 mM Glucose-6-Phosphate Sodium.

Reagent 4 (G6PD Lysis): Ready to use 0.2% Saponin aqueous solution.

Hemolysate was prepared according to the manufacturer’s manual for the Udilipse kit by adding a well-mixed 100 microliters of whole blood to 400 microliters of 0.2% Saponin. The G6PD activity was measured within one hour of hemolysate preparation. The activity of the G6PD enzyme is measured by the rate of NADPH formation, which is measured spectrophotometrically by means of the increase in extinction at 340, 334 or 365 nanometer. As the Udilipse kit only measures G6PD activity, hemoglobin was measured spectrophotometrically on the same sample by CELL-DYN Sapphire (Abbott Diagnostics, Santa Clara, CA, USA). Then, the values of G6PD activity and hemoglobin were entered and stored in Cerner Lab Information System Software (Cerner Corporation, Kansas City, MO, USA), which is already programmed to express the G6PD activity in U/g Hb by dividing the obtained G6PD activity by the obtained hemoglobin values.

Reference sample group

The IFCC-CLSI recommends estimating RLs from a healthy population [20, 25]. In our case, the healthy population would be the population free from G6PD mutations. The Al-Ahsa area is composed of an oasis part inhabited by Arab Muslims of urban descent and a desert part inhabited by Arab Muslims of Bedouin descent. The overall prevalence of G6PD deficiency in the Al-Ahsa area is 23% in males and 13% in females [30]. The G6PD-Mediterranean (WHO class II) constitutes 84% of G6PD mutations in this population and the G6PD-A- (WHO class III) represents 5.8% [2, 28, 31]. G6PD deficiency is confined to the oasis part of the area, as malaria was much more prevalent in the oasis than in the desert [3234]. A recent retrospective Ahsai study found no single severe G6PD deficient case among 236 neonates from one Arabic Bedouin tribe subjected to cord blood screening by FST [29]. Therefore, this tribe can function as the reference population for estimating the LRLs of G6PD [35]. Henceforth, this tribe will be referred to as the reference tribe.

The reference sample group was identified from the delivery room log book via the mothers’ names. It consisted of all consecutive neonates born to mothers from the reference tribe between March, 2010 and June, 2012 that satisfied all the following posteriori inclusion criteria: 1) full-term (37 weeks of gestation); 2) roomed in with mothers and were not admitted to the special care nursery; 3) no phototherapy treatment during the neonatal period; 4) negative direct antiglobulin test; and 5) fathers of female neonates were from the reference tribe. For male neonates it does not matter whether the father is from the reference tribe or not, but it does matter for female neonates as the G6PD is inherited as X-linked recessive. Thus, only copies of birth notices of female neonates were reviewed to ascertain that both parents were from the reference tribe. At KAH, birth notices that include the full name of both parents are issued upon home discharges and copies of these birth notices are stored in hard and electronic medical records of neonates. We selected only full-term neonates as we did not expect to have an adequate number of preterm neonates of the reference sample group during the study period and as G6PD values have been shown to be higher in preterm neonates <34 weeks of gestation [15, 36, 37]. Gestational age was calculated according to the best obstetric estimate at KAH based on the first or second trimester obstetric ultrasound and/or the last menstrual period, and on the Ballard score when the best obstetric estimate was uncertain [38, 39]. The phototherapy threshold used in this hospital has been published previously and is more conservative than those used by the American Academy of Pediatrics [29]. All the study data were extracted from electronic health care records (QuadraMed CPR 5.0.9, Reston, VA, USA). This study was exempted from review by the Institutional Review Board.

Statistical analysis

We performed statistical analyses as endorsed by the IFCC- CLSI [24, 25]. The IFCC-CLSI has endorsed the use of the 2.5th percentile of the values as the LRL. Further, this body endorses the use of a minimum of 146 reference individuals in each gender partition to calculate the 95% confidence interval (CI) of the 2.5th percentile [25]. We estimated the 2.5th percentiles for males, females, and the combined group.

The IFCC-CLSI suggests detecting outliers by Tukey’s boxplot for unskewed data and by method of Horn for skewed data. The method of horn uses Tukey’s boxplot on Box-Cox transformed data [40]. Although the method of Horn is implanted in RefVal 4.11, [41] we also estimated the best lambda and performed the Box-Cox transformation by using already available SPSS syntax [42]. The two Tukey’s inner fences are the 25th percentile minus 1.5 interquartile range (IQR) and the 75th percentile plus 1.5 IQR. The two Tukey’s outer fences are the 25th percentile minus 3.0 IQR and the 75th percentile plus 3.0 IQR. Values outside the inner fences but inside the outer fences are considered mild outliers and values outside the outer fences are considered extreme outliers.

Data extraction and entry of the outliers were double-checked. It was difficult to ascertain whether these outliers were due to G6PD mutations or to pre-analytic errors as these values were not cross-referenced with molecular testing. Assuming that the reference tribe is in Hardy-Weinberg equilibrium, the Hardy-Weinberg equation was used as a surrogate to determine whether these outliers represent G6PD mutations [43, 44]. The asymptotic Pearson’s chi-square goodness-of-fit test with 2 degrees of freedom was calculated to test the departure from Hardy-Weinberg proportions [45]. Distributions of G6PD values with and without outliers were presented as histograms with superimposed best-fitting normal distribution curves. Non-integer G6PD values were rounded down to the nearest integer at the boundaries of histogram bins. The type of distribution of the G6PD values was assessed by visual examination of the histograms, comparison between measures of central tendency (mean, median, and mode), skewness, kurtosis, the Anderson-Darling test, and the Shapiro-Wilk test [46, 47]. A distribution was considered normal if it had a two-sided P value > .05 for skewness coefficient, kurtosis coefficient, and the Anderson-Darling or Shapiro-Wilk tests [46, 47]. The IFCC recommends testing for normal distribution using coefficients of skewness and kurtosis and the Anderson-Darling test [46]. The Shapiro-Wilk test, which has not been evaluated by the IFCC, was also used as it has recently been found to be more powerful than the Anderson-Darling test [47]. For comprehensiveness, z-scores of both skewness and kurtosis ware calculated as [Skewness (or Kurtosis)/SD of Skewness (or Kurtosis)] and z-scores between ±1.96 were considered statistically insignificant [48].

The IFCC-CLSI recommends excluding values from unhealthy individuals and estimating the LRLs from healthy individuals only [20, 25]. Thus, LRLs were estimated both in presence and absence of the detected outliers as these outliers might be due to G6PD mutations. As laboratory professionals diverge on the best method to estimate LRL, LRLs and their 95% CIs were estimated by three methods when appropriate: parametric, non-parametric bootstrap based on 500 bootstrap samples and non-parametric rank- based [24, 25, 41, 49]. Lower reference limits were also estimated based on methods of previous studies in order to compare between LRLs based on those methods and those based on method of the IFCC-CLSI. The Harris-Boyd method was used to assess the reliability of gender combined LRLs [50, 51]. The combined gender LRLs were considered to be reliable when the larger SD divided by the smaller SD <1.5 and when the normal deviate z value was less than the critical z value of 5. Additionally, the proportion criteria were used for non-normal distributions as it is more accurate than the Harris-Boyd method for such distributions. (Lahti A, 2004), [52] The combined gender LRLs were considered to be reliable when proportion of G6PD values less than the combined LRL did not exceed 4.1% in any gender subgroups [52].

The binomial test was used to compare the observed frequencies of male and female neonates. The two-sample t-test assuming equal variances or Mann–Whitney U test was used to compare differences between male and female neonates’ continuous variables when appropriate. The two-sample t-test assuming unequal variances was used to compare means of this study with those of previous studies. The one-sample Wilcoxon signed rank test was used to compare medians of this study with hypothetical values equal to medians of previous studies. A two-sided P value < .05 was considered statistically significant. Data analysis was performed using the RefVal 4.11, IBM SPSS Statistics 20 (Chicago, IL, USA), and OpenEpi 2.2.1 programs.

Results

Figure 1 depicts the process of study selection. A total of 463 term neonates (241 males and 222 females) of the reference tribe were born during the study period. Of these, 62 neonates (13.4%) did not meet the posteriori inclusion criteria (32 males and 30 females). The 2.5th percentile of the G6PD values of the excluded neonates was 9.6 U/g Hb (range: 7.3-19.6). Of the 401 neonates (86.6%) that met the posteriori inclusion criteria (209 males and 192 females), 2 males and 4 females had no quantitative G6PD testing on their cord blood as the reagents were not available at the time of birth. Thus, 395 term neonates with a similar proportion of males (n = 207) and females (n =188) left for analysis (P = .37). Mean (SD) gestational age of male and female neonates were similar [39.6 (1.3) versus 39.6 (1.4) weeks, P = .68). Mean (SD) birth weight of male and female neonates were similar [3315 (460) versus 3243 (442) grams, P = .12).

Figure 1
figure 1

Flow chart describing neonates’ selection.

The combined gender reference sample group had a mean (SD) birth weight of 3281 (452) grams and a mean gestational age of 39.6 (1.3) weeks.

The outliers

The G6PD values among the males had no significant skewness (0.07). The SPSS syntax indicated that the original data had the lowest skewness. Thus, Tukey’s boxplot was performed on their original values (Figure 2). Both the RefVal and SPSS programs detected exactly the same eight (3.9%) G6PD values as outliers among males (Figure 2).

Figure 2
figure 2

Tukey’s boxplot for outlier detection in male neonates. The open circles represent mild outliers which are the G6PD values fall between the Tukey’s inner and outer fences. The star represents extreme outlier which is the G6PD value fall outside the Tukey’s outer fences.

The G6PD values among the females had significant negative skewness (-0.93), thus, Tukey’s boxplot was performed on their Box-Cox transformed values (Method of Horn). The SPSS program indicated that a lambda of 1.84 is the best. The Box-Cox transformed G6PD values had a skewness of -0.003. Both the RefVal and SPSS programs detected exactly the same12 (6.4%) G5PD values as outliers among the females (Figure 3).

Figure 3
figure 3

Tukey's boxplot on box-cox transformed G6PD values of female neonates (method of Horn for outlier detection). The open circles represent mild outliers which are the G6PD values fall between the Tukey’s inner and outer fences.

Some of the outliers were separated by a gap from the rest of the G6PD values (see below). Assuming that the frequency of mutated G6PD alleles among males is 0.039 (number of detect outliers among males/total number of males), the Hardy-Weinberg equation revealed that the expected number of both heterozygous and homozygous females would be 14, nearly the same as number of the detected outliers (P = .87).

Estimation of LRLs in presence of the outliers

The histograms (Figures 4 and 5), measures of central tendency, and all statistical tests of normality (Table 1) suggested that the G6PD values of males, females and the combined group exhibited non-normal distributions. The G6PD values of males had one gap in each ends. The gap in lower end bounded by 9.10 and 7.90 U/g Hb and the gap in upper end bounded by 22.80 and 20.80 U/g Hb (Figure 4). The G6PD values of females had one gap in the lower end only which bounded by 9.60 and 7.50 U/g Hb (Figure 4).

Figure 4
figure 4

Distributions of quantitative cord G6PD values in male and females in presence of the outliers with superimposed best-fitting normal distribution curves.

Figure 5
figure 5

Distributions of quantitative cord G6PD values in the combined gender group in presence of the outliers with superimposed best-fitting normal distribution curve

Table 1 Characteristics of G6PD values and lower reference limits (2.5th percentiles) in presence of outliers

Table 1 depicts that G6PD levels in male and female neonates were similar (P = .50). It depicts the estimated LRLs in presence of the outliers. The RefVal indicated that even the two-stage transformation of the combined group did not yield a normal distribution. As a result, it recommended using the bootstrap method for the combined LRLs. The combined LRLs were reliable as the proportion of G6PD values less than the common LRL (10.0 U/g Hb) did not exceed 4.1% among male or female neonates, the ratio of SDs < 1.5 and the normal deviate z test < the critical z value of 5 (Table 1).

Estimation of LRLs in absence of the outliers

After the outliers were excluded, LRLs were estimated from cord blood samples of 375 neonates (199 males and 176 females). The histograms (Figures 6 and 7) and measures of central tendency (Table 2) suggested that the G6PD values of males, females and the combined group all exhibited normal distributions. This was confirmed by the non-statistically significant coefficients of skewness and kurtosis and by the Anderson-Darling and Shapiro-Wilk tests (Table 2). As expected, the observed gaps in presence of the outliers disappeared.

Figure 6
figure 6

Distributions of quantitative cord G6PD values in male and females in absence of the outliers with superimposed best-fitting normal distribution curves.

Figure 7
figure 7

Distributions of quantitative cord G6PD values in the combined gender group in absence of the outliers with superimposed best-fitting normal distribution curve.

Table 2 Characteristics of G6PD values and lower reference limits (2.5th percentiles) in absence of the outliers

Table 2 depicts that G6PD levels in male and female neonates were similar (P = .62). Estimated parametric LRLs had the tightest 95% CI, followed by bootstrap and non-parametric LRLs (Table 2). The combined LRLs were reliable as the ratio of SDs < 1.5 and the normal deviate z test < the critical z value of 5 (Table 2). As a result, the suggested LRL is 11.0 U/g Hb.

Table 3 depicts LRL estimations based on methods of previous studies. They were nearly the same, lower, or higher than those estimations based on statistical analyses of the IFCC-CLSI (Table 3). The observed mean and median G6PD values were nearly the same, lower, or higher than those of 8 previous studies (Table 3).

Table 3 Summary of 8 previous studies on the lower reference limit (LRL)/lower decision limit (LDL) of cord G6PD activity and our LRL/LDL estimations based on methods of those studies

Discussion

The present study estimated the LRLs of quantitative cord G6PD using standard statistical methods endorsed by the IFCC-CLSI for establishing RIs [24, 25]. The study used a homogenous reference sample group of healthy term neonates expected to have a low prevalence of G6PD deficiency. The results showed that 11.0 U/g Hb is a reliable combined parametric LRL for both genders. This value may help to identify G6PD-deficient heterozygous female neonates [28, 53]. This LRL can be transferred to other laboratories using a similar G6PD measurement method, satisfying the conditions for transference validation.

Reference limits are the cornerstones of interpretation of quantitative cord G6PD activity measurements. Standardization of statistical analysis for RI estimation is critical and is enabled by the IFCC-CLSI’s endorsement of a standard statistical analysis for RI estimation [24, 25]. Yet, the present study is the first to implement this standard method of statistical analysis to estimate one-sided RIs of quantitative cord G6PD. It demonstrated clearly that different methods are yielding different LRL estimations as most of our estimations based on methods of previous studies were inconsistent with our estimations based on method of the IFCC-CLSI (Table 3). Thus, standardized statistical analysis and reporting of studies on LRL estimations should be encouraged.

Our observed gaps in G6PD values of males confirmed the findings of others [8, 13, 16, 54]. Apparently, presence of these gaps depends on type of G6PD mutation and age of RBCs [55]. An old Italian study found that 98-100% of RBC of 115 normal schoolboys (G6PD-B) were G6PD (+), 98-100% of RBCs of 45 schoolboys bearing G6PD-Meditrianian were G6PD (-), and six schoolboys bearing G6PD-Seattle-like had a mosaic population of G6PD (+) and (-) RBCs [54]. For that reason, gaps were observed between G6PD values of G6PD-Meditrianian and normal boys while no gap was observed between G6PD-Seattle-like and normal boys [54]. Having said that, our observed gaps developed as their outermost borders and afterward values might be due to pre-analytic errors.

Our observed gap in the cord G6PD values of females has not been reported before. It is rather contradicting that G6PD values of females are continuum as heterozygotes may express a spectrum of values; normal, intermediate, and deficient [16, 56]. The intermediate values will bridge the gap between normal and deficient values. That is why it is difficult to diagnosis heterozygotes without a family history or molecular testing [5659]. Thus, this observed gap developed by chance or because their outermost borders and afterward values might be due to pre-analytic errors. However, excessive skewing of X-chromosome inactivation (allele ratios 3:1) in heterozygotes might account partly for our observed gap in G6PD values of females [60]. The excessive skewing has been reported to be in favor of RBCs bearing normal G6PD gene [54]. As a result, G6PD values for heterozygotes would be within RI of non-G6PD deficient females. On other hand, the skewing could be in favor of abnormal RBCs, and then G6PD values for heterozygotes would be within RI of homozygous G6PD-deficient females [60, 61]. Reported incidence of excessive skewing of X-chromosome in cord blood of healthy, term, female neonates is 9- 24% [62, 63]. Nevertheless, as the observed gap could be due to pre-analytic errors or chance and has no plausible explanation, further confirmatory and exploratory studies are required before firm conclusions can be drawn.

Comprehensive normal distribution testing in the present study showed that the distribution of quantitative cord G6PD exhibited a normal distribution only in absence of the outliers. Previous studies that have run limited normal distribution testing diverge on the normal distribution of G6PD values. The G6PD values exhibited a normal distribution in the study by Boo et al. but a non-normal distribution in the study by Fok et al. [7, 8]. The fact that only full-term neonates were included in both the present study and that by Boo et al. might explain the normal distribution in those populations, while the inclusion of both preterm and term neonates in the study by Fok et al. might account for the non-normal distribution, as G6PD activity has been shown to vary with gestational age (Table 3) [15, 36, 37]. Some researchers based their LRL and LDL estimations on the mean without addressing normal distribution testing [13]. Taken together, these observations emphasize the importance of conducting comprehensive normal distribution testing or estimating non-parametric RIs that do not make any assumptions about the data distribution [26].

Outliers are known to skew statistical tests based on sample means and variances, making it strongly advisable to detect outliers before performing these tests. Thus, it was not surprising that presence of these outliers widened the 95% CIs and distorted the homogeneity and the normal distributions, particularly for the combined gender group where even the two-stage transformation did not yield a normal distribution. The IFCC-CLSI emphasizes retaining outliers unless they are due to human or laboratory errors or representing unhealthy individuals. Retrospective nature of study precludes ascertaining whether or not these outliers were due to human or laboratory errors. However, there were three evidences suggesting that the outliers in the present study might be due to G6PD mutations. First, the frequencies of these outliers were in Hardy-Weinberg proportions. Second, presence of gaps between some of these outliers and the rest of presumed normal G6PD values. Third, a recent study from Jeddah area of western Saudi Arabia that found all the 20 female neonates with a cord G6PD value 6.6 U/g Hb subjected to molecular testing were G6PD- deficient [64]. However, the Jeddah study did not subject female neonates with a cord G6PD value > 6.6 U/g Hb to molecular testing [64].

We are suggesting using estimated LRLs in absence of the outliers for three reasons. First, the outliers of the present study might be due to pre-analytic errors. Second, they might be due to G6PD mutations. Third, all estimated LRLs were higher in absence than in presence of the outliers. Increasing cut-off point of LRLs will increase their sensitivity at the expense of their specificity. This is quite acceptable as mislabeling G6PD-deficient neonates as normal is intolerable whereas mislabeling G6PD-normal neonates as deficient is tolerable [65]. The parametric method is preferable to bootstrap and non-parametric methods for normal distributions [49]. This study showed that parametric LRLs had tight 95% CIs. The reference sample group in the present study satisfied the assumptions of the Harris-Boyd method as the G6PD values exhibited a normal distribution and the proportions of males and females were similar [50]. The Harris-Boyd method showed that the LRL of the combined genders is reliable. The standard statistical tests, the two-sample t-test assuming equal variances and Mann–Whitney U test also showed that G6PD levels were similar in male and female neonates. This is logic and consistent with the fact that the LRL was estimated from non-G6PD deficient neonates, as normal males and females should have the same G6PD values [58]. Consequently, we suggested the parametric LRL of the combined gender group (11.0 U/g Hg) as the LRL for cord G6PD of term neonates. Future studies on the LRLs for preterm neonates using the IFCC-CLSI are warranted as G6PD activity has been shown to vary across gestational age categories [15, 36, 37].

Some laboratory professionals have advocated estimating the 5.0th percentile as the LRL when a one-sided RI is assumed [66]. We elected the 2.5th percentile as the LRL in compliance with the IFCC-CLSI guidelines. However, the 5.0th percentiles of the G6PD values in this study were similar to the 2.5th percentiles (data not shown).

The IFCC-CLSI allow transferal of a reference limit from one laboratory to another when pre-analytical factors and testing methods are comparable in both laboratories and when one of their 3 validation methods is satisfied [25]. One of these validation methods is the N-20 reference sample group method, in which 20 samples from carefully assessed reference individuals are tested. When no more than 2 of these 20 samples fall outside the transferred reference limit, it is statistically valid to accept the transferred reference limit. Various commercial G6PD kits are available (Table 3), but as few studies have addressed their comparability, [65] future studies using the CLSI guidelines for methodological comparison are warranted [67].

Observed mean and median G6PD values were statistically significantly different from those reported in previous studies. The discrepancy could be due to differences in the timing of these studies, the ethnicities of the reference sample groups, gestational ages, G6PD kits and their reagents, and statistical strategies for LRL estimation (Table 3). However, not all of these differences are clinically significant, as trivial differences were sometimes statistically significant due to large sample sizes (Table 3).

Controversy exists as to whether G6PD activity should be expressed as U/g Hb or per number of RBCs. The WHO and some researchers have endorsed expressing the G6PD activity per number of RBCs to per U/g Hb, as the hemoglobin content may fluctuate independently of the G6PD activity, particularly when RBCs are hypochromic [6, 12]. Other researchers have compared the two expressions using neonatal peripheral blood and found that they are perfectly correlated [35]. However, the correlation between the two expressions of cord G6PD activity still needs to be determined.

The present study has four limitations that should be noted. First, retrospective nature of the study precluded assuring that the grandmothers were form the reference tribe as their names were not printed in the birth notices. Nevertheless, we do not assume that prevalence of neonates of grandmothers who were not from the reference tribe was high enough to alter prevalence of G6PD deficiency in the reference tribe significantly. The intermarriages between the reference tribe and other tribes are not a familiar tradition. The intermarriages between Bedouin descents and urban descents rarely occurred. For instance, in the present study, 3 out of 222 (1.4%) female neonates were product of intermarriages between the reference tribe and other tribes that all were Bedouin descents. Second, reported LRLs were not cross-referenced with molecular testing, so a future study using molecular testing to verify these LRLs is warranted. Third, the G6PD activity was only expressed as U/g Hb, precluding comparison with other studies that expressed G6PD activity per number of RBCs [12]. Fourth, methods of Horn and Harris-Boyd have some limitations that have been addressed elsewhere; yet, these are the best available methods for assessing outliers and common LRL, respectively [51, 68, 69]. However, even the standard statistical tests, the two-sample t-test assuming equal variances and Mann–Whitney U test also showed that G6PD levels were similar in male and female neonates.

Conclusions

The present study estimated the LRLs of cord G6PD activity using the standard method of statistical analysis endorsed by the IFCC-CLSI. We believe that we have included healthy neonates as we used stringent inclusion criteria and we included a reference group that has a low prevalence of G6PD deficiency. The results showed that the LRL of cord G6PD activity for full-term neonates of both genders is 11.0 U/g Hb and that this value is transferable to other settings when the IFCC-CLSI requirements of transference are satisfied. Further studies are warranted on LRLs for preterm neonates using the statistical standards of the IFCC-CLSI, cross-referencing the LRL with molecular testing, comparing methods of expressing cord G6PD activity, and comparing various commercial kits. Standardized statistical analysis and reporting of studies on RI estimation should be encouraged.

Abbreviations

KAH:

King Abdulaziz Hospital

G6PD:

Glucose-6-phosphate dehydrogenase

LRL:

Lower reference limit

LDL:

Lower decision limit

RI:

Reference interval

IFCC:

International federation of clinical chemistry

CLSI:

Clinical laboratory standard institute

WHO:

World Health Organization

u/gHb:

Units/gram hemoglobin

RBCs:

Red blood cells

FST:

Fluorescence spot test

°C:

Degree celsius

EDTA:

Ethylenediaminetetraacetic acid

NADP:

Nicotinamide adenine dinucleotide phosphate

IQR:

Interquartile range

95% CI:

95% Confidence interval.

References

  1. Cappellini MD, Fiorelli G: Glucose-6-phosphate dehydrogenase deficiency. Lancet. 2008, 371 (9606): 64-74. 10.1016/S0140-6736(08)60073-2.

    Article  CAS  PubMed  Google Scholar 

  2. Minucci A, Moradkhani K, Hwang MJ, Zuppi C, Giardina B, Capoluongo E: Glucose-6-phosphate dehydrogenase (G6PD) mutations database: review of the “old” and update of the new mutations. Blood Cells Mol Dis. 2012, 48 (3): 154-165. 10.1016/j.bcmd.2012.01.001.

    Article  CAS  PubMed  Google Scholar 

  3. Yoshida A, Beutler E, Motulsky AG: Human glucose-6-phosphate dehydrogenase variants. Bull World Health Organ. 1971, 45 (2): 243-253.

    CAS  PubMed  PubMed Central  Google Scholar 

  4. Kaplan M, Hammerman C: Glucose-6-phosphate dehydrogenase deficiency: a hidden risk for kernicterus. Semin Perinatol. 2004, 28 (5): 356-364. 10.1053/j.semperi.2004.09.001.

    Article  PubMed  Google Scholar 

  5. Beutler E, Gaetani G, Kaloustain V, Niwa S, Pannich V, Sodeinde O: Glucose-6-phosphate dehydrogenase deficiency. WHO working group. Bull World Health Organ. 1989, 67 (6): 601-611.

    Google Scholar 

  6. Betke K, Brewer GJ, Kirkman HN, Luzzato L, Motulsky AG, Ramot B, Siniscalco M: Standardization of procedures for the study of glucose-6-phosphate dehydrogenase. Report of a WHO scientific group. World Health Organ Tech Rep Ser. 1967, 366: 1-53.

    Google Scholar 

  7. Boo NY, Ainoon BO, Ooi LH, Cheong SK, Haliza BM: Glucose-6-phosphate dehydrogenase enzyme activity of normal term Malaysian neonates of different ethnic origins. J Paediatr Child Health. 1994, 30 (3): 273-274. 10.1111/j.1440-1754.1994.tb00633.x.

    Article  CAS  PubMed  Google Scholar 

  8. Fok TF, Lau SP, Fung KP: Cord blood G-6-PD activity by quantitative enzyme assay and fluorescent spot test in Chinese neonates. Aust Paediatr J. 1985, 21 (1): 23-25.

    CAS  PubMed  Google Scholar 

  9. Reclos GJ, Schulpis KH, Gavrili S, Vlachos G: Evaluation of glucose-6-phosphate dehydrogenase activity in two different ethnic groups using a kit employing the haemoglobin normalization procedure. Clin Biochem. 2003, 36 (5): 393-395. 10.1016/S0009-9120(03)00034-1.

    Article  CAS  PubMed  Google Scholar 

  10. Ainoon O, Alawiyah A, Yu YH, Cheong SK, Hamidah NH, Boo NY, Zaleha M: Semiquantitative screening test for G6PD deficiency detects severe deficiency but misses a substantial proportion of partially-deficient females. Southeast Asian J Trop Med Public Health. 2003, 34 (2): 405-414.

    CAS  PubMed  Google Scholar 

  11. Azma RZ, Hidayati N, Farisah NR, Hamidah NH, Ainoon O: G6PD enzyme activity in normal term Malaysian neonates and adults using a OSMMR2000-D kit with Hb normalization. Southeast Asian J Trop Med Public Health. 2010, 41 (4): 982-988.

    CAS  PubMed  Google Scholar 

  12. Ardati KO, Bajakian KM, Mohammed AM, Coe EL: Glucose-6-phosphate dehydrogenase phenotypes in Bahrain: quantitative analysis and alectrophoresis characterization. Saudi Med J. 1995, 16 (2): 102-104.

    Google Scholar 

  13. Kaplan M, Hoyer JD, Herschel M, Hammerman C, Stevenson DK: Glucose-6-phosphate dehydrogenase activity in term and near-term, male African American neonates. Clin Chim Acta. 2005, 355 (1–2): 113-117.

    Article  CAS  PubMed  Google Scholar 

  14. Obasa TO, Adesiyun OO, Mokuolu OA, Ojuawo AI: Comparative analysis of glucose-6-phosphate dehydrogenase levels in pre-term and term babies delivered at University of Ilorin Teaching Hospital. Pediatric reports. 2012, 4 (1): e7-

    Article  PubMed  PubMed Central  Google Scholar 

  15. Riskin A, Gery N, Kugelman A, Hemo M, Spevak I, Bader D: Glucose-6-phosphate dehydrogenase deficiency and borderline deficiency: association with neonatal hyperbilirubinemia. J Pediatr. 2012, 161 (2): 191-196 e191. 10.1016/j.jpeds.2012.02.018.

    Article  CAS  PubMed  Google Scholar 

  16. Algur N, Avraham I, Hammerman C, Kaplan M: Quantitative neonatal glucose-6-phosphate dehydrogenase screening: distribution, reference values, and classification by phenotype. J Pediatr. 2012, 161 (2): 197-200. 10.1016/j.jpeds.2012.02.045.

    Article  CAS  PubMed  Google Scholar 

  17. Kaplan M, Algur N, Hammerman C: Intermediate values of glucose-6-phosphate dehydrogenase. J Pediatr. 2012, 161 (3): 571-

    Article  PubMed  Google Scholar 

  18. Riskin A, Kugelman A, Bader D: Reply. J Pediatr. 2012, 161 (3): 571-572.

    Article  Google Scholar 

  19. Solberg HE: International federation of clinical chemistry (IFCC), scientific committee, clinical section, expert panel on theory of reference values, and international committee for standardization in haematology (ICSH), standing committee on reference values. Approved recommendation (1986) on the theory of reference values. Part 1. The concept of reference values. J Clin Chem Clin Biochem. 1987, 25 (5): 337-342.

    CAS  PubMed  Google Scholar 

  20. PetitClerc C, Wilding P: International Federation of Clinical Chemistry (IFCC), Scientific Committee, Clinical Section. The theory of reference values. Part 2. Selection of individuals for the production of reference values. J Clin Chem Clin Biochem. 1984, 22 (2): 203-208.

    CAS  PubMed  Google Scholar 

  21. Solberg HE, PetitClerc C: International Federation of Clinical Chemistry (IFCC), Scientific Committee, Clinical Section, Expert Panel on Theory of Reference Values. Approved recommendation (1988) on the theory of reference values. Part 3. Preparation of individuals and collection of specimens for the production of reference values. J Clin Chem Clin Biochem. 1988, 26 (9): 593-598.

    CAS  PubMed  Google Scholar 

  22. Solberg HE, Stamm D: International Federation of Clinical Chemistry, Scientific Division: approved recommendation on the theory of reference values. Part 4. Control of analytical variation in the production, transfer and application of reference values. Eur J Clin Chem Clin Biochem. 1991, 29 (8): 531-535.

    CAS  PubMed  Google Scholar 

  23. Dybkaer R: International federation of clinical chemistry (IFCC)1),2) the theory of reference values. Part 6. Presentation of observed values related to reference values. J Clin Chem Clin Biochem. 1982, 20 (11): 841-845.

    CAS  PubMed  Google Scholar 

  24. Solberg HE: International Federation of Clinical Chemistry (IFCC), Scientific Committee, Clinical Section, Expert Panel on Theory of Reference Values (EPTRV), and International Committee for Standardization in Haematology (ICSH), Standing Committee on Reference Values. Approved Recommendation (1987) on the Theory of Reference Values. Part 5. Statistical treatment of collected reference values. Determination of reference limits. J Clin Chem Clin Biochem. 1987, 25 (9): 645-656.

    CAS  Google Scholar 

  25. Horowitz GL, Altaie S, Boyd JC, Ceriotti F, Garg U, Horn P, Pesce A, Sine HE, Zakowski J: Defining, establishing, and verifying reference intervals in the clinical laboratory; approved guideline vol. 28. 2008, Wayne, PA: Clinical and Laboratory Standards Institute (CLSI), (Corrected 2010): 72, Third

    Google Scholar 

  26. Horowitz GL: Estimating reference intervals. Am J Clin Pathol. 2010, 133 (2): 175-177. 10.1309/AJCPQ4N7BRZQVHAL.

    Article  PubMed  Google Scholar 

  27. Ioannidis JP: Why most published research findings are false. PLoS medicine. 2005, 2 (8): e124-10.1371/journal.pmed.0020124.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Minucci A, Giardina B, Zuppi C, Capoluongo E: Glucose-6-phosphate dehydrogenase laboratory assay: how, when, and why?. IUBMB Life. 2009, 61 (1): 27-34. 10.1002/iub.137.

    Article  CAS  PubMed  Google Scholar 

  29. Al-Abdi SY, Mousa TA, Al-Aamri MA, Ul-Rahman NG, Abou-Mehrem AI: Hyperbilirubinemia in glucose-6-phosphate dehydrogenase-deficient male newborns in Al-Ahsa Saudi Arabia. Saudi Med J. 2010, 31 (2): 175-179.

    PubMed  Google Scholar 

  30. Warsy AS, El-Hazmi MA: G6PD deficiency, distribution and variants in Saudi Arabia: an overview. Ann Saudi Med. 2001, 21 (3–4): 174-177.

    CAS  PubMed  Google Scholar 

  31. Al-Ali AK, Al-Mustafa ZH, Al-Madan M, Qaw F, Al-Ateeq S: Molecular characterization of glucose-6-phosphate dehydrogenase deficiency in the Eastern Province of Saudi Arabia. Clin Chem Lab Med. 2002, 40 (8): 814-816.

    Article  CAS  PubMed  Google Scholar 

  32. Gelpi AP: Glucose-6-phosphate dehydrogenase deficiency in Saudi Arabia: a survey. Blood. 1965, 25 (4): 486-493.

    CAS  PubMed  Google Scholar 

  33. Gelpi AP: Glucose-6-phosphate dehydrogenase deficiency, the sickling trait, and malaria in Saudi Arab children. J Pediatr. 1967, 71 (1): 138-146. 10.1016/S0022-3476(67)80246-4.

    Article  CAS  PubMed  Google Scholar 

  34. Gelpi AP: Glucose-6-phosphate dehydrogenase deficiency in Saudi Arabia. Bull World Health Organ. 1967, 37 (4): 539-546.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Kaplan M, Leiter C, Hammerman C, Rudensky B: Enzymatic activity in glucose-6-phosphate dehydrogenase-normal and -deficient neonates measured with a commercial kit. Clin Chem. 1995, 41 (11): 1665-1667.

    CAS  PubMed  Google Scholar 

  36. Mesner O, Hammerman C, Goldschmidt D, Rudensky B, Bader D, Kaplan M: Glucose-6-phosphate dehydrogenase activity in male premature and term neonates. Arch Dis Child Fetal Neonatal Ed. 2004, 89 (6): F555-F557. 10.1136/adc.2004.049148.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Ko CH, Wong RP, Ng PC, Li K, Chui KM, Yuen PM, Fok TF: Oxidative challenge and glucose-6-phosphate dehydrogenase activity of preterm and term neonatal red blood cells. Neonatology. 2009, 96 (2): 96-101. 10.1159/000205720.

    Article  CAS  PubMed  Google Scholar 

  38. Al-Abdi SY, Al-Aamri MA, Dabelah KI, Mousa TA, Al-Rahman NG, Bukhamsien HI: Associations between spontaneous preterm birth and maternal-newborn ABO blood phenotype pairs. Saudi Med J. 2012, 33 (6): 660-664.

    PubMed  Google Scholar 

  39. Ballard JL, Khoury JC, Wedig K, Wang L, Eilers-Walsman BL, Lipp R: New Ballard Score, expanded to include extremely premature infants. J Pediatr. 1991, 119 (3): 417-423. 10.1016/S0022-3476(05)82056-6.

    Article  CAS  PubMed  Google Scholar 

  40. Horn PS, Feng L, Li Y, Pesce AJ: Effect of outliers and nonhealthy individuals on reference interval estimation. Clin Chem. 2001, 47 (12): 2137-2145.

    CAS  PubMed  Google Scholar 

  41. Solberg HE: The IFCC recommendation on estimation of reference intervals. The RefVal program. Clin Chem Lab Med. 2004, 42 (7): 710-714.

    Article  CAS  PubMed  Google Scholar 

  42. Osborne JW: Improving your data transformations: applying the Box-Cox transformation. Practical Assessment, Research & Evaluation. 2010, 15 (12): 1-9.

    Google Scholar 

  43. Dudek RW: BRS Genetics, 1st edn. 2010, China: Wolters Kluwer Health: Lippincott Williams & Wilkins

    Google Scholar 

  44. Edwards AWF: G. H. Hardy (1908) and Hardy–Weinberg Equilibrium. Genetics. 2008, 179 (3): 1143-1150. 10.1534/genetics.104.92940.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Wang J, Shete S: Testing departure from Hardy-Weinberg proportions. Statistical Human Genetics-Methods in Molecular Biology. Volume 850. Edited by: Elston RC, Satagopan JM, Sun S. 2012, New York, NY, USA: Humana Press, 77-102. 1

    Google Scholar 

  46. Solberg HE: Statistical treatment of reference values in laboratory medicine: testing the goodness-of-fit of an observed distribution to the Gaussian distribution. Scand J Clin Lab Invest Suppl. 1986, 184: 125-132.

    CAS  PubMed  Google Scholar 

  47. Razali NM, Wah YB: Power comparisons of Shapiro-Wilk, Kolmogorov-Smirnov, Lilliefors and Anderson-Darling tests. Journal of Statistical Modeling and Analytics. 2011, 2 (1): 21-33.

    Google Scholar 

  48. Cramer D, Howitt DL: The SAGE Dictionary of Statistics: A Practical Resource for Students in the Social Sciences. 2004, London, UK: Sage Publications Ltd, 1

    Book  Google Scholar 

  49. Pavlov IY, Wilson AR, Delgado JC: Resampling approach for determination of the method for reference interval calculation in clinical laboratory practice. Clin Vaccine Immunol. 2010, 17 (8): 1217-1222. 10.1128/CVI.00112-10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Harris EK, Boyd JC: On dividing reference data into subgroups to produce separate reference ranges. Clin Chem. 1990, 36 (2): 265-270.

    CAS  PubMed  Google Scholar 

  51. Lahti A: Are the common reference intervals truly common? Case studies on stratifying biochemical reference data by countries using two partitioning methods. Scand J Clin Lab Invest. 2004, 64 (4): 407-430. 10.1080/00365510410006027.

    Article  CAS  PubMed  Google Scholar 

  52. Lahti A, Petersen PH, Boyd JC, Rustad P, Laake P, Solberg HE: Partitioning of nongaussian-distributed biochemical reference data into subgroups. Clin Chem. 2004, 50 (5): 891-900. 10.1373/clinchem.2003.027953.

    Article  CAS  PubMed  Google Scholar 

  53. Watchko JF: Screening for glucose-6-phosphate dehydrogenase deficiency in newborns-practical considerations. J Pediatr. 2012, 161 (2): 179-180. 10.1016/j.jpeds.2012.04.001.

    Article  PubMed  Google Scholar 

  54. Rinaldi A, Filippi G, Siniscalco M: Variability of red cell phenotypes between and within individuals in an unbiased sample of 77 heterozygotes for G6PD deficiency in Sardinia. Am J Hum Genet. 1976, 28 (5): 496-505.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. Luzzatto L, Poggi V: Glucose-6-phosphate dehydrogenase deficiency. Nathan and Oski’s hematology of infancy and childhood. Edited by: Orkin SH, Nathan DG, Ginsburg D, Look AT, Fisher DE SELIV. 2008, Canada: Saunders Elsevier, 883-907. 7

    Google Scholar 

  56. Kaplan M, Hammerman C: Neonatal screening for glucose-6-phosphate dehydrogenase deficiency: biochemical versus genetic technologies. Semin Perinatol. 2011, 35 (3): 155-161. 10.1053/j.semperi.2011.02.010.

    Article  PubMed  Google Scholar 

  57. Luzzatto L: Glucose 6-phosphate dehydrogenase deficiency: from genotype to phenotype. Haematologica. 2006, 91 (10): 1303-1306.

    CAS  PubMed  Google Scholar 

  58. Beutler E: Glucose-6-phosphate dehydrogenase deficiency: a historical perspective. Blood. 2008, 111 (1): 16-24. 10.1182/blood-2007-04-077412.

    Article  CAS  PubMed  Google Scholar 

  59. May J, Meyer CG, Großterlinden L, Ademowo OG, Mockenhaupt FP, Olumese PE, Falusi AG, Luzzatto L, Bienzle U: Red cell glucose-6-phosphate dehydrogenase status and pyruvate kinase activity in a Nigerian population. Trop Med Int Health. 2000, 5 (2): 119-123. 10.1046/j.1365-3156.2000.00529.x.

    Article  CAS  PubMed  Google Scholar 

  60. Burton C, Kaczmarski R: Glucose 6-phosphate dehydrogenase deficiency and oxidative hemolysis. Clinical Studies in Medical Biochemistry. Volume 1. Edited by: Glew RH, Rosenthal MD. 2007, New York: Oxford University Press, 123-133. Third

    Google Scholar 

  61. Hsia YE, Miyakawa F, Baltazar J, Ching NS, Yuen J, Westwood B, Beutler E: Frequency of glucose-6-phosphate dehydrogenase (G6PD) mutations in Chinese, Filipinos, and Laotians from Hawaii. Hum Genet. 1993, 92 (5): 470-476. 10.1007/BF00216453.

    Article  CAS  PubMed  Google Scholar 

  62. Busque L, Mio R, Mattioli J, Brais E, Blais N, Lalonde Y, Maragh M, Gilliland D: Nonrandom X-inactivation patterns in normal females: lyonization ratios vary with age. Blood. 1996, 88 (1): 59-65.

    CAS  PubMed  Google Scholar 

  63. El Kassar N, Hetet G, Brière J, Grandchamp B: X-Chromosome inactivation in healthy females: incidence of excessive lyonization with age and comparison of assays involving DNA methylation and transcript polymorphisms. Clin Chem. 1998, 44 (1): 61-67.

    CAS  PubMed  Google Scholar 

  64. Jarullah J, AlJaouni S, Sharma MC, Busha MSJ, Kamal MA: Detection of glucose-6-phosphate dehydrogenase deficiency in heterozygous Saudi female neonates. Enzyme Engineering. 2012, 1 (105): 1-3.

    Google Scholar 

  65. Kaplan M, Leiter C, Hammerman C, Rudensky B: Comparison of commercial screening tests for glucose-6-phosphate dehydrogenase deficiency in the neonatal period. Clin Chem. 1997, 43 (7): 1236-1227.

    CAS  PubMed  Google Scholar 

  66. Boyd JC: Defining laboratory reference values and decision limits: populations, intervals, and interpretations. Asian journal of andrology. 2010, 12 (1): 83-90. 10.1038/aja.2009.9.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Krouwer JS, Tholen DW, Garber CC, Goldschmidt HMJ, Kroll MH, Linnet K, Meier K, Robinowitz M, Kennedy JW: Method Comparison and Bias Estimation Using Patient Samples; Approved Guideline. vol. 30. 2010, Wayne, PA: Clinical and Laboratory Standards Institute (CLSI), Second Edition (Interim Revision)

    Google Scholar 

  68. Solberg HE, Lahti A: Detection of outliers in reference distributions: performance of Horn’s algorithm. Clin Chem. 2005, 51 (12): 2326-2332. 10.1373/clinchem.2005.058339.

    Article  CAS  PubMed  Google Scholar 

  69. Henny J: The IFCC recommendations for determining reference intervals: strengths and limitations. LaboratoriumsMedizin. 2009, 33 (2): 45-51.

    Article  CAS  Google Scholar 

  70. Reclos GJ, Hatzidakis CJ, Schulpis KH: Glucose-6-phosphate dehydrogenase deficiency neonatal screening: preliminary evidence that a high percentage of partially deficient female neonates are missed during routine screening. J Med Screen. 2000, 7 (1): 46-51. 10.1136/jms.7.1.46.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Financial/funding support

The language editing of this manuscript by American Journal Experts was financed by the King Abdullah International Medical Research Center.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Sameer Yaseen Al-Abdi.

Additional information

Competing interests

We do not have any conflicts of interest to disclose.

Authors’ contributions

SYA conceptualized and designed the study, double checked and cleaned data, performed the statistical analysis, drafted the initial manuscript, and approved the final manuscript as submitted. He had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. ASS contributed substantially to the designing the study, data collection, interpretation of this work, reviewed and revised, and approved the final manuscript as submitted. FLD contributed substantially to the designing the study, data collection, interpretation of this work, reviewed and revised, and approved the final manuscript as submitted. AAS contributed substantially to the designing the study, data collection, interpretation of this work, reviewed and revised, and approved the final manuscript as submitted. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Cite this article

Al-Abdi, S.Y., Alsaigh, A.S., Aldawoud, F.L. et al. Lower reference limits of quantitative cord glucose-6-phosphate dehydrogenase estimated from healthy term neonates according to the clinical and laboratory standards institute guidelines: a cross sectional retrospective study. BMC Pediatr 13, 137 (2013). https://doi.org/10.1186/1471-2431-13-137

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2431-13-137

Keywords