Literature DB >> 35517435

A novel catalytic kinetic method for the determination of mercury(ii) in water samples.

Abhinav Agarwal1, Amit Kumar Verma1, Masafumi Yoshida2, Radhey Mohan Naik1, Surendra Prasad3.   

Abstract

Mercury(ii) ions act as catalyst in the substitution of cyanide ion in hexacyanoruthenate(ii) by pyrazine (Pz) in an acidic medium. This property of Hg(ii) has been utilized for its determination in aqueous solutions. The progress of reaction was followed spectrophotometrically by measuring the increase in absorbance of the yellow colour product, [Ru(CN)5Pz]3- at 370 nm (λ max, ε = 4.2 × 103 M-1 s-1) under the optimized reaction conditions; 5.0 × 10-5 M [Ru(CN)6 4-], 7.5 × 10-4 M [Pz], pH 4.00 ± 0.02, ionic strength (I) = 0.05 M (KCl) and temp. 45.0 ± 0.1 °C. The proposed method is based on the fixed time procedure under optimum reaction conditions. The linear regression (calibration) equations between the absorbance at fixed times (t = 15, 20 and 25 min) and [Hg(ii)] were established in the range of 1.0 to 30.0 × 10-6 M. The detection limit was found to be 1.5 × 10-7 M of Hg(ii). The effect of various foreign ions on the proposed method was also studied and discussed. The method was applied for the determination of Hg(ii) in different wastewater samples. The present method is simple, rapid and sensitive for the determination of Hg(ii) in trace amount in the environmental samples. This journal is © The Royal Society of Chemistry.

Entities:  

Year:  2020        PMID: 35517435      PMCID: PMC9055178          DOI: 10.1039/d0ra03487h

Source DB:  PubMed          Journal:  RSC Adv        ISSN: 2046-2069            Impact factor:   3.361


Introduction

Mercury (Hg) occurs naturally in the earth's crust.[1] Its concentration increases in the environment and water bodies due to the natural processes but mostly accelerated by the anthropogenic activities from pesticides, paints, batteries, industrial waste and land application of industrial or domestic sludge.[1,2] Thus, it gradually concentrates and moves away through the surface run-off transporting to the aquatic systems as rivers, lakes, seas and oceans.[2,3] Antagonistically to the organic contaminants, heavy metals such as Hg does not undergo microbial or chemical degradation and thus persist for a long time in the environment.[3,4] Therefore, heavy metals including Hg are the most hazardous pollutants to the environment because of the rapid industrialization and urbanization.[4] Thus, there has been continuous considerable emphasis on the Hg analysis in different type of samples.[5] Among the various heavy metals, particularly Hg has been listed as a priority pollutant by the international environmental and health agencies because of its persistence, bioaccumulation and toxicity (PBT) in the environment.[6] It not only causes serious health problems but also poses a great challenge to the environmental protection.[7,8] Among other pollutants, Hg is of major concern in aquatic environments. Thus, Hg is considered as a global and recalcitrant pollutant due to its biogeochemical properties and its toxicity that can affect the health of human and ecosystems.[4] Though metallic Hg is an insoluble element, it is easily oxidized to the soluble ionic form Hg2+ in the freshwater reservoirs and constitute a serious threat through the process of biomagnification.[9,10] Thus, mercury is of major concern in aquatic environments.[11] The consumption of contaminated fish, sea mammals and ground water are the prominent environmental sources of Hg exposure in humans.[11,12] It forms quite stable complexes with sulfhydryl (–SH) groups in the human body forming mercaptides having mobility through the tissues,[13] which lead to the inactivation of numerous enzyme reactions, amino acids and sulfur containing antioxidants such as N-acetyl cysteine (NAC), alpha-lipoic acid (ALA), and glutathione (GSH) and makes it quite toxic to the human being (cf. following general reaction).[14]2R–SH + Hg The accumulation of Hg in the body is associated with hazardous health effects, such as gastrointestinal and nervous system disorders, respiratory and acute renal failures, hypertension, coronary heart disease (CHD), and cardiovascular disease (CVD).[15] Thus, even a small amount of Hg can have adverse effects on the human health.[16] One of the most stable forms of Hg is the highly toxic, water-soluble, divalent mercuric ion (Hg2+), which is widely distributed in soil, air and aquatic environments as a result of illegal waste release from various anthropogenic and industrial activities. Hg2+ can cause serious environmental pollution and result in permanent damage to the biological organisms due to its acute and/or chronic toxicity.[17] Thus, the health and environmental issues associated with Hg2+ has stimulated researchers to develop inexpensive, effective and reliable methods for the detection and determination of Hg2+ sensitively and selectively.[5,7,8,17,18] Attempts have been made to determine Hg in the environmental,[18] geological,[19] food[20,21] and biomedical[22,23] samples. Majority of the analytical methods proposed for the determination of Hg are based on the sophisticated techniques viz. cold vapor atomic absorption spectrometry (CV-AAS),[23-25] inductively coupled plasma mass spectrometry (ICP-MS),[26,27] atomic fluorescence spectrometry (AFS),[28,29] anodic stripping voltammetry (ASV),[30-32] high performance liquid chromatography (HPLC),[33,34]etc. All these techniques require expensive instrumentation, pre-concentration for the enrichment of analyte and its separation from the other constituents of the sample to minimize the matrix effect prior to their application. In addition, skilled operators, high cost, long time for running per sample, etc., have been other disadvantages. To overcome these problems, attempts have been made to develop analytical methods based on the kinetics and catalysis,[5,18,35-38] formation of colored complexes of Hg(ii) with other compounds,[39,40]etc. We have also been interested in developing analytical methods for the determination as well as removal of various analytes and toxic species of environmental, biological and medicinal interest.[5,18,35-38,41,42] Keeping the above explained backgrounds in mind, in the present communication, a successful attempt has been made to develop a novel, simple and precise catalytic kinetic method (CKM) for the trace level determination of Hg(ii) in the aqueous samples which is based on the Hg(ii) catalyzed substitution of (CN−) in the [Ru(CN)6]4− by pyrazine (Pz). The kinetics and mechanism of this reaction has already been reported by us[43] and has now been exploited as an indicator reaction system for the CKM in the determination of Hg(ii). The developed CKM has successfully been utilized for the determination of Hg(ii) in the wastewater samples.

Experimental

Reagents and instrumentation

All the reagents used were of analytical grade and deionized distilled water (DDW) was used throughout the study for preparation of all the solutions. The stock solutions (1.0 × 10−2 M) of K4[Ru(CN)6]·3H2O, Pz and HgCl2 were prepared by dissolving their appropriate amounts in DDW. K4[Ru(CN)6]·3H2O solution was kept in a dark amber colour volumetric flask to prevent its photodecomposition. The working solutions of these reagents were prepared by their appropriate dilution from the respective stock solutions as required. Potassium hydrogen phthalate (PHP)–HCl buffer of pH 4.00 ± 0.02 was prepared according to the literature reported method.[44] Kinetic measurements and recording of various spectral scans were carried out in 10 mm matched quartz cuvettes using a Shimadzu UV-240 double beam spectrophotometer equipped with a self-designed thermostatic cell compartment. A remi ultra-cryostat was used to maintain the temperature of the reaction system. All the pH measurements were made on a Toshniwal digital pH meter model CL46. A certified ‘A’ grade volumetric apparatus were used throughout the work.

Procedure

All the required solutions were first placed in the thermostat maintained at 45.0 ± 0.1 °C for 30 min prior to their use to attain thermal equilibrium. 2.0 mL each of the solutions of Pz, buffer, HgCl2 and K4[Ru(CN)6] was pipetted out and mixed in a 10 mL volumetric flask, which was also placed in the same thermostat, shaken quickly and transferred immediately into a 10 mm spectrophotometric cuvette placed in the cell compartment of the spectrophotometer at 45.0 ± 0.1 °C. The progress of the reaction was monitored spectrophotometrically by measuring the increase in absorbance due to the formation of the product [Ru(CN)5Pz]3− at its λmax of 370 nm. A fixed time procedure was used to record the absorbance as a function of the concentration of Hg(ii).The cuvettes were cleaned with acetone after few kinetic runs in order to remove the deposited intense yellow-colored [Ru(CN)5Pz]3− complex.

Results & discussion

The indicator reaction

The detailed kinetics and mechanism of Hg(ii) catalyzed substitution of CN− in the [Ru(CN)6]4− by Pz have been studied and reported by us earlier.[43] The reaction product, [Ru(CN)5Pz]3− complex has strong absorption band at 370 nm and molar extinction coefficient ε = 4.2 × 103 M−1 s−1. The UV-visible spectra of the reactants and the product shown in Fig. 1 clearly shows no interference between the absorbance of the product and reactants i.e. [Ru(CN)6]4− and Pz. Hence, the progress of the reaction was easily followed spectrophotometrically by measuring the increase in absorbance at 370 nm i.e. λmax of [Ru(CN)5Pz]3− without making any corrections for the absorbance due to the reactants [Ru(CN)6]4− and Pz.
Fig. 1

UV-visible absorption spectra of reactants and products under the conditions: (A) [Pz] = 2.5 × 10−6 M, (B) [Ru(CN)64−] = 1.25 × 10−5 M, (C) product [Ru(CN)5Pz3−] formed by mixing [Ru(CN)64−] = 5.0 × 10−5 M, [Pz] = 7.5 × 10−4 M, [Hg(ii)] = 1.5 × 10−4 M, I = 0.05 M (KCl), pH = 4.00 ± 0.02 and temp. = 45.0 ± 0.1 °C, (D) [Hg(ii)] = 1.0 × 10−4 M.

The rate law and its analytical application

Based on the kinetic and mechanistic studies,[43] the overall rate of the indicator reaction comprises of the sum of the rate of the uncatalyzed as well as catalyzed reactions and is given by eqn (1). If K[Ru(CN)64−] ≫ 1, eqn (1) reduces to eqn (2) as follows: In eqn (2), first term on the right hand side corresponds to the rate of the uncatalyzed path and involves the rate constant for the slow decomposition of the activated complex (formed between [Ru(CN)6]4−, Hg2+ and H2O) and some other concentration terms.[43]Eqn (2) clearly shows a direct correlation between the rate of the indicator reaction and the concentration of Hg(ii) and envisages for the accurate determination of Hg(ii). In the present study, the rate of the uncatalyzed reaction is almost negligible. Thus, the catalytic kinetic method (CKM) for the determination of Hg(ii) based on the [Ru(CN)64−]–Pz indicator reaction will be more accurate in comparison to the other CKM involving the [Fe(CN)64−] – ligand indicator reactions.[35-38] The increase in absorbance (A) at different time t (t = 15, 20 and 25 min) after mixing the reactants was recorded as a measure of the initial rate under the optimized reaction condition as 5.0 × 10−5 M [Ru(CN)54−], 7.5 × 10−4 M [Pz], pH 4.00 ± 0.02, ionic strength (I) 0.05 M (KCl) and temperature 45.0 ± 0.1 °C,[43] using fixed time procedure. The calibration curves were obtained by plotting the values of A against different Hg(ii) concentrations where the corresponding linear regression equations correlating A to the [Hg(ii)] obtained are shown in eqn (3)–(5).

Quantitative determination of Hg(ii)

In order to validate the analytical applicability of the proposed CKM, recovery experiments were performed in various water samples spiked with Hg(ii) and the results obtained for the three fixed times along with the standard deviations and the percentage errors are shown in Table 1. From Table 1, it is clear that the calibration curve corresponding to the A15 is a closer measure to the initial rate with the minimal percentage error as compared to A20 and A25. Time less than 15 min and more than 25 min were not taken to minimize the experimental error. Based on this observation, the A15 calibration curve is recommended for the Hg(ii) determination in the aqueous samples.

Evaluation of percentage error in the determination of Hg(ii) under conditions: [Ru(CN)64−] = 5.0 × 10−5 M, [Pz] = 7.5 × 10−4 M, pH = 4.00 ± 0.02, I = 0.05 M (KCl), temp. = 45.0 ± 1.0 °C

[Hg2+] × 106 M (taken) A 15 A 20 A 25
[Hg2+] × 106 M (found)b ± SDaError (%)[Hg2+] × 106 M (found)b ± SDaError (%)[Hg2+] × 106 M (found)b ± SDaError (%)
1.000.94 ± 0.08−6.000.94 ± 0.02−6.001.06 ± 0.05+6.00
2.002.25 ± 0.06+12.502.00 ± 0.060.001.88 ± 0.03+6.00
3.003.00 ± 0.080.003.05 ± 0.02+1.673.06 ± 0.06+2.00
4.004.02 ± 0.02+0.504.05 ± 0.01+0.804.06 ± 0.05+1.50
5.005.05 ± 0.05+1.005.05 ± 0.02+1.005.03 ± 0.09+6.00
10.0010.00 ± 0.070.0010.05 ± 0.06+0.5010.05 ± 0.06+0.50
15.0015.25 ± 0.02+1.7015.25 ± 0.03+1.7015.13 ± 0.04+0.87
20.0020.00 ± 0.030.0020.02 ± 0.03+0.1020.25 ± 0.04+1.25
30.0029.80 ± 0.05−0.6729.50 ± 0.04−1.7030.25 ± 0.02+0.83

The ±SD values represent the % relative standard deviation of the mean for three determinations.

Mean of three determinations.

The ±SD values represent the % relative standard deviation of the mean for three determinations. Mean of three determinations. The detection limit (Xd), defined as three times standard deviation of the blank, was evaluated using Tanaka's method[45] using eqn (6), where Xb is the average of the blank values, (SD)b is the standard deviation of blank and equal to R/d2, R is blank kmax − blank kmin and 1/d2 is a factor for obtaining (SD)b from the range R of n replicates whose value is recommended as 0.5908 for the use in three blank measurements. The detection limit for the present method corresponding to the A15 calibration curve under the optimized experimental conditions was found to be 1.5 × 10−7 M.

Interference study

The general precision and accuracy for the determination of Hg(ii) by the proposed CKM was tested in the presence of several cations, anions and complexing agents. For this purpose, the recovery i.e. the determination of Hg(ii) was performed using a fixed concentration of Hg(ii), 4.0 × 10−6 M, in the presence of different known concentrations of individual cations, anions and complexing agents. The concentration of Hg(ii) was determined using the A15 calibration eqn (3). The recovery results on the determination of 4.0 × 10−6 M Hg(ii) in the presence of different individual cations, anions and complexing agents along with their concentrations have been presented in Table 2. Based on the results of Hg(ii) recovered, it was confirmed that Mn2+, Ni2+, Ag+, Sn2+, Pb2+, Mg2+, Ca2+, Cd2+, Cu2+, Al3+, Fe3+, Br−, I−, NO3−, CO32−, C2O42−, SO42− ions did not interfere in the determination of Hg(ii) but Li+, Co2+, Zn2+, NTA, IDA, EDTA interfered significantly at their concentrations reported in Table 2. Metal hexacyanoferrates are very efficient sorbents for the recovery of alkali metal ions and especially Li+, which leads to the formation of the corresponding lithium hexacyanoruthenate(ii) complex as its Fe(ii) counterpart in the solution.[46] Hence, a low recovery of Hg(ii) is very much expected in the presence of Li+ ions in the aqueous medium. The interference by NTA, IDA, EDTA may be attributed to the possible ligand substitution between the monodentate cyanide ligand in [Ru(CN)6]4− and the ligands viz. NTA, IDA, EDTA. Interference by Co2+, Cu2+, Zn2+ is probably due to the complex formation between these metal ions and [Ru(CN)6]4− or Pz.

Effect of different foreign ions on the determination of Hg(ii) under the conditions: [Ru(CN)64−] = 5.0 × 10−5 M, [Pz] = 7.5 × 10−4 M, [Hg2+] = 4.0 × 10−6 M, pH = 4.00 ± 0.02, I = 0.05 M (KCl), temp. = 45.0 ± 1.0 °C

Foreign ion[Foreign ion] (added), M[Hg2+] (found), MError (%)Inference
Mn2+2 × 10−44.1× 10−6+2.50No inference
Ni2+4 × 10−43.9 × 10−6−2.50No inference
Ag+8 × 10−44.0 × 10−60.00No inference
Sn2+8 × 10−44.1 × 10−6+2.50No inference
Pb2+4 × 10−43.8 × 10−6−5.00No inference
Mg2+4 × 10−44.1 × 10−6+2.50No inference
Ca2+4 × 10−43.9 × 10−6−2.50No inference
Cd2+4 × 10−44.0 × 10−60.00No inference
Cu2+8 × 10−43.8 × 10−6−5.00No inference
Al3+4 × 10−44.0 × 10−60.00No inference
Fe3+8 × 10−43.9 × 10−6−2.50No inference
Li+8 × 10−51.2 × 10−6−70.00Interfered
Co2+2 × 10−51.6 × 10−6−60.00Interfered
Zn2+4 × 10−51.9 × 10−6−52.50Interfered
Br4 × 10−44.1 × 10−6+2.50No inference
I8 × 10−43.9 × 10−6−2.50No inference
NO38 × 10−44.0 × 10−60.00No inference
CO32−8 × 10−43.9 × 10−6−2.50No inference
C2O42−4 × 10−44.0 × 10−60.00No inference
SO42−4 × 10−44.1 × 10−6+2.50No inference
NTA2 × 10−42.4 × 10−6−40.00Interfered
IDA2 × 10−42.2 × 10−6−45.00Interfered
EDTA2 × 10−41.9 × 10−6−52.50Interfered

Validation and analytical application of CKM

In order to validate the proposed method, five water spiked synthetic mixtures (SMs) were prepared, which contained other metal ions along with the catalyst Hg(ii). This was followed by the determination of Hg(ii) in five different synthetic mixtures (SM-1 to SM-5) using the proposed CKM. The results were further confirmed by atomic absorption spectrophotometry (AAS), as shown in Table 3. The results obtained by CKM are in close agreement with those determined by AAS (Table 3).

Determination of Hg(ii) in synthetic mixtures under optimized reaction conditions: [Ru(CN)64−] = 5.0 × 10−5 M, [Pz] = 7.5 × 10−4 M, [Hg2+] = 4.0 × 10−6 M, pH = 4.00 ± 0.02, I = 0.05 M (KCl), temp. = 45.0 ± 1.0 °C and its comparison atomic absorption spectrometry (AAS) method

Synthetic mixtures (SM)Composition of SM (ng mL−1)[Hg2+] foundb ± SDa by CKM (ng mL−1)Recovery (%)[Hg2+] foundb ± SDa by AAS (ng mL−1)
SM-1Hg 1086.0 + Ca 4202.01086.2 ± 1.1100.021086.1 ± 4.3
Cu 5006.0 + Cd 4808.0
SM-2Hg 1086.0 + Ag 1206.51085.9 ± 1.299.991085.9 ± 5.2
Mg 4024.0 + Ni 1435.0
SM-3Hg 1086.0 + Mn 3264.01086.6 ± 1.6100.051086.1 ± 3.6
Ba 2244.0 + Cr 206.0
SM-4Hg 1086.0 + Pd 2206.51085.2 ± 1.899.931086.3 ± 3.6
Al 1188.0 + Ag 608.0
SM-5Hg 1086.0 + Ca 1224.01086.1 ± 1.1100.011085.9 ± 5.5
Mg 1108.0 + Fe 906.0

The ±SD values represent the % relative standard deviation of the mean for three determinations.

Mean of three determinations.

The ±SD values represent the % relative standard deviation of the mean for three determinations. Mean of three determinations. After validation, the method was successfully applied for the determination of Hg(ii) in the wastewater. For this, six wastewater samples (WWS-1 to WWS-6) were collected from our laboratory on different days and Hg(ii) was determined using the proposed CKM. The results obtained for different WWS are shown in Table 4. The results obtained by CKM were also validated by AAS and were in excellent agreement with those determined by CKM (Table 3). The maximum error in the proposed CKM with respect to AAS is less than 5%, which is well accepted in the environmental analysis.

Determination of Hg(ii) in wastewater samples (WWS) under optimized reaction conditions: [Ru(CN)64−] = 5.0 × 10−5 M, [Pz] = 7.5 × 10−4 M, pH = 4.00 ± 0.02, I = 0.05 M (KCl), temp. = 45.0 ± 1.0 °C and its comparison with the (AAS) method

Wastewater sample (WWS)[Hg2+] foundb ± SDa by CKM (ng mL−1)[Hg2+] foundb ± SDa by AAS (ng mL−1)Error in CKM versus AAS (%)
WWS-1224.06 ± 2.50228.70 ± 3.26+2.07
WWS-2329.64 ± 4.35331.65 ± 2.35+1.44
WWS-3244.16 ± 3.84237.45 ± 4.50−2.75
WWS-4402.25 ± 4.36408.24 ± 3.45+1.50
WWS-5398.17 ± 5.50390.75 ± 3.60−1.86
WWS-6389.07 ± 5.50407.08 ± 4.60+4.63

The ±SD values represent % relative standard deviation of the mean for three determinations.

Mean of three determinations.

The ±SD values represent % relative standard deviation of the mean for three determinations. Mean of three determinations.

Conclusions

Though there are many sophisticated methods available for the Hg(ii) determination but the present CKM based on the indicator reaction between [Ru(CN)6]4−–Pz offers a number of advantages over them, such as the unanalyzed reaction rate is negligible under specified reaction conditions and only 15 min is required for the Hg(ii) determination without using any costly solvents and instruments. In addition, the method can be successfully applied for the determination of Hg(ii) in the presence of a number of cations, anions and complexing agents with very good sensitivity. In conclusion, the proposed CKM is quite sensitive, quick, economical and superior to few known methods including CKM.[35-38] It can selectively be applied for the Hg(ii) determination at micro-level in environmental water samples.

Conflicts of interest

The authors are declare that there is not any conflict of interest including any financial, personal or other relationships with other people or organizations.
  16 in total

Review 1.  The toxicology of mercury--current exposures and clinical manifestations.

Authors:  Thomas W Clarkson; Laszlo Magos; Gary J Myers
Journal:  N Engl J Med       Date:  2003-10-30       Impact factor: 91.245

2.  Chemical and ecological effects of contaminated tunnel wash water runoff to a small Norwegian stream.

Authors:  Sondre Meland; Reidar Borgstrøm; Lene Sørlie Heier; Bjørn Olav Rosseland; Oddvar Lindholm; Brit Salbu
Journal:  Sci Total Environ       Date:  2010-06-12       Impact factor: 7.963

3.  Determination of mercury in water and fish samples by cold vapor atomic absorption spectrometry after solid phase extraction on agar modified with 2-mercaptobenzimidazole.

Authors:  N Pourreza; K Ghanemi
Journal:  J Hazard Mater       Date:  2008-04-20       Impact factor: 10.588

4.  A simple and low cost dual-wavelength β-correction spectrophotometric determination and speciation of mercury(II) in water using chromogenic reagent 4-(2-thiazolylazo) resorcinol.

Authors:  A H Al-Bagawi; W Ahmad; Z M Saigl; H Alwael; E A Al-Harbi; M S El-Shahawi
Journal:  Spectrochim Acta A Mol Biomol Spectrosc       Date:  2017-07-01       Impact factor: 4.098

5.  Mercury contamination in the sludge of drinking water treatment plants dumping into a reservoir in Rio de Janeiro, Brazil.

Authors:  Julio Cesar Wasserman; Letícia de Oliveira Silva; Gabriela Cugler de Pontes; Evaldo de Paiva Lima
Journal:  Environ Sci Pollut Res Int       Date:  2018-08-11       Impact factor: 4.223

6.  Biological methylation of mercury in aquatic organisms.

Authors:  S Jensen; A Jernelöv
Journal:  Nature       Date:  1969-08-16       Impact factor: 49.962

7.  A novel enhanced fluorescence method based on multifunctional carbon dots for specific detection of Hg2+ in complex samples.

Authors:  Ying Xu; Yao Fan; Lei Zhang; Qin Wang; Haiyan Fu; Yuanbin She
Journal:  Spectrochim Acta A Mol Biomol Spectrosc       Date:  2019-05-11       Impact factor: 4.098

8.  A simple method for methylmercury, inorganic mercury and ethylmercury determination in plasma samples by high performance liquid chromatography-cold-vapor-inductively coupled plasma mass spectrometry.

Authors:  Samuel S de Souza; Andres Dobal Campiglia; Fernando Barbosa
Journal:  Anal Chim Acta       Date:  2012-11-29       Impact factor: 6.558

9.  Nitrogen-doped carbon quantum dots: facile synthesis and application as a "turn-off" fluorescent probe for detection of Hg2+ ions.

Authors:  Ruizhong Zhang; Wei Chen
Journal:  Biosens Bioelectron       Date:  2013-12-10       Impact factor: 10.618

Review 10.  A review on the distribution of Hg in the environment and its human health impacts.

Authors:  Ki-Hyun Kim; Ehsanul Kabir; Shamin Ara Jahan
Journal:  J Hazard Mater       Date:  2015-11-21       Impact factor: 10.588

View more

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