Literature DB >> 35368883

Optimization of COD, nitrate-N and phosphorus removal from hatchery wastewater with acclimatized mixed culture.

Norazwina Zainol1,2, Kamaliah Abdul Samad3, Che Asilah Ilyana Che Jazlan1, Nurul Ain Razahazizi1.   

Abstract

The goal of this study is to optimize the condition of the pollutant removal process by using acclimatized mixed culture (AMC) in the treatment of contaminated waste from the hatchery industry. The removal of chemical oxygen demand (COD), nitrate-N, and total phosphorus was optimized using a central composite design and the Response Surface Methodology (RSM) with two parameters: AMC content and retention time (days). Each factor had a range value of 15%-35% AMC content and a retention time of 3-5 days, with COD, nitrate-N, and total phosphorus removal as responses. Prior to experimentation, the synthetic wastewater was prepared, and the mixed cultures were acclimatized. In 13 runs, the experiment was carried out in accordance with the setup created by the Design-Expert software. The sample was tested for COD, nitrate-N, and total phosphorus using a Hach spectrophotometer. The findings show a strong relationship between predicted and experimental COD, nitrate-N, and total phosphorus removal values. At optimum conditions of 29% AMC content and 4 days of retention time, removal of COD, nitrate-N, and total phosphorus was observed to be 28%, 80% and 36%, respectively. The discovery also revealed that maximum values of removal of 62% COD, 94% nitrate-N, and 46% total phosphorus could be obtained under various optimum conditions. The study shows that, the acclimatized mixed culture (AMC) can be used as a potential biological wastewater treatment as well as a natural removal of COD, nitrate-N, and total phosphorus.
© 2022 The Authors. Published by Elsevier Ltd.

Entities:  

Keywords:  Acclimatized mixed culture; Chemical oxygen demand; Hatchery wastewater; Nitrate-N; Phosphorus

Year:  2022        PMID: 35368883      PMCID: PMC8971312          DOI: 10.1016/j.heliyon.2022.e09217

Source DB:  PubMed          Journal:  Heliyon        ISSN: 2405-8440


Introduction

Hatcheries are critical components of the aquaculture system because they are designed and operated specifically for the purpose of raising larvae economically [1]. Yet there are many challenges faced by the industry around the world which among that is finding a way to provide a hatchery environment that is best for larval growth without negatively impacting production economics [2]. The poor quality of water is the most common stressor present in the hatchery [2]. The high amounts of wastewater containing compounds including total nitrogen and total phosphorus were the main obstacles for hatchery wastewater [3]. The main source of nitrogen and phosphorus is derived from feed application, where a large fraction of it will remain as uneaten food and metabolic waste in the pond [4]. Shrimps are one of the most fundamental commodities of the global fishery trade [5] that has experienced significant expansion and contributes significantly to worldwide aquaculture productivity. However, the rise of shrimp aquaculture has had unforeseen repercussions, such as wastewater management challenges and other concerns related to the environmental impact of the effluent. Most of the nutrients used in shrimp aquaculture are obtained from outside sources such as on-farm formulated feeds and commercially available fertilizers. Shrimp feed that is not consumed or digested accumulates at the bottom of the pond and is destroyed by microbes [6]. The biological oxygen demand (BOD) of the pond water is increased as a result of this activity. To increase the primary productivity of shrimp ponds, fertilizers are added [7]. Yet, shrimp feeds and fertilizers have higher quantities of nitrogen and phosphorus than the water used in the cultivation process. Nitrogen and phosphorus are the primary constituents that cause eutrophication in aquaculture wastewater, which leads to the destruction of the ecosystem. Therefore, proper effluent treatment prior to discharge into water bodies is critical, because the sustainability of shrimp hatcheries and appropriate handling are dependent on the availability of a clean water source in continuity [8]. The need to abide by the current National Pollution Discharge Elimination System (NPDES) authorising measures has inducted the thought of substitutive treatment techniques for hatchery wastewater. The potential pollution problems for regular hatchery effluent include oxygen depletion, nutrient enrichment, taste and odour in cases where the receiving water flow is low. Conventional wastewater treatment used in aqua culturing systems relies on a variety of biological processes [3]. In the hatchery industry, mixed culture is one of the methods that are used in biological treatment to treat the wastewater. The advantages of applying mixed culture include reducing the levels of ammonia, nitrite, nitrate, and phosphorus, increasing the concentrations of dissolved oxygen, and promoting the decomposition of organic matter [9]. Additionally, mixed culture is capable of absorbing soluble organic compounds directly, which consequently contributed to the removal of chemical oxygen demand removal (COD) [10]. Nutrients may be removed chemically, physically, or biologically in conventional wastewater treatment. Of all the methods, the biological method seems to be the most efficient due to its cost-efficiency, as less energy consumption and maintenance requirements are needed [11]. The application of microbials in wastewater treatment has been widely used. The anaerobic-anoxic-aerobic method (A2O method) is one of the most effective methods for removing nitrogen and phosphorus from wastewater, which the reaction tank is made up of three tanks: anaerobic, anoxic, and aerobic. Many previous studies used the method to remove nitrogen and phosphorus from domestic wastewater [12, 13] and slaughterhouse wastewater [14]. However, a large amount of space is required to construct the system, making it unsuitable for small scale studies. Therefore, using microbes that have been treated directly in wastewater makes this study more feasible without requiring a large amount of experimental space. The use of microbes to treat hatchery wastewater has received relatively little attention, particularly in the application of mixed culture. Many previous studies have highlighted the numerous benefits of mixed culture over pure culture. Priyadharshini and Kumar [15] discovered that the removal efficiencies for COD mixed culture were 63% compared to 57% for COD removal in pure Bacillus sp. culture. While, a study conducted by Tangahu et al. [16] on Batik wastewater treatment using the intermittent method indicated that mixed culture of Scirpus grossus and Iris pseudacorus removed up to 89% and 97% of COD and BOD, respectively, as compared to a single culture. In other studies, Masser et al. [17] showed that a mixed culture of four species of bacteria removed 69% of total phosphorus from rubber wastewater, which was much greater than when each bacterium was treated alone. Hence, this study attempts to highlight the application of mixed cultures instead of pure culture for the treatment of hatchery wastewater. The goals of this research were to characterize the wastewater and to identify the optimum condition of the biological treatment of contaminated wastewater using acclimatized mixed culture. The optimum condition obtained has the potential to significantly remove nutrients and COD from wastewater by optimizing it using a central composite design (CCD).

Material and methods

Collection of samples

All samples were taken from a pond located in Setiu, Terengganu, Malaysia. The pond water was collected in order to create synthetic wastewater with a similar composition to the pond wastewater, and the pond sediment was collected to serve as the mixed culture [18]. Store-bought stones with diameters ranging from 1.5 to 2.0 cm were chosen.

Preparation of synthetic wastewater

Synthetic wastewater was formulated in accordance with hatchery wastewater characterization to achieve a COD concentration of 250 mg/L, a nitrate-N concentration of 10 mg/L, and a total phosphorus concentration of 10 mg/L. It was critical to generate synthetic wastewater with precise COD, nitrate-N, and phosphorus concentrations in order to standardize their contents throughout the studies and minimize experimental error. To prepare 10 L synthetic wastewater, 7g of fish food pellet solution (CNS) was mixed with 10 L tap water and autoclaved at 121 °C for 15 min. 250 mL of prepared synthetic wastewater was used to feed mixed cultures daily for 14 days.

Acclimatization of mixed cultures

Acclimatization was accomplished to ensure that the COD content in the sample could be adequately adapted by the mixed culture. A bioreactor was filled with 1250 mL of mixed culture and 3750 mL of nutrient stock [16]. Next, the total solids and COD concentrations were identified, which were 29294 mg/L and 285.67 mg/L, respectively. This procedure was carried out by administering 250 mL of nutritious stock daily for two weeks.

Experimental setup for optimization

Table 1 shows the factors and ranges of the CCD experimental setup. The experiments were developed using the central composite design (CCD) to determine the optimum process parameters of two factors, namely acclimatized mixed culture (AMC) content and retention time, using Response Surface Methodology (RSM). The experimental setup created by Design-Expert software is shown in Table 2. The experiment was conducted in a conical flask accordingly, with a total of 13 experimental runs. Conical flasks with varying concentrations of AMC and synthetic wastewater were filled to a total volume of 200 mL. Samples were kept at room temperature at various retention times. Samples were collected. Chemical oxygen demand (COD), nitrate-N, and total phosphorus levels were determined in samples. The outputs were inserted into the software and analysed with a 95% confidence level using analysis of variance (ANOVA).
Table 1

The range of factor values used in central composite design.

Factor-10+1
A: Acclimatized Mixed Culture (AMC) content1520253035
B: Retention time (days)33.544.55
Table 2

Experimental data of CCD.

StdFactors
Responses
ABCODNitrate-NTotal P
1203.5026980
2303.50146936
3204.50305617
4304.5050416
5154.0024620
6354.0038.8945
7253.0019980
8255.0062621.6
9254.00287440
10254.00277546
11254.0031.57842
12254.00307235
13254.0031.57738

A: Acclimatized mixed culture (AMC); B: Retention time.

The range of factor values used in central composite design. Experimental data of CCD. A: Acclimatized mixed culture (AMC); B: Retention time.

Sample analysis

The analysis of the sample was carried out immediately after collection in accordance with the retention time (Table 2). Chemical oxygen demand (COD), nitrate-N, and total phosphorus were the three responses that were determined.

COD analysis

COD test was performed to identify the initial and final readings used to calculate the removal of COD from the wastewater. In order to calculate the COD, 0th-day sample was used to detect the initial value of COD for all 10 containers in the 10th day. The COD vial was mixed with 2 mL of the 0th-day sample. Blank was prepared by mixing 2 mL of deionized water into the COD vial. Both the sample and blank vials were heated in a DRB200 reactor for 120 min at 150 °C before being tested. The test for the blank and sample was repeated three times to get an average reading using a DR900 Hach spectrophotometer.

Nitrate-N analysis

To calculate nitrate-N, the sample and the blank were prepared by adding 10 ml of the sample. One NitraVer 5 nitrate reagent pillow was added into the sample bottle. The DR900 Hach spectrophotometer was used to test both the sample and the blank.

Total phosphorus analysis

Molybdovanadate with Acid Persulfate Digestion Method (Method 10127) was used to measure total phosphorus. Deionized water was added to the wastewater and acclimatized mixed culture (AMC) mixtures to dilute to dilution factors of five and ten, respectively. The Hach spectrophotometer was then used to determine the total removal of phosphorus.

Results and discussions

Fitting the model

Table 2 displays the experimental results where, among the center points of 25% AMC and 4 days of retention time, run 11 exhibits the highest percentage removal for COD and nitrate-N. Meanwhile, run 10 had the greatest total phosphorus removal. It is also observed that different parameters contribute to the removal of pollutants differently. The highest removal of COD was observed at 62%, 98% for nitrate-N and 46% for total phosphorus. The analysis of variance demonstrates that COD and total phosphorus are adequately represented by a quadratic polynomial model after being fitted by different models (linear, two-factorial, quadratic and cubic). According to the analyses, the cubic polynomial model was determined to be the best for nitrate-N. The adjusted R2 of the quadratic model for COD and total phosphorus (0.9577 and 0.8907, respectively) was the highest among the linear, two factorial and cubic models. In the meantime, the cubic model had the highest adjusted R2 when compared to the linear and two factorial models for nitrate-N (0.7366). The quadratic polynomial model was insignificant for nitrate-N.

Statistical analysis for COD

Analysis of variance (ANOVA)

The ANOVA (F-test) and p-value of COD, which model coefficient is being assessed, and the parameters are examined for significance, are described in Table 3. It is also used to indicate the interaction strength between factors. The confidence level portrayed by the ANOVA table was considerably higher than 95%, while the p-value was much lower than 0.0001. A p-value less than 0.05 indicated that the model was significant, indicating that it was appropriate for this experiment. Both main effects and their interactions were significant based on the p-value of less than 0.05. Meanwhile, the “lack of fit” of this model is not significant with a p-value of 0.2053, indicating the quadratic model is adequate. The R2 of 0.9753 and adjusted R2 of 0.9577 in which higher than 0.80 signified that the proposed model represents the experimental data well. As mentioned by Lee et al. [19], the R2 for a good fit of a model should be at least 0.80, hence this regression model described the process effectively. The predicted against actual data on COD removal is presented in Figure 1a. A linear distribution was seen, indicating that the model fit the data well. The predicted values were in close agreement to the observed values for COD removal. For COD removal, the predicted values were in close agreement to the actual values.
Table 3

ANOVA for COD removal response surface model.

SourceSum of squaresMean squareF valueP-valueProb > FLevel of significance
Model1864.23372.8555.35<0.0001Significant
A-AMC Content117.81117.8117.490.0041
B-Retention time1323.001323.00196.42<0.0001
AB256.00256.0038.010.0005
A23.713.710.550.4824
B2164.23164.2324.380.0017
Residual47.156.74
Lack of Fit30.4510.152.430.2053Not significant
Cor Total1911.38
R20.9753
Radj0.9577
Figure 1

Predicted against actual plot of (a) COD, (b) nitrate-N, and (c) total phosphorus removal.

ANOVA for COD removal response surface model. Predicted against actual plot of (a) COD, (b) nitrate-N, and (c) total phosphorus removal.

Effect of AMC content and retention time on COD removal

Figure 2 depicts the influence of AMC content and retention time on COD removal. Figure 2a shows that increasing the AMC content from 20% to 30% resulted in an increase in COD removal. Similarly, increasing the retention time from 3.5 days to 4.5 days resulted in an increase in COD removal, as shown in Figure 2b. Oljira et al. [20] had demonstrated that the removal of COD in the brewery effluents increases with the increase in incubation period. Moreover, instead of a single culture application, the combination of three different bacteria species showed a remarkable increment in the removal of COD. The treatment with the consortium of Aeromonas sp., Pseudomonas sp., and Bacillus sp. gave rise to a 93.25% removal efficiency of COD, as opposed to the individual treatments of 76.78%, 79.61%, and 82.67% removal, respectively. The organic material presents in the wastewater operates as a substrate for the aerobic microbial metabolism, which explains the reduction of COD concentration. The removal efficiency increased as the incubation period lengthened. Their results were consistent with this study, which found that the rate of COD removal increased over time and that the bacterial content is critical to the process of improving removal.
Figure 2

Effect of AMC content on, % (a) COD, (c) nitrate-N, and (e) total phosphorus removal, and effect of retention time, days on (b) COD, (d) nitrate-N, and (f) total phosphorus removal.

Effect of AMC content on, % (a) COD, (c) nitrate-N, and (e) total phosphorus removal, and effect of retention time, days on (b) COD, (d) nitrate-N, and (f) total phosphorus removal. The influence of AMC content and retention time on COD removal is shown in Figure 3. The plot demonstrated that the retention time had an effect on COD removal for any AMC content. The amount of COD removal increased as the retention time was increased.
Figure 3

Contour plots of (a) COD, (c) nitrate-N, and (e) total phosphorus removal, and 3D response surface plots of (b) COD, (d) nitrate-N, and (f) total phosphorus removal as function of AMC content and retention time.

Contour plots of (a) COD, (c) nitrate-N, and (e) total phosphorus removal, and 3D response surface plots of (b) COD, (d) nitrate-N, and (f) total phosphorus removal as function of AMC content and retention time.

Statistical analysis for Nitrate-N

The ANOVA for nitrate-N was presented in Table 4. It was found that the model had a p-value of 0.0354, which was less than 0.05, indicating that it was statistically significant. The p-value for both main effects was lower than 0.05 and therefore was significance; meanwhile, the interaction effect for nitrate-N was insignificant, proved by its p-value of 0.4485, which was higher than 0.05. Meanwhile, the “lack of fit” of this model is significant with a p-value of 0.0015, indicating the quadratic model is not adequate. The R2 value of 0.8903 explains that the estimated model suits the experimental data adequately. The predicted versus actual plot on nitrate-N removal was presented in Figure 1b. The linear distribution showed that the model was fitted well, and that the values of predicted and observed were in good agreement.
Table 4

ANOVA for nitrate-N removal.

SourceSum of squaresMean squareF valueP valueProb > FLevel of significance
Model2942.54420.365.800.0354Significant
A-AMC Content512.00512.007.060.0451
B-Retention time648.00648.008.930.0305
AB49.0049.000.680.4485
A25.935.930.0820.7864
B223.3123.310.320.5953
Residual362.6972.54
Lack of Fit339.89339.8959.630.0015Significant
Cor Total3305.23
R20.8903
Radj0.7366
ANOVA for nitrate-N removal.

Effect of independent and interactive parameters on nitrate-N removal

Figure 2 depicts the influence of two independent parameters on nitrate-N removal. As shown in Figure 2c, the removal of nitrate-N was increased with the increase in AMC content, while an increase in retention time caused the decrease in nitrate-N removal, as depicted in Figure 2d. These suggest that three days were sufficient for the microbes to remove nitrate-N, as the maximum removal of nitrate-N was observed on day 3. Nitrate was required for the growth of the mixed culture. Increased retention time could promote the growth of bacteria that digest all impurities, allowing the wastewater to be cleansed [21]. Nitrate removal via the denitrification process requires an ample amount of oxygen [22, 23]. Therefore, oxygen was reduced as a result of the longer retention time, which eventually diminished the removal of nitrate. Lu et al. [24] in their research found that 0.7 h was the optimum hydraulic retention time instead of 12 h for the removal of nitrogen. Furthermore, the nitrite oxidation activity was significantly increased with a decreased hydraulic retention time from 30 to 5 h. This is due to the increase in nitrite-oxidizing bacteria, where the short hydraulic retention time favoured the relative growth of these bacteria, particularly the fast-growing Nitrobacter sp [25]. As pointed out by Oljira et al. [20], the mixed combination of Aeromonas sp., Pseudomonas sp., and Bacillus sp., had shown a significant contribution to the removal of total nitrogen, as compared to their single treatment. The synergistic effect on the removal of the pollutants by these three isolates might explain this condition. The contour plot and 3D response surface were presented in Figures 3c and 3d, representing the effect of interactions between parameters on nitrate-N removal. The optimum removal of nitrate-N was up to 98% at a condition of 20% AMC content and 3.5 days of retention time. An increase in retention time was seen to unnecessarily affect the removal of nitrate-N. Also, it was remarked that the removal of nitrate-N at any AMC content was affected by the retention time.

Statistical analysis for total phosphorus

Table 5 presents the summary of ANOVA for total phosphorus removal. The p-value of less than 0.05, which is 0.0005 was an indication of the significance of the model. Both the main effects of AMC content and retention time displayed an insignificant effect with p-values of 0.1518 and 0.6660, respectively. However, the interaction effect was statistically significant at 0.0073. Meanwhile, the “lack of fit” of this model is not significant with a p-value of 0.1062, indicating the quadratic model is adequate. Both R2 and adjusted R2 of 0.9362 and 0.8907, respectively, demonstrated the estimated model adequately fits the experimental data. Figure 1c portrays the predicted against the actual plot for total phosphorus removal. Again, as seen in all responses, the linear distribution was observed, indicating that the model was well-fitted. The values of predicted and observed were in close agreement.
Table 5

ANOVA for total phosphorus removal.

SourceSum of squaresMean squareF valueP-valueProb > FLevel of significance
Model4054.03810.8120.550.0005Significant
A-AMC Content102.08102.082.590.1518
B-Retention time8.008.000.200.6660
AB552.25552.2514.000.0073
A22081.672081.6752.760.0002
B22271.452271.4557.570.0001
Residual276.1839.45
Lack of Fit207.3869.134.020.1062Not significant
Cor Total4330.21
R20.9362
Radj0.8907
ANOVA for total phosphorus removal.

Effect of AMC content and retention time on total phosphorus removal

The influence of AMC content and retention time on total phosphorus is depicted in Figure 2. There was an increase in the removal of phosphorus with the increase in AMC content until a maximum value of 46% at the centre point and was decreased after reaching the maximum value (Figure 2e). A comparable patent was also seen for retention time, represented by Figure 2f where the removal of total phosphorus was increased with the increase of retention time and achieved its maximum value at the centre point. Naili et al. [26] mentioned that the mixed culture of Alcaligenes denitrificans, Moraxella lacunata, Pseudomonas aeruginosa, and Acinetobacter junii had contributed significantly to the removal of phosphate in the wastewater and that the synergistic relationship between the co-cultured species had accelerated the removal process. The influence of interaction between AMC content and retention time was graphically illustrated by Figures 3e and 3f using contour plot and 3D response, respectively. The presence of a significant interaction between the parameters was indicated by an elliptical contour plot. It was noticed from the figure that the optimum removal of total phosphorus was obtained at 25% AMC content and at 4 days of retention time.

Optimum condition

Table 6 displays the optimum conditions suggested by the Design-Expert software obtained from the analysis of RSM. It was suggested that at AMC content of 29% with a retention time of 4 days, as much as 28%, 80%, and 36% of COD, nitrate-N and total phosphorus removal could be obtained. At 25% of AMC content and 4 days of retention time, the highest COD, nitrate-N and total phosphorus measured were 31.5%, 78% and 42% respectively. However, to achieve maximum removal for each response, the conditions were seen varied among responses. The removal of COD was up to 62% in a condition of 25% AMC content and 5 days of retention time. For nitrate-N, at 25% AMC content and 3 days of retention time had given rise to the removal value of up to 94%. To reach the maximal removal of 46% for total phosphorus, the optimal condition was at 25% of AMC content and 4 days of retention time. All responses demonstrated the same AMC content of 25% for the maximum removal with varied retention times.
Table 6

Suggested optimum conditions.

FactorValue
A: AMC content29 %
B: Retention time4 days
COD removal28 %
Nitrate-N removal80 %
Total phosphorus removal36 %
Suggested optimum conditions.

Conclusion

This study focused on the application of acclimatized mixed culture (AMC) in removing chemical oxygen demand (COD), nitrate-N and total phosphorus. The conditions of 29% AMC content and a retention time of 4 days were optimum for the removal process as determined by central composite design, with a maximum removal of 31.5% COD, 78% nitrate-N, and 46% total phosphorus obtained. While separating conditions to achieve the maximum removal value resulted in 62%, 94%, and 46% removal of COD, nitrate-N, and total phosphorus, respectively. The results demonstrate that acclimatized mixed culture can naturally remove COD, nitrate-N and total phosphorus.

Declarations

Author contribution statement

Norazwina Zainol: Conceived and designed the experiments; Contributed reagents, materials, analysis tools or data. Kamaliah Abdul Samad: Analyzed and interpreted the data; Wrote the paper. Che Asilah Ilyana Che Jazlan: Conceived and designed the experiments; Performed the experiments; Analyzed and interpreted the data; Wrote the paper. Nurul Ain Razahazizi: Conceived and designed the experiments; Performed the experiments.

Funding statement

This work was supported by Universiti Malaysia Pahang (RDU210341).

Data availability statement

The data that has been used is confidential.

Declaration of interests statement

The authors declare no conflict of interest.

Additional information

No additional information is available for this paper.
  5 in total

Review 1.  Enhancing microalgal photosynthesis and productivity in wastewater treatment high rate algal ponds for biofuel production.

Authors:  Donna L Sutherland; Clive Howard-Williams; Matthew H Turnbull; Paul A Broady; Rupert J Craggs
Journal:  Bioresour Technol       Date:  2014-10-24       Impact factor: 9.642

2.  Denitrifying phosphorus removal and impact of nitrite accumulation on phosphorus removal in a continuous anaerobic-anoxic-aerobic (A2O) process treating domestic wastewater.

Authors:  Wei Zeng; Lei Li; Ying-ying Yang; Xiang-dong Wang; Yong-zhen Peng
Journal:  Enzyme Microb Technol       Date:  2010-11-04       Impact factor: 3.493

3.  Nitrogen removal from poultry slaughterhouse wastewater in anaerobic-anoxic-aerobic combined reactor: Integrated effect of recirculation rate and hydraulic retention time.

Authors:  Carla Limberger Lopes; Tatiane Martins de Assis; Fernando Hermes Passig; Adriana Neres de Lima Model; Juliana Bortoli Rodrigues Mees; Francisco J Cervantes; Jackeline T Gotardo; Simone Damasceno Gomes
Journal:  J Environ Manage       Date:  2021-11-30       Impact factor: 6.789

4.  The effect of hydraulic retention time on ammonia and nitrate bio-removal over nitrite process.

Authors:  Yaofeng Lu; Zhibin Ding; Kexin Gao; Zijun Duo; Shaoyi Xu; Xiaolan Zhao; Wei Li; Ping Zheng
Journal:  Environ Technol       Date:  2018-10-12       Impact factor: 3.247

5.  Potential Applications of Some Indigenous Bacteria Isolated from Polluted Areas in the Treatment of Brewery Effluents.

Authors:  Temesgen Oljira; Diriba Muleta; Mulissa Jida
Journal:  Biotechnol Res Int       Date:  2018-02-11
  5 in total

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