Literature DB >> 28144260

The rs2228145 polymorphism in the interleukin-6 receptor and its association with long-term prognosis after myocardial infarction in a pilot study.

Anna Szpakowicz1, Witold Pepinski2, Ewa Waszkiewicz1, Małgorzata Skawronska2, Anna Niemcunowicz-Janica2, Wlodzimierz J Musial1, Karol A Kaminski1.   

Abstract

INTRODUCTION: Interleukin-6 (IL-6) is a cytokine with a complex function that is described as both pro- and anti-inflammatory. One factor that influences its function is the rs2228145 A/C single nucleotide polymorphism (SNP) of the IL-6 receptor (IL6R) gene. C allele carriers have a decreased inflammatory response and decreased prevalence of ischemic heart disease. The aim of the study was to investigate the association of the rs2228145 SNP of the IL6R gene with long-term total mortality in patients with ST-elevation myocardial infarction (STEMI) treated invasively.
MATERIAL AND METHODS: We analyzed the data of consecutive patients with ST elevation myocardial infarction (STEMI) treated with primary percutaneous coronary intervention (PCI). Genotyping was performed with the TaqMan method. The analyzed end-point was total long-term mortality (median: 2875 days).
RESULTS: The registry comprised 553 patients (mean age: 62.4 ±11.9 years; 25.6% females, n = 142; TIMI 3 obtained in 91.7% of patients, n = 507). No significant differences in baseline characteristics were found between the genotypes. During long-term follow-up 171 (30.9%) patients died. There was non-significantly higher mortality in the rs2228145 AA homozygotes compared to C allele carriers (OR = 1.34, 95% CI: 0.93-1.93, p = 0.1).
CONCLUSIONS: The rs2228145 polymorphism of IL6R was not significantly associated with long-term mortality after STEMI. However, AA homozygotes (high-risk genotype for ischemic heart disease) showed a trend towards adverse outcome compared to C allele carriers. The observed trend is promising, but it requires independent replication studies.

Entities:  

Keywords:  IL-6 receptor; ST elevation myocardial infarction; acute coronary syndrome; interleukin-6

Year:  2016        PMID: 28144260      PMCID: PMC5206359          DOI: 10.5114/aoms.2016.58636

Source DB:  PubMed          Journal:  Arch Med Sci        ISSN: 1734-1922            Impact factor:   3.318


Introduction

Interleukin-6 (IL-6) is a cytokine with a very complex function that is described as both pro- and anti-inflammatory. Therefore, it is not resolved yet whether it is a causative factor for atherosclerosis and heart failure, just a risk marker, or even a compensatory agent (like B-type natriuretic peptide in heart failure). An increased concentration of IL-6 is associated with higher prevalence of acute coronary syndrome [1, 2], adverse subsequent prognosis in this group of patients [3-5], and increased risk for development of heart failure [6, 7]. Interleukin-6 is secreted by myocardium during acute myocardial infarction and following reperfusion [8]. The multidirectional influence of IL-6 on the cardiovascular system can be explained by its very complex interaction with receptors: the IL-6 receptor (IL6R, glycoprotein 80, CD126) and glycoprotein 130 (gp130, IL6STIL-6 signal transducer, CD130). Both of them can exist in a soluble form (sIL6R and sgp130) or as membrane bound receptors. A specific IL-6 binding subunit, IL-6R, is present only on a few cell types (hepatocytes, leukocytes), whereas the subunit common for all IL-6 family cytokines, gp130 (necessary to activate intracellular transduction cascades), is expressed on all types of cells. The signaling pathways of IL-6 are different depending on the primary form of IL-6 (free or bound to soluble receptors), adequate effector membrane receptors, and stimulated cell lines. Another factor that influences IL-6 function is the rs2228145 A/C single nucleotide polymorphism (SNP) in the intron of an IL6R gene. The C allele is associated with increased sIL6R concentrations [9-11], decreased IL6R membrane expression on monocytes and CD4+ lymphocytes, and decreased response to IL-6 (decreased phosphorylation of STAT1 and STAT3) [10]. This phenomenon can potentially be explained by the mechanism of increased cleavage of the membrane receptor to plasma [12]. C allele carriers have a decreased inflammatory response as well as decreased prevalence of atrial fibrillation, diabetes, metabolic syndrome, aortic aneurysm and rheumatoid arthritis [10, 13–15]. In patients with aortic stenosis the C allele was associated with lower transvalvular gradients [16]. A recent meta-analysis of 82 studies confirmed that the rs2228145 SNP is also associated with all forms of coronary heart disease (CHD) in patients of European descent [9]. This finding was not confirmed in Han Chinese with premature CHD (age < 55 years in men or < 65 years in women) [17]. In this analysis the SNP also failed to associate with severity of atherosclerotic lesions or clinical phenotypes of CHD [17]. The study, however, was limited by the small number of patients included (n = 187 cases). The association between the rs2228145 polymorphism and CHD in the European population has already been reliably confirmed. Still, there are no studies focusing on long-term prognosis after myocardial infarction, which is a specific form of CHD. Therefore, the aim of the study was to investigate the association of the rs2228145 SNP of the IL6R gene with 5-year total mortality in patients with ST-elevation myocardial infarction (STEMI) treated invasively.

Material and methods

This was a retrospective observational study performed in a real-life registry of patients with STEMI. As described previously [18-20], the registry included consecutive patients with STEMI admitted to our department in the years 2001–2005 and treated invasively within 12 h from symptoms onset. On the day of admission, after obtaining informed written consent, blood samples were collected. Afterwards, clinical data were retrieved retrospectively from hospital documentation. In this way we constructed a real-world registry, with no exclusion criteria except for lack of consent. In this analysis, however, we focused on long-term outcome and therefore selected only those individuals who survived the first 48 h after hospital admission. ST-elevation myocardial infarction was diagnosed based on a history of typical chest pain, ECG changes (ST-segment elevation or a new left bundle branch block) and a rise in cardiac necrosis markers. Pharmacological treatment was consistent with contemporary guidelines. The database comprised patients’ history, data from physical examination on admission, laboratory tests, echocardiography, and coronary angiography. The GRACE risk score was calculated retrospectively, based on data collected on admission [21]. The control group comprised 51 adult men and 50 adult women, who took part in paternity testing. Detailed clinical characteristics were not available, but such a random group should be highly representative in terms of the genetic background for our region. Blood samples for genotype testing were collected in EDTA tubes and stored at –20°C. Commercial kits were used for DNA extraction (Blood Mini, A&A Biotechnology). Genotypes were determined with TaqMan SNP Genotyping Assay on the ABI 7500 real-time PCR platform (Applied Biosystems). Ten percent of the samples were genotyped in duplicate. Long-term follow-up was performed (maximum: 4167, minimum: 2558 days). The analyzed end-point was all-cause mortality. Survival status was retrieved from the local population registry run by a Government Office. The research was performed in accordance with the ethical standards laid down in the 1964 Declaration of Helsinki. Informed written consent was obtained from all the subjects before their inclusion in the study. The protocol was approved by the Ethics Review Board of the Medical University of Bialystok.

Statistical analysis

Statistical analysis was performed with Statistica 9.0 and ARLEQUIN v.3.0 (Hardy-Weinberg equilibrium) software. The 3 genotypes of the rs2228145 SNP were compared in terms of clinical parameters (χ2 or Kruskal-Wallis ANOVA test, as appropriate) and mortality rates (χ2 test). Based on visual analysis of Kaplan-Meier survival curves, patients with CC and AC genotypes were combined into one group for additional analysis. Probability of survival was compared with the log-rank test. For univariate and multivariate analyses, logistic regression was used. Variables from univariate analysis significantly associated with mortality (except for the GRACE risk score) were entered into a primary model of multivariate analysis and further eliminated in a backward stepwise manner. A two-sided p value < 0.05 was considered statistically significant.

Results

We enrolled in the registry 652 subjects, 9 of whom were lost to follow-up (1.4%). Genotyping was successfully performed in 553 patients, which constituted our study group. No discrepancies were observed in the results of samples genotyped in duplicate. The control group comprised 101 subjects who underwent paternity testing. The percentages of specific genotypes of the rs2228145 SNP and both study and control groups are shown in Table I.
Table I

Percentages of specific genotypes of the rs2228145 polymorphism and associated mortality rates

GenotypeAAACCC
Study group (n = 553):
 Number (%)243 (43.9)242 (43.8)68 (12.3)
 Long-term mortality, n (%)34.6 (84)*28.5 (69)26.5 (18)
Control group (n = 101):
 Number (%)35 (34.65)51 (50.5)15 (14.85)

p = 0.1 compared to C allele carriers

χ2 test. Percentages of specific genotypes in the study group vs. control group: p = 0.21, χ2 test.

Percentages of specific genotypes of the rs2228145 polymorphism and associated mortality rates p = 0.1 compared to C allele carriers χ2 test. Percentages of specific genotypes in the study group vs. control group: p = 0.21, χ2 test. There was a trend towards higher prevalence of the C allele (protective one) in the control group, p = 0.21, χ2 test. Both groups were in Hardy-Weinberg equilibrium (p-values 0.52 and 0.6, consecutively). Mean age in the study group was 62.4 ±11.9 years, 25.6% of the group were female (n = 142), and TIMI 3 after PCI was obtained in 91.7% of patients (n = 507). Detailed clinical characteristics of the study group are presented in Table II. No significant clinical differences were found between the genotypes.
Table II

Baseline characteristics of the study group based on rs2228145 genotype (risk allele for myocardial infarction – A)

CharacteristicOverall population N = 553rs2228145 AA homozygotes N = 243rs2228145 AC heterozygotes N = 242rs2228145 CC homozygotes N = 68 P-value
Age [years]62.4 (11.9)63.1 (11.7)61.8 (12.6)62.4 (10.0)0.44
Female gender (%)25.6 (n = 142)25.1 (n = 61)24.0 (n = 58)33.8 (n = 23)0.24
Hypertension (%)55.1 (n = 305)56.4 (n = 137)52.9 (n = 128)58.8 (n = 40)0.60
Type 2 diabetes (%)22 (n = 122)23.9 (n = 58)22.3 (n = 54)14.7 (n = 10)0.27
Previous myocardial infarction (%)10.1 (n = 56)9.0 (n = 22)11.6 (n = 28)8.8 (n = 6)0.61
Systolic blood pressure [mm Hg]138.5 (28.7)139.6 (27.9)137.6 (28.9)138.1 (30.9)0.45
Heart rate [beats/min]75.5 (17.8)76.7 (18.9)75.4 (17.3)72.1 (14.6)0.17
Killip class III or IV (%)6.5 (n = 36)7.8 (n = 19)6.6 (n = 16)1.5 (n = 1)0.17
ST-elevation in anterior leads (%)39.0 (n = 216)42.4 (n = 103)36.8 (n = 89)35.3 (n = 24)0.39
TIMI flow grade 3 after procedure (%)91.7 (n = 507)93.4 (n = 227)90.5 (n = 219)89.7 (n = 61)0.41
Stent implantation (%)76.8 (n = 425)76.5 (n = 186)78.9 (n = 191)70.6 (n = 48)0.35
No. of vessels with significant stenosis1.71 (0.8)1.72 (0.82)1.71 (0.79)1.69 (0.79)0.93
Creatinine [mg/dl]1.036 (0.41)1.06 (0.48)1.0 (0.34)1.03 (0.32)0.12
Total cholesterol [mg/dl]195.8 (42.3)191.7 (45.6)198.8 (40.2)199.8 (36.2)0.06
LDL cholesterol [mg/dl]128.4 (37.8)124.9 (40.6)131.2 (36.2)131.1 (32.4)0.09
HDL cholesterol [mg/dl]43.8 (13.2)44.0 (14.6)43.6 (11.9)43.7 (12.1)0.85
Triglycerides [mg/dl]123.0 (66.3)120.7 (64.5)122.3 (61.0)133.7 (88.0)0.77
Hemoglobin [g/dl]13.5 (7.0)13.5 (6.7)13.6 (8.2)13.1 (1.7)0.88
Ejection fraction (%)46.0 (9.5)45.3 (9.6)46.6 (9.7)46.5 (8.8)0.19
GRACE risk score149.8 (34.8)151.9 (33.7)148.8 (35.5)146 (28.0)0.23
Baseline characteristics of the study group based on rs2228145 genotype (risk allele for myocardial infarction – A) A long-term follow-up was performed (median of censored observations 2875 days, maximum 4167, minimum 2558). During observation 171 (30.9%) patients died. Mortality rates for the rs2228145 genotypes are presented in Table I. There was a trend for lower mortality in C allele carriers compared to AA homozygotes (28.1% vs. 34.6%, p = 0.1, χ2 test). Figure 1 shows Kaplan-Meier survival curves for the rs2228145 genotypes and long-term observation.
Figure 1

Kaplan-Meier survival curves for the rs2228145 genotypes and long-term observation. P = 0.38 (analysis of 3 genotypes); p = 0.14 (AA homozygotes and C allele carriers), log-rank test

Kaplan-Meier survival curves for the rs2228145 genotypes and long-term observation. P = 0.38 (analysis of 3 genotypes); p = 0.14 (AA homozygotes and C allele carriers), log-rank test The trend towards adverse outcome in AA homozygotes compared to C allele carriers was not statistically significant (p = 0.14, log-rank test). In univariate analysis (logistic regression, Table III) the rs2228145 AA genotype was not significantly associated with long-term outcome (OR = 1.34, 95% CI: 0.93–1.93, p = 0.1). In multivariate analysis the variables independently associated with long-term outcome were age, Killip class on admission, and ejection fraction (Table III).
Table III

Univariate and multivariate analysis for long-term mortality

VariableHazard ratio95% CI P-value
Univariate analysis:
 Age [years]1.071.05–1.088< 0.0001
 Heart rate [beats/min]1.0121.002–1.0220.018
 Systolic blood pressure [mm Hg]0.921.0003–0.990.91
 Killip class2.071.6–2.65< 0.0001
 Type 2 diabetes1.91.2–2.80.003
 Previous myocardial infarction2.01.1–3.50.017
 Anterior myocardial infarction1.30.91–1.890.14
 TIMI 3 flow after PCI0.520.28–0.940.03
 No. of vessels with significant stenosis1.180.94–1.470.14
 Ejection fraction (%)0.950.93–0.96< 0.0001
 Creatinine [mg/dl]2.21.3–3.70.003
 GRACE risk score1.0211.015–1.027< 0.0001
 rs2228124 AA genotype1.340.93–1.930.1
Multivariate analysis:
 Age [years]1.061.04–1.08< 0.0001
 Killip class1.671.26–2.20.0002
 Ejection fraction (%)0.970.94–0.990.004
Univariate and multivariate analysis for long-term mortality

Discussion

We did not prove an association between the rs2228145 polymorphism and long-term mortality after STEMI. However, AA homozygotes showed a trend towards adverse outcome compared to C allele carriers. In the early phase of follow-up, Kaplan-Meier survival curves overlapped, and they started to diverge in the third year of observation. This confirms that the reported effect was independent of potential differences in the acute phase of myocardial infarction. Our findings are consistent with previous papers defining adenosine as a risk allele for ischemic heart disease and cytosine as a protective one [9]. In a recent study performed in patients with inflammatory polyarthritis (the Norfolk Arthritis Register) the rs2228145 polymorphism was associated neither with all-cause nor with cardiovascular mortality [22]. The investigated cohort, however, cannot be directly compared to our homogenous group of patients with STEMI, who are at notably high risk of death. Inflammatory polyarthritis is a disorder of a very heterogeneous etiology. It includes various forms of systemic rheumatic illnesses (e.g. rheumatoid arthritis, lupus), inflammatory osteoarthritis, infectious and postinfectious arthritis, seronegative spondyloarthritides, and other systemic disorders (e.g. sarcoidosis, malignancies). All those illnesses affect multiple organs and also increase the risk of cardiovascular death, but not as specifically as a history of myocardial infarction (only 54% of deaths in the Norfolk Arthritis Register were reported as cardiovascular). The next important difference is the large dispersion in time of observation: in the Norfolk Arthritis Register, patients were enrolled in the years 1990–2011. During that long time, 23% of them died, which is comparable to the 5-year mortality of patients with STEMI [23]. Finally, patients from that study were not clinically verified for previous cardiovascular disease, and cardiovascular deaths were defined solely based on death certificates, which may be unreliable. The C allele carriers have increased sIL6R concentrations [9-11], decreased IL6R membrane expression, and decreased response to IL-6 [10]. Under normal conditions, a majority of the IL-6 in the plasma forms a complex with sIL6R. The sIL6R-IL-6 complex binds membrane gp130, which is available on most cell types in the body and is responsible for further signal transduction. This phenomenon is called IL-6 trans-signaling. Alternatively, the complex can be bound by a soluble form of gp130 that acts as a trans-signaling inhibitor. The free form of IL-6, which is in minority, binds membrane IL6R that acts further via membrane gp130. The membrane IL6R is present only on selected cell types: hepatocytes, macrophages/monocytes, neutrophils, and some subpopulations of lymphocytes. Its availability on cardiomyocytes has been discussed. To sum up, an IL-6 complex with sIL6R interacts with a majority of cells in the body (if not inhibited by sgp130) and free IL-6 can stimulate only specific cells with membrane IL6R. In C allele carriers, the balance between the soluble and membrane form of IL6R is shifted, probably due to increased cleavage of the membrane receptor to plasma [12] and further change in the effector cells. This results in a phenotype of decreased prevalence of cardiovascular disease, including ischemic heart disease [9]. Such an effect after myocardial infarction can potentially extend to a decreased number of subsequent cardiovascular events and mortality. The function of IL-6 is very complex: both pro- and anti-inflammatory effects are described. Classically, IL-6 increases production of acute phase proteins in liver and expression of adhesive molecules in endothelium. These mechanisms promote atherosclerosis. On the other hand, IL-6 enhances production of the IL-10 and IL-1 receptor antagonist and has a cytoprotective effect for cardiomyocytes [24]. An increase in local expression of IL-6 has been observed in ischemia-reperfusion models, experimental [25, 26] as well as clinical [27, 28]. Previous studies in our material supported the hypothesis that during scheduled percutaneous coronary interventions an increase in IL-6 plasma concentration has a global rather than a local character, as there were no significant gradients between the aorta and coronary sinus [29]. All the reports make interpretation of the role of IL6R and its associated SNP more complex and difficult; however, there is a general opinion that the effects of IL-6 trans-signaling are predominantly pro-inflammatory [30]. Despite all the advances in treatment of acute myocardial infarction, long-term prognosis in this group of patients remains poor. Therefore, every effort is being made to improve not only therapy but also risk stratification. Novel parameters helpful in the assessment of prognosis are continuously searched for [31-33]. High-risk patients require more intensive ambulatory monitoring and more aggressive treatment. Basic clinical parameters will always have the greatest value in this field. However, due to decreasing costs and increasing availability, genetic tests are still a promising option that could supplement the traditional approach. Several limitations of the study need to be acknowledged. First, the number of patients enrolled in the research was not sufficient to prove a significant correlation between the genotype and survival. Further research is required in this field. Next, the effect of the genotype was revealed after a very long-term follow-up. In our case, the mean time of observation was almost 8 years, and the first effect started in the 3rd year of observation. The longer the observation, the less specific the effect – in the case of this analysis, the less specific for cardiac mortality. Moreover, we had no possibility to analyze plasma samples for IL-6, IL6R, or CRP concentrations, which would have improved the quality of the study. Finally, the research was performed retrospectively (although the data were collected in a prospective manner). In conclusion, The rs2228145 polymorphism of the IL6R was not significantly associated with long-term mortality after STEMI. However, AA homozygotes (a high-risk genotype for ischemic heart disease) showed a trend towards an adverse outcome compared to C allele carriers. Further research in this field is required, because the observed effect size is promising for clinical utilization.
  33 in total

Review 1.  The role of interleukin-6 in the failing heart.

Authors:  K C Wollert; H Drexler
Journal:  Heart Fail Rev       Date:  2001-03       Impact factor: 4.214

2.  Relationship between interleukin 6 and mortality in patients with unstable coronary artery disease: effects of an early invasive or noninvasive strategy.

Authors:  E Lindmark; E Diderholm; L Wallentin; A Siegbahn
Journal:  JAMA       Date:  2001-11-07       Impact factor: 56.272

3.  Acute cardiac inflammatory responses to postischemic reperfusion during cardiopulmonary bypass.

Authors:  S Zahler; P Massoudy; H Hartl; C Hähnel; H Meisner; B F Becker
Journal:  Cardiovasc Res       Date:  1999-03       Impact factor: 10.787

4.  The rs9982601 polymorphism of the region between the SLC5A3/MRPS6 and KCNE2 genes associated with a prevalence of myocardial infarction and subsequent long-term mortality.

Authors:  Anna Szpakowicz; Marek Kiliszek; Witold Pepiński; Ewa Waszkiewicz; Maria Franaszczyk; Małgorzata Skawrońska; Sławomir Dobrzycki; Anna Niemcunowicz-Janica; Rafał Ploski; Grzegorz Opolski; Włodzimierz J Musiał; Karol A Kamiński
Journal:  Pol Arch Med Wewn       Date:  2015-02-20

5.  Variation in the interleukin-6 receptor gene associates with type 2 diabetes in Danish whites.

Authors:  Yasmin H Hamid; Søren A Urhammer; Dorit P Jensen; Charlotte Glümer; Knut Borch-Johnsen; Torben Jørgensen; Torben Hansen; Oluf Pedersen
Journal:  Diabetes       Date:  2004-12       Impact factor: 9.461

6.  Cardiac cytokine expression is upregulated in the acute phase after myocardial infarction. Experimental studies in rats.

Authors:  Alexander Deten; Hans Christian Volz; Wilfried Briest; Heinz-Gerd Zimmer
Journal:  Cardiovasc Res       Date:  2002-08-01       Impact factor: 10.787

7.  Proinflammatory cytokine levels in patients with depressed left ventricular ejection fraction: a report from the Studies of Left Ventricular Dysfunction (SOLVD).

Authors:  G Torre-Amione; S Kapadia; C Benedict; H Oral; J B Young; D L Mann
Journal:  J Am Coll Cardiol       Date:  1996-04       Impact factor: 24.094

8.  Inflammatory markers and risk of heart failure in elderly subjects without prior myocardial infarction: the Framingham Heart Study.

Authors:  Ramachandran S Vasan; Lisa M Sullivan; Ronenn Roubenoff; Charles A Dinarello; Tamara Harris; Emelia J Benjamin; Douglas B Sawyer; Daniel Levy; Peter W F Wilson; Ralph B D'Agostino
Journal:  Circulation       Date:  2003-03-25       Impact factor: 29.690

9.  The IL6R variation Asp(358)Ala is a potential modifier of lung function in subjects with asthma.

Authors:  Gregory A Hawkins; Mac B Robinson; Annette T Hastie; Xingnan Li; Huashi Li; Wendy C Moore; Timothy D Howard; William W Busse; Serpil C Erzurum; Sally E Wenzel; Stephen P Peters; Deborah A Meyers; Eugene R Bleecker
Journal:  J Allergy Clin Immunol       Date:  2012-05-01       Impact factor: 14.290

10.  Functional IL6R 358Ala allele impairs classical IL-6 receptor signaling and influences risk of diverse inflammatory diseases.

Authors:  Ricardo C Ferreira; Daniel F Freitag; Antony J Cutler; Joanna M M Howson; Daniel B Rainbow; Deborah J Smyth; Stephen Kaptoge; Pamela Clarke; Charlotte Boreham; Richard M Coulson; Marcin L Pekalski; Wei-Min Chen; Suna Onengut-Gumuscu; Stephen S Rich; Adam S Butterworth; Anders Malarstig; John Danesh; John A Todd
Journal:  PLoS Genet       Date:  2013-04-04       Impact factor: 5.917

View more
  2 in total

1.  Ectopic upregulation of membrane-bound IL6R drives vascular remodeling in pulmonary arterial hypertension.

Authors:  Yuichi Tamura; Carole Phan; Ly Tu; Morane Le Hiress; Raphaël Thuillet; Etienne-Marie Jutant; Elie Fadel; Laurent Savale; Alice Huertas; Marc Humbert; Christophe Guignabert
Journal:  J Clin Invest       Date:  2018-04-09       Impact factor: 14.808

2.  Transmembrane protease serine 2 (TMPRSS2) rs75603675, comorbidity, and sex are the primary predictors of COVID-19 severity.

Authors:  Gonzalo Villapalos-García; Pablo Zubiaur; Rebeca Rivas-Durán; Pilar Campos-Norte; Cristina Arévalo-Román; Marta Fernández-Rico; Lucio García-Fraile Fraile; Paula Fernández-Campos; Paula Soria-Chacartegui; Sara Fernández de Córdoba-Oñate; Pablo Delgado-Wicke; Elena Fernández-Ruiz; Isidoro González-Álvaro; Jesús Sanz; Francisco Abad-Santos; Ignacio de Los Santos
Journal:  Life Sci Alliance       Date:  2022-05-30
  2 in total

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