Risk Factors and Predictive Model for Dermatomyositis Associated with Rapidly Progressive Interstitial Lung Disease

Introduction

The main manifestations of dermatomyositis (DM) are varying degrees of inflammation of the muscles, skin and lungs, and are closely related to the severity and prognosis of the disease.1 Interstitial lung disease (ILD) and malignancy are significant complications that determine the prognosis of patients with DM.2 Rapidly progressive ILD (RP-ILD) is a subtype of DM-ILD, which is rapidly progressive, challenging to treat, and has a high mortality rate compared to DM-ILD, and is critical in the clinical management of DM. A large cohort study of myositis in the United States showed that the 1-, 5- and 10-year survival rates of myositis patients with ILD were 97, 91 and 81%, respectively.3 The 6-month survival rate of patients with RP-ILD is 28%-66% after several aggressive treatments.4,5 Therefore, early diagnosis of RP-ILD can improve the efficiency of screening and diagnosis and provide meaningful guidance for carrying out early and aggressive treatment.6

Positive anti-melanoma differentiation-associated gene 5 (MDA5) antibody is the one antibody that is known to be highly associated with RP-ILD.7 About half of the patients with MDA5-positive DM develop RP-ILD.8,9 The relationship between other myositis-associated antibodies and RP-ILD is unknown. In addition to this, Myofascia-dominant involvement, elevated CRP, decreased lymphocyte count, elevated Soluble CD206, elevated interleukin-15 (IL-15), elevated carcinoembryonic antigen, and prolonged activated partial thromboplastin time (aPTT) ratio are also risk factors for the development of RP-ILD in DM patients.10–15 High plasma ferritin levels, generalized worsening of pulmonary infiltrates, ground glass shadow, and elevated serum neopterin levels are significantly associated with reduced survival in DM patients with RP-ILD.16,17 Studies on RP-ILD associated with juvenile DM, on the other hand, found that IL-18, KL-6, ferritin and anti-MDA5 antibodies were associated with RP-ILD.18 However, most of the risk factors reported in these studies are single clinical features or laboratory indicators, and it would be a better method to use multiple clinical features or laboratory indicators to develop a disease prediction model. The purpose of this study was to collect clinical characteristics and laboratory data from DM patients, to identify risk factors for RP-ILD in DM patients, and to develop a valid predictive model for RP-ILD in DM patients.

Materials and Methods Study Cohort and Patients

All 257 patients with DM included in this study were from January 2017 to December 2020 inpatients of the Affiliated Huai’an No.1 People’s Hospital of Nanjing Medical University, and 27 patients with DM in the validation cohort were from January 2021 to August 2021 inpatients of the Affiliated Huai’an No.1 People’s Hospital of Nanjing Medical University. The diagnosis of DM [including clinically amyopathic DM (CADM)] was based on the Bohan and Peter classification criteria. In this cohort, myositis-specific autoantibodies (MSAs) were measured by immunoblotting (EUROLINE, EUROIMMUN AG, Germany) by the same central laboratory in all patients. The medical records of all patients were retrospectively analyzed. Demographic data, including sex, age at onset and disease duration; clinical features, including muscle weakness, Gottron’s sign, heliotrope rash, V rush, shawl sign, periungueal erythema, arthritis, mechanic’s hand, skin ulcers (clinical features come from patient’s symptom record and are confirmed by doctor’s physical examination record); laboratory characteristics such as glutamate transaminase (ALT), aspartate transaminase (AST), lactate dehydrogenase (LDH), creatine kinase (CK), erythrocyte sedimentation rate (ESR), C-reactive protein (CRP), serum ferritin (SF), antinuclear antibody (ANA), anti-Ro-52 antibody, anti-aminoacyl-tRNA synthetase (ARS) antibodies, and anti-MDA5 antibody were extracted from the electronic records of patients at the beginning of hospitalization. The study complies with the Declaration of Helsinki, and the ethics committee of the Affiliated Huai’an No.1 People’s Hospital of Nanjing Medical University approved the study protocol. The requirement for informed consent was waived due to the study’s retrospective nature and anonymous data.

Inclusion and Exclusion Criteria

The diagnosis of ILD is based on the clinical signs and characteristics of CT scans. Diagnosis of ILD (notable ILD findings such as hairy glass attenuation, solid, reticular and honeycomb) is based on respiratory symptoms (dry cough and dyspnea on exertion), physical examination (such as velcro rales at the base of the lungs) and high-resolution CT (HRCT) findings, excluding interstitial changes due to infection and medications. RP-ILD is defined as the presence of one of the following four conditions within one month: 1) Progressive worsening of dyspnea or cough symptoms with a significant decrease in quality of life; 2) A decrease in lung function, including forceful lung volume (FVC) of more than 10% or a decrease in pulmonary diffusion capacity for carbon monoxide (DLCO) of more than 15%; 3) A persistent increase in the degree of interstitial pneumonia on high-resolution CT (HRCT) of the chest; 4) Arterial blood gas analysis demonstrating respiratory failure or a decrease in partial pressure of oxygen more significant than 10 mmHg, independently determined by a senior physician in the rheumatology and respiratory departments. Exclusion criteria include (1) patients with other autoimmune diseases, including rheumatoid arthritis, systemic lupus erythematosus and systemic sclerosis; (2) patients with lung cancer and chronic lung disease, including environmental, occupational and chronic obstructive pulmonary disease.

Vacancy Value Processing

Patients with more missing data (≥ 3 variables) are directly deleted; for patients with less missing data (< 3 variables), the KNN (K-Nearest Neighbor) algorithm is used to fill in the missing values.

Logistic Regression

The R package “rms” was used to perform univariate logistic regression (LR) for all independent variables, and the R package “caret” was used to perform multivariate LR for independent variables with significant differences, and the “step” function was used to achieve stepwise regression to filter the characteristic variables to develop the predictive model.

Machine Learning Approach

To develop predictive models of RP-ILD, three machine learning algorithms were applied: least absolute shrinkage and selection operator (LASSO), random forest (RF), and extreme gradient boosting (XGBoost). These machine algorithms (LASSO, RF, and XGBoost) were applied using the R package glmnet,19 ranger20 and xgboost,21 respectively. Metabolomics data was divided into two parts: the training set (70%) and the test set (30%). We performed a repeated K-fold cross-validation (repeats = 5, K = 5) on the models in order to obtain an unbiased estimate for their performance. Using the R package pROC, the receiver operating characteristic (ROC) curve was prepared.22

Results

A total of 284 patients with DM were included in this study, 257 of which were used for model development and internal validation studies and 27 for external validation. The main characteristics of the participants are shown in Table 1 and Supplementary Table 1.

Table 1 Clinical Characteristic of DM with RP-ILD and DM Without RP-ILD

Logistic Regression

Since the difference between the results in the univariate LR analysis does not reflect the effect of the factor on the outcome event very realistically, we used the variables with P < 0.2 in the results of the univariate LR analysis as the first echelon of candidate variables23 and screened out nine independent variables, namely age at onset, disease duration, muscle weakness, Gottron’s sign, LDH, ESR, CRP, anti-Ro-52 antibody, and anti-MDA5 antibody (Table 2). Among them, disease duration is a protective factor for RP-ILD, and ESR, CRP, anti-Ro-52 antibody and anti-MDA5 antibody are risk factors for RP-ILD. Previous studies have shown that when risk factors are screened using logistic regression, the regression results are more stable and more consistent with the original results when the events per variable (EPV) is greater than or equal to 10.24,25 There were 9 variables screened by univariate LR in this study; therefore, the number of EPVs required needed to exceed 90, and the number of EPVs required to select variables with P < 0.5 needed to exceed 50, and the number of DM patients with RP-ILD occurred in 41 cases; therefore, it was not suitable for further analysis using multivariate LR.

Table 2 Univariate Logistic Regression Analysis of Clinical and Laboratory Characteristics

Machine Learning Approach

To develop predictive models for RP-ILD, we adopted three machine learning methods. The AUC, sensitivity, and specificity of LASSO test set were 0.669, 0.417 and 0.922, respectively (Figure 1A); the AUC, sensitivity, and specificity of RF test set were 0.542, 0.083 and 1.000, respectively (Figure 1B); the AUC, sensitivity, and specificity of XGBoost test set were 0.663, 0.881 and 0.444, respectively (Figure 1C); and the performance summary of RF, LASSO and XGBoost on test set is shown in Table 3. On the basis of the first round of machine learning, we screened out the important variables in the models. There were 13 important variables in LASSO, as shown in Figure 1D. Figure 1E showed the importance of each variable in RF. There are 17 important variables in XGBoost, as shown in Figure 1F. The top 10 important variables of the 3 models were intersected to obtain common important variables. As shown in the Venn diagram (Figure 1G), there were 5 common important variables, namely disease duration, LDH, CRP, anti-Ro-52 antibody and anti-MDA5 antibody. The second round of machine learning is based on these five common variables. The AUC, sensitivity, and specificity of LASSO test set were 0.661, 0.417 and 0.906, respectively (Figure 2A, Table 3); the AUC, sensitivity, and specificity of RF test set were 0.667, 0.333 and 1.000, respectively (Figure 2B); the AUC, sensitivity, and specificity of XGBoost test set were 0.867, 0.900 and 0.833, respectively (Figure 2C); and the performance summary of RF, LASSO and XGBoost on test set is shown in Table 4.

Table 3 Performance Summary of Different Machine Learning Models Based on 23 Variables

Table 4 Performance Summary of Different Machine Learning Models Based on 5 Variables

Figure 1 Different machine learning methods on 23 variables. (A) ROC curve of LASSO on the test set. (B) ROC curve of RF on the test set. (C) ROC curve of XGBoost on the test set. (D) Important variables screening of LASSO. (E) Important variable screening of RF. (F) Important variable screening of XGBoost. (G) Venn diagram of machine learning methods.

Figure 2 ROC curve of prediction models based on 5 variables. (A) ROC curve of LASSO on the test set. (B) ROC curve of RF on the test set. (C) ROC curve of XGBoost on the test set. (D) ROC curve of LASSO on the validation set. (E) ROC curve of RF on the validation set. (F) ROC curve of XGBoost on the validation set.

Machine Learning Models Validation

We further validated the performance of these three models on validation set, and the results showed that, the AUC, sensitivity, and specificity of LASSO were 0.764, 0.800 and 0.727, respectively (Figure 2D); the AUC, sensitivity, and specificity of RF were 0.664, 0.600 and 0.727, respectively (Figure 2E); the AUC, sensitivity, and specificity of XGBoost were 0.909, 1.000 and 0.818, respectively (Figure 2F); and the performance summary of RF, LASSO and XGBoost on validation set is shown in Table 4. Based on the results of the models, XGBoost is the optimal model in this study. To make the prediction model more convenient and easy to use, we used the R package “shiny” to create a web version of the predictive model tool based on XGBoost at https://morrosun.shinyapps.io/XGBoost/.

Discussion

The mechanism by which RP-ILD occurs is unclear, and multiple cytokines are involved in the pathogenesis of RP-ILD, such as interferon (IFN)-γ, IFN-α, IFN-inducible protein-10 (IP-10), IL-6, IL-8, IL-10, IL-15, and tumor necrosis factor-α (TNF-α).26,27 There is no standard treatment protocol for DM-related RP-ILD, but aggressive combinations with high-dose glucocorticoids, calcium-modulated neurophosphatase inhibitors and intravenous cyclophosphamide have been proposed.17,28 Plasma exchange, tocilizumab and rituximab have been used to treat refractory RP-ILD.29–31 Treatments such as veno-venous extracorporeal membrane lung oxygenation or even lung transplantation should be considered in life-threatening cases.17,28

Early identification of RP-ILD is vital in the management of patients with DM, and a considerable number of studies have reported multiple risk factors associated with RP-ILD. Most of the risk factors reported in these studies are single clinical features or laboratory indicators, and it is better to use multiple clinical features or laboratory indicators to develop a disease prediction model. Disease prediction models are mainly various types of condition scoring systems. These tools are mainly developed based on LR or Cox regression models, and have been applied to the prediction of many diseases.32–34 Multivariable models can better assess disease severity, predict outcomes, and guide individualized treatment than univariate models. However, in this study, there are many independent variables, and the number of DM patients is relatively small, so it is impossible to establish a prediction model by LR. Therefore, we used three different machine learning methods to develop prediction models for RP-ILD, and their AUC values were 0.669, 0.542 and 0.663, respectively. All three models have low AUC values.

After further screening of the variables according to their importance to each model, as shown in the results of the venn diagram, we obtained five most important variables, namely disease duration, LDH, CRP, anti-Ro-52 antibody and anti-MDA5 antibody. Univariate LR results showed that anti-MDA5 antibody was the strongest risk factor for RP-ILD (OR = 9.17, P < 0.001), followed by anti-Ro-52 antibody (OR = 5.27, P < 0.001). Several studies have confirmed that anti-Ro-52 antibody positivity is an independent risk factor for myositis and even in ILD associated with Sjogren’s syndrome.35–38 RP-ILD occurs more frequently in anti-ARS antibody positive patients who coexist with anti-Ro-52 antibody than in patients without anti-Ro-52 antibody.39 In addition, the positive rate of MSAs was high in patients with isolated anti-Ro-52-ILD.40 The effect of anti-Ro-52 antibody positivity on RP-ILD in DM patients has not been specifically studied. Our study showed that anti-Ro-52 antibody positivity is also a vital risk factor for RP-ILD in patients with DM. The AUCs of the prediction models developed based on the five variables were 0.661, 0.667 and 0.867, respectively, with XGBoost performing the best. After validation using an independent validation set, it is still XGBoost that has the highest AUC value. Therefore, we chose fewer variables to construct the XGBoost-based prediction model. We created a web-based version of the prediction tool on this basis, which is better able to predict the risk of RP-ILD in DM patients at early disease stages.

There are still many deficiencies in this research. First of all, this is a single-center cohort study with a small sample; Second, this is a retrospective study that only included hospitalized patients, and there may be deviations in patient selection, and some patients’ medical records are incomplete and some data are missing. Thirdly, all patients in this study are Chinese, and the prediction effect on other races needs to be further verified. Fourthly, in addition to the laboratory data we collected, there are many laboratory indicators that are highly correlated with RP-ILD, such as cytokines, soluble CD 206, Krebs von den Lungen-6 and serum neopterin levels.12,16,18 However, this study is a retrospective study. The vast majority of patients lack these indicators. We will continue this work in the future, providing a better reference for the early diagnosis and treatment of RP-ILD.

Conclusions

In conclusion, anti-MDA5 antibody and anti-Ro-52 antibody are vital risk factors for RP-ILD in DM. The prediction model constructed using XGBoost can be used for risk identification and early intervention in DM patients with RP-ILD and practical application.

Disclosure

The authors report no conflicts of interest in this work.

References

1. O’Connell K, LaChance A. Dermatomyositis. N Engl J Med. 2021;384(25):2437. doi:10.1056/NEJMicm2033425

2. Waldman R, DeWane M, Lu J. Dermatomyositis: diagnosis and treatment. J Am Acad Dermatol. 2020;82(2):283–296. doi:10.1016/j.jaad.2019.05.105

3. Johnson C, Pinal-Fernandez I, Parikh R, et al. Assessment of Mortality in Autoimmune Myositis With and Without Associated Interstitial Lung Disease. Lung. 2016;194(5):733–737. doi:10.1007/s00408-016-9896-x

4. Nakashima R, Hosono Y, Mimori T. Clinical significance and new detection system of autoantibodies in myositis with interstitial lung disease. Lupus. 2016;25(8):925–933. doi:10.1177/0961203316651748

5. Koga T, Fujikawa K, Horai Y, et al. The diagnostic utility of anti-melanoma differentiation-associated gene 5 antibody testing for predicting the prognosis of Japanese patients with DM. Rheumatology. 2012;51(7):1278–1284. doi:10.1093/rheumatology/ker518

6. Jablonski R, Bhorade S, Strek M, et al. Recognition and Management of Myositis-Associated Rapidly Progressive Interstitial Lung Disease. Chest. 2020;158(1):252–263. doi:10.1016/j.chest.2020.01.033

7. Li J, Liu Y, Li Y, et al. Associations between anti-melanoma differentiation-associated gene 5 antibody and demographics, clinical characteristics and laboratory results of patients with dermatomyositis: a systematic meta-analysis. J Dermatol. 2018;45(1):46–52. doi:10.1111/1346-8138.14092

8. Motegi S, Sekiguchi A, Toki S, et al. Clinical features and poor prognostic factors of anti-melanoma differentiation-associated gene 5 antibody-positive dermatomyositis with rapid progressive interstitial lung disease</span>. Eur J Dermatol. 2019;29(5):511–517. doi:10.1684/ejd.2019.3634

9. Wu W, Guo L, Fu Y, et al. Interstitial Lung Disease in Anti-MDA5 Positive Dermatomyositis. Clin Rev Allerg Immu. 2021;60:5.

10. Karino K, Kono M, Kono M, et al. Myofascia-dominant involvement on whole-body MRI as a risk factor for rapidly progressive interstitial lung disease in dermatomyositis. Rheumatology. 2020;59(7):1734–1742. doi:10.1093/rheumatology/kez642

11. Xu Y, Yang C, Li Y, et al. Predictive factors of rapidly progressive-interstitial lung disease in patients with clinically amyopathic dermatomyositis. Clin Rheumatol. 2016;35(1):113–116. doi:10.1007/s10067-015-3139-z

12. Shen Y, Zhang Y, Huang Z, et al. Increased Levels of Soluble CD206 Associated with Rapidly Progressive Interstitial Lung Disease in Patients with Dermatomyositis. Mediators Inflamm. 2020;2020:7948095. doi:10.1155/2020/7948095

13. Sagawa T, Kida T, Inaba T, et al. Utility of Coagulation Markers for the Prediction of Rapidly Progressive Interstitial Lung Disease in Patients with Dermatomyositis. Lung. 2019;197(4):437–442. doi:10.1007/s00408-019-00245-0

14. Shimizu T, Koga T, Furukawa K, et al. IL-15 is a biomarker involved in the development of rapidly progressive interstitial lung disease complicated with polymyositis/dermatomyositis. J Intern Med. 2021;289(2):206–220. doi:10.1111/joim.13154

15. Zhu D, Qiao J, Tang S, et al. Elevated carcinoembryonic antigen predicts rapidly progressive interstitial lung disease in clinically amyopathic dermatomyositis. Rheumatology. 2021;60(8):3896–3903. doi:10.1093/rheumatology/keaa819

16. Peng Q, Zhang Y, Liang L, et al. A high level of serum neopterin is associated with rapidly progressive interstitial lung disease and reduced survival in dermatomyositis. Clin Exp Immunol. 2020;199(3):314–325. doi:10.1111/cei.13404

17. Selva-O’Callaghan A, Romero-Bueno F, Trallero-Araguás E, et al. Pharmacologic Treatment of Anti-MDA5 Rapidly Progressive Interstitial Lung Disease. Curr Treatm Opt Rheumatol. 2021;1:1–15.

18. Kobayashi N, Takezaki S, Kobayashi I, et al. Clinical and laboratory features of fatal rapidly progressive interstitial lung disease associated with juvenile dermatomyositis. Rheumatology. 2015;54(5):784–791. doi:10.1093/rheumatology/keu385

19. Friedman J, Hastie T, Tibshirani R. Regularization paths for generalized linear models via coordinate descent. J Stat Softw. 2010;33(1):1. doi:10.18637/jss.v033.i01

20. Wright MN, Ziegler A. ranger: a fast implementation of random forests for high dimensional data in C++ and R. J Stat Softw. 2017;77(1):1–17. doi:10.18637/jss.v077.i01

21. Chen T, Guestrin C XGBoost: a Scalable Tree Boosting System. 22nd SIGKDD Conference on Knowledge Discovery and Data Mining; 2016.

22. Robin X, Turck N, Hainard A, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinform. 2011;12(1):1–8. doi:10.1186/1471-2105-12-77

23. Kang S, Cho Y, Park G, et al. Predictors for functionally significant in-stent restenosis: an integrated analysis using coronary angiography, IVUS, and myocardial perfusion imaging. JACC Cardiovasc Imaging. 2013;6(11):1183–1190. doi:10.1016/j.jcmg.2013.09.006

24. Harrell F, Lee K, Mark D. Multivariable prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors. Stat Med. 1996;15(4):361–387. doi:10.1002/(SICI)1097-0258(19960229)15:4<361::AID-SIM168>3.0.CO;2-4

25. Peduzzi P, Concato J, Feinstein A, et al. Importance of events per independent variable in proportional hazards regression analysis. II. Accuracy and precision of regression estimates. J Clin Epidemiol. 1995;48(12):1503–1510. doi:10.1016/0895-4356(95)00048-8

26. Ishikawa Y, Iwata S, Hanami K, et al. Relevance of interferon-gamma in pathogenesis of life-threatening rapidly progressive interstitial lung disease in patients with dermatomyositis. Arthritis Res Ther. 2018;20(1):240. doi:10.1186/s13075-018-1737-2

27. Gono T, Kaneko H, Kawaguchi Y, et al. Cytokine profiles in polymyositis and dermatomyositis complicated by rapidly progressive or chronic interstitial lung disease. Rheumatology. 2014;53(12):2196–2203. doi:10.1093/rheumatology/keu258

28. Romero-Bueno F, Diaz Del Campo P, Trallero-Araguás E, et al. Recommendations for the treatment of anti-melanoma differentiation-associated gene 5-positive dermatomyositis-associated rapidly progressive interstitial lung disease. Semin Arthritis Rheum. 2020;50(4):776–790. doi:10.1016/j.semarthrit.2020.03.007

29. Takahashi R, Yoshida T, Morimoto K, et al. Successful Treatment of Anti-MDA5 Antibody-Positive Dermatomyositis-Associated Rapidly Progressive-Interstitial Lung Disease by Plasma Exchange: two Case Reports. Clin Med Insights Case Rep. 2021;14:11795476211036322. doi:10.1177/11795476211036322

30. Zhang X, Zhou S, Wu C, et al. Tocilizumab for refractory rapidly progressive interstitial lung disease related to anti-MDA5-positive dermatomyositis. Rheumatology. 2021;60(7):e227–e228. doi:10.1093/rheumatology/keaa906

31. So H, Wong V, Lao V, et al. Rituximab for refractory rapidly progressive interstitial lung disease related to anti-MDA5 antibody-positive amyopathic dermatomyositis. Clin Rheumatol. 2018;37(7):1983–1989. doi:10.1007/s10067-018-4122-2

32. Eddama M, Fragkos K, Renshaw S, et al. Logistic regression model to predict acute uncomplicated and complicated appendicitis. Ann R Coll Surg Engl. 2019;101(2):107–118. doi:10.1308/rcsann.2018.0152

33. Nusinovici S, Tham Y, Chak Yan M, et al. Logistic regression was as good as machine learning for predicting major chronic diseases. J Clin Epidemiol. 2020;122:56–69. doi:10.1016/j.jclinepi.2020.03.002

34. Cadrin-Tourigny J, Bosman L, Nozza A, et al. A new prediction model for ventricular arrhythmias in arrhythmogenic right ventricular cardiomyopathy. Eur Heart J. 2019;40(23):1850–1858. doi:10.1093/eurheartj/ehz103

35. Xing X, Li A, Li C. Anti-Ro52 antibody is an independent risk factor for interstitial lung disease in dermatomyositis. Respir Med. 2020;172:106134. doi:10.1016/j.rmed.2020.106134

36. Xu A, Ye Y, Fu Q, et al. Prognostic values of anti-Ro52 antibodies in anti-MDA5-positive clinically amyopathic dermatomyositis associated with interstitial lung disease. Rheumatology. 2021;60(7):3343–3351. doi:10.1093/rheumatology/keaa786

37. Sabbagh S, Pinal-Fernandez I, Kishi T, et al. Anti-Ro52 autoantibodies are associated with interstitial lung disease and more severe disease in patients with juvenile myositis. Ann Rheum Dis. 2019;78(7):988–995. doi:10.1136/annrheumdis-2018-215004

38. Buvry C, Cassagnes L, Tekath M, et al. Anti-Ro52 antibodies are a risk factor for interstitial lung disease in primary Sjögren syndrome. Respir Med. 2020;163:105895. doi:10.1016/j.rmed.2020.105895

39. Shi J, Li S, Yang H, et al. Clinical Profiles and Prognosis of Patients with Distinct Antisynthetase Autoantibodies. J Rheumatol. 2017;44(7):1051–1057. doi:10.3899/jrheum.161480

40. Shao C, Sun Y, Huang H, et al. Myositis specific antibodies are associated with isolated anti-Ro-52 associated interstitial lung disease. Rheumatology. 2022;61(3):1083–1091. doi:10.1093/rheumatology/keab488

留言 (0)

沒有登入
gif