Literature DB >> 34021692

Risk for second bladder and rectal malignancies from cervical cancer irradiation.

Michalis Mazonakis1, Efrosyni Lyraraki2, Maria Tolia2, John Damilakis1.   

Abstract

The objective of this study was to estimate the risk of developing second malignancies to partially in-field organs from volumetric modulated arc therapy (VMAT) of cervical cancer and to compare the above risks with those from the conventional three-dimensional conformal radiotherapy (3D-CRT). Seventeen consecutive patients with uterine cervix carcinoma were selected. VMAT and 3D-CRT plans were generated with 6 and 10 MV photons, respectively. The prescribed tumor dose was 45 Gy given in 25 fractions. Differential dose-volume histogram data from the treatment plans were obtained for the partially in-field organs such as bladder and rectum. These data were used to estimate the patient-specific lifetime attributable risk (LAR) for bladder and rectal cancer induction with a non-linear model based on a mixture of plateau and bell-shaped dose-response relationships. The estimated risks per 10000 people were compared with the baseline risks for unexposed population. The patient-specific rectal cancer risk estimates from VMAT were significantly lower than those from 3D-CRT (P = 0.0144). The LARs for developing bladder malignancies from VMAT were significantly high compared to those from conventional irradiation (P = 0.0003). The mean difference between the patient-specific LARs for radiation-induced bladder and rectal malignancies as derived from 3D-CRT and VMAT plans was 6.6% and 2.0%, respectively. The average LAR for developing bladder and rectal malignant diseases due to VMAT was 9.2 × 10-4 and 43.7 × 10-4 , respectively. The corresponding risks following 3D-CRT were 8.6 × 10-4 and 44.6 × 10-4 . These average risks showed that pelvic irradiation increases the baseline probability for cancer induction by 12.6-19.1%. The differences in the second cancer risks associated with the VMAT and 3D-CRT for cervical cancer were found to be small. Both treatment techniques resulted in considerable increased probabilities for developing bladder and rectal malignancies relative to those of unirradiated population.
© 2021 The Authors. Journal of Applied Clinical Medical Physics published by Wiley Periodicals LLC on behalf of American Association of Physicists in Medicine.

Entities:  

Keywords:  3D-CRT; VMAT; cervical cancer; second malignancies

Mesh:

Year:  2021        PMID: 34021692      PMCID: PMC8292701          DOI: 10.1002/acm2.13274

Source DB:  PubMed          Journal:  J Appl Clin Med Phys        ISSN: 1526-9914            Impact factor:   2.102


INTRODUCTION

Cervical cancer is the fourth most frequent malignant disease and, also, the fourth main cause of cancer mortality in females globally. The 5‐year survival rate of this type of cancer is 66.1% and it reaches to 91.8% for localized disease stage. Radiation therapy plays a major role in the management of the uterine cervix carcinoma. Previous studies reported that the therapeutic irradiation of this gynecological cancer is associated with an increased probability for developing second malignancies to surrounding heavily exposed sites such as urinary bladder and rectum. , , The above results were derived from patients irradiated before 2001 without advanced treatment delivery and planning approaches. The use of intensity modulated techniques is currently recommended for external‐beam radiotherapy of cervical cancer. The intensity modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT) may lead to a significant reduction of the acute toxicity compared to that from conventional pelvic irradiation. However, the application of IMRT and/or VMAT in clinical practice always creates a concern about the second cancer risk magnitude due to the long delivery times and the large volume of exposed normal tissues in respect to those related to three‐dimensional conformal radiation therapy (3D‐CRT). To the best of our knowledge, there is only one report dealing with carcinogenesis from IMRT for cervical cancer. They found that IMRT results in a higher cumulative cancer risk than 3D‐CRT. The VMAT for cervical cancer has gained popularity the last decade due to the decrease of the treatment delivery time and monitor unit (MU) usage in comparison with IMRT. , , No attempts have been made to assess the probability for second cancer development in patients subjected to VMAT for gynecological carcinomas. This study was conducted to estimate the risks for developing second malignancies to critical partially in‐field organs such as bladder and rectum following VMAT for cervical cancer and to compare the aforementioned risks with those resulting from conventional 3D‐CRT.

MATERIALS AND METHODS

CT scanning and contouring

Seventeen consecutive patients with localized carcinoma of the uterine cervix were included in this work. The patient’s age ranged from 40 to 68 yrs old with a mean age of 53.1 ± 9.2 yrs (Fig. 1). The study participants underwent a planning computed tomography (CT) examination in supine treatment position. The CT slice thickness was 5 mm without a gap. All CT datasets were transferred to Monaco system (Elekta AB, Stockholm, Sweden) for contouring and planning. The RTOG‐0418 guidelines were used to define the targets and organs at risk (OARs) on the CT scans acquired with a full urinary bladder and an empty rectum. The clinical target volume (CTV) comprised the upper 3.0 cm of vagina and paravaginal soft tissue lateral to the vagina. The internal, external and common iliac lymph nodes, and presacral nodal regions were also included within the CTV. A 7‐mm margin was applied to uniformly expand the CTV to the planning target volume (PTV). The bowel, rectum, urinary bladder, and femoral heads were manually delineated on the CT scans as the organs‐at‐risk (OARs). The study protocol was approved by the ethics committee of our institution. A written informed consent was obtained by the participants.
Fig. 1

Age of the cervical cancer patients during external‐beam radiation therapy.

Age of the cervical cancer patients during external‐beam radiation therapy.

3D‐CRT and VMAT planning

For each patient, treatment plans were generated with 3D‐CRT and VMAT for delivery on a medical linear accelerator (Infinity, Elekta AB, Stockholm, Sweden) emitting 6 and 10 MV photons. The PTV prescribed dose was 45 Gy given in 25 daily fractions of 1.8 Gy. The 3D‐CRT plans involved 10 MV photon beams. The beam arrangement consisted of an anteroposterior, a posteroantrior and two lateral opposed treatment fields at the following gantry angles: of 0°, 90°, 180°, and 270°. The beam weight was optimized for each separate plan in order to achieve a homogeneous dose distribution within the PTV and reduce the radiation dose to surrounding structures. The planning aim was to cover at least 95% of the PTV with the 95% of the prescribed dose. The VMAT plans were created with 6 MV photons. Two coplanar full arcs rotating in clockwise and anticlockwise directions were used. The segment width was set to 1.0 cm and the increment of gantry was 30°. All VMAT plans were designed to satisfy the dose constraints introduced by the RTOG 0418 for modulated radiotherapy. These constraints are summarized in Table 1. Furthermore, the maximum PTV dose was kept below 49.5 Gy whereas less than 1% of the target volume received a dose smaller than 93% of the prescribed dose. No more than 1% or 1 cm3 (whichever is lower) of the healthy tissues outside PTV absorbed a radiation dose exceeding 110% of the tumor dose.
Table 1

Dosimetric constraints of the planning target volume (PTV) and surrounding critical organs applied for VMAT planning.

StructureVi a
PTVV45 ≥97%
BowelV40 <30%
RectumV30 <60%
BladderV45 <35%
Femoral headsV30 ≤15%

The Vi corresponds to the structure volume absorbing a radiation dose of i Gy.

Dosimetric constraints of the planning target volume (PTV) and surrounding critical organs applied for VMAT planning. The Vi corresponds to the structure volume absorbing a radiation dose of i Gy.

Estimation of the patient‐specific cancer risk

Pelvic radiation therapy for gynecological carcinomas may result in second bladder and rectal malignant diseases. , , , The above structures were characterized as partially in‐field organs in accordance with a previous study of Howell et al. The above OARs were adjacent to the target volume and they received an inhomogeneous dose distribution in the 3D‐CRT and VMAT plans of all study participants. Parts of both rectum and urinary bladder absorbed radiation doses up to that of the PTV. Hall reported that the risk of radiation carcinogenesis is linearly related to the absorbed dose for a dose range of 0.1–2.5 Gy. Dasu and Toma‐Dasu also stated that the linear relationship exists up to doses of 1–2 Gy. The extrapolation of the linear‐no‐threshold approach to doses exceeding the above levels is not recommended. There are several non‐linear models in the literature which ignore either the fractionation or the cell proliferation effects occurring during radiotherapy. Schneider et al. introduced a mechanistic model providing site‐specific dose response relations using data obtained by A‐bomb survivors and irradiated Hodgkin’s disease patients. This mechanistic model is based on a mixture of plateau and bell‐shaped dose response relationships. The model‐based cancer risk estimates account for the tumor dose fractionation and the interfraction repair of the exposed tissues. The above non‐linear mechanistic model was applied to estimate the cancer risk to bladder and rectum from radiation therapy for cervical cancer. The probability of carcinogenesis to critical organs has been previously assessed with this model for patients irradiated for carcinomas , , and benign disorders. The application of the mechanistic model required the knowledge of the organ equivalent dose (OED) of each OAR. Differential dose‐volume histograms were used to calculate the OED of the bladder or rectum from the 3D‐CRT and VMAT plans of each patient as follows: where is the total organ volume, is the organ volume receiving a dose equal to , is a factor associated with the organ‐dependent repopulation and is the cell killing factor. The was calculated using the formula: where and are the linear quatradic parameters, is the PTV dose of 1.8 Gy delivered in each fraction and is the total target dose of 45 Gy. The patient‐dependent excess absolute risk (EAR) for the appearance of radiation‐induced bladder or rectal malignancies due to 3D‐CRT or VMAT for cervical cancer was estimated as follows: where is the initial slope at low doses of the dose‐response curve for a Western population, and are the age modifying parameters for each organ of interest, is the age of the female patient during radiation therapy and is the attained age. The lifetime attributable risk (LAR) for second cancer induction was estimated with the formula: where the was taken equal to 80 years, is a cancer‐risk free period of 5 years for the development of radiation‐induced solid tumors and is the probability of a healthy female to survive from to obtained by the United States life tables. The parameters, , , , and for the bladder were taken from the literature and they were 0.06, 0.219 Gy‐1, −0.024, 2.38, and 3.8/(104 PY Gy), respectively. The corresponding parameters for the rectum were 0.56, 0.033 Gy‐1, −0.056, 6.9 and 0.73/(104 PY Gy). , The LAR estimates were expressed as the risk per 10000 people.

Estimation of the average cancer risk

The average OEDav of bladder and rectum was found from the 3D‐CRT and VMAT plans of all patients. The calculated OEDav was employed to estimate the respective average lifetime risk (LARav) for the appearance of radiotherapy‐induced bladder and rectal malignancies using the equations of the previous subsection. The LARav was estimated for a typical 50‐year‐old patient at the time of irradiation and an attained age of 80 yrs. The LARav was combined with the baseline risk (BR) to estimate the average relative risk (RRav) of carcinogenesis as follows: Based on SEER data, the BR of a 50‐year‐old healthy female to be diagnosed with bladder and rectal malignancies in 30 yrs is 0.68% and 2.33%, respectively.

Statistics

The patient‐specific second cancer risk estimates related to 3D‐CRT for cervical carcinoma were compared with the assessments derived from the VMAT plans by using a Wilcoxon signed‐rank test. Differences at the level of P‐value of 0.05 or less were considered as significant. The Bland‐Altman statistical test was also applied to determine the percentage mean difference (MD) between the organ‐specific cancer risks estimated by the two treatment delivery techniques. The 95% confidence intervals between the risks derived from the two techniques were set to MD±1.96 SD, where SD is the standard deviation of the differences. Statistical analysis was performed using the GraphPad Prism software (Graph Pad Software Inc., CA, USA).

RESULTS

Patient‐specific cancer risk estimates

The OED calculations are presented in Figs. 2 and 3. The OED range of the bladder from 3D‐CRT and VMAT was 17.1–18.1 cGy and 18.3–20.1 cGy, respectively (Fig. 2). The corresponding range for rectum was 831.0–980.7 cGy and 813.9–977.9 cGy (Fig. 3). The lifetime risks for the development of radiotherapy‐induced malignancies are presented in Table 2. The LAR range for bladder cancer induction associated with 3D‐CRT and VMAT for cervical cancer was found to be 2.3 × 10‐4 to 13.2 × 10‐4 and 2.5 × 10‐4 to 13.7 × 10‐4, respectively. The corresponding LARs for rectal malignancies were 10.3 × 10‐4 to 82.8 × 10‐4 and 9.9 × 10‐4 to 78.3 × 10‐4.
Fig. 2

Organ equivalent dose calculations for bladder derived from the 3D‐CRT and VMAT plans of the cervical cancer patients.

Fig. 3

Organ equivalent dose calculations for rectum derived from the 3D‐CRT and VMAT plans of the cervical cancer patients.

Table 2

Lifetime attributable risk (LAR) estimates per 10000 people for the development of second bladder and rectal malignancies derived from 3D‐CRT and VMAT plans of patients with cervical cancer.

PatientLAR for bladder cancer (×10‐4)LAR for rectal cancer (×10‐4)
3D‐CRTVMAT3D‐CRTVMAT
14.14.220.319.6
211.812.263.763.9
36.37.127.927.1
412.613.069.369.9
52.32.510.39.9
62.62.811.411.3
713.213.782.878.3
810.612.163.559.9
93.94.117.517.1
109.410.049.450.3
114.24.319.418.7
127.17.733.833.2
1310.511.156.655.6
148.38.745.245.1
158.79.849.449.3
164.95.223.022.4
179.09.849.247.7
Organ equivalent dose calculations for bladder derived from the 3D‐CRT and VMAT plans of the cervical cancer patients. Organ equivalent dose calculations for rectum derived from the 3D‐CRT and VMAT plans of the cervical cancer patients. Lifetime attributable risk (LAR) estimates per 10000 people for the development of second bladder and rectal malignancies derived from 3D‐CRT and VMAT plans of patients with cervical cancer. The LARs for second cancer induction from 3D‐CRT were significantly different from those related to VMAT for cervical cancer (Bladder: P = 0.0003; Rectum: P = 0.0144). The use of 3D‐CRT resulted in lower bladder cancer risks than VMAT for all patients examined (Table 2). Conventional treatment led to an increased rectal cancer risk compared to that from VMAT in 14 of 17 study participants (Table 2). The Bland‐Altman scatter plots are shown in Fig. 4. The MD between the bladder cancer risk estimates associated with 3D‐CRT and VMAT was −6.6 ± 3.4% with 95% confidence intervals of 0.1% to −13.3%. The corresponding MD in the assessment of the probability for rectal cancer development was 2.0 ± 2.1%. The 95% limits of agreement were equal to −2.1% and 6.1%.
Fig. 4

Scatter plots presenting the differences between the patient‐specific risks for developing bladder cancer (a), and rectal cancer (b) as estimated by the 3D‐CRT and VMAT plans of patients with primary carcinoma of the uterine cervix against the average risk value. The dotted lines illustrate the 95% confidence intervals and the solid line is the mean% risk difference.

Scatter plots presenting the differences between the patient‐specific risks for developing bladder cancer (a), and rectal cancer (b) as estimated by the 3D‐CRT and VMAT plans of patients with primary carcinoma of the uterine cervix against the average risk value. The dotted lines illustrate the 95% confidence intervals and the solid line is the mean% risk difference.

Average cancer risk estimates

The OEDav of bladder attributable to 3D‐CRT and VMAT for cervical carcinoma was 17.7 and 18.9 cGy, respectively. The corresponding quantity for the rectum was 909.2 and 890.4 cGy. The average risks are presented in Table 3. The LARav for the induction of second malignancies ranged from 8.6 × 10‐4 to 44.6 × 10‐4 by the organ of interest and the treatment delivery technique. The RRav related to bladder and rectal cancer development was up to 1.135 and 1.191, respectively.
Table 3

Average lifetime attributable risk (LARav) and average relative risk (RRav) for the development of second bladder and rectal malignancies following 3D‐CRT and VMAT of a 50‐year‐old cervical cancer patient.

Organ‐at‐riskTechniqueLARav (×10‐4)RRav
Bladder3D‐CRT8.61.126
VMAT9.21.135
Rectum3D‐CRT44.61.191
VMAT43.71.187
Average lifetime attributable risk (LARav) and average relative risk (RRav) for the development of second bladder and rectal malignancies following 3D‐CRT and VMAT of a 50‐year‐old cervical cancer patient.

DISCUSSION

The unavoidable exposure of healthy tissues to ionizing radiation during external‐beam radiotherapy of primary neoplasms may elevate the risk for subsequent carcinogenesis. The appearance of second tumors at distant sites from the primarily irradiated area is relatively rare. Dorr and Hermann found that almost 50% of the second tumors appear in the margin area of the treatment volume and less than 10% within the radiation fields. Welte et al. showed that 69% of the total number of second malignancies are presented at the field margins or inside the primarily irradiated region. The current study was focused on the assessment of the second cancer risk to the partially in‐field bladder and rectum due to radiation therapy for uterine cervix carcinoma. These organs, which are characterized by the high susceptibility for radiation carcinogenesis, were partly exposed to primary radiation due to their adjacent location to the treatment volume. The probabilities for the appearance of second malignancies were estimated for 3D‐CRT with 10 MV photons and 6 MV VMAT of cervical carcinoma. The above photon energies are usually employed in our department for non‐modulated and modulated treatment of primary pelvic tumors. The neutron contribution to the exposure of critical sites due to irradiation with 10 MV X‐rays was considered as minimal. The second bladder and rectal cancer risks varied considerably by the patient’s age during treatment. Based on the risks shown in Table 2, the LAR for developing second bladder malignancies due to irradiation of a 40‐year‐old female (patient no. 7) was 5.5 times higher than that for patient no. 5 who was 68‐year‐old at the time of treatment. The corresponding ratio related to the second rectal cancer was about 8.0. The use of 3D‐CRT resulted in significantly lower patient‐specific LARs for bladder cancer induction than those from VMAT. The opposite result was found for the rectal cancer risk where VMAT significantly reduced this risk compared to 3D‐CRT. Despite the above statistical significance, the differences between the second cancer risks estimated by the two treatment techniques were small for all patients. Bland‐Altman analysis showed that the MD for the bladder and rectal cancer risk derived from pelvic VMAT and 3D‐CRT is only 6.6% and 2.0%, respectively. The RRav for developing a second malignancy to bladder or rectum after radiation therapy for cervical cancer was assessed for a typical 50‐year‐old female. The RRav for rectal cancer induction due to VMAT was 1.187. This result revealed that the application of the above treatment technique may elevate the lifetime rectal cancer risk by 18.7% in comparison to the probability of unirradiated 50‐year‐old females. The corresponding increase associated with the conventional 3D‐CRT was 19.1%. The use of VMAT and 3D‐CRT was also found to give an increased probability for developing bladder malignancies in respect to unexposed people by 13.5% and 12.6%, respectively. The above baseline cancer risk increases are considerable and they cannot be ignored by physicians who follow‐up the irradiated patients. Limited data have been published about the theoretical risk of carcinogenesis to adjacent heavily exposed organs attributable to radiation therapy for cervical cancer. , The cancer risk to the entire colon was estimated in these reports. Based on our results, the lifetime bladder cancer risks in females treated during the 5th decade of life with 3D‐CRT were 8.7 × 10‐4 to 13.2 × 10‐4. These probabilities are similar to the lifetime risks of (12–14) × 10‐4 for 40‐year‐old patients undergoing four‐field box irradiation of gynecological carcinomas with 6 MV photons. Zwaheln et al. showed that the 6MV‐IMRT leads to a bladder cancer risk increase of 11.4% compared to 3D‐CRT when a linear‐exponential model is used, whereas, the data analysis with a plateau model gave a decrease of 0.4%. Their results are comparable with the low differences presented in our study based on the application of a non‐linear mechanistic model. We found that VMAT increases the probability for developing bladder malignancies in respect to 3D‐CRT by 2.3%–13.2% in 95% of the irradiated patients. The presented OEDs of bladder and rectum and the relevant cancer risks might contain uncertainties arising from the analysis of the DVHs derived from a commercial treatment planning software. Parts of the above organs‐at‐risk were located close but outside the treatment volume. Inaccuracies may exist in the out‐of‐field dose calculations generated by treatment planning systems. The results of this work may be also limited by errors in the definition of the organ‐specific parameters of the mechanistic model applied for cancer risk estimation. The quantity was taken equal to 3 Gy for all tissues under investigation in accordance with the recommendation of Schneider et al. Reported experience has shown an insignificant change of the breast cancer risk with values from 1 to 5 Gy. It has to be mentioned that the estimated probabilities for radiation‐induced bladder and rectal malignancies presented here solely referred to therapeutic doses. The use of image guidance procedures in clinical practice may result in small absorbed doses to these critical organs compared to those from radiation therapy. Further research is required to evaluate whether these imaging doses have an impact on the second cancer risk magnitude.

CONCLUSIONS

The use of VMAT for cervical cancer significantly reduced the probability for developing second rectal malignancies than 3D‐CRT. The bladder cancer risk from VMAT was significantly high compared to that from the conventional irradiation. However, the absolute differences in the patient‐specific probabilities for the appearance of second malignancies due to VMAT and 3D‐CRT were found to be small. Both delivery techniques led to noticeable elevated second cancer risks compared to those for unirradiated population.

CONFLICT OF INTEREST

The authors declare that there is no conflict of interest regarding the publication of this article.
  31 in total

1.  Out-of-field photon and neutron dose equivalents from step-and-shoot intensity-modulated radiation therapy.

Authors:  Stephen F Kry; Mohammad Salehpour; David S Followill; Marilyn Stovall; Deborah A Kuban; R Allen White; Isaac I Rosen
Journal:  Int J Radiat Oncol Biol Phys       Date:  2005-07-15       Impact factor: 7.038

2.  Statistical methods for assessing agreement between two methods of clinical measurement.

Authors:  J M Bland; D G Altman
Journal:  Lancet       Date:  1986-02-08       Impact factor: 79.321

3.  United States Life Tables, 2017.

Authors:  Elizabeth Arias
Journal:  Natl Vital Stat Rep       Date:  2019-06

Review 4.  Models for the risk of secondary cancers from radiation therapy.

Authors:  Alexandru Dasu; Iuliana Toma-Dasu
Journal:  Phys Med       Date:  2017-02-24       Impact factor: 2.685

5.  A phase II study of intensity modulated radiation therapy to the pelvis for postoperative patients with endometrial carcinoma: radiation therapy oncology group trial 0418.

Authors:  Anuja Jhingran; Kathryn Winter; Lorraine Portelance; Brigitte Miller; Mohammad Salehpour; Rakesh Gaur; Luis Souhami; William Small; Lawrence Berk; David Gaffney
Journal:  Int J Radiat Oncol Biol Phys       Date:  2012-04-28       Impact factor: 7.038

6.  AAPM TG 158: Measurement and calculation of doses outside the treated volume from external-beam radiation therapy.

Authors:  Stephen F Kry; Bryan Bednarz; Rebecca M Howell; Larry Dauer; David Followill; Eric Klein; Harald Paganetti; Brian Wang; Cheng-Shie Wuu; X George Xu
Journal:  Med Phys       Date:  2017-08-20       Impact factor: 4.071

7.  The European Society of Gynaecological Oncology/European Society for Radiotherapy and Oncology/European Society of Pathology guidelines for the management of patients with cervical cancer.

Authors:  David Cibula; Richard Pötter; François Planchamp; Elisabeth Avall-Lundqvist; Daniela Fischerova; Christine Haie Meder; Christhardt Köhler; Fabio Landoni; Sigurd Lax; Jacob Christian Lindegaard; Umesh Mahantshetty; Patrice Mathevet; W Glenn McCluggage; Mary McCormack; Raj Naik; Remi Nout; Sandro Pignata; Jordi Ponce; Denis Querleu; Francesco Raspagliesi; Alexandros Rodolakis; Karl Tamussino; Pauline Wimberger; Maria Rosaria Raspollini
Journal:  Radiother Oncol       Date:  2018-05-01       Impact factor: 6.280

8.  Effect of intensity-modulated pelvic radiotherapy on second cancer risk in the postoperative treatment of endometrial and cervical cancer.

Authors:  Daniel R Zwahlen; Jeremy D Ruben; Phillip Jones; Frank Gagliardi; Jeremy L Millar; Uwe Schneider
Journal:  Int J Radiat Oncol Biol Phys       Date:  2009-06-01       Impact factor: 7.038

9.  Radiation dose and cancer risk to out-of-field and partially in-field organs from radiotherapy for symptomatic vertebral hemangiomas.

Authors:  Michalis Mazonakis; Antonis Tzedakis; Efrossyni Lyraraki; John Damilakis
Journal:  Med Phys       Date:  2016-04       Impact factor: 4.071

10.  Site-specific dose-response relationships for cancer induction from the combined Japanese A-bomb and Hodgkin cohorts for doses relevant to radiotherapy.

Authors:  Uwe Schneider; Marcin Sumila; Judith Robotka
Journal:  Theor Biol Med Model       Date:  2011-07-26       Impact factor: 2.432

View more

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