Literature DB >> 34104089

A clinical study on plasma biomarkers for deciding the use of adjuvant corticosteroid therapy in bronchopulmonary dysplasia of premature infants.

Haiyan Zhu1, Yian Tian2, Huaiping Cheng1, Yafei Zheng1, Wei Wang1, Tianping Bao1, Rong Wu3, Zhaofang Tian1.   

Abstract

Objective: The study was designed to investigate some plasma markers which help us to decide the use of adjuvant corticosteroid therapy in bronchopulmonary dysplasia (BPD) of premature infants.
Methods: Thirty BPD infants were treated by dexamethasone. Among these cases, dexamethasone was significant effective in 10 cases, and no significant effective in 20 cases. These patients were divided into two groups as the significant effect (SE) group (n=10) and the non-significant effect (NE) group (n=20) according to the curative effect of dexamethasone. Fifteen non-BPD infants with gestational age and gender matching were selected as the control group. Plasma samples before and after dexamethasone treatment were collected from three infants chosen randomly from SEG for the data-independent acquisition (DIA) analysis. ELISA was further used to detect the levels of differential proteins LRP1 and S100A8 in all individuals, including SE, NE and control groups.
Results: DIA analysis results showed that after dexamethasone treatment, there were a total of 52 plasma proteins that showed significant differences, of which 43 proteins were down-regulated and 9 proteins were up-regulated. LRP1 and S100A8 were two plasma proteins that were significantly changed after dexamethasone treatment. Compared with the control group, plasma LRP1 was significantly increased in BPD. Interestingly, the plasma concentration of LRP1 in the NE group was significantly higher than that in the SE group. S100A8, as an indicator of plasma inflammation, was significantly higher in BPD than the control group. Unlike LRP1, there was no significantly difference between the SE and NE group (P=0.279) before dexamethasone treatment.
Conclusion: Elevated plasma LRP1 and S100A8 in BPD infants are two indicators that correlated with the efficacy of dexamethasone, and might be used as biomarkers for deciding the use of adjuvant corticosteroids therapy in the BPD. © The author(s).

Entities:  

Keywords:  biomarker; bronchopulmonary dysplasia; corticosteroids; proteome

Mesh:

Substances:

Year:  2021        PMID: 34104089      PMCID: PMC8176188          DOI: 10.7150/ijms.58650

Source DB:  PubMed          Journal:  Int J Med Sci        ISSN: 1449-1907            Impact factor:   3.738


Introduction

Bronchopulmonary dysplasia (BPD) is the result of a complex process in which several prenatal and/or postnatal factors interfere with lower respiratory tract development, leading to a severe, lifelong disease. In the underlined disease the lungs and bronchi are damaged which destroys tissues in the alveoli of the lungs. Although its definition, epidemiology, pathophysiology, prevention, and management have continued to evolve. The prematurity of newborns is the major risk factor for BPD. However, the severity of the BPD differs from infant to infant 1. Each year, over 10000 new BPD cases are reported in the US 2. In China, the incidence of BPD in ultra-premature or ultra-low birth weight infants (ELBW) varied in 12.5%-48.1% 3,4. Many non-pharmacologic and pharmacologic therapies have been suggested for lowering the lung's injury and tend to enhance the recovery of patients 5 i.e., protective ventilation strategies, surfactant supplementation, optimal oxygen saturation goals, and the antenatal corticosteroids 6. The presently used drug therapies for BPD treatment are as following, caffeine 7, bronchodilator 8, surfactant 9, diuretics 10, inositol 11, corticosteroids12, viral immunization 13 and cardiac medications 14 etc. However, these treatments are partially effective; the majority of the interventions tried so far have not proven to be beneficial in rigorous meta-analyses of eligible studies. And still, there is a lack of significant therapy for curing the underlined disease. In BPD the corticosteroids are used to reduce the inflammatory processes that lead to the pathogenesis of the disease. Corticosteroids have a highly potent inhibitory activity against inflammation and have an important contribution to the prevention of BPD. With the risks for Perforated Bowel, hypertrophic cardiomyopathy, major neurosensory disability, and cerebral palsy (CP), the “early” dexamethasone therapy has not been recommended, and the “later” dexamethasone therapy can be used selectively 15. The statement policy of an American Academy of Pediatrics was published in 2010 regarding the use of postnatal corticosteroids and admitted that the available data were not enough to support the corticosteroids recommendations 16. Since the BPD infants have different therapeutic response to dexamethasone therapy, so it has an important clinical value to find biomarkers that can evaluate the effectiveness of the BPD affected infants for dexamethasone therapy. DIA is a non-standard quantitative technique employed for the identification of differentially expressed proteins. This technique is well suited for the evaluation of lung injury in children to identify key biomarkers 17. The purpose of the study is to find differentially expressed plasma proteins of BPD infants responsive to dexamethasone therapy by DIA technique to promote effective treatment of BPD infants.

Materials and Methods

Patients

The existing work was conducted in the Neonatal Intensive Care Unit of The Affiliated Huaian No.1 People's Hospital of Nanjing Medical University. At 36 weeks of postmenstrual age from Jan 2019 to Mar 2020, premature babies with gestational period < 32 weeks were detected with BPD. The diagnosed was based on “New consensus on BPD definitions and diagnostic criteria, developed in 2018 NICHD Seminar” (abbreviated as 2018 NICHD) 18. The following inclusion criterion is applicable in the subjects: clinical diagnosis as BPD grade Ⅲ, the details were as follows. Treatment with more than 21% oxygen for minimum of 28 days, required for 30% oxygen and/or positive pressure or discharge from the hospital, whichever comes first. Any one of the following exclusion criteria was applicable in the current study: I) There existed some congenital malformations such as digestive tract, respiratory tract, central nervous system, and heart; II) There existed ventilator-associated pneumonia (VAP); III) any glucocorticoid e.g., dexamethasone, hydrocortisone, budesonide, etc. administrated before the research. IV) Clinical signs of infection were existed. Totally, thirty BPD infants were enrolled in this study. All eligible cases were treated with a 10-day dexamethasone course after diagnosing with BPD. The details are as follows: Infants received a 10-day tapering course of dexamethasone sodium phosphate (0.15 mg/kg per day for 3 days, 0.10 mg/kg per day for 3 days, 0.05 mg/kg per day for 2days, and 0.02 mg/kg per day for 2 days; total of 0.89 mg/kg over 10 days 19. Plasma samples were retained before and after dexamethasone therapy, respectively. The curative effect was judged immediately after 10-day course of dexamethasone. According to the curative effect of dexamethasone, these BPD infants were divided into two groups as the significant effect (SE) group (n = 10) and the non-significant effect (NE) group (n = 20). SE was defined as patients with no need of respiratory support and oxygen dependence. NE was defined as lowering the inhaled oxygen concentration and/or the decrease in positive pressure support, or no change in inhaled oxygen concentration and/or no improvement in respiratory support. In addition, 15 non-BPD infants with gestational age and gender matching were selected as the control group. The approval for this clinical research was provided by the ethics committee of The Affiliated Huaian No.1 People's Hospital of Nanjing Medical University, and the relevant informed consent was signed from the parents before the trial. The registration number of the Chinese Clinical Trail Registry is ChiCTR2000033796 (http://www.chictr.org.cn/historyversionpub.aspx?regno=ChiCTR200003796).

DA analysis

Plasma samples before and after dexamethasone treatment were collected from three infants chosen randomly from SEG for the data-independent acquisition (DIA) analysis (Guangzhou Jidio Biotechnology Co. Ltd, Guangzou, China). Plasma samples were transferred into lysis buffer, and the proteins in lysates were concentrated, followed by proteins digesting via sequence-grade modified trypsin using the method of BCA Protein Assay Kit (Guangzhou Jidio Biotechnology Co. Ltd, Guangzhou, China). The digested peptides were dried to prepare for the analysis of nano-HPLC-MS/MS. Unprocessed data obtained from DIA were analyzed and sort out via Spectronaut X (Biognosys AG, Switzerland) with default parameters. After Student's t-test, various expressed proteins (having Q value < 0.05 and Absolute AVG log2 ratio > 0.58) were filtered.

Enzyme-linked immunosorbet assay (ELISA)

To confirm the differentially expressed proteins identified by DIA, we focused on LDL receptor-related protein 1 (LRP1) and S100 calcium-binding protein A8 (S100A8). Venous blood was centrifuged at 1,000 × g to obtain a supernatant. The levels of LRP1 and S100A8 in the supernatant were identified via the ELISA kit (abcam, Cambridge, England) for the control group and the BPD group (including pre-and post- dexamethasone administration), following the given instructions of the manufacturer.

Statistical analysis

GraphPad Prism 8.0 (GraphPad Software Inc.) was used for statistical processing. The difference was evaluated by paired Student t test (2-group comparisons) or ANOVA followed by the post hoc Tukey's Multiple Comparison Test (multi-group comparisons) as appropriate. Data are presented as the mean ± SD. P < 0.05 was considered statistically significant.

Results

The characteristics of BPD infants

A total of 30 eligible infants were included in this study, of which 10 cases showed significant effect (SE), while 20 non-significant effects (NE) after dexamethasone treatment. The two groups' patients had similar perinatal characteristics, as depicted in Table 1. The difference was that three patients in the NE group had ventilator associated pneumonia, while none in the SE group. Three BPD patients were randomly selected from the SE group. The perinatal characteristics were presented in Table 2. The perinatal characteristics of the selected three BPD patients for DIA analyses were representative.
Table 1

Sociodemographic and prerandomization characteristics of study infants [mean ±SD or n (%)]

Perinatal characteristicsSE group (n=10)NE group (n=20)X2/tP value
Gestational age (weeks)31.2 ± 1.131.5 ± 1.40.0920.761
Birthweight (grams)1325.0 ± 249.01365.0 ± 262.00.0530.816
Male sex [n (%)]6 (60)13 (65)0.6010.421
Maternal dexamethasone [n (%)]9 (90)18 (90)>0.05
Vaginal delivery [n (%)]6 (60)11 (55)0.0930.758
Apgar score at 5 min7 ± 27 ± 2>0.05
Surfactant doses administered [n (%)]10 (100)20 (100)>0.05
BPD grade (III) [n (%)]10 (100)20 (100)>0.05
Respiratory distress syndrome [n (%)]10 (100)20 (100)>0.05
Pulmonary hypertension (PH) [n (%)]2 (20)5 (25)0.1440.739
Patent ductus arteriosus (PDA) [n (%)]8 (80)17 (85)0.6290.416
Ventilator associated pneumonia [n (%)]0 (0)3 (15)<0.01
Atrial septal defect (ASD) [n (%)]8 (80)18 (90)0.8380.341
Ventricle septal defect (VSD) [n (%)]8 (80)17 (85)1.8760.170
Table 2

Sociodemographic and prerandomization characteristics of dexamethasone sensitive group infants (n = 10)

Perinatal characteristicsDIA groupNo-DIA group
n37
Gestational age, weeks30.0 ± 1.030.6 ± 1.5
Birthweight, grams1350 ± 1601365 ± 180
Male sex, n (%)2 (66.7)4 (57.1)
Maternal dexamethasone, n (%)3 (100.0)6 (85.7)
Vaginal delivery, n (%)2 (66.7)4 (57.1)
Apgar score at 5 min7 ± 27 ± 2
Respiratory distress syndrome, n (%)3 (100.0)7 (100.0)
Surfactant doses administered, n (%)3 (100.0)7 (100.0)
BPD grade (III), n (%)3 (100.0)7 (100.0)

There were no significant differences between groups. Data are mean ± SD or n (%).

Differentially expressed proteins before and after dexamethasone administration in SE group

There were 52 proteins displayed significant expression changes before and after dexamethasone treatment. In this cohort, the upregulation of 9 proteins and downregulation of 43 proteins were found, as indicated in Table 3. The volcanic map directly showed the location of various proteins between the groups that were compared, as depicted in Fig. 1. Hierarchical cluster analysis and Z-score homogenization were used to map the differentially expressed proteins in groups that were compared (each column represents a sample and each row represents a protein). The protein expression in different samples is expressed in different colors. The red color represents the elevated expression, while the blue color showed decreased expression level, as depicted in Fig. 2.
Table 3

The information table of significant difference protein between after and before treatment of corticosteroids (top 10 proteins)

idSymbollog2fcP value
ENSG00000123384LRP1-2.7641527050.006
ENSG00000127083OMD-1.7132228680.013
ENSG00000129538RNASE1-1.756946770.019
ENSG00000143546S100A8-2.0125649870.043
ENSG00000011052NME-13.437609290.0459
ENSG00000167815PRDX2-1.756656930.165
ENSG00000090013BLVRB-2.1065225020.183
ENSG00000285043AC093512.2-2.216570550.211
ENSG00000111341MGP-2.7878951940.222
ENSG00000041982TNC-2.697200620.366
Figure 1

Distinct pulmonary proteomic patterns in T1 group and T2 group. Red represents up-regulated protein, blue represents down regulated protein, and black has no difference (T1 Pre-dexamethasone administration; T2 Post-dexamethasone administration).

Figure 2

Heat map of differentially expressed proteins (T1 Pre- dexamethasone administration; T2 Pos-t dexamethasone administration).

Bioinformatics analyses of the differentially expressed proteins

The analysis of GO (The Gene Ontology) indicated that the differentially expressed proteins participated in the cofactor metabolic process (GO : 0051186) and the catabolic process (GO : 0009056). The most significantly affected molecular functions were catalytic activity, binding, and structural activity relationship (Supplementary Online Fig. S1). Pathway analysis by KEGG suggested that the most significantly affected cascades were purine metabolism, pyrimidine metabolism, drug metabolism-other enzymes (Supplementary Online Fig. S2). We further map genes to each term of the DO database (http://disease-ontology.org/). The results of DO enrichment were shown in Supplementary Online Fig. S3. The DO items that pulmonary hypertension and persistent fetal circulation syndrome might relate to the pathogenesis of BPD. Then, we map the gene to each term of the react database (https://react.org/). The underlined results were depicted in Supplementary Online Fig. S4, and the top twelve Reactome items were significantly enriched in the gene (p < 0.05).

Validation of changed LRP1 and S100A8levels in the plasma samples

ELISA analysis showed that LRP1 were significantly different among the three groups before dexamethasone administration (F = 30.42, P < 0.001). Compared with the control group, plasma LRP1 was significantly increased in BPD. Interestingly, the plasma concentration of LRP1 in the NE group was significantly higher than that in the SE group (Fig. 3A). After dexamethasone administration than, the level of LRP1 were differentially downregulated either in the SE group or in the NE group (Fig. 3B). Likewise, S100A8 was significantly different among the three groups before dexamethasone administration (F = 58.61, P < 0.001). The level of S100A8 was significantly higher in both the SE and NE group than that in the control group (P < 0.05). However, there was no significantly difference between the SE and NE group (P = 0.279) before dexamethasone treatment (Fig. 3C). The level of S100A8 was differentially downregulated in both the SE group and NE group after dexamethasone administration than that before dexamethasone administration. (P < 0.05) (Fig. 3D).
Figure 3

The levels of LRP1 and S100A8 in the plasma of different groups. *P<0.05, ***P<0.001 (SE-the significant effect group, NE-the non-significant effect group, control-non-BPD group).

Discussion

In this study, we used DIA technology to analyze the plasma proteins of BPD infants before and after dexamethasone treatment. We found that both LRP1 and S00A8 decreased significantly after dexamethasone treatment. If the baseline LRP1 is high, it suggests that dexamethasone may not show significant effectiveness. This provides a basis for clinical selection of adjuvant corticosteroids to treat BPD. Dexamethasone had been widely used since the 1980s because of its potent activity against inflammation. The underlined corticosteroid can significantly change the course of BPD and decrease the BPD severity 20. However, the adverse effects (such as cerebral palsy and major neurosensory disability, etc.) caused by glucocorticoids were gradually found from late 1990 21.Therefore, in 2002, the American Academy of Pediatrics had suggested that the use of glucocorticoids in BPD should be more cautious 22. Then the therapeutic use of corticosteroids in BPD had reduced effectively, but the rate of BPD had increased 23. In recent years, glucocorticoid use in BPD infants has been increased. A European study revealed that 13.9% of premature infants (aged less than 30 weeks) received corticosteroid treatment after birth 24, however, in Japan the use of systemic steroid was published by 56% of units in 2005, 57% in 2010, and 68% of units in 2015, respectively 25. There is still no consensus on the use of corticosteroids in BPD 26, 27, 28. In our study, 10 BPD infants were found to be significantly positive response to dexamethasone treatment, and the efficacy of the other 20 patients is not very satisfactory. The results suggested that some sensitive biomarkers are needed to choose the use of adjuvant corticosteroid therapy in BPD. Protein mass spectrometry had been applied to some experimental studies regarding BPD 29, 30, and clinical research by DIA indicated that calcium-associated proteins can be used as sensitive markers for BPD progress 31. In this study, the results revealed that a total of 52 proteins were differentially expressed before and after dexamethasone treatment. In these differentially expressed proteins, the decreased expression of 43 proteins and elevated expression of 9 proteins were determined. The underlined differentially expressed proteins were further evaluated by literature search and the LRP1 and S100A8 were found to be closely correlated with lung diseases and/or glucocorticoids. In this study, the plasma LRP1 level in the two BPD subgroups was significantly increased than that in the control group before dexamethasone administration. We also found that, the levels of LRP1 in the SE group were lower than those in the NE group before dexamethasone treatment. These results suggested that LRP1 might be either a diagnostic index or a target of dexamethasone therapy for BPD. These results suggest that a high baseline level of LRP1 may indicate that dexamethasone is not significant effective. This will provide a basis for clinical selection of dexamethasone to treat BPD. LRP1 was first discovered in 1988, and its amino acid sequence was found to be similar to the LDL receptor. LRP1 is a cell surface receptor and regulates cellular and molecular mechanisms that drive the pathophysiological inflammatory reactions and reorganization of tissues in many organs 32. In current decades, LRP1 has been evolved as a significant regulator of the inflammatory reactions. LRP1 negatively regulates adaptive immune responses in the site of HDM-induced eosinophilic airway inflammation 33. LRP has a close association with dexamethasone. The expression of LRP was up-regulated by dexamethasone, and down-regulated by lipopolysaccharide (LPS), gamma interferon (IFN-gamma) or a combination of both. LRP was less sensitive to dexamethasone in activated astrocytes than in microglia 34. Dexamethasone increases cell-surface LRP activity in HepG2 cells by increasing the steady state mRNA levels and suggest that post-transcriptional mechanisms play a role in controlling LRP mRNA levels 35. Calprotectin S100A8 (S100 calcium binding protein A8, S100A8) is an important member of calcium binding protein family. The evaluation of its role in the pathogenesis of respiratory diseases or its usefulness as a biomarker for the appropriate diagnosis and prognosis of lung diseases have only gained attention in recent years 36. S100A8/A9 is an essential factor for neutrophil recruitment to lungs, and S100A8 promotes acute lung injury via Toll-like receptor 4-dependent activation of AECs 37; S100A8/A9 serum levels along with chemokines are useful in distinguishing between active tuberculosis and asymptomatic Mycobacterium tuberculosis-infected latent individuals, and targeting S100A8/A9 pathways as host-directed therapy for TB 38; Calprotectin is significantly elevated in the sputum of Bronchiolitis obliterans patients and reflects ongoing neutrophilic inflammation 39; serum S100A8/S100A9 was negatively associated with pulmonary function in acute exacerbation COPD(AE-COPD) patients, indicating that the serum S100A8/S100A9 heterodimer may be involved in the progression of AE-COPD, and may be a relevant serum biomarker in the diagnosis for AE-COPD 40. Elevated S100A8/S100A9 expression causes glucocorticoid resistance in MLL-rearranged infant acute lymphoblastic leukemia 41. We found that S100A8 levels were significantly different between BPD group and the control group. Unlike LRP1, there was no difference in the level of S100A8 between the SE and NE groups. In conclusion, we determined some differentially expressed plasma proteins of the dexamethasone significant effect infants with BPD via the DIA technique, among which, LRP1 and S100A8 might be the biomarkers for BPD diagnosis and deciding the use of adjuvant corticosteroids therapy in the BPD infants. This study is only a single center study and the sample number of is still small. In order to support the conclusions of this study, a large sample and multicenter study is necessary. Supplementary figures. Click here for additional data file.
  41 in total

1.  Elevated S100A8/S100A9 expression causes glucocorticoid resistance in MLL-rearranged infant acute lymphoblastic leukemia.

Authors:  J A P Spijkers-Hagelstein; P Schneider; E Hulleman; J de Boer; O Williams; R Pieters; R W Stam
Journal:  Leukemia       Date:  2012-01-27       Impact factor: 11.528

2.  Describing the Increase in Preterm Births in the United States, 2014-2016.

Authors:  Joyce A Martin; Michelle J K Osterman
Journal:  NCHS Data Brief       Date:  2018-06

3.  Management of bronchopulmonary dysplasia in Japan: A 10-year nationwide survey.

Authors:  Fuyu Miyake; Masato Ito; Hirotaka Minami; Masanori Tamura; Fumihiko Namba
Journal:  Pediatr Neonatol       Date:  2019-11-22       Impact factor: 2.083

4.  S100A8/A9 regulates CD11b expression and neutrophil recruitment during chronic tuberculosis.

Authors:  Ninecia R Scott; Rosemary V Swanson; Noor Al-Hammadi; Racquel Domingo-Gonzalez; Javier Rangel-Moreno; Belinda A Kriel; Allison N Bucsan; Shibali Das; Mushtaq Ahmed; Smriti Mehra; Puthayalai Treerat; Alfredo Cruz-Lagunas; Luis Jimenez-Alvarez; Marcela Muñoz-Torrico; Karen Bobadilla-Lozoya; Thomas Vogl; Gerhard Walzl; Nelita du Plessis; Deepak Kaushal; Thomas J Scriba; Joaquín Zúñiga; Shabaana A Khader
Journal:  J Clin Invest       Date:  2020-06-01       Impact factor: 14.808

Review 5.  Surfactant instillation in spontaneously breathing preterm infants: a systematic review and meta-analysis.

Authors:  Vincent Rigo; Caroline Lefebvre; Isabelle Broux
Journal:  Eur J Pediatr       Date:  2016-09-27       Impact factor: 3.183

6.  Controlled trial of dexamethasone therapy in infants with bronchopulmonary dysplasia.

Authors:  M C Mammel; T P Green; D E Johnson; T R Thompson
Journal:  Lancet       Date:  1983-06-18       Impact factor: 79.321

7.  Trends in Care Practices, Morbidity, and Mortality of Extremely Preterm Neonates, 1993-2012.

Authors:  Barbara J Stoll; Nellie I Hansen; Edward F Bell; Michele C Walsh; Waldemar A Carlo; Seetha Shankaran; Abbot R Laptook; Pablo J Sánchez; Krisa P Van Meurs; Myra Wyckoff; Abhik Das; Ellen C Hale; M Bethany Ball; Nancy S Newman; Kurt Schibler; Brenda B Poindexter; Kathleen A Kennedy; C Michael Cotten; Kristi L Watterberg; Carl T D'Angio; Sara B DeMauro; William E Truog; Uday Devaskar; Rosemary D Higgins
Journal:  JAMA       Date:  2015-09-08       Impact factor: 56.272

8.  Neonatal mortality and morbidity among infants between 24 to 31 complete weeks: a multicenter survey in China from 2013 to 2014.

Authors:  XiangYong Kong; FengDan Xu; Rong Wu; Hui Wu; Rong Ju; XiaoLin Zhao; XiaoMei Tong; HongYan Lv; YanJie Ding; Fang Liu; Ping Xu; WeiPeng Liu; HongBin Cheng; TieQiang Chen; ShuJuan Zeng; WenZheng Jia; ZhanKui Li; HuiXian Qiu; Jin Wang; ZhiChun Feng
Journal:  BMC Pediatr       Date:  2016-11-03       Impact factor: 2.125

9.  Calcium signaling-related proteins are associated with broncho-pulmonary dysplasia progression.

Authors:  Cinzia Magagnotti; Piero Giuseppe Matassa; Angela Bachi; Valentina Vendettuoli; Isabella Fermo; Maria Rosa Colnaghi; Rose Mary Carletti; Domenica Mercadante; Elena Fattore; Fabio Mosca; Annapaola Andolfo
Journal:  J Proteomics       Date:  2013-10-16       Impact factor: 4.044

Review 10.  The proteomics of lung injury in childhood: challenges and opportunities.

Authors:  Prue M Pereira-Fantini; David G Tingay
Journal:  Clin Proteomics       Date:  2016-02-29       Impact factor: 3.988

View more
  1 in total

1.  Role of the LRP1-pPyk2-MMP9 pathway in hyperoxia-induced lung injury in neonatal rats.

Authors:  Ya-Fei Zheng; Hai-Yan Zhu; Wei Wang; Jing-Jing Hu; Tian-Ping Bao; Zhao-Fang Tian
Journal:  Zhongguo Dang Dai Er Ke Za Zhi       Date:  2021-12-15
  1 in total

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