Literature DB >> 21912501

A comparison of three-dimensional and two-dimensional cephalometric evaluations of children with cleft lip and palate.

Ozlem Tulunoglu1, Elcin Esenlik, Ayse Gulsen, Ibrahim Tulunoglu.   

Abstract

OBJECTIVES: The aim of this retrospective study was to compare the consistency of orthodontic measurement performed on cephalometric films and 3D CT images of cleft lip and palate (CLP) patients.
METHODS: The study was conducted with 2D radiographs and 3D CT images of 9 boys and 6 girls aged 7-12 with CLP. 3D reconstructions were performed using MIMICS software.
RESULTS: Frontal analysis found statistical differences for all parameters except occlusal plane tilt (OcP-tilt) and McNamara analysis found statistical differences in 2D and 3D measurements for all parameters except ANS-Me and Co-Gn; Steiner analysis found statistical differences for all parameters except SND, SNB and Max1-SN. Intra-group variability in measurements was also very low for all parameters for both 2D and 3D images.
CONCLUSIONS: Study results indicate significant differences between measurements taken from 2D and 3D images in patients with cleft lip and palate.

Entities:  

Keywords:  3D reconstruction; CLP; Cephalometry; Computerized tomography

Year:  2011        PMID: 21912501      PMCID: PMC3170029     

Source DB:  PubMed          Journal:  Eur J Dent


INTRODUCTION

Conventional cephalometry is an inexpensive and essential method for evaluating orthodontic patients.1,2 However, it is difficult to evaluate patients with major craniofacial deformities, in particular facial asymmetry, using conventional cephalometry.1 For this reason, three-dimensional computed tomography (3D-CT) imaging is widely used in the diagnosis and surgical treatment planning of patients with craniofacial deformities.3–5 3D-CT data is especially helpful in evaluating patients with asymmetry because they do not suffer from the magnification- and distortion-related problems inherent in 2D projections. With the general trend away from traditional 2D cephalometric analysis towards new 3D techniques, it is helpful to compare the use of different applications in various settings.6,7 A study conducted by Kragskov et al1 comparing lateral and frontal cephalograms and CT scans of human dry skulls found lateral cephalograms to be more reliable than 3D-CT, with less than 1 mm of variation between observers for most points, compared to about 2 mm for 3D-CT images. Differences between 3D-CT and frontal cephalograms were less obvious. They also stated that, frontal cephalometrics showed significantly less inter-observer variation, and for standard lateral and frontal cephalometric points, there is no evidence that 3D-CT is more reliable than the conventional cephalometric methods in normal skull.1 However, conventional cephalometrics are inadequate for assessing patients with severe asymmetric craniofacial syndromes; therefore, 3D-CT cephalometrics are indicated in such patients.1,8 Adams et al8 found 3D evaluation to be more precise than 2D evaluation. Their results indicate that, when the actual distance is measured on a human skull in its true dimensions of 3D space, is more precise and 4 to 5 times more accurate than the 2D approach.8 Various authors have stated that in patients with clefts, 3D imaging provides a better understanding of bone and soft-tissue defects than 2D imaging.1,8–10 Most patients requiring routine orthodontic treatment do not require successive CT scans. However, predicting growth is very important in patients with severe soft and bony tissue defects, especially in cases where repeated surgical intervention or reconstruction is envisaged.11 In such cases, patients needed to be followed closely so that surgeons can track the changes that occur between examinations in order to ensure accurate planning and successful intervention.12,13 The MIMICS 10.02 (Materialise, Belgium) software program has advanced features that perform 3D reconstruction of CT data as well as most types of anthropometric analysis. There are several studies on patients with severe craniofacial defects using software allowing 3D reconstruction of CT scan data and anthropometric analyses, reporting conflicting findings on the advantages and disadvantages of the usage of software versus cephalometric analyses performed on conventional two dimensional lateral cephalography.1,5,8,13,14 This retrospective study aimed to compare Mc-Namara, Steiner and Frontal analyses performed using conventional 2D lateral and frontal cephalometric radiograms and 3D CT images in a population of children aged 7–12 years with cleft lip and palate.

MATERIALS AND METHODS

Study participants were selected from among the 118 patients treated for CLP at the Medical Faculty’s Department of Reconstructive and Plastic Surgery between 2005–2009. Of these, 21 patients had lateral and frontal cephalometric films and 3D CT images in the department archives; however, 6 patients were excluded because of the time lapse (3–6 months) between when the cephalograms and 3D CT images were obtained. In total, orthodontic analysis was performed using radiograms and cranial CT images of 15 children with CLP (9 boys, 6 girls; age range: 7–12 years). CT scans (Discovery, General Electric’s, USA) were obtained using a strict, standardized protocol, with the patient in a horizontal position and the patient’s head positioned with the help of a head positioner. Three-dimensional reconstruction and orthodontic measurement were performed using the MIMICS software program (MIMICS 10.02, Materialise, Belgium). Standard lateral and frontal cephalograms were taken at the same magnification with the same equipment using a cephalostat incorporated into a conventional x-ray device (Proline 2002, Planmeca OY, Finland). Cephalometric films were traced using a protractor with an accuracy of 0.5 mm and 0.5°. All landmark identification and analysis was performed by the same orthodontist, and measurements were repeated after 15 days to control for examiner error. For some symmetrical parameters (Co-Gn, Co-A, SN-Go-Gn, SpPo-GoMe) mean values of left- and right-side measurements from 3D images were calculated in order to compare 2D and 3D measurements. McNamara, Steiner and Frontal analyses landmarks included in Mimics software program were measured. Measurements included both vertical and antero-posterior components of the craniofacial form.

Statistical analysis

All statistical analysis was carried out using the statistical software package SPSS version 11.5 (SPSS Inc., Ill, USA). Mann Whitney U tests was performed to identify differences between the 2D and 3D groups, with a level of significant set at P<.05. Intra-class correlation (ICC) coefficients were calculated to determine the reliability of measurements.

RESULTS

ICC coefficients were very high for both 2D (0.94–0.99) and 3D measurements (0.88–0.99). Results of Frontal, McNamara and Steiner analyses are given, respectively, in Tables 1, 2 and 3. Frontal analysis found all parameters except for OcPtilt to be significantly higher in the 2D group when compared to the 3D group. McNamara analysis found no differences between ANS-Me and Co-Gn measurements in the 2D and 3D groups; however, Co-A was significantly higher (P<.01) and Maxillo-Mandibular difference (max-mand) and SpPo-GoMe were significantly lower (P<.05 and P<.001, respectively) in the 3D group when compared to the 2D group. Steiner analysis found no significant differences between 2D and 3D groups for the parameters SNB, SND and Max1-SN; however, ANB, Mand1-NB, Max1-NA, SNGoGN and SNA were significantly greater and interincisal angle (1–1) significantly smaller in the 3D group when compared to the 2D group.
Table 1

Comparison of 2D and 3D cephalometric evaluation for Frontal analysis.

ParametersMean± Std. DevMedianMinimumMaximumtest
B3L-B3R (mm) CEPH27.41±3.5628.0023.0034.00P<.012*
B3L-B3R (mm) 3D25.43±2.0024.7222.9028.87
B6L-B6R (mm) CEPH54.00±6.9653.0042.5064.50P<.01**
B6L-B6R (mm) 3D47.33±3.2646.3743.7452.98
ImL-ImR (mm) CEPH82.91±6.8784.0072.0092.00P<.03*
ImL-ImR (mm) 3D78.03±4.5678.6569.9583.61
JL-JR (mm) CEPH65.00±11.6970.0044.5076.50P<.02*
JL-JR (mm) 3D58.96±7.0459.0442.5467.68
Nasal width (mm) CEPH29.32±3.4729.0024.0034.50P<.000***
Nasal width (mm) 3D23.433.4822.4819.7931.73
OcP-tilt (mm) CEPH1.14±0.841.000.002.50P<.107ns
OcP-tilt (mm) 3D1.84±0.681.951.023.12
ZA-AZ (mm) CEPH128.73±10.66133.00111.00141.00P<.000***
ZA-AZ (mm) 3D114.42±5.69116.79103.00120.62
ZL-MoL (mm) CEPH70.41±7.1266.5060.5082.00P<.000***
ZL-MoL (mm) 3D60.83±4.4060.7251.7567.32
ZR-MoR (mm) CEPH70.73±7.3667.0060.0084.00P<.000***
ZR-MoR (mm) 3D60.52±4.2261.2654.0468.34
ZL-ZR (mm) CEPH93.59±9.6588.0080.00110.00P<.022*
ZL-ZR (mm) 3D88.27±4.1989.2982.1994.56
Table 2

Comparison of 2D and 3D cephalometric evaluation for Mc Namara analysis.

ParametersMean± Std. DevMedianMinimumMaximumtest
ANS-Me (mm) CEPH58.18±7.3757.0047.0072.00P<.54ns
ANS-Me (mm) 3D56.82±7.1056.5447.9169.05
Co-A (mm) CEPH81.41±7.2582.5071.0091.00P<.01**
Co-A (mm) 3D87.21±4.5189.1580.2693.60
Co-Gn (mm) CEPH104.14±6.75105.5093.00113.00P<.93ns
Co-Gn (mm) 3D104.09±4.37103.0999.09115.23
Max-Mand (mm) CEPH22.23±6.6222.0012.0031.00P<.02*
Max-Mand (mm) 3D16.87±3.8716.7211.4623.02
SpPo-GoMe (°) CEPH69.77±7.6970.0056.0080.00P<.000***
SpPo-GoMe (°) 3D48.82±7.3847.4537.2165.82
Table 3

Comparison of 2D and 3D cephalometric evaluation for Steiner analysis.

ParametersMean± Std. DevMedianMinimumMaximumtest
SNA (°) CEPH75.86±6.4476.5064.0085.50P<.007**
SNA (°) 3D77.17±5.6177.1667.6386.67
SNB (°) CEPH73.27±4.0975.0064.0077.00P<.062ns
SNB (°) 3D74.52±3.1475.1468.4978.22
ANB (°) CEPH2.59±4.722.00−7.009.50P<.036*
ANB (°) 3D2.73±4.802.72−7.5711.62
SND (°) CEPH70.68±3.9972.0062.0074.00P<.11ns
SND (°) 3D70.22±4.1171.7559.7773.63
Max1-NA (°) CEPH0.55±18.446.00−42.0024.00P<.006**
Max1-NA (°) 3D20.72±7.9717.7311.8137.63
Max1-SN (°) CEPH75.91±17.4282.0035.0091.00P<.286ns
Max1-SN (°) 3D80.90±9.0884.4263.1489.56
1–1 (°) CEPH162.05±15.39156.50148.00192.00P<.008**
1–1 (°) 3D140.56±18.91140.40109.37170.41
Mand1-NB (°) CEPH17.18±6.2517.007.0030.00P<.005**
Mand1-NB (°) 3D21.62±7.4023.458.1030.94
SN-Go-Gn (°) CEPH35.64±4.6735.0030.0045.00P<.000***
SN-Go-Gn (°) 3D48.27±3.3146.7743.5954.62

DISCUSSION

CT allows for the measurement of real lengths and angles that is not available with conventional 2D cephalometry. Furthermore, 3D CT makes it possible to assess age-related changes and post-surgical changes in patients with facial deformity, regardless of differences in patient position between examinations.1,15 Studies have shown that the additional spatial information provided by 3D analysis makes it possible to obtain a more accurate assessment of the formation, position, characteristics and dimensions of facial soft tissue as well as the eruption of permanent teeth in all patients using 3D images rather than standard 2D methods.4,16 Craniofacial CT analysis has also been reported to be an excellent method for quantifying and analyzing surface and deep craniofacial structures, especially in cleft palate patients.17 While CT is not routinely recommended, it is helpful in those rare cases where plain films may not be reliable enough to assist in making clinical decisions.18 Errors related to projection, landmark identification and measurement technique can affect the accuracy of measurements taken from frontal cephalograms.19 Probably the most important problem of the PA cephalogram is differential magnification; however, correction of differential magnification is neither widely understood nor routinely practiced.20 Lateral cephalometric films have been reported to have similar magnification problems.21 In line with earlier studies22 2D PA cephalogram measurements were statistically larger than 3D measurements for all parameters (with the exception of OcP-tilt), with the most noticeable difference (14 mm) in the parameter ZA-AZ. These differences are likely the result of magnification error. Vlijmen et al22 found statistically significant differences between conventional frontal cephalometric radiographs and CBCT-constructed 3D models for all measurements except antegonial notch and left lateral orbital margin. The authors explained the discrepancies as resulting from differences in how the landmarks are viewed in 2D images (on a single tomographic plane) and 3D images (on multiple planes). For this reason, comparison of 2D and 3D images should be undertaken with great caution. Reproducing head posture and errors in identifying landmarks have also been shown to be problematic with posteroanterior cepahlograms.23–25 Major et al25 found that rotation of the vertical axis did not affect the relationship of landmarks to the best horizontal line, but did affect their relationship to the best vertical line. In a study conducted with dry skulls, van Vlijmen et al26 showed that patient position, i.e., head tilt/rotation, can significantly influence measurements in frontal cephalometric analysis, and Malkoç et al27 similarly reported that head rotation is likely to affect both linear and angular measurements in lateral cephalometric analysis. In the present study, most parameters (SN-GoGn, ANB, SNA, Max1-NA, Mand1-NB and interincisal angle) varied significantly between the lateral cephalograms and the 3D images. The higher values in the 2D group may be the result of head rotation and magnification error with the lateral cephalograms. Moreover, the complex structure of CLP makes it difficult to define this region in cephalometric films. In a study of 6 human dry skulls comparing cephalometric radiographs and CT imaging processed using Mimics, Varghese et al14 found linear measurements on spiral CT images to be comparable to anatomical measurements and more reliable than cephalometric measurements, which were acceptable for the midsagittal region in the anteroposterior plane, but which varied significantly from anatomical and CT measurements of most parameters. Another study on dry skulls by Major et al28 regarding the posteroanterior cephalometric films. Given that landmark identification is easier on dry skulls devoid of soft tissue, it is likely that the results of the above-mentioned studies represent the minimum error and for this reason differences of up to four times greater have been reported for some measurements in the presence of soft tissue.30 The clinical significance of measurement accuracy varies depending upon the level of accuracy required. For example, if CT data is used for pre-surgical planning, information needs to be sufficiently accurate to prevent surgical interventions that result in insufficient jaw lengthening or shortening that might require additional surgery. Despite the difficulty in accurately identifying landmarks in patients with severe maxillary or mandibulary defects, this information is critical for all orthognathic surgery patients, especially CLP patients.31 Leonardia et al32 have stated that it is generally easier to identify landmarks lying on a sharp curve or at the intersection of two curves than landmarks located on flat or broad curves. The same authors noted that points located in areas of high contrast are easier to identify than points located in areas of low contrast and that the mastoid, latero-orbitale and antegonion landmarks are easier to identify than dental landmarks.32 The statistically significant differences noted by Steiner analyses of SNGoGN and SNA measurements using 2D and 3D images in the present study may be related to the nature of the point Sella, one of the most important reference points. Although generally considered to be the center of a cephalometric view, Sella Turcica is, in fact, not a planar structure localized on the mid-cranial sagittal plane. This represents a significant difficulty in the assessment of the point Sella on conventional cephalometric film on three-dimensional reconstruction of computerized tomography scans. In line with this, Yitschaky et al33 found measurements involving Sella Turcica to vary significantly between 2D and 3D images, with the 3D location of the mid-sella point to be more precise using CT. The distance from the condylion to Point A (Co-A) represents the position of the maxilla in relation to the mandibular condyles. Even if we consider that the crest of the contour of maxillary anterior concavity in the lateral cephalogram can also be determined likewise as a point in the 3D reconstruction, the mandibular condyles are at a great distance from the sagittal mid-cranial plane. However, in lateral cephalograms, this distance is on a 2D plane they are neither in the in vivo situation nor in 3D reconstruction. Considering that it is nearly impossible to separately measure Co-A for the right and left condyles using cephalometric films, the present study calculated a mean value for the separate measurements obtained using 3D images and used this mean value to obtain an accurate comparison between 2D and 3D measurements. When considering the reference points in 3D, we tried to reach to the nearest points corresponding to their 2D cephalometric landmarks in order to maximize the comparability of 2D and 3D measurements. It appears likely that 3D evaluation represents a more accurate method for monitoring craniofacial development than 2D lateral cephalometric films, especially among children with craniofacial deformities, where facial asymmetry is generally more pronounced. Maeda et al10 suggested that 3D CT might have a higher sensitivity than PA cephalometric film imaging in the diagnosis of maxillary asymmetry. For example, on PA cephalograms, it was not possible to evaluate anteroposterior asymmetry, for example, the distance between a landmark and the coronal reference plane, and the authors suggested that because 3D-CT permits more detailed evaluation, it has the potential to replace conventional cephalographic classification of facial asymmetry.10 With regard to intra- and inter-observer reliability, the degree of reported variation in landmark identification varies widely among studies. Major et al28 stated that landmarks with identification errors greater than 1.5 mm should probably be avoided and that landmarks with identification errors greater than 2.5 mm are definitely inappropriate for PA cephalograms. In the present study intra-observer variability was very high; however, the measurements should be analyzed carefully due to the differences found between the groups for clinical evaluation. Each individual landmark is known to have its own characteristic noncircular envelope of error,29,34,35 and, as Lou et al29 point out, differences in measurement of short linear distances will be reflected as greater percentages of error than similar absolute differences in measurement of longer distances. Similarly, for angular measurements, the absolute difference in measurement of an acute angle appears as a greater percent of error than that of a more obtuse angle. In addition to landmark shape, measurement accuracy is also affected by defects and asymmetries. For example, in CLP patients, locating the deepest point of the maxilla (point A) is very difficult. Among the measurements obtained in the present study, the greatest errors were found for SNA angle and Co-A length. Imaging of the pre-maxillary region in 3D imaging is very different from imaging of 2D cephalometric films. In general, superimposition makes it difficult to accurately identify upper incisor position on cephalometric films, and in the present study, excessive inclination and rotation in CLP patients added to the difficulty of accurately determining the position of the upper incisors, which explains the statistically significant differences in maxillary incisor measurements between the 2D and 3D groups. Other factors involved in landmark identification accuracy include clinical experience,28 the clarity of the definition used to describe the landmark, image quality, the geometry (straight or ovoid) of the object to be identified and the contrast between adjacent objects.36 In the present study, for instance, the deepest point of the maxilla (jugular point) used in frontal analysis appeared as a fossa in the 3D image, but as a point in the cephalometric films. Similarly, due to the 3D nature of the meatus acusticus externus, the point Porion appeared in different locations in the 2D and 3D images. This would explain the great variability between 2D and 3D measurements of the parameter SpPo-GoMe in the present study (Table 2). Another reason for this variability could be the fact that SpPo-GoMe is located in different tomographic planes, as discussed above. Troulis et al37 stated that 3D-CT driven software has great potential for use in the diagnosis and analysis of congenital and acquired craniomaxillofacial skeletal deformities; for the calculation of vectors of surgical movement, localization of osteotomies and placement of distraction devices as part of surgical treatment planning; and, possibly, as a basis for a surgical navigation system. The CT images taken in the present study were used by plastic surgeons for treatment planning. If standard and repeatable reference points can be obtained, 3D cephalometric evaluation can provide accurate data on craniofacial measurements, as other studies have suggested.38,39 However, major concerns continue to exist regarding their widespread usage.1,2,5,12 First, CTs expose patients to considerably higher amounts of radiation than conventional cephalograms. This limits the number of CT scans that can be taken over short intervals to monitor craniofacial development, especially in growing children. Second, CT scans requires very expensive equipment, making the cost of an individual CT scan much more expensive than a lateral cephalogram. For this reason, CBCT seems likely to replace medical CT for use in orthodontics and plastic surgery in the near future. Perhaps most important, there is still a need for the establishment of widely accepted reference points and norms appropriate for 3D evaluation in order to facilitate the consistency of inter-method evaluation.

CONCLUSIONS

This study found significant differences between measurements made on 2D and 3D images taken of CLP patients. Although there are some previous studies reporting that cephalometric and three-dimensional evaluations show a considerably well correlation, but one should be very careful when interpreting the results. This can be due to the fact that the subjects of this study are consisted of children with cranio-facial anomalies while all the norms in conventional cephalometric analyses are considered on normal subjects. Within the limitations of this study, it can be concluded that further studies are needed in order to more accurately determine the anatomical landmarks on the three dimensional reproduction and thus to establish new standards for 3D evaluations of CLP patients.
  39 in total

1.  Effect of head rotation on posteroanterior cephalometric radiographs.

Authors:  Young-Jooh Yoon; Dong-Hoon Kim; Pil-Sik Yu; Heung-Joong Kim; Eui-Hwan Choi; Kwang-Won Kim
Journal:  Angle Orthod       Date:  2002-02       Impact factor: 2.079

2.  Craniofacial measurements based on 3D-CT volume rendering: implications for clinical applications.

Authors:  M G P Cavalcanti; S S Rocha; M W Vannier
Journal:  Dentomaxillofac Radiol       Date:  2004-05       Impact factor: 2.419

3.  Comparison between traditional 2-dimensional cephalometry and a 3-dimensional approach on human dry skulls.

Authors:  Gregory L Adams; Stuart A Gansky; Arthur J Miller; William E Harrell; David C Hatcher
Journal:  Am J Orthod Dentofacial Orthop       Date:  2004-10       Impact factor: 2.650

4.  The effect of head rotation on cephalometric radiographs.

Authors:  Siddik Malkoc; Zafer Sari; Serdar Usumez; Alp Erdin Koyuturk
Journal:  Eur J Orthod       Date:  2005-06       Impact factor: 3.075

5.  Three-dimensional cephalometry: spiral multi-slice vs cone-beam computed tomography.

Authors:  Gwen R J Swennen; Filip Schutyser
Journal:  Am J Orthod Dentofacial Orthop       Date:  2006-09       Impact factor: 2.650

6.  Three-dimensional facial morphology following surgical repair of unilateral cleft lip and palate in patients after nasoalveolar molding.

Authors:  G D Singh; D Levy-Bercowski; M A Yáñez; P E Santiago
Journal:  Orthod Craniofac Res       Date:  2007-08       Impact factor: 1.826

Review 7.  Accuracy of measurements and reliability of landmark identification with computed tomography (CT) techniques in the maxillofacial area: a systematic review.

Authors:  Leo Lou; Manuel O Lagravere; Sharon Compton; Paul W Major; Carlos Flores-Mir
Journal:  Oral Surg Oral Med Oral Pathol Oral Radiol Endod       Date:  2006-10-27

8.  Craniofacial skeletal measurements based on computed tomography: Part I. Accuracy and reproducibility.

Authors:  A A Waitzman; J C Posnick; D C Armstrong; G E Pron
Journal:  Cleft Palate Craniofac J       Date:  1992-03

9.  In vivo comparison of conventional and cone beam CT synthesized cephalograms.

Authors:  Vandana Kumar; John Ludlow; Lucia Helena Soares Cevidanes; André Mol
Journal:  Angle Orthod       Date:  2008-09       Impact factor: 2.079

10.  Comparison of common hard tissue cephalometric measurements between computed tomography 3D reconstruction and conventional 2D cephalometric images.

Authors:  Oded Yitschaky; Meir Redlich; Yossi Abed; Marina Faerman; Nardy Casap; Nurith Hiller
Journal:  Angle Orthod       Date:  2011-01       Impact factor: 2.079

View more
  4 in total

1.  Accuracy of 3D cephalometric measurements based on an automatic knowledge-based landmark detection algorithm.

Authors:  Abhishek Gupta; Om Prakash Kharbanda; Viren Sardana; Rajiv Balachandran; Harish Kumar Sardana
Journal:  Int J Comput Assist Radiol Surg       Date:  2015-12-24       Impact factor: 2.924

Review 2.  Three-dimensional imaging methods for quantitative analysis of facial soft tissues and skeletal morphology in patients with orofacial clefts: a systematic review.

Authors:  Mette A R Kuijpers; Yu-Ting Chiu; Rania M Nada; Carine E L Carels; Piotr S Fudalej
Journal:  PLoS One       Date:  2014-04-07       Impact factor: 3.240

3.  Dento-skeletal characteristics of cleft patients with missing teeth.

Authors:  Waeil Batwa; Hala A Almarhoon; Khalid A Almoammar; Nasser Alqahtani; Sahar F Albarakati; Thikriat Al-Jewair
Journal:  Clin Cosmet Investig Dent       Date:  2018-11-08

4.  Three-Dimensional Analysis of Posterior Mandibular Displacement in Rats.

Authors:  Ioannis Lyros; Efstratios Ferdianakis; Demetrios Halazonetis; Theodoros Lykogeorgos; Antigoni Alexiou; Konstantina-Eleni Alexiou; Maria Georgaki; Emmanouil Vardas; Zafeiroula Yfanti; Apostolos I Tsolakis
Journal:  Vet Sci       Date:  2022-03-20
  4 in total

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