Determination of metabolites for NBS by MS/MS
Screening for MMAs and PA was performed with a cutoff of C3 acylcarnitine ≥ 4.3 µmol/L or C3/C2 ratio ≥ 0.2, and both parameters were equally important for the initial NBS by MS/MS. From 2015 to 2020, 53 screen-positive individuals benefited from the expanded NBS program in Xuzhou, China, including 45 individuals with MMAs and 8 individuals with PA. Notably, 1 of 45 MMA individuals presented C3 and C3/C2 values below the established thresholds in the recall test and thus was technically a false negative in the NBS by MS/MS. This individual was subjected to second-tier metabolic analyses, and a significantly increased MMA concentration was identified by UPLC-MS/MS, which ultimately resulted in an isolated MMA diagnosis. To investigate MS/MS metabolic patterns, we determined the levels of three C3-associated NBS metabolic analytes in confirmed individuals and healthy infants (Fig. 1). Not unexpectedly, the concentrations of C3 showed a significant difference between the MMA/PA individuals and the control group individuals (p = 3.1e-29), but there was no significant difference between MMA individuals and PA individuals (p = 0.50). Significant differences were also identified for the C3/C2 ratio (p = 8.01e-21) between MMA/PA individuals and control group individuals. Met showed significant differences between MMA/PA subtypes, with relatively decreased Met levels in individuals with remethylation defects (cblC type) compared to those in individuals with propionyl-CoA carboxylase deficiency. However, Met levels were not significantly different between the isolated MMA group and the cblC type group (Fig. 1C). Unfortunately, there were almost no differences in other NBS analytes between neonates with MMA/PA and healthy infants, indicating that the other analytes were metabolically similar.
Biochemical and molecular genetic detection of MMA and PA individuals
Of the 53 confirmed individuals who were analyzed, multiple individuals carried compound heterozygous variants in pathogenic genes (Fig. 1D), resulting in a low percentage of homozygotes with only the cblC type in Xuzhou, China. DNA sequence analysis of 86 genes for NBS detected more variants in confirmed MMA/PA cases, in whom five associated genes (MUT, MMAB, MMACHC, PCCA, and PCCB) were identified. Forty-five MMA individuals (5 of whom had isolated MMA and 40 of which had cblC deficiency) and eight PA individuals had two P/LP variants or one P/LP variant and a variant of unknown significance (VUS) in a pathogenic gene. Eight different known MUT variants (c.323G > A, c.441T > A, c.581C > T, c.1106G > A, c.1219A > T, c.1741C > T, c.1677-1G > A and c.1880A > G) were identified in the 5 isolated MMA individuals from 5 unrelated families, and all the variants were located in exons, except for one splicing type variant. In addition, eleven variants were detected in the PCCA gene, and four variants were detected in the PCCB gene in 8 PA individuals, as shown in Supplementary Table 1. Indeed, not only isolated MMA individuals but also PA individuals were confirmed to carry compound heterozygous variants. We observed a strikingly high frequency of cblC-type individuals in our cohort (75.47%, n = 40/53), in which 80 independent cblC alleles were affected. The sequencing results confirmed the compound heterozygous or homozygous variants in the MMACHC gene in 40 infants, of which 7 were homozygous for c.609G > A (71.4%, 5/7), c.658_660del (14.3%, 1/7) and c.80A > G (14.3%, 1/7) (Supplemantal Table 1). To calculate metabolite levels, individuals homozygous for the MMACHC gene (patients 3, 9, 17, 30, 31, 33 and 40) were placed into one group for further analyses, with the other individuals placed in a second group. The mean C3/C2 ratio of the homozygous group obtained from the initially measured NBS data was 0.926, which was significantly different from that of the compound heterozygous group (C3/C2 ratio = 0.483). Correspondingly, the Met values in the homozygous group with MMACHC variants ranged from 3.34 to 11.01 µmol/L, with a mean value of 7.03 µmol/L, which was near the minimum value of the normal range of Met and lower than that in the compound heterozygous group. However, no significant difference in C3 values for individuals with cblC deficiency was found based on the biochemical data (Fig. 1).
Evaluation of NBS for MMAs and PA based on season, gestational age and birth weight
The levels of specific metabolites detectable (e.g., C2) by MS/MS for NBS are known to have seasonal variation. We assessed the seasonal profile of MMA/PA screening parameters (Fig. 2A/C/E) and found a significantly increased percentage of false negative results in newborns caused by low concentrations of C2 and Met in the summer than in other seasons from 2015 to 2020 (p < 0.05) (Fig. 2B). The relatively higher temperature and humidity in summer reduced the positive predictive value (PPV) from 3.86–1.83% with numerous false positive cases in Xuzhou, China. A previous study indicated that metabolic profiles may vary by gestational age (GA) or birth weight. We further compared GA between the MMA/PA-positive (MMA/PA-TP) group, MMA/PA-false-positive (MMA/PA-N) group and healthy control group (Fig. 2D). Within the MMA/PA-N group, there was a significantly higher proportion of preterm (GA ≤ 37 weeks) than full-term (> 37 weeks) births after analysis with the Mann-Whitney U test. An additional analysis of neonate birth weight (normal: 2,500-4,000 grams) revealed a relatively lower birth weight for MMA screening-positive newborns (both true positives and false positive newborns) than that for the healthy group (Fig. 2E).
To evaluate the sensitivity of MS/MS analytes for MMA/PA screening, receiver operating characteristic (ROC) curves were utilized to distinguish true and false positive MMA/PA individuals. The optimal screening efficiency was determined using a sensitivity of 0.957 and a specificity of 0.972 under the curve area with 0.996 for parameter C3 (Fig. 3A and B). However, the parameter C3/C2 had a better sensitivity and specificity than C3 with ROC analysis. The mean decrease in the accuracy (MDA) index (Fig. 3D) suggested that the C3/C2 ratio was a significant parameter for MMA/PA screening. However, other metabolic analytes, including the C0 and Met levels, showed similar rankings by MS/MS. Of the above several covariates analyzed, seasonal variation played a fundamental role in primary NBS. Ranking analysis of the MS/MS metabolites showed significant differences based on phenotypic subtype. Individuals with cobalamin deficiency might have lower Met levels than those with other phenotypes (Fig. 1C). Moreover, the C2 and Met levels declined sharply during the summer due to the higher air temperature in Xuzhou, and parameter C2 was crucial for identifying MMA/PA in the initial MS/MS screening (Fig. 3E/F).
Generalization of a novel second-tier approach: determination of cutoffs for MMA, MCA and HCY
From January 2019 to March 2020, 263 samples with increased C3 concentrations at the initial NBS test were further assessed by novel second-tier UPLC-MS/MS screening. As the dataset was not symmetrical, a total of 1,830 negative subjects were used to calculate the cutoff value based on the percentile distribution of metabolites. The statistical results of the percentile distribution of the negative controls are shown in Table 1. Notably, the value "4.54" of the MMA metabolite and the value "12.64" of the HCY metabolite were two suspicious outliers. According to the NCCLS C28-A3 file, D/R (MMA)= (4.54–3.96)/(5.55-0.00) = 0.10, while D/R (HCY)= (12.64–10.03)/(17.14–0.9) = 0.16; however, both D/R ratios were ≤ 1/3, and there was no outlier in any of the three metabolite indexes. Previous literature and clinical research have indicated that sex does not affect the three analytes that were measured with UPLC-MS/MS. Based on the distribution of the histogram chart, the 99.5th percentile and the 99.7th percentile were considered the cutoff values for the second-tier test (Table 1). When using the 99.5th percentile as the cutoff value, the positive rate of MMA/PA detection was estimated to be 0.45%, 0.40% and 0.76%, respectively. When the 99.7th percentile was reached, the positive rate of detection for each disorder was reduced, and the number of recall cases decreased from 43 to 29 (Table 2). A comprehensive analysis of several published laboratory parameters was performed, and as shown in Table 3, the 99.7th percentile was similar to the cutoff values established by the Mayo Clinic, USA, and the 99.5th percentile is used as the cutoff value in Nanjing, China. A strict test of the generalizability of the novel second-tier test should require confirmed MMA/PA individuals to be identified by this screening method. Considering the significantly increased concentrations of MMA, MCA and HCY in patients, the cutoff values for MMA, MCA and HCY for the second-tier test were 3.13, 0.29 and 9.77 µmol/L (the 99.7th percentile), respectively.
Here, we developed a following nonderivatized approach using UPLC-MS/MS instead of primal derivatized testing. Despite there being only 207 negative controls for nonderivatized detection for the second-tier test, the new cutoffs for MMA and HCY showed significantly increased performance (Supplementary Table 2 and Supplementary Figure II). For comparison of the distribution characteristics of MMA, MCA and HCY values before and after modification, the Mann-Whitney test was performed. As shown in Fig. 2, there were significant differences among the three indexes before and after experimental improvement. First, we compared the 10th, 25th, 75th and 90th percentile MMA, MCA and HCY values. The 75th percentile values of MMA and HCY with primal derivatized detection were similar to the 25th percentile values with the nonderivatized approach, and the difference between the 10th and 90th percentiles was greater than 30%; there were only small differences in the MCA values between the 10th and 90th percentiles for the second-tier test. Taken together, these data show that the distribution of MCA values from the derivatized and nonderivatized tests was similar, which allowed us to conclude that the previous cutoff values could be generalized. In a specific test of the generalizability of the cutoff values for MMA and HCY, however, the previous cutoff values were not equal to the novel nonderivatized test values.
Comparison of the second-tier test with NBS analysis for newborns suspected to have MMA/PA disorders
Thirteen individuals with confirmed MMA/PA disorders were identified by second-tier screening from January 2019 to March 2020, and their disease was confirmed by urinary GC/MS and DNA sequencing. As shown in Fig. 1C, the individuals were identified with a one-step second-tier test before genetic analysis. The C3 level, C3/C2 ratio and Met level were analyzed in patients grouped based on the initial MS/MS NBS results: the isolated MMA group (n = 1), the Cbl group (n = 6), the PA group (n = 2) and the healthy group (n = 1,830). As shown in Fig. 1A, C3, Met and C3/C2 were generally increased in the three groups compared to the healthy group, while the C3/Met ratio was not significantly different. Unfortunately, none of the differences between the three disease groups were significant. However, the second-tier UPLC-MS/MS test could directly distinguish these groups, and further statistical analysis of the parameters revealed significant differences (Fig. 4B, C). The mass spectrum peaks (shown in Fig. 4) revealed that a solely increased MMA peak likely indicated isolated MMA, an increased MMA peak combined with an increased HCY peak likely indicated Cbl deficiency, an increased MCA peak likely indicated PA, and other patterns likely indicated none of these disorders. Patients who fall within the uncertain category in initial NBS would be differentially screened using a second-tier test, as confirmation via genetic analysis and clinical follow-up is required, which means that this second-tier test can reduce the initial recall rate of patients suspected of having these disorders on NBS. The consideration of MMA, MCA and HCY as markers for MMA/PA would have eliminated 88.97% (234/263) of unnecessary referrals and recall tests.