Literature DB >> 31446631

Development and accuracy of an artificial intelligence algorithm for acne grading from smartphone photographs.

Sophie Seité1, Amir Khammari2, Michael Benzaquen3, Dominique Moyal1, Brigitte Dréno2.   

Abstract

We developed an artificial intelligence algorithm (AIA) for smartphones to determine the severity of facial acne using the GEA scale and to identify different types of acne lesion (comedonal, inflammatory) and postinflammatory hyperpigmentation (PIHP) or residual hyperpigmentation. Overall, 5972 images (face, right and left profiles) obtained with smartphones (IOS and/or Android) from 1072 acne patients were collected. Three trained dermatologists assessed the acne severity for each patient. One acne severity grade per patient (grade given by the majority of the three dermatologists from the two sets of three images) was used to train the algorithm. Acne lesion identification was performed from a subgroup of 348 images using a tagging tool; tagged images served to train the algorithm. The algorithm evolved and was adjusted for sensibility, specificity and correlation using new images. The correlation between the GEA grade and the quantification and qualification of acne lesions both by the AIA and the experts for each image were evaluated for all AIA versions. At final version 6, the GEA grading provided by AIA reached 68% and was similar to that provided by the dermatologists. Between version 4 and version 6, AIA improved precision results multiplied by 1.5 for inflammatory lesions, 2.5 for non-inflammatory lesions and by 2 for PIHP; recall was improved by 2.6, 1.6 and 2.7. The weighted average of precision and recall or F1 score was 84% for inflammatory lesions, 61% for non-inflammatory lesions and 72% for PIHP.
© 2019 The Authors. Experimental Dermatology published by John Wiley & Sons Ltd.

Entities:  

Keywords:  GEA scale; Global Acne Severity scale; Smartphone; acne; algorithm; artificial intelligence; lesion count

Mesh:

Year:  2019        PMID: 31446631      PMCID: PMC6972662          DOI: 10.1111/exd.14022

Source DB:  PubMed          Journal:  Exp Dermatol        ISSN: 0906-6705            Impact factor:   3.960


INTRODUCTION

Acne is a chronic multi‐factorial inflammatory disease of the pilosebaceous follicle that may affect the patient's quality of life. Acne evolves through flare‐ups and affects more than 85% of adolescents. It often continues into adulthood.1 Over the past years, adult acne and acne relapse (such as acne reappearance at regular intervals) have increased in frequency and recent data indicate that, in France, the total number of working days lost due to acne relapses add up to a total of 350 000 days per year.2 Moreover, with a decreasing number of dermatologists over the last decades, getting appointments with dermatologists has become more and more difficult for patients, resulting in irregular patient follow‐up, poor treatment outcome and a positive correlation between the duration of the disease prior to treatment and the development of scars.3, 4, 5 This highlights the need for an appropriate and timely management of acne patients to avoid the development of residual hyperpigmentation and scars or a negative impact on the patient's quality of life. Therefore, the correct grading and identification of acne lesions is the first step in a successful treatment of acne, allowing the most suitable therapeutic approach to be chosen. To date, health agencies recommend both the inflammatory and non‐inflammatory acne lesion count using acne grading scales such as the Investigator Global Assessment (IGA) for the USA, with up to five ordinal grades (0‐4, ie clear, almost clear, mild, moderate, severe and very severe) or the Global Acne Severity Scale (GEA) for Europe, with a scale ranging from 0 to 5 (no lesions, virtually no lesions, mild, moderate, severe and very severe).6, 7 Computational methods, including the automatic counting of lesions and/or evaluation of the severity grading, similar to that used by doctors in daily clinical practice, are increasing in interest. However, to date, these methods have still not been validated.8, 9, 10 A recent study conducted in the USA assessed the severity of acne using Artificial Intelligence (AI) and the IGA scale.8 AI was able to classify acne lesions according to an IGA ordinal scale with high accuracy, no human intervention and with no need for lesion counting. Moreover, several smartphone applications have been proposed as diagnostic self‐monitoring tools.11 In this article, we present a validated AI algorithm (AIA) which, using smartphones, allows the grading of acne severity based on the GEA scale. It provides a fast determination of the severity of acne and identifies the different lesion types, thus supporting the early therapeutic management of acne patients.

METHODS

The AIA was developed in six steps as described in Figure 1.
Figure 1

Development and validation process of artificial intelligence algorithm. Internal data testing was performed using Krippendorff's alpha and Cohen's kappa. Clinical tests were performed using an interclass correlation coefficient and the Cicchetti interpretation table. The F1 score (0‐1) equalled the weighted average of precision and recall. PIHP, postinflammatory hyperpigmentation

Development and validation process of artificial intelligence algorithm. Internal data testing was performed using Krippendorff's alpha and Cohen's kappa. Clinical tests were performed using an interclass correlation coefficient and the Cicchetti interpretation table. The F1 score (0‐1) equalled the weighted average of precision and recall. PIHP, postinflammatory hyperpigmentation

Step 1—Data collection

A total of 5972 images from the face, and from right and left profiles of 1072 acne patients from different racial groups having signed an informed consent, were collected from France, South Africa, China and India. All patients had to have different acne severity and images were taken using, if possible, at the same time, two smartphones equipped with IOS or Android systems, as described in Table 1.
Table 1

Overall patient demographics and acne severity (1072 acne patients—5972 images)

 Male (%)Female (%)Age ± SD (y)IOS (n)AN (n)GEA 0 (n)GEA 1 (n)GEA 2 (n)GEA 3 (n)GEA 4 (n)GEA 5 (n)
Caucasian326724.1 ± 9.01430153742248160106132
African366422.2 ± 8.28828827125129321 
Asian307028.0 ± 11.742942966747221 
Latin386221.7 ± 7.13939  274 
Indien316921.4 ± 3.7153153   15333
Total356523.9 ± 9.22933303955440338182525

Abbreviations: AN, Android; GEA, Group of experts in acne; IOS, Apple system; SD, standard deviation.

Overall patient demographics and acne severity (1072 acne patients—5972 images) Abbreviations: AN, Android; GEA, Group of experts in acne; IOS, Apple system; SD, standard deviation.

Step 2—Grading of data and statistical analysis

Three trained dermatologists with expertise in acne graded each patient's acne severity using the European GEA scale and the three image views as described inFigure 1.6 For each patient, the final GEA grade chosen was that confirmed by at least 2 out of the 3 dermatologists (Table 1). In addition, inflammatory and non‐inflammatory lesions and PIHP were tagged by one of the three dermatologists on 348 images from 117 acne patients of the three main racial groups, and from GEA 0 to 4 for their identification and recognition (Figure 1 and Table 2).
Table 2

Patient demographics and acne severity used for tagging by a dermatologist (117 acne patients—348 images)

 Male (%)Female (%)Age ± SD (y)IOS (n)AN (n)GEA 0 (n)GEA 1 (n)GEA 2 (n)GEA 3 (n)GEA 4 (n)GEA 5 (n)
Caucasian (n = 46)485219.7 ± 7.2904509161830
African (n = 36)227819.9 ± 6.048600721800
Asian (n = 35)267424.1 ± 12.55451310111010
Total (n = 117)336721.6 ± 9.3192156326483640

Abbreviations: GEA, Group of experts in acne; SD, standard deviation.

Patient demographics and acne severity used for tagging by a dermatologist (117 acne patients—348 images) Abbreviations: GEA, Group of experts in acne; SD, standard deviation.

Step 3—Training Of The AIA

A total of 4958 images (from the 5972 pool) corresponding to 903 acne participants were used to develop the algorithm (Table 3). Each set of three images associated with one GEA grade and the images with tagged lesions and PIHP were used to train the algorithm. As the difficulty to obtain images from patients with severe acne resulted in a low number of Grade 4 and 5, participants with an acne grade of ‘4’ or ‘5’ were pooled and graded ‘4+’.
Table 3

Patient demographics and acne severity to develop the algorithm (903 acne patients—4958 images)

 Male (%)Female (%)Age ± SD (y)IOS (n)AN (n)GEA 0 (n)GEA 1 (n)GEA 2 (n)GEA 3 (n)GEA 4 (n)
Caucasian336724.5 ± 9.412231329422231388315
African336722.9 ± 8.67237237115100181
Asian307028.0 ± 11.742942966747221
Latin336728.3 ± 7.899    3
Indian217923.1 ± 5.84242    15
Total336724.6 ± 9.7242625325540528512335

Abbreviations: AN, Android; GEA, Group of experts in acne; IOS, Apple system; SD, standard deviation.

Patient demographics and acne severity to develop the algorithm (903 acne patients—4958 images) Abbreviations: AN, Android; GEA, Group of experts in acne; IOS, Apple system; SD, standard deviation. From these images, image classification and segmentation techniques based on deep learning provided by Perfect Mobile Corp (New Taipei City, Taiwan) were applied to develop the algorithm. Image classification techniques were used for classifying the input image into five different groups (GEA 0, 1, 2, 3, 4+) in order to define the final GEA score. In addition, image detection/segmentation techniques were used to localize the different lesion types and PIHP on the target areas and results were submitted to a dermatologist for correction and improvement.

Step 4—Testing of the AIA versions

The algorithm was tested on a set of 1014 new images from 169 different acne patients than in step 2.3 (Table 4).
Table 4

Patient demographics and acne severity for internal testing of the algorithm (169 acne patients—1014 images)

 Male (%)Female (%)Age ± SD (y)IOS (n)AN (n)GEA 0 (n)GEA 1 (n)GEA 2 (n)GEA 3 (n)GEA 4 (n)
Caucasian297122.0 ± 5.9210210 262123 
African475318.7 ± 5.3159159 112913 
Asian          
Latin406019.7 ± 5.83030  271
Indian336721.0 ± 2.4108108   1521
Total435720.6 ± 5.3507507037525822

Abbreviations: AN, Android; GEA, Group of experts in acne; IOS, Apple system; SD, standard deviation.

Patient demographics and acne severity for internal testing of the algorithm (169 acne patients—1014 images) Abbreviations: AN, Android; GEA, Group of experts in acne; IOS, Apple system; SD, standard deviation.

Step 5—Clinical testing of AIA in acne patients

Given good internal test results, the 4th version of the algorithm was tested through a clinical evaluation performed in 53 acne patients from the department of Dermatology Nantes France with a severity grade ranging from ‘1’ to ‘3’ on the GEA scale 6 (Table 5). Face‐to‐face evaluations and from photographs of GEA grades and lesion counts by the dermatologists and by the AIA were compared.
Table 5

Patient demographics and acne severity in the clinical test (53 acne patients—159 images)

Male (%)Female (%)Age ± SD (y)MinMaxIOS Photograph 1 (n)IOS Photograph 2 (n)GEA0 (n)GEA1 (n)GEA2 (n)GEA3 (n)GEA4 (n)GEA5 (n)
28% (n = 15)72% (n = 38)21.8 ± 7.3134715915903113900

Abbreviations: GEA, Group of experts in acne; IOS, Apple system; SD, standard deviation.

Patient demographics and acne severity in the clinical test (53 acne patients—159 images) Abbreviations: GEA, Group of experts in acne; IOS, Apple system; SD, standard deviation.

Step 6—Improvement and validation of the algorithm

Any mistake in the GEA evaluation and location of the different lesion types and PIHP was indicated to the developer in order to improve the AIA, and the AIA was to be tested again.

Statistics

Inter‐dermatologist reproducibility was checked using Krippendorff's alpha and Cohen's kappa statistical tests.12 Intra‐class correlation coefficient or ICC and a descriptive statistic test were used to measure the reproducibility of numerical measurements made by the different examiners. Similar tests were used to compare lesion counts by dermatologists and the algorithm. Furthermore, for their location, precision, defined as the correct number and type of lesions identified by the AIA, and recall, defined as the number and type of lesions correctly identified by the dermatologists and by the AIA, were evaluated. Tags made by dermatologists were compared with algorithm tags by identifying true (TP)‐ and false (FP)‐positive tags and true (TN)‐ and false (FN)‐negative tags and precision (TP/ (TP + FP)) and recall (TP/(TP + FN)). The weighted average of precision and recall or F1 score corresponding to the ratio between 2 (precision x recall)/(precision + recall) was calculated for all the images. Coefficients were interpreted using the Cicchetti interpretation table (ICC ≥ 0.75 = Excellent; [0.60‐0.75] = Good; [0.40‐0.60] = Fair and <0.40 = poor).13

RESULTS

Full face and right and left profile images of 1072 acne patients from different racial groups (10 to 58 years old) were collected and graded for their severity by three dermatologists (Table 1). The weighted kappa for the evaluation performed by the three dermatologists using the GEA grade for the three main racial groups of different phototype exceeded 0.6, corresponding to a substantial correlation. When comparing the GEA grade obtained for the same patient, a substantial agreement for African Black (for IOS [0.69 (95%: −0.64 to 0.75; P < .0001)] and for Android [0.65 (95%: −0.6 to 0.71; P < .0001)]), a moderate agreement for Caucasian White (for IOS [0.59 (95%: −0.54 to 0.64; P < .0001)] and Android [0.52 (95%: −0.47 to 0.57; P < .0001)]) and a moderate/fair agreement for Asian participants (moderate for IOS [0.52 (95%: −0.43 to 0.6; P < .0001)] and fair for Android [0.38 (95%: −0.3 to 0.45; P < .0001)]) was obtained. In addition, inflammatory and non‐inflammatory lesions and PIHP were tagged on 348 photographs from 117 patients by a dermatologist (Table 2). Results indicate that the highest accuracy among the three dermatologists was obtained for Black Africans using ‘IOS images’ (‘Substantial’—Alpha = 0.69, P < .0001) and the lowest agreement was obtained for Chinese, using ‘Android images’ (‘Fair’—Alpha = 0.38, P < .0001). Thus, agreement depends on the type of device using IOS or Android, even if they were equipped with cameras of similar quality, as the difference in colour rendition influenced on the acne severity assessment (more red for IOS than Android). Moreover, results obtained for Asians have to be considered with caution, due to the low number of participants. A total of 4958 images corresponding to 903 acne participants with the GEA grade given by the majority of the 3 dermatologists from the two sets of three images (Table 3) associated with a total of 2939 tagged inflammatory lesions, 7603 tagged non‐inflammatory lesions and 5702 tagged PIHP were used to develop the algorithm. During the development, tags performed by the algorithm were submitted to the same dermatologist for correction; 1014 new images from 169 acne patients (Table 4) were used for testing using the GEA scale and counting of inflammatory lesions, non‐inflammatory lesions and PIHP. The clinical evaluation conducted in the 53 acne patients (Table 5) by three dermatologists showed a correspondence between the 4th version of the AIA and the different assessments made by the dermatologists using the GEA scale reaching 48%; conversely, it was 59% between dermatologists using images and 67% during face‐to‐face assessments. The intra‐class correlation coefficient (ICC) for the face‐to‐face counting by the three dermatologists varied from 0.68 (Good) to 0.8 (excellent) for inflammatory lesions, from 0.47 (fair) to 0.89 (excellent) for non‐inflammatory lesions and from 0.16 (poor) to 0.55 (fair) for PIHP. In comparison, for the algorithm, the ICC obtained for images counting were lower and varied from 0.51 (fair) to 0.74 (good) for inflammatory lesions, from 0.2 (poor) to 0.51 (fair) for non‐inflammatory lesions and from 0.1 (poor) to 0.57 (fair) for PIHP. During this validation step, we also evaluated the reproducibility of the algorithm through two sets of three images of the same participants and correspondence rate of 74% was noted. After improvement of the AIA, when testing final version 6, the GEA grading provided by the algorithm reached 68% and was similar to that provided by the dermatologists. Furthermore, the algorithm provided improved precision results multiplied by 1.5 for inflammatory lesions, 2.5 for non‐inflammatory lesions and 2 for PIHP; recall was improved by 2.6, 1.6 and 2.7, respectively, as shown in Table 6. The F1 score calculated at the end of the process on the global image basis for the sixth version considered accurate was 84% for inflammatory lesions, 61% for non‐inflammatory lesions and 72% for PIHP.
Table 6

Precision and recall provided by the different version of the algorithm

 PrecisionRecall
AIA V4 vs dermatologist (%)AIA V5 vs dermatologist (%)AIA V6 vs dermatologist (%)AIA V4 vs dermatologist (%)AIA V5 vs dermatologist (%)AIA V6 vs dermatologist (%)
Comedonal lesions192448143237
Inflammatory lesions475972436473
PIHP253249224160

Abbreviations: AIA, Artificial intelligence algorithm; PIHP, postinflammatory hyperpigmentation.

Precision and recall provided by the different version of the algorithm Abbreviations: AIA, Artificial intelligence algorithm; PIHP, postinflammatory hyperpigmentation.

DISCUSSION

Even though acne is a very common disease, only little time or no time at all is spent on medical and dermatology education.14 For this reason, it is important for clinicians to regularly update their practice to reflect current diagnostic and treatment standards. Moreover, clinical assessment of acne severity does not always correlate with the patient's perception of its severity.15, 16, 17, 18, 19, 20 Because the correct grading and identification of acne lesions is the first step in a successful treatment of acne, this AIA may help dermatologists and other medical practitioners to choose the most suitable treatment options according to their patientsacne.15 We have developed an AIA able to assess acne severity and to identify and to count acne lesions from three images (face and right and left profiles) of acne patients using smartphones with IOS and Android systems, which both are used worldwide, thus also allowing global use of our AIA. The present AIA was able to classify acne patients according to the severity of their disease based on the GEA scale, a scale that is widely used in Europe to diagnose acne severity and which has been validated through direct clinical and image grading.6 Moreover, we have studied the suitability of the GEA scale for images of Black African and Asian (Chinese) in which the problem of PIHP is more important than in the Caucasian population. Research studies use a multitude of scales to grade acne, and there is still no standard evaluation tool.21 The present AIA may be an interesting and easy‐to‐use tool for standardized clinical investigations in acne. Moreover, this AIA, when using geo‐localization, may also allow to assess the impact of external factors that influence the occurrence, duration and severity of acne such as pollution and climatic condition.22 Non‐adherence to acne treatment is a common problem worldwide.3, 4 Acne requires continued and regular treatment, with topical and oral therapies, but the use of both is currently associated with a high level of non‐adherence.23 Adherence may be improved by early management, including an increased frequency of visits and an improved more frequent evaluation of improvement and by simplifying treatment regimens. Effective cosmetic management alongside pharmaceutical treatment was shown to increase adherence.15, 24 Other work has already been done to develop an automatic classification system of acne patients, focusing on the automatization of acne lesion counts or for the severity measurement using the IGA scale, widely used in the USA.8, 9, 10 However, to our knowledge, this is the first time that an AIA for smartphones has been developed that assesses acne severity based on the GEA scale and that allows the identification and counting of acne lesions with reliability comparable to that of a trained dermatologist. Moreover, the present algorithm may be suitable for the three main racial groups, suggesting a high generalizability. The main limitation of our work is the fact that the algorithm has been built on a defined aim and variables. GEA evaluation and lesion tagging were performed by trained dermatologists on low‐quality smartphone images, using visible light exposure without following any precise protocol. Therefore, automatization of these evaluations with any smartphone working for three racial groups may reduce the variability. This is critical because of the higher reproducibility of results generated by the algorithm compared with those provided by the dermatologists. Despite these limitations, we believe that the present artificial intelligence algorithm is a promising tool for an easy and reproducible evaluation of acne severity and lesion identification, allowing in the near future to adapt acne treatment to each patient as best as possible.

CONFLICT OF INTERESTS

S. Seité and Dominique Moyal are employees of La Roche‐Posay, France. Brigitte Dréno received consultancy honoraria from La Roche‐Posay, Galderma, Pierre Fabre. Amir Khammari received consultancy honoraria from La Roche‐Posay. Michael Benzaquen received consultancy honoraria from La Roche‐Posay.

AUTHOR CONTRIBUTIONS

SS, DM and BD have participated in the conception and design of the studies, AK and MB in the acquisition and analysis of data, and all authors in the preparation and critical review of the manuscript. All read and approved the final submitted version.
  23 in total

1.  Scarred for life?

Authors:  A M Layton; D Seukeran; W J Cunliffe
Journal:  Dermatology       Date:  1997       Impact factor: 5.366

2.  [Recent data on epidemiology of acne].

Authors:  B Dréno
Journal:  Ann Dermatol Venereol       Date:  2010-11       Impact factor: 0.777

3.  Simulation study of a panel of reliability indicators applied to paired measurements.

Authors:  Carl M Russell; David F Williamson; John J Bartko; Edwin L Bradley
Journal:  Am J Hum Biol       Date:  1994       Impact factor: 1.937

4.  Artificial Intelligence for the Objective Evaluation of Acne Investigator Global Assessment.

Authors:  Antonella Melina; Nhan Ngo Dinh; Benedetta Tafuri; Giusy Schipani; Steven Nisticò; Carlo Cosentino; Francesco Amato; Diane Thiboutot; Andrea Cherubini
Journal:  J Drugs Dermatol       Date:  2018-09-01       Impact factor: 2.114

Review 5.  Practical management of acne for clinicians: An international consensus from the Global Alliance to Improve Outcomes in Acne.

Authors:  Diane M Thiboutot; Brigitte Dréno; Abdullah Abanmi; Andrew F Alexis; Elena Araviiskaia; Maria Isabel Barona Cabal; Vincenzo Bettoli; Flordeliz Casintahan; Steven Chow; Adilson da Costa; Tam El Ouazzani; Chee-Leok Goh; Harald P M Gollnick; Minerva Gomez; Nobukazu Hayashi; Maria Isabel Herane; Juan Honeyman; Sewon Kang; Lajos Kemeny; Raj Kubba; Julien Lambert; Alison M Layton; James J Leyden; Jose Luis López-Estebaranz; Nopadon Noppakun; Falk Ochsendorf; Cristina Oprica; Beatriz Orozco; Montserrat Perez; Jaime Piquero-Martin; Jo-Ann See; Dae Hun Suh; Jerry Tan; Vicente Torres Lozada; Patricia Troielli; Leihong Flora Xiang
Journal:  J Am Acad Dermatol       Date:  2017-11-08       Impact factor: 11.527

Review 6.  Going digital: a narrative overview of the effects, quality and utility of mobile apps in chronic disease self-management.

Authors:  Ian A Scott; Paul Scuffham; Deepali Gupta; Tanya M Harch; John Borchi; Brent Richards
Journal:  Aust Health Rev       Date:  2020-02       Impact factor: 1.990

7.  Creation of a tool to assess adherence to treatments for acne.

Authors:  Henry Pawin; Claire Beylot; Martine Chivot; Michel Faure; Florence Poli; Jean Revuz; Brigitte Dréno
Journal:  Dermatology       Date:  2008-10-22       Impact factor: 5.366

8.  Global Burden of Skin Disease: Inequities and Innovations.

Authors:  Divya Seth; Khatiya Cheldize; Danielle Brown; Esther F Freeman
Journal:  Curr Dermatol Rep       Date:  2017-08-07

9.  Interrater reliability: the kappa statistic.

Authors:  Mary L McHugh
Journal:  Biochem Med (Zagreb)       Date:  2012       Impact factor: 2.313

Review 10.  The influence of exposome on acne.

Authors:  B Dréno; V Bettoli; E Araviiskaia; M Sanchez Viera; A Bouloc
Journal:  J Eur Acad Dermatol Venereol       Date:  2018-02-15       Impact factor: 6.166

View more
  4 in total

Review 1.  Lack of Transparency and Potential Bias in Artificial Intelligence Data Sets and Algorithms: A Scoping Review.

Authors:  Roxana Daneshjou; Mary P Smith; Mary D Sun; Veronica Rotemberg; James Zou
Journal:  JAMA Dermatol       Date:  2021-11-01       Impact factor: 11.816

2.  Automatic Acne Object Detection and Acne Severity Grading Using Smartphone Images and Artificial Intelligence.

Authors:  Quan Thanh Huynh; Phuc Hoang Nguyen; Hieu Xuan Le; Lua Thi Ngo; Nhu-Thuy Trinh; Mai Thi-Thanh Tran; Hoan Tam Nguyen; Nga Thi Vu; Anh Tam Nguyen; Kazuma Suda; Kazuhiro Tsuji; Tsuyoshi Ishii; Trung Xuan Ngo; Hoan Thanh Ngo
Journal:  Diagnostics (Basel)       Date:  2022-08-03

Review 3.  COVID-19 and artificial intelligence: Experts and dermatologists perspective.

Authors:  Yaser Goldust; Farah Sameem; Samia Mearaj; Atula Gupta; Anant Patil; Mohamad Goldust
Journal:  J Cosmet Dermatol       Date:  2022-08-17       Impact factor: 2.189

Review 4.  Artificial Intelligence Applications in Dermatology: Where Do We Stand?

Authors:  Arieh Gomolin; Elena Netchiporouk; Robert Gniadecki; Ivan V Litvinov
Journal:  Front Med (Lausanne)       Date:  2020-03-31
  4 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.