2. Screening status for LSDs genes
Among the 22,687 individuals screened, 21,296 tested negatives for LSD gene screening, 1,361 were carriers (24 carriers with two different disease variants), and 30 had initial positive results (Fig. 1). Among the initial positive cases, Fabry disease had the highest proportion (50.00%) and the highest initial positive rate (0.066%, 1/1,512), followed by Krabbe disease (26.67% and 0.035%, 1/2,836), glycogen storage disease type II (GSD II) (13.33% and 0.018%, 1/5,672), Niemann‒Pick disease (NPD) (6.66% and 0.009%, 1/11,344) and mucopolysaccharidosis type II (MPS II) (3.33% and 0.004%, 1/22,687). The overall initial positive rate for LSDs was approximately 0.13% (1/756) (Table 2).
Table 2
Results from the NBGS of the LSDs
Disease | Negative screening result | Positive screening result | Positive screening rate (1/n) | Diagnosed | Prevalence (1/n) | Excluded (Carrier) | Follow-up | FPR | PPV |
Fabry | 22672 | 15 | 1/1 512 | 5 | 1/4537 | 0 | 10 | NA | NA |
Krabbe | 22679 | 8 | 1/2 836 | 8 | 1/2 836 | 0 | 0 | 0% | 100.00% |
GSD II | 22683 | 4 | 1/5 672 | 1 | 1/22 687 | 3 | 0 | 0.0132% | 25.00% |
NPD | 22685 | 2 | 1/11 344 | 1 | 1/22 687 | 1 | 0 | 0.0044% | 50.00% |
MPS II | 22686 | 1 | 1/22 687 | 0 | NA | 0 | 1 | NA | NA |
All | 22657 | 30 | 1/756 | 15 | 1/1 512 | 4 | 11 | 0.0177% | 50.00% |
FPR, false positive rate = false positives/(false positives + true negatives) = newborns without the disease but with positive NBGS results in the related gene/all newborns without the disease; PPV, positive predictive value = true positives/(true positives + false positives) = patients with the disease, and positive NBGS results in related genes/all newborns with positive NBGS results in the disease-related gene. NA, not applicable.
NBGS identified a total of 1,367 pathogenic variants for LSDs carried by 1,344 newborns. The overall carrier rate was approximately 6.03% (1/17), with Krabbe disease having the highest carrier rate (2.56%, 1/39), followed by NPD (1.05%, 1/95) and GSD II (0.86%, 1/116). Among Krabbe disease carriers, the c.1901T > C variant is highly prevalent, accounting for up to 86.38%. The c.955C > G variant is more prevalent at 65.19% in carriers of the NPD-A/B type, and the most prevalent variant observed in the GSD II type is c.2132_2133delinsGG, accounting for approximately 23.08%. (Table S2).
3. Enzyme activity detection for LSDs.
The lysosomal enzyme assay kit is limited to detecting GALC, GAA, ASM, IDUA, and GLA. Lysosomal enzyme activity was detected in a total of 900 samples. Among these, 290 healthy newborns with negative results in traditional newborn screening were utilized to establish high-risk cutoff values (0.2 or 0.3 multiples of the median)7, 25 and 0.5–99.5 percentile reference ranges for LSDs enzyme activity (Table S3). The histogram and statistical data illustrating the distribution of enzymatic activity in healthy newborns are displayed in Figure S1. No pathogenic variants in the LSDs genes were detected in these samples. A cutoff value lower than the high-risk cutoff value with relevant gene variants confirmed the diagnosis.
There were 29 genomic screening positive samples for LSDs (excluding one case of MPS II without corresponding enzyme activity detection) (Fig. 1 and Table 3) and 38 positive samples for other diseases (Table S4). Additionally, 543 samples were LSDs carriers (220 Krabbe carriers, 153 glycogen GSD II carriers, 117 NPD-A/B carriers, 45 mucopolysaccharidosis type I (MPS I) carriers, and 8 individuals with more than one gene carrier). The characteristics of these newborns are shown in Table 1. The enzyme activity results for individuals with negative, carrier, and positive genomic screening results are illustrated in Fig. 2, which shows a distinct decrease in enzyme activity in positive samples compared with that in negative and carrier samples. Additionally, the carriers presented significantly lower enzyme activity than did the negative samples. There were no statistically significant differences in gestational age or birth weight (Figure S2).
Table 3
Genomic screening positive samples for LSDs.
No. | Disease | Gene | Genotype | Location | Gender | Pathogenetic | Enzyme activity | Status |
1 | Fabry disease | GLA | c.640-801G > A/. | IVS4/. | Male | P | 1.91 | Diagnosed |
2 | | | c.640-801G > A/. | IVS4/. | Male | P | 1.34 | Diagnosed |
3 | | | c.640-801G > A/. | IVS4/. | Male | P | 1.41 | Diagnosed |
4 | | | c.1067G > A/. | EX7E/. | Male | LP | 1.57 | Diagnosed |
5 | | | c.1067G > A/. | EX7E/. | Male | LP | 1.03 | Diagnosed |
6 | | | c.911G > C/. | EX6/. | Male | LP | 4.76 | Follow-up |
7 | | | c.911G > C/. | EX6/. | Male | LP | 4.27 | Follow-up |
8 | | | c.640-801G > A/- | IVS4/- | Female | P | 3.38 | Follow-up |
9 | | | c.640-801G > A/- | IVS4/- | Female | P | 4.32 | Follow-up |
10 | | | c.640-801G > A/- | IVS4/- | Female | P | 8.62 | Follow-up |
11 | | | c.640-801G > A/- | IVS4/- | Female | P | 3.43 | Follow-up |
12 | | | c.911G > C/- | EX6/- | Female | LP | 9.75 | Follow-up |
13 | | | c.911G > C/- | EX6/- | Female | LP | 8.58 | Follow-up |
14 | | | c.1072_1074delGAG/- | EX7E/- | Female | LP | 4.41 | Follow-up |
15 | | | c.593T > C/- | EX4/- | Female | VUS | 2.62 | Follow-up |
16 | Krabbe | GALC | c.1901T > C/ c.1912G > A | EX16/EX17E | Male | LP/LP | 0.74 | Diagnosed |
17 | | | c.1901T > C/ c.2041G > A | EX16/EX17E | Male | LP/LP | 0.21 | Diagnosed |
18 | | | c.1901T > C/ c.1901T > C | EX16/ EX16 | Male | LP/LP | 0.29 | Diagnosed |
19 | | | c.1901T > C/ c.1901T > C | EX16/ EX16 | Male | LP/LP | 0.34 | Diagnosed |
20 | | | c.1901T > C/ c.1901T > C | EX16/ EX16 | Female | LP/LP | 0.17 | Diagnosed |
21 | | | c.1901T > C/ c.1901T > C | EX16/ EX16 | Female | LP/LP | 0.24 | Diagnosed |
22 | | | c.1901T > C/ c.1901T > C | EX16/ EX16 | Female | LP/LP | 0.31 | Diagnosed |
23 | | | c.1901T > C/ c.1901T > C | EX16/ EX16 | Female | LP/LP | 0.95 | Diagnosed |
24 | GSD II | GAA | c.2237G > A, c.503G > A /- | EX16, EX2/- | Male | P/LP | 3.51 | Excluded |
25 | | | c.2132_2133delinsGG, c.1669A > T/- | EX11/EX20E | Male | LP/P | 3.80 | Excluded |
26 | | | c.1634C > T/ c.2815_2816delGT | EX15, EX12/- | Male | LP/LP | 0.63 | Excluded* |
27 | | | c.2132_2133delinsGG, c.1669A > T/- | EX15, EX12/- | Female | LP/LP | 0.36 | Diagnosed |
28 | NPD-A/B | SMPD1 | c.995C > G/ c.995C > G | EX2/EX2 | Female | LP/LP | 0.26 | Diagnosed |
29 | NPD-C | NPC1 | c.1351G > A, c.3734_3735delCT/- | EX9/EX24 | Female | LP/LP | 14.74 | Excluded |
30 | MPS II | IDS | c.817C > T/. | EX6/. | Male | LP/. | - | Follow-up |
*, The pseudo-defective alleles c.1726G > A and c.2065G > A were detected.
4. Prevalence and genetic profiles of LSDs.
Among the 30 newborns with positive results from the genomic screening for lysosomal diseases, 4 were excluded after familial verification and enzyme activity testing. Fifteen newborns were diagnosed, and 11 newborns were followed-up (Table 2 and Table 3). The combined birth incidence rate, which included Fabry disease, Krabbe disease, GSD II, NPD and MPS II, was 1/1,512, and the diagnosis rate was 0.066%.
There were 15 newborns identified as potential Fabry patients, with an overall incidence rate of 1/1,512. Among them, seven were male, with a male incidence rate of approximately 1/1,590, and eight were female, with a female incidence rate of approximately 1/1,391. Enzyme activity detection of GLA was detected in only five male patients, yielding an incidence rate of approximately 1/2,226, which is lower than the rate determined by the NBGS. The primary GLA variant was c.640-801G > A (7/15), followed by c.911G > C (4/15). Eight newborns were diagnosed with Krabbe disease, with an incidence rate of 1/2,782. The primary genotypes included homozygous variations of c.1901T > C (3/4) and compound heterozygous variations with c.1901T > C (1/4). There is one potential patient each for GSD II, NPD, and MPS II, leading to an incidence rate of 1/22,256 for each of these three diseases.
5. NBGS is more efficient than enzyme activity detection for LSDs.
After our statistical investigation, the false positive rates (FPRs) and positive predictive values (PPVs) for NBGSs of LSDs such as Krabbe disease, GSD II, and NPD were obtained, as shown in Table 2. Enzyme activity results in a lack of representativeness in female patients with Fabry disease26, 27, preventing the evaluation of PPV and FRP. Compared with previous studies, NBGS shows better potential than enzyme activity screening for lysosomal storage disorders, especially in Asian regions such as Japan, Taiwan, and China (Table 4)6, 7, 10, 28–35.
Table 4
FRP and PPV in NBGS + Enzyme activity detection group
| NBGS | Enzyme activity screening |
Disease | Screened newborns | FPR | PPV | Region | Screened newborns | FPR | PPV | Reference |
Krabbe | 22,687 | 0% | 100% | Shanghai of China | 50,108 | 0.046% | 28.13% | 9 |
| | | | New York of USA | 550,000 | 0.004% | 8.00% | 24 |
GSD II | 22,687 | 0.0132% | 25% | Shanghai of China | 50,108 | 0.584% | 1.01% | 9 |
| | | | Taiwan of China | 191,786 | 0.44% | 1.83% | 25 |
| | | | Taiwan of China | 402,281 | 1.03% | 0.62% | 26 |
| | | | Taiwan of China | 473,738 | 0.47% | 1.25% | 27 |
| | | | Taiwan of China | 132,538 | 0.82% | 0.37% | 28 |
| | | | Japan | 297,387 | 0.033% | 7.34% | 29 |
| | | | Japan | 103,204 | 0.215% | 1.33% | 30 |
| | | | North East Italy | 44,411 | 0.007% | 25.00% | 22 |
| | | | Missouri of USA | 43,702 | 0.023% | 44.44% | 7 |
| | | | Hungary | 40,024 | 0.14% | 14.06% | 31 |
| | | | Austria | 34,736 | 0.003% | 80.00% | 32 |
NPD-A/B | 22,687 | 0% | 100% | Shanghai of China | 50,108 | 0.002% | 83.33% | 9 |
| | | Austria | 34,736 | 0.003% | 0% | 32 |
| | | Hungary | 40,024 | 0.007% | 40.00% | 31 |
FPR, false positive rate = false positives / (false positives + true negatives); PPV, positive predictive value = true positives / (true positives + false positives).
Furthermore, when the FRP and PPV were compared between the genomic screening and enzyme activity screening methods for Krabbe, GSD II, and NPD in the NBGS + enzyme activity detection group, it was observed that the genomic screening FRP for LSDs was significantly lower than that of the enzyme activity detection (P < 0.05), whereas the PPV was significantly greater than that of the enzyme activity detection method (P < 0.05) (Table 5). This finding indicates that, compared with relying solely on lysosomal enzyme activity screening, incorporating NBGS has a more proactive role in detecting LSDs.
Table 5
FRP and PPV in NBGS + Enzyme activity detection group
|
Disease
|
Krabbe
|
GSD II
|
NPD-A/B
|
NBGS
|
Negative
|
290
|
290
|
290
|
Carrier
|
220
|
153
|
117
|
True positive
|
8
|
1
|
1
|
False positive
|
0
|
3
|
1
|
FPR
|
0%
|
0.67%
|
0.24%
|
PPV
|
100.00%
|
25.00%
|
50.00%
|
Enzyme activity detection
|
Negative
|
495
|
431
|
403
|
True positive
|
6
|
1
|
1
|
False positive
|
7
|
11
|
3
|
NBGS Negative
|
2
|
2
|
2
|
NBGS Carrier
|
5
|
9
|
1
|
False negative
|
2
|
0
|
0
|
FPR
|
1.37%
|
2.48%
|
0.74%
|
FPR (exclude NBGS carrier)
|
0.009%
|
0.009%
|
0.009%
|
PPV
|
46.15%
|
8.33%
|
25.00%
|
PPV (exclude NBGS carrier)
|
75.00%
|
33.33%
|
33.33%
|
P-value
|
FPR1 (NBGS vs Enzyme activity detection)
|
< 0.05
|
FPR2 (NBGS vs Enzyme activity detection excluded NBGS carrier)
|
NS
|
PPV1 (NBGS vs Enzyme activity detection)
|
< 0.05
|
PPV2 (NBGS Negative vs Enzyme activity detection excluded NBGS carrier)
|
NS
|
FPR, false positive rate = false positives / (false positives + true negatives); PPV, positive predictive value = true positives / (true positives + false positives).