Abstract

Dupilumab (DUP) is the first biological agent used treating atopic dermatitis (AD). Notwithstanding its high cost, the type of patient group for which the drug is effective remains unclear. In this retrospective study, we aimed to identify novel and reliable biomarkers which can be measured before DUP administration and to predict the efficacy of DUP. Serum samples from 19 patients with AD treated with DUP were analysed by metabolome analysis using gas chromatography–mass spectrometry. Total 148 metabolites were detected, and the relative values of the metabolites were compared between the patient group that achieved 75% improvement in Eczema Area and Severity Index 16 weeks after administration of DUP (high responders: HR; n = 11) and that did not (low responders: LR; n = 8). The HR and LR groups had significant differences in the relative values of the eight metabolites (lactic acid, alanine, glyceric acid, fumaric acid, nonanoic acid, ribose, sorbitol, and ornithine), with ribose emerging as the best. Furthermore, we evaluated the serum concentrations of ribose and found that ribose may be a useful metabolite biomarker for predicting the efficacy of DUP in AD.

1. Introduction

Atopic dermatitis (AD) is characterised by an itchy eczema with recurrent exacerbations and relapses. The Eczema Area and Severity Index (EASI) score is used to measure the severity of AD, with higher scores reflecting worse severity. EASI-75 responders are participants who have achieved 75% improvement in EASI score from the baseline. Dupilumab (DUP) is a human monoclonal antibody that blocks the signalling of interleukin-4 and interleukin-13 and has been used for the treatment of AD in Japan from April 2018. Although DUP treatment is effective in several patients with AD in the clinical practice, a few do not respond adequately to it. In the CHRONOS trial [1], 69% patients achieved EASI-75 16 weeks after DUP administration; however, the remaining 31% did not, and the factors predicting the efficacy of DUP were unknown [1]. Another clinical trial also reported that DUP had similar effects regardless of eosinophil count, thymus and activation-regulated chemokine (TARC) and IgE at baseline [2], and there are currently no reliable predictors. The high cost of DUP treatment is a major concern in the choice of treatment for AD because patients typically pay 30% of their medical expenses in Japan. Therefore, identifying robust biomarkers to predict the efficacy of DUP prior to treatment is essential. Recently, metabolome analysis has been attracting attention and applied to the search for biomarkers for diagnosis and treatment, especially in the field of cancer. However, there have been few studies using metabolome analysis in the field of dermatology. In this retrospective study, we aimed to identify the biomarkers that predict the efficacy of DUP in patients with AD using metabolome analysis.

2. Materials and Methods

2.1. Study Participants and Data Collection

This study included 19 male patients with typical features of AD out of 109 patients with AD who commenced DUP at our hospital between April 2018 and July 2020. Patients with prurigo-type AD, associated with predominantly excoriated nodular lichenified lesions [3], were excluded because the surface area of the skin rash was small and the EASI score cannot be accurately measured. This study complied with the guidelines of the Declaration of Helsinki and was approved by the Ethics Committee of Hyogo Medical University (permission number 0212). Written informed consent was obtained from all the participants prior to their enrolment in the study. The patients satisfied all the following criteria: (1) they were under appropriate treatment, with topical steroids (strong class or higher) or topical calcineurin inhibitors at least for the last 6 months, (2) Investigator’s Global Assessment score was ≥3, (3) EASI score was ≥16 or EASI score of head and neck was ≥2.4, (4) body surface area score was ≥10%. The patients received 600 mg of DUP as a loading dose and 300 mg every 2 weeks thereafter, in addition to topical steroids or calcineurin inhibitors during DUP treatment. The topical applications and emollients were used according to the proactive therapy recommended by the Japanese AD guidelines [4]. Clinical and laboratory data were collected at baseline (before DUP treatment). Disease severity was determined using the EASI and Dermatology Life Quality Index (DLQI) at baseline and 16 weeks after DUP administration. Serum samples for metabolome analysis were obtained in the morning under fasting conditions at baseline and 16 weeks. In this study, the primary goal of AD treatment was to achieve EASI-75. We defined patients who achieved and did not achieve EASI-75 at 16 weeks as high responders (HR) and low responders (LR), respectively. In addition, we defined patients who achieved and did not achieve EASI-90 at 16 weeks as super high responders (SHR) and non-super-high responders (non-SHR), respectively. A secondary goal of AD treatment was to achieve a 4-point or greater improvement in the DLQI from baseline. Patients who achieved 4-point or greater improvement at 16 weeks were defined as DLQI-high responders (DLQI-HR) and those who did not as DLQI-low responders (DLQI-LR).

3. Materials

2-Isopropylmalic acid (2-IPMA), methoxyamine hydrochloride, and D-ribose were purchased from Sigma–Aldrich Japan (Tokyo, Japan). N-Methyl-N-trimethylsilyltrifluoroacetamide (MSTFA) was purchased from GL Sciences, Inc. (Tokyo, Japan). A standard alkane series mixture (C7–C33) was obtained from Restek Co. (PA, USA). [U-13C5]-labelled D-ribose was purchased from Cambridge Isotope Laboratories, Inc. (Andover, MA, USA).

3.1. Pretreatment for Gas Chromatography–Mass Spectrometry (GC/MS)-Based Metabolome Analysis

For extraction of metabolites from the serum, 50 µL of each sample was mixed with 250 µL of methanol and 10 µL of 0.5 mg/mL 2-IPMA, and then shaken at 1,200 rpm for 30 min at 37°C before being centrifuged at 20,600 × g for 5 min at 4°C. Next, 225 µL of the supernatant obtained from each mixture was added to a clean tube, centrifuged, and freeze-dried overnight. For oximation, 60 µL of 20 mg/mL methoxyamine hydrochloride dissolved in pyridine was added into the tube, and then sonicated for 20 min before being shaken at 1,200 rpm for 90 min at 30°C. The mixture was centrifuged at 20,600 × g for 5 min at 20°C, and then 40 µL of the resultant supernatant was subjected to GC/MS as described in the “GC/MS analysis” section.

3.2. GC/MS Analysis

GC/MS analysis was performed using an AOC-6000 (Shimadzu Co., Kyoto, Japan) and GCMS-TQ8040 (Shimadzu Co.) equipped with a BPX-5 capillary column (internal diameter: 30 m × 0.25 mm; film thickness: 0.25 µm; SEG, Victoria, Australia). In the AOC-6000, 20 µL of MSTFA were added to the sample supernatant, and then the mixture was incubated at 750 rpm for 30 min at 37°C before 1.0 µL of the derivatised solution was injected into the GCMS-TQ8040. GC/MS analysis was performed using the Smart Metabolites Database Ver.2 (Shimadzu, Co.), which contains information regarding the GC analytical conditions, multiple reaction monitoring (MRM) parameters, and the retention index employed for metabolite measurement. To correct the retention time, the Automatic Adjustment of Retention Time (AART) function of the GCMSsolution software (Shimadzu Co.) and a standard alkane series mixture (C7–C33) were used. Peak identification was performed automatically and confirmed manually based on the specific precursor and productions, and the retention time.

3.3. Quantitative Analysis of D-Ribose

D-ribose was quantitatively analysed using the internal standard method with D-ribose and the corresponding stable isotope, [U-13C5]-labelled D-ribose. In brief, 50 µL of each serum sample was mixed with 270 µL of methanol containing [U-13C5] -labelled D-ribose. As pre-treatment for the GC/MS analysis, the mixture was treated according to the protocols described in the “Pretreatment for GC/MS-based metabolome analysis” section. The analysis was performed as described in the “GC/MS analysis” section. The serum concertation (µM) of D-ribose was calculated based on the multipoint calibration curves derived from the peak area (y-axis; native/isotope) and concentration (x-axis; native/isotope) ratios.

3.4. Statistical Analysis

The normality of the data distribution was verified using Shapiro–Wilk test. Student’s t-test and Mann–Whitney U test were used to compare means for normally and not-normally distributed values, respectively. Student’s t-test was used to compare the relative values of metabolites at baseline analysed using metabolome analysis. Based on the metabolites that exhibited significant differences, principal component analysis (PCA) was performed to differentiate between the HR and LR samples. For the metabolites displaying significant differences, logistic regression and receiver operating characteristic (ROC) analyses were performed to estimate the area under the curve (AUC). The serum concentration of the metabolite with the highest AUC was then provided for logistic regression and ROC curve analyses to estimate its performance as a predictive biomarker. The optimal cut-off value was determined using the Youden index. Changes in the EASI at 16 weeks from baseline were calculated as percentages. The change in EASI was compared based on the predicted probability using Student’s t-test. The relative values of the metabolites at baseline were also compared between the SHR and non-SHR groups and between the DLQI-HR and DLQI-LR groups using Student’s t-test and Mann–Whitney U test. Changes in the relative values of metabolites at 16 weeks from baseline were calculated as percentages. Spearman’s R correlation was used to measure the association between changes in the EASI and relative values of the metabolites. Statistics were calculated using GraphPad Prism 9 and JMP Pro 16.

4. Results

Table 1 presents the characteristics of the patients. Eleven and eight patients were included in the HR and LR groups, respectively. In total, 148 metabolites were detected in the serum samples using metabolome analysis. Among these, the relative values of eight metabolites (lactic acid, alanine, glyceric acid, fumaric acid, nonanoic acid, ribose, sorbitol, and ornithine) exhibited significant differences at baseline between the HR and LR groups by the Student’s t-test (Table 2, Figure 1; , respectively). PCA revealed that HR and LR could be distinguished based on the relative values of the eight metabolites, and ribose was the metabolite that characterised HR (Figure 2). Logistic regression and ROC analyses were performed using the relative values of the eight metabolites. The AUC of ribose (0.920) was the highest among the eight metabolites (Table 3). To examine the results of ribose in greater detail, ribose was quantitatively analysed, and the calculated serum concentration was used for logistic regression and ROC analyses (Table 4). The AUC value was 0.898, and the predicted probability cut-off value of 0.855 in the ROC curve yielded a sensitivity of 73% and a specificity of 100% (Figure 3). There was a significant difference in the change of EASI at 16 weeks from baseline between patients with predicted probabilities of ≥0.855 and <0.855 by the Student’s t-test (Figure 4; ). Four and 15 patients were included in the SHR and non-SHR groups, respectively. Among the eight metabolites, the relative values of two metabolites (glyceric acid and fumaric acid, but not ribose; , 0.0145, 0.5955) displayed significant differences between the SHR and non-SHR groups by the Student’s t-test and Mann–Whitney U test. Fifteen and four patients were included in the DLQI-HR and DLQI-LR groups, respectively. The relative values of all eight metabolites did not significantly differ between DLQI-HR and DLQI-LR groups by the Student’s t-test and Mann–Whitney U test. All four patients in the DLQI-LR group were included in the HR group and had a relatively low-baseline DLQI of less than 7. No correlation was found between the change in EASI and the relative values of ribose at 16 weeks from baseline (r = −0.0525).

5. Discussion

In this study, eight metabolites exhibited significant differences between HR and LR groups before the administration of DUP, and logistic regression and ROC analyses suggested that ribose could be the best candidate biomarker. The serum concentration of ribose was measured to apply the results to the clinical practice, and logistic regression and ROC analyses were reperformed to evaluate its ability as a biomarker. The results indicate that ribose could be a useful biomarker for predicting the efficacy of DUP. The reason for the efficacy of DUP in patients with high levels of ribose is unclear; however, ribose is associated with poly (ADP) ribose polymerase (PARP), which is involved in numerous cellular processes such as DNA repair, protein turnover, inflammatory regulation, ageing, and metabolic regulation [5]. Additionally, PARP may play a critical role in AD [6, 7]. A few studies have reported the use of metabolomics or lipidomics for AD; however, ribose has not yet been discussed [814]. The values of ribose in healthy subjects have not been clearly defined in previous studies, and they were not measured in this study. Comparing ribose levels in healthy subjects and AD patients may help us understand the role of ribose in AD and is an important topic for future research. Using metabolomics and lipidomics, Zhang et al. [15] demonstrated that DUP treatment for patients with AD affected the concentrations of specific metabolites and altered specific metabolic pathways associated with them, including glycerophospholipid metabolism, valine, leucine and isoleucine biosynthesis, citrate cycle, arachidonic acid metabolism, pyrimidine metabolism, and sphingolipid metabolism. [15]. Ribose is not included in these metabolic pathways, but may be slightly related to pyrimidine metabolism because it is a nucleotide component. Further investigation into the relationship between ribose and the mechanism of DUP treatment for patients with AD is required.

Metabolome analysis, an omics technology, is the comprehensive study of low-molecular-weight metabolites in the body, including saccharides, amino acids, and fatty acids. Metabolome data cover information based on environmental factors as well as information derived from genetic factors because metabolites are the downstream products of cellular processes, including deoxyribonucleic acid, ribonucleic acid, and proteins. Therefore, the metabolome, which is closer to the phenotype, has been recognised as useful for the detailed evaluation of final phenotypes, and metabolome analysis has been widely applied in biomarker research [16].

Nettis et al. [17] considered complete responders as patients who achieved all three domains (EASI-75, peak pruritus Numerical Rating Scale (PP-NRS) ≥ 4 point improvement, and DLQI ≥ 4 point improvement). Predictors of complete responders after 16 weeks of DUP treatment were lower age, low baseline PP-NRS score, early onset AD, history of allergic conjunctivitis, and complete responders after 4 weeks of DUP treatment. Ariëns et al. [18] reported that the baseline EASI was significantly higher among patients who achieved all the previously described three domains compared with those who did not after 16 weeks of DUP treatment. Ferrucci et al. [19] reported that the predictors of patients who achieved EASI-75 after 4 weeks of treatment were early AD onset and the absence of hypereosinophilia (>500 × 103/L), while no predictors achieved EASI-75 after 16 weeks of treatment. Linder et al. [20] defined patients who achieved a SCORing Atopic Dermatitis (SCORAD) score of 75 after 16 weeks as optimal responders, and their predictors were lower weight and less frequent exposure to tobacco and alcohol. Olesen et al. [21] found that the predictors of patients who achieved EASI-75 were female sex, lower age, lower baseline EASI, and lower baseline IgE after 4 weeks of treatment, and only female patients after 12 weeks. Lee et al. [22] found that the predictors for patients who received DUP treatment every 4 weeks and achieved EASI-75 after 16 weeks were low baseline eosinophil count and low baseline LDH level. Thus, although several researchers have reported factors predicting the efficacy of DUP, the results vary and consistent conclusions have not yet been reached.

The non-correlation between the changes in EASI and in the relative values of ribose suggests that ribose does not reflect disease activity during DUP treatment. Discrepancy was found between the changes in the EASI and DLQI in this study. The reason may be that patients with severe to moderate disease in this study had fewer subjective symptoms owing to chronic symptoms of AD, and a few in the HR group had a low-baseline DLQI.

This study has a few limitations. First, the number of patients was too small, and the biomarkers suggested in this study were exploratory. Second, this study was conducted using data from a single centre in Japan, and the results may differ according to race, ethnicity, or region. Third, only male patients were selected in this study to eliminate differences in metabolites between sexes. It is unclear whether similar results would be obtained in females.

6. Conclusions

Ribose may be a useful biomarker for predicting the effects of DUP. Our results provide clues for establishing personalised treatment regime for AD. Currently, several metabolites, including ribose, cannot be measured using routine blood tests. However, the means of measuring these metabolites will be established if their abilities as biomarkers are proven. The results of this study are exploratory, and a greater number of cases need to be investigated in the near future.

Data Availability

The data presented in this study can be made available by the corresponding author/s upon reasonable request.

Conflicts of Interest

YI has received honoraria for lectures from Sanofi and Maruho. SN belongs to the Department of Omics Medicine, Hyogo Medical University, which is funded by Shimadzu Co. The remaining authors declare no conflict of interest regarding the publication of this article.

Authors’ Contributions

YI, KY, and SN contributed in the conceptualisation; SM and MM contributed in the data curation; SM, YI, and SN contributed in the formal analysis; YI contributed in the funding acquisition; SM, YI, and MN contributed in the investigation; SM, YI, and SN contributed in the methodology; YI and KY contributed in the project administration; YI and NM contributed in the resources; NK contributed in the supervision; SM and YI contributed in the validation; SM and YI contributed in the visualisation; SM contributed in the writing–original draft preparation; YI, NK, and SN contributed in the writing–review and editing.

Acknowledgments

This work was supported in part by the JSPS KAKENHI 21K08337 for YI and 21K08313 for KY. In this study, AOC-6000, GCMS-TQ8040, and Smart Metabolites Database Ver.2 from Shimadzu Co. were used as research tools for metabolome analysis. We would like to thank Editage (http://www.editage.jp) for English language editing. This work was also supported by the Hyogo Medical University Diversity Grant for Research Promotion under MEXT Funds for the Development of Human Resources in Science and Technology, “Initiative for Realizing Diversity in the Research Environment (Characteristic-Compatible Type).”