Literature DB >> 30627245

Reduction of radiation dose for coronary computed tomography angiography using prospective electrocardiography-triggered high-pitch acquisition in clinical routine.

Narumol Chaosuwannakit1, Pattarapong Makarawate2.   

Abstract

PURPOSE: To evaluate the image quality, radiation exposure, and means of application in a group of patients who underwent coronary computed tomography angiography (CCTA) performed with low-dose prospective electrocardiography (ECG)-triggered acquisition in which a standard sequence was added if the low-dose sequence did not allow reliable exclusion of coronary stenosis with respect to image quality.
MATERIAL AND METHODS: The present study was approved by the Ethics Committee of the Faculty of Medicine, and informed consent was obtained from all patients. The authors performed a retrospective review of 256 consecutive patients referred for CCTA using dual-source CT scanner (Definition FLASH, Siemens, Germany). CCTA was performed using prospective ECG-triggered high-pitch acquisition. In patients with higher heart rates (> 65 bpm) or in whom irregular heart rates were noted prior to the scan, a subsequent CCTA was performed immediately (double flash protocol). The effective radiation dose was calculated for each patient. All images were evaluated by two independent observers for quality on a four-point scale with 1 being non-diagnostic image quality and 4 being excellent.
RESULTS: Mean effective whole-body dose of CCTA was 1.6 ± 0.4 mSv (range, 0.4-5.4) for the entire cardiac examination and 0.9 ± 0.3 mSv (range, 0.4-2.8) for individual prospective ECG-triggered high-pitch CCTAs. In 27 of these patients with higher heart rates or occasional premature ventricular contractions or atrial fibrillation, subsequent CCTAs were performed immediately. The average image quality score was good to excellent with less than 1% unevaluable coronary segments. The double flash protocol resulted in a fully diagnostic CCTA in all cases.
CONCLUSIONS: The prospective ECG-triggered high-pitch CCTA technique is feasible and promising in clinical routine with good to excellent image quality and minimal radiation dose. The double flash protocol might become a more robust tool in patients with higher heart rates or arrhythmia.

Entities:  

Keywords:  computed tomography; coronary CTA; image quality; radiation exposure

Year:  2018        PMID: 30627245      PMCID: PMC6323584          DOI: 10.5114/pjr.2018.76791

Source DB:  PubMed          Journal:  Pol J Radiol        ISSN: 1733-134X


Introduction

Coronary computed tomography angiography (CCTA) has increasingly gained importance as a noninvasive, fast, and accurate study for diagnosing coronary artery disease (CAD) [1-7]. The relatively high radiation exposure in CCTA compared with invasive conventional coronary angiography remains a challenge [8,9]. Hence, not only careful selection of patients suitable for CCTA examination, but also reduction of the patient radiation exposure in CCTA without compromising diagnostic accuracy should be aimed for. CCTA performed with 64-detector CT without the use of electrocardiography (ECG) pulsing typically results in radiation doses ranging from 15 to 21 mSv [5] and 9 mSv with the use of ECG pulsing [10,11]. It is well known that the radiation dose in cardiac imaging can substantially be reduced by the use of prospective ECG-triggering data acquisitions [12]. The method of prospective ECG-triggering (PT) has been evaluated for image quality of the coronary arteries and diagnostic accuracy as well as effective radiation dose in several studies [13-19]. PT is a technique used with CCTA, which uses forward-looking prediction of R-wave timing, step-and-shoot non-spiral acquisition during imaging, and unique cone beam reconstruction [20]. By contrast, CCTA with standard retrospective ECG-triggering (RT) uses backward-looking measurement of R-wave timing, spiral scanning during table motion, and more traditional cone beam reconstruction. In addition, with PT the X-ray beam is turned on for only a short portion of diastole, and it is turned off during the rest of the R-R cycle, whereas with RT the X-ray beam is turned on throughout the R-R interval. Initial results showed that PT may be able to provide sufficient image quality with low radiation dose [13-16]. The latest generation of dual-source CT (DSCT) scanners has introduced a new scan mode – a prospective ECG-triggered helical data acquisition with very high pitch values of more than 3.0. This technique enables acquisition of the entire volumetric dataset of the heart within a fraction of a single cardiac cycle. The high pitch allows acquisition with very low radiation exposure (< 1 mSv) in patients with low, regular heart rates [17,18]. Accurate and consistent visualisation of the coronary arteries at high image quality is critical for coronary CTA. Hence, the purpose of this study was to describe the results of prospective ECG-triggered high-pitch mode for coronary CTA in clinical routine practice in a tertiary care centre with respect to image quality and radiation exposure.

Material and methods

Patients

This retrospective study included 256 patients (130 males and 126 females) referred for CCTA between February 2013 and February 2014. The indications for CCTA were an abnormal, equivocal, or non-diagnostic stress test, chest pain, evaluation of cardiomegaly and congestive heart failure, as well as the evaluation of cardiac aetiology of syncope. Patients with an intermediate probability of coronary artery disease (CAD) were also referred for a CCTA as a first test. The above are considered appropriate indications for CCTA, based on the criteria of the American College of Cardiology (ACC) [19-21] and the recent American Heart Association Scientific Statement on Cardiac CT [19-21]. Exclusion criteria for CCTA included the presence of renal failure, pregnancy, and a history of allergic reaction to iodine-containing contrast agents. The present study was approved by the Ethics Committee of the Faculty of Medicine. Written, informed consent was not required because of the retrospective nature of the investigation.

Cardiac computed tomography angiography protocol

CCTA examinations were performed on a DSCT scanner (Definition FLASH, Siemens Healthcare, Forchheim, Germany). The system was equipped with two X-ray tubes and two corresponding detectors mounted on a single gantry with an angular offset of 90° [22]. The CCTA used automatic tube current modulation in x, y, z directions (Care Dose 4D, Siemens Healthcare). The scanner technology enabled a prospectively ECG-triggered high-pitch (3.4) spiral acquisition (FLASH Spiral Cardio, Siemens Healthcare, Forchheim, Germany). The CCTA scan parameters were as follows: two X-ray sources, detector collimation 128 × 0.6 mm with double sampling by rapid alteration of the focal spot in the longitudinal direction (Z-flying focal spot) [23], rotation time 0.28 seconds, and tube voltage 100 or 120 kV. Image acquisition was performed during inspiratory breath-hold. To familiarise the patient with the protocol, breath-holding was practiced before the examination. A contrast agent bolus of 50-80 ml was injected with a mean flow rate of 5 ml/s followed by a 50 ml saline flush. Circulation time from the antecubital fossa to the aortic root was determined in each patient with administration of a timed bolus of 20 ml of iodinated contrast. The entire volume of the heart was covered during one breath-hold in approximately five seconds with simultaneous recording of the ECG trace. Patients were scanned in the supine position twice, first without contrast medium to calculate the calcium score and secondly after contrast medium injection. Studies were acquired in the cranio-caudal direction from the level of the carina to just below the diaphragm. In 27 out of the 256 patients a second CCTA was acquired immediately during the same contrast injection timed (“double flash protocol”). In these cases, the delay time between the two acquisitions was set to four seconds, the minimal allowable time on this system. Therefore, the total scan time was longer and the total contrast volume was increased by 25 to 35 cc compared to single high-pitch CCTA acquisition ([4 seconds + ~1second scan] × flow rate). This protocol was used in patients with higher heart rates (> 65 bpm) or in whom irregular heart rate was noted prior to the scan. Neither nitroglycerine nor beta-blockers were given to the patient.

Cardiac computed tomography angiography image reconstruction

All CCTA datasets were transferred to a dedicated workstation. To evaluate the coronary arteries, the images were reconstructed with a small FOV (120-190 mm), which was restricted to the heart region. The images were reviewed in the curve multiplanar reconstruction, axial, coronal, and sagittal planes, using a mediastinal window (width: 450, level: 35), lung window (width: 1500, level: –700), and bone window (width: 1500, level: 450) for all examinations.

Coronary computed tomography angiography image analysis

CCTA image analysis was performed by two independent experienced cardiovascular and thoracic radiologists (with a respective nine and seven years of experience in examining cardiovascular and thoracic CT scans) on a dedicated post processing workstation. Image quality was analysed for every segment according to a four-point scale (Table 1) [24,25] (1 = non-diagnostic image quality; 2 = acceptable image quality with artefacts but not compromising diagnostic vessel assessment; 3 = good image quality with minimal artefacts but fully evaluable coronary vessel structures; 4 = excellent image quality without artefacts). Typical examples for each grade (1-4) are given in Figures 1 and 2. The two readers rated the image quality independently with a subsequent consent read.
Table 1

Four-point rating scale for assessment of image quality. Scores of 3 and 4 correspond to a non-diagnostic image quality that does not allow assessment of the entire coronary artery for the presence of a flow-limiting stenosis

ScoreCharacteristics
1Non-diagnostic image quality for multiple segments within one coronary artery considered unevaluable
2Acceptable image quality with artefacts but not compromising diagnostic vessel assessment
3Good image quality with minimal artefacts but fully evaluable coronary vessels
4Excellent image quality without artefacts
Figure 1

Example of excellent image quality using prospective electrocardiography (ECG)-triggered high-pitch (pitch 3.4) acquisition of a 65-year-old male patient with atypical chest pain. Images show a curved multiplanar reconstruction of the right coronary artery (A) and a 3D volume-rendered reconstruction of coronary arteries (B) demonstrating excellent image quality. Average heart rate is 59 bpm, estimated radiation dose is 0.5 mSv, and amount of contrast agent is 50 ml

Figure 2

The figure demonstrates typical images for each grade (1-4) regarding the right coronary artery (RCA). A) Excellent image quality (grade 4). B) Good image quality (grade 3) with minimal artefacts but fully evaluable coronary vessel (dash arrows). C) Moderate image quality (grade 2) due to beam hardening artefacts of pacemaker leads (black arrows). D) Severe blurring artefacts (white arrows) cause non-diagnostic image quality (grade 1) of the proximal and mid RCA

Four-point rating scale for assessment of image quality. Scores of 3 and 4 correspond to a non-diagnostic image quality that does not allow assessment of the entire coronary artery for the presence of a flow-limiting stenosis Example of excellent image quality using prospective electrocardiography (ECG)-triggered high-pitch (pitch 3.4) acquisition of a 65-year-old male patient with atypical chest pain. Images show a curved multiplanar reconstruction of the right coronary artery (A) and a 3D volume-rendered reconstruction of coronary arteries (B) demonstrating excellent image quality. Average heart rate is 59 bpm, estimated radiation dose is 0.5 mSv, and amount of contrast agent is 50 ml The figure demonstrates typical images for each grade (1-4) regarding the right coronary artery (RCA). A) Excellent image quality (grade 4). B) Good image quality (grade 3) with minimal artefacts but fully evaluable coronary vessel (dash arrows). C) Moderate image quality (grade 2) due to beam hardening artefacts of pacemaker leads (black arrows). D) Severe blurring artefacts (white arrows) cause non-diagnostic image quality (grade 1) of the proximal and mid RCA

Measurement of the radiation exposure

The dose-length product (DLP) displayed by the CT unit was recorded for each CCTA. The effective dose of CCTA was estimated by a method proposed by the Fleischner Society [26]. The effective dose is derived from the product of the DLP and a conversion coefficient for the anatomical region examined which is 0.014 mSv mGy–1 cm–1 for the chest [27]. However, it should be noted that a uniform conversion coefficient for all images is not entirely accurate because it does not account for the different conditions in the examinations.

Statistical analysis

Statistical analyses were performed using SPSS software version 16 (SPSS, Inc., Chicago, IL, USA). Continuous data were expressed as mean ± SD, and categorical variables were given as percentages. Comparisons of the patient characteristics, the CT parameters, and the dose estimates were performed using the t-test for normally distributed data, and using the χ2 test if not. Intraclass correlation coefficient (ICC) was used for interobserver agreement of subjective image quality. The coefficient represents concordance, where 1 is perfect agreement and 0 is no agreement at all. A significance level of p < 0.05 was considered a statistically significant result, and all reported p-values were two-sided.

Results

A total of 256 patients (130 men, 126 women) were included in this study over a period of 12 months. Clinical characteristics of the patients are summarised in Table 2. A total of 283 CCTAs in 256 patients were evaluated. In 27 patients two prospectively triggered high-pitch acquisitions were performed. A second CCTA was acquired immediately during the same contrast injection timed (“double flash protocol”). This protocol was used in patients with higher heart rates (> 65 bpm) or in whom irregular heart rates were noted prior to the scan.
Table 2

Patient characteristics of the single prospective electrocardiography (ECG)-triggered high-pitch coronary computed tomography angiography (SP) and double flash protocol (DF)

CharacteristicSPDFp-value
Age (years)*54 ± 10 (25-81)58 ± 15 (45-81)0.15
Men#130 (50.8%)130 (50.8 %)1.0
Height (cm)*161.8 ± 13158.9 ± 150.25
Weight (kg)*59.8 ± 1252.5 ± 100.16
Body mass index (kg/m²)*23.1 ± 6 (15-35)21.2 ± 8 (18-35)0.41
Heart rate (bpm)*58 ± 9 (50-64)75 ± 11 (65-100)0.02*

Data are mean ± standard deviation, numbers in parentheses are range.

Data are numbers of patients, and numbers in parentheses are percentages.

p value < 0.05

bmp – beats per minute

Patient characteristics of the single prospective electrocardiography (ECG)-triggered high-pitch coronary computed tomography angiography (SP) and double flash protocol (DF) Data are mean ± standard deviation, numbers in parentheses are range. Data are numbers of patients, and numbers in parentheses are percentages. p value < 0.05 bmp – beats per minute Patients in the single prospectively ECG-triggered high-pitch group (SP) and those in the double flash protocol (DF) group did not differ significantly in age, sex, or mean body mass index. However, there was a significant difference between the groups for mean heart rate (58 ± 9 beats/min for the SP group, 75 ± 11 for the DF group; p = 0.02) (Table 2). The distribution of calcium score of all patients are shown in Table 3.
Table 3

Distribution of Agaston calcium score

Total calcium scoreN = 256 (%)
024 (9.3)
1-10091 (35.6)
101-20079 (30.9)
201-30030 (11.7)
301-40018 (7.0)
> 40014 (5.5)
Distribution of Agaston calcium score In SP group, the mean DLP was 63.7 ± 21 mGy . cm, resulting in an effective radiation dose per examination of 0.9 ± 0.3 mSv. In the DF group, mean DLP was 114.3 ± 28 mGy . cm, resulting in an effective radiation dose per examination of 1.6 ± 0.4 mSv. The SP group showed a significantly lower DLP and consequently a significantly lower mean effective radiation dose compared with the RT group (p = 0.015). Mean effective whole-body dose of a cardiac scan ranged from 0.3 to 3.4 mSv. In total, 3679 coronary artery segments, 849 coronary artery vessels, and 283 patients were evaluated for image quality. Among segments 2538 (69%) had an image quality score of 4 (“excellent”), 699 segments (19 %) a score of 3 (“good”), 418 segments (11%) a score of 2 (“moderate”), and 24 segments (< 1%) were scored as “unevaluable”. Mean rating score for all patients and segments was 3.7 (“good” to “excellent”). An example is given in Figure 2. Unevaluable segments were observed in 20 patients. The predominant reasons given by the observers were extensive cardiac motion artefact, noise, and streak artefacts due to extensive calcifications. However, all of these patients belonged to the subgroup that was examined by the “double flash protocol”. In these patients, all segments were evaluable after taking both scans into account (Figure 3). Twenty of these 24 segments were small distal branches with little clinical relevance in the individual patients. The two radiologists demonstrated a good agreement regarding subjective image quality (ICC-coefficient of 0.92).
Figure 3

Non-diagnostic image quality (A; white arrow) compared with excellent image quality (B; black arrow) of right coronary artery (RCA) in the same patient after immediate subsequent coronary computed tomography angiography were performed (double flash protocol)

Non-diagnostic image quality (A; white arrow) compared with excellent image quality (B; black arrow) of right coronary artery (RCA) in the same patient after immediate subsequent coronary computed tomography angiography were performed (double flash protocol)

Discussion

The use of CCTA with RT results in good image quality and few non-evaluable coronary artery segments; however, the radiation dose to patients is relatively high [27-29]. The second-generation DSCT prospective ECG-triggered high-pitch CCTA offers a unique opportunity to achieve maximal dose reduction while preserving fully diagnostic image quality. DSCT is necessary for this technique because rapid table motion leads to gaps in the trajectory of the first detector. A second detector is needed to fill these gaps. Radiation exposure is very low using this technique since slice overlap (and therefore redundant exposure) is avoided. Using a pitch value of 3.4, data acquisition of the entire z-axis of the adult heart within a fraction of one cardiac cycle (approximately 260 ms) is possible. Whereas in conventional prospective triggered axial acquisition protocols (“step and shoot”) some redundant radiation exposure is still applied at slab interfaces, redundant radiation exposure is nearly eliminated in a high-pitch protocol and occurs only at the beginning and end of the scanned volume [30,31]. The aim of the present study was to describe the results of prospective ECG-triggered high-pitch mode for coronary CTA in clinical routine practice in a tertiary care centre, with respect to image quality and radiation exposure. However, dose reduction must not compromise image quality. Non-diagnostic scans because of ‘too little radiation’ are a worst-case scenario because they expose the patient to radiation without the hoped for benefit. Several studies have recently described the use of this novel technique [17,18,32,33]. Previous studies have shown that prospective ECG-triggered high-pitch spiral acquisition for cardiac CTA is feasible and provides excellent image quality at a consistent dose below 1.0 mSv in carefully selected patients with a low and regular heart rate (< 60 bpm) and a body weight of less than 100 kg [16,18]. We retrospectively analysed 256 patients in which this technique was applied in clinical routine at a tertiary academic centre. In contrast to previous studies, we included all patients scanned on our service regardless of body habitus and with much less rigorous heart rate and rhythm restrictions. Our data show that CCTA using this technique is technically feasible in a less strict clinical routine setting, resulting in an average “good to excellent” image quality at a reasonably low mean radiation exposure. Our studies demonstrated mean radiation dose of 1.6 mSv for all patients and 0.9 mSv for a single-acquisition CCTA, which is in the same range as reported in the previous publications [17,18,32,33]. The percentage of unevaluable segments and subjective image quality analysis was lower in the present study (0.008%) compared to the data presented by Lell et al. (0.6%) [32]. Nevertheless, the majority of non-diagnostic segments were small, distal, and of little clinical relevance. A wide range of heart rates (range 50-100), excessive calcifications of the coronary arteries in bypass graft cases, which lead to “unevaluable” segments due to blooming and streak artefacts. Compared with studies reporting on DSCT of coronary arteries with standard axial “step and shoot” technique, the prospective high-pitch mode technique enables another dose reduction of approximately 40-50% [13-16]. Earls et al. [16] report a mean effective dose of 2.8 mSv (DLP 170, conversion factor 0.017) for a prospective ECG-triggered single-source coronary CTA, whereas an effective dose of 1.2 mSv (DLP 85.7, conversing factor 0.014) was found in the study presented herein. Because the radiation dose of a prospective ECG-triggered high-pitch CCTA is very low, a protocol consisting of two scans, one immediately after another (“double flash mode”), might be an opportunity to obtain diagnostic image quality despite a patient’s elevated heart rate or occasional PVCs. In comparison to increased “padding” in prospectively triggered single-source CT, which allows assessment of the heart in additional phases within the same cardiac cycle [34], the “double flash protocol” performed in this study concentrates on the same phase (60%) at different cardiac cycles. In the present study, this protocol led to diagnostic image quality of all coronary segments after interpreting both acquisitions. The patients in whom several coronary segments would have been unevaluable after a single scan were diagnostic due to this protocol. Therefore, this technique might obviate a second retrospective ECG-gated cardiac scan with a second injection of contrast, but it still results in low total radiation exposure. Despite promising initial results, our study has potential limitations. First, our results represent a single-centre experience, the generalisability of the present results is limited. Secondly, a retrospective analysis of clinical data, a systemical bias cannot be excluded. Thirdly, a systematic comparison of CCTA and invasive coronary angiography was not performed, so diagnostic accuracy versus a gold standard is not available. Lastly, radiation doses were calculated and not directly measured in our analysis.

Conclusions

In the present study, our findings suggest that CCTA performed with prospective ECG-triggered high-pitch is a feasible and promising technique. When applied in clinical routine, evaluation of the coronary arteries is possible at good to excellent image quality with very low radiation exposure. The “double flash protocol” may be a promising technique for patients with higher heart rates > 65 bpm or with irregular heart rates, who might otherwise not receive a diagnosis with a single acquisition in this mode. Further studies are necessary to validate these findings.

Conflict of interest

The authors report no conflict of interest.
  34 in total

1.  Image quality, motion artifacts, and reconstruction timing of 64-slice coronary computed tomography angiography with 0.33-second rotation speed.

Authors:  Bernd J Wintersperger; Konstantin Nikolaou; Franz von Ziegler; Thorsten Johnson; Carsten Rist; Alexander Leber; Thomas Flohr; Andreas Knez; Maximilian F Reiser; Christoph R Becker
Journal:  Invest Radiol       Date:  2006-05       Impact factor: 6.016

Review 2.  Radiation exposure at chest CT: a statement of the Fleischner Society.

Authors:  John R Mayo; John Aldrich; Nestor L Muller
Journal:  Radiology       Date:  2003-07       Impact factor: 11.105

3.  Accuracy of 64-MDCT in the diagnosis of ischemic heart disease.

Authors:  Konstantin Nikolaou; Andreas Knez; Carsten Rist; Bernd J Wintersperger; Alexander Leber; Thorsten Johnson; Maximilian F Reiser; Christoph R Becker
Journal:  AJR Am J Roentgenol       Date:  2006-07       Impact factor: 3.959

4.  Quantitative parameters of image quality in 64-slice computed tomography angiography of the coronary arteries.

Authors:  Maros Ferencik; Cesar H Nomura; Pal Maurovich-Horvat; Udo Hoffmann; Antonio J Pena; Ricardo C Cury; Suhny Abbara; Koen Nieman; Umaima Fatima; Stephan Achenbach; Thomas J Brady
Journal:  Eur J Radiol       Date:  2006-01-24       Impact factor: 3.528

5.  High-resolution spiral computed tomography coronary angiography in patients referred for diagnostic conventional coronary angiography.

Authors:  Nico R Mollet; Filippo Cademartiri; Carlos A G van Mieghem; Giuseppe Runza; Eugène P McFadden; Timo Baks; Patrick W Serruys; Gabriel P Krestin; Pim J de Feyter
Journal:  Circulation       Date:  2005-10-03       Impact factor: 29.690

6.  Visual and automatic grading of coronary artery stenoses with 64-slice CT angiography in reference to invasive angiography.

Authors:  Stephanie Busch; Thorsten R C Johnson; Konstantin Nikolaou; Franz von Ziegler; Andreas Knez; Maximilian F Reiser; Christoph R Becker
Journal:  Eur Radiol       Date:  2006-12-16       Impact factor: 5.315

7.  Accuracy of MSCT coronary angiography with 64-slice technology: first experience.

Authors:  Sebastian Leschka; Hatem Alkadhi; André Plass; Lotus Desbiolles; Jürg Grünenfelder; Borut Marincek; Simon Wildermuth
Journal:  Eur Heart J       Date:  2005-04-19       Impact factor: 29.983

8.  Clinical value of MDCT in the diagnosis of coronary artery disease in patients with a low pretest likelihood of significant disease.

Authors:  Konstantin Nikolaou; Carsten Rist; Bernd J Wintersperger; Tobias F Jakobs; Roland van Gessel; Miles A Kirchin; Andreas Knez; Franz von Ziegler; Maximilian F Reiser; Christoph R Becker
Journal:  AJR Am J Roentgenol       Date:  2006-06       Impact factor: 3.959

Review 9.  Diagnostic performance of multislice spiral computed tomography of coronary arteries as compared with conventional invasive coronary angiography: a meta-analysis.

Authors:  Michèle Hamon; Giuseppe G L Biondi-Zoccai; Patrizia Malagutti; Pierfrancesco Agostoni; Rémy Morello; Marco Valgimigli; Martial Hamon
Journal:  J Am Coll Cardiol       Date:  2006-09-26       Impact factor: 24.094

10.  Radiation dose estimates from cardiac multislice computed tomography in daily practice: impact of different scanning protocols on effective dose estimates.

Authors:  Jörg Hausleiter; Tanja Meyer; Martin Hadamitzky; Ester Huber; Maria Zankl; Stefan Martinoff; Adnan Kastrati; Albert Schömig
Journal:  Circulation       Date:  2006-03-06       Impact factor: 29.690

View more
  1 in total

1.  Comparison of radiation dose and its correlates between coronary computed tomography angiography and invasive coronary angiography in Northeastern Thailand.

Authors:  Phatraporn Aupongkaroon; Pattarapong Makarawate; Narumol Chaosuwannakit
Journal:  Egypt Heart J       Date:  2022-01-25
  1 in total

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