Literature DB >> 26190684

Peripheral refraction and higher-order aberrations with cycloplegia and fogging lenses using the BHVI-EyeMapper.

Ravi Chandra Bakaraju1, Cathleen Fedtke2, Klaus Ehrmann3, Darrin Falk2, Varghese Thomas2, Brien Anthony Holden3.   

Abstract

PURPOSE: To determine if a fogging lens ameliorates accommodative effects driven by the closed-view design of the BHVI-EyeMapper (EM) instrument. We compared cycloplegic refraction and higher-order aberration measurements of the EM with those obtained with a fogging lens.
METHODS: Twenty-six, young, participants (15F, 25±5 years, range: 18-35 years, SE: +0.25 D and -3.50 D) with good ocular health were recruited. Five independent measurements of on- and off-axis refraction and higher-order aberrations were recorded across the horizontal visual field, under two conditions: non-cycloplegic measurements with +1.00 D fogging lens and cycloplegia, always in the same sequence. The contralateral eye was occluded during the measurements. Two drops of 1% Tropicamide delivered within 5 min facilitated cycloplegic measurements. All participants were refracted 30 min after installation of the second drop.
RESULTS: Mean spherical equivalent measures of the non-cycloplegic condition were significantly more myopic than their cycloplegic counterparts (p<0.05); approximately by 0.50 D centrally, increasing to 1.00 D towards the periphery. The horizontal astigmatic component, J180, demonstrated small but statistically significant differences between the test conditions. Differences were predominant for eccentricities greater than 30°, in both nasal and temporal meridians. The oblique astigmatic component, J45, was not significantly different between the test conditions. The primary spherical aberration coefficient C(4, 0) was significantly less positive for the non-cycloplegic state than its cycloplegic counterpart. This result held true across the entire horizontal visual field. The horizontal coma and trefoil coefficients C(3, 1) and C(3, 3) were not significantly different between the two conditions.
CONCLUSIONS: The use of +1.00 D fogging lens without cycloplegia did not provide complete relaxation of accommodation. The discrepancies between cycloplegic and non-cycloplegic EM measurements were found to be more pronounced for peripheral field angles than central measures, for both M and J180 components.
Copyright © 2015 Spanish General Council of Optometry. Published by Elsevier Espana. All rights reserved.

Entities:  

Keywords:  Accommodation; Acomodación; Peripheral refraction and higher-order aberrations; Refracción periférica y aberraciones de alto orden

Mesh:

Substances:

Year:  2015        PMID: 26190684      PMCID: PMC4705317          DOI: 10.1016/j.optom.2015.06.003

Source DB:  PubMed          Journal:  J Optom        ISSN: 1989-1342


Introduction

The actual mechanism of myopia genesis and its rate of progression are not clearly understood. Nevertheless, accumulated evidence from animal data suggests an involvement of the retinal image shell, as guided by peripheral refraction and off-axis higher-order aberrations (HOA).1, 2, 3 There is some clinical evidence in humans supporting the role of peripheral refraction in myopia progression.4, 5, 6, 7, 8, 9, 10 The emerging paradigm of the contribution of peripheral refraction to myopia development has triggered considerable interest in the peripheral optics of the eye.11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25 In the majority of these studies conventional, commercial instruments or procedures, have been modified to facilitate peripheral optics measurement, requiring repeated turning of the head or eye and tedious instrument re-alignment. For this reason, techniques that facilitate expedited measurements of peripheral refraction and/or aberrations have gained much attention in the last few years.26, 27, 28 In response to the desire for automated peripheral refraction, the Brien Holden Vision Institute, Sydney, Australia developed the BHVI-EyeMapper (EM), a closed-view global aberrometer (Fig. 1). The EM performs fast (<0.5 s) refraction scans across 100° of the visual field in 10° steps, along horizontal, vertical and oblique visual field meridians, without the need for eye or head turn. Complete details of the working principles of the instrument can be found elsewhere. The EM has been previously validated against a conventional aberrometer (COAS-HD®, Wavefront Sciences, USA) in conjunction with Dynamic Stimulation Aberrometry (Optana®, Germany) and an open field autorefractor (Shin-Nippon NVision K5001®, Grand Seiko, Japan) for refraction measures obtained at distance and various near vergences (Bakaraju RC, et al. IOVS. 2012;53:ARVO E-Abstract 1354). Overall, the EM being influenced by its closed-view instrument design produces slightly ‘more negative’ results, about −0.50 D more myopic for distance refraction measures than those obtained with the Shin-Nippon autorefractor. However, results were in modest agreement with those obtained with COAS-HD aberrometer, another closed-view on-axis aberrometer. In the present study, we determined trends and relative differences in the on- and off axis refraction and HOA obtained using a fogging lens versus cycloplegia. Finally we sought to determine if, driven by its closed view instrument design, the fogging lens ameliorates accommodative effects on the refraction measurements.
Figure 1

(a) The BHVI-EyeMapper (EM), Brien Holden Vision Institute, Sydney, Australia and (b) three-dimensional view describing the instrument design.

Methods

All the research procedures conformed to the regulations guided by the declaration of Helsinki. A local ethics committee (Bellberry, Adelaide, South Australia) approved the research protocol and the study was registered with the Australia New Zealand Clinical Trials Registry (#ACTRN12611001103954). Participants were recruited in response to an advertisement circulated via email to staff and eligible candidates on the patient database of the Brien Holden Vision Institute, Sydney, Australia. All participants were screened for suitability using routine ocular examination, which included auto-refraction (Nidek Tonoref II, Japan), subjective refraction, visual acuity and general ocular health examination with a slit-lamp biomicroscope and direct ophthalmoscope. A total of 26 ocular healthy participants (15 female) in the age range of 18–35 years (25 ± 5 years) with best-corrected visual acuity ≥20/20 in each test eye and a spherical equivalent refractive error between +0.25 D and −3.50 D were recruited for the study. The EM is a closed-view aberrometer, designed to rapidly measure refractive errors and HOA of the eye over a wide range of viewing angles. The special feature of the EM is the use of numerous fixed optical components to steer the illumination and reflection beam to peripheral angles. A deflection system consisting of 33 mirrors and one scanning mirror was arranged to maintain equal optical path length while providing rapid scanning of the illumination beam across an angular range of up to ±50°, measuring peripheral refractive errors and HOA in discrete 10° steps along the selected meridian, including horizontal, vertical and some oblique meridians. One measurement scan (i.e. 11 positions) takes 0.45 s. The EM's working distance is 100 mm. The instrument is equipped with an internal back-illuminated fixation target mounted on a translation stage to facilitate accommodative response measurements. The software adjustable linear stage can be manoeuvred on a continuous scale to create a range of object vergences from +1.00 D (fogging) to −5.00 D (accommodative target). Five independent wavefront measurements on the unaided eyes of each participant facilitated on- and off-axis, refraction and HOA, across the horizontal visual field under two test conditions. The measurements were always recorded in the following sequence: (1) no-cycloplegia with +1.00 D fogging; and (2) under cycloplegia with Tropicamide 1%. The fellow eye was always occluded with an eye patch. Two drops of 1% w/v Tropicamide within 5 min facilitated all cycloplegic measurements. For the cycloplegic test case, all participants were measured on the EM 30 min after the installation of the last drop. Under cycloplegia, visual acuity, auto-refraction (Nidek, Japan) and subjective refraction were assessed again. For all EM measurements, the illumination of the testing room was kept below 10 lux, to ensure the wavefront measures under ‘no-cycloplegia’ were greater than 4 mm pupil diameter. A Maltese cross with 8-arms designed to be telecentric in object space served as the fixation target. Each wavefront measurement provided a Hartmann–Shack image for further post-processing, the commercial camera software package HASO (Imagine Optics, France) was used to derive the spot locations. Once the spot locations were obtained, custom mathematical routines written in CVI/Lab windows were deployed to obtain slopes and fit them to Zernike polynomials derivatives. Refraction measures were derived as the Fourier vector components (M, J180 and J45) from the Zernike coefficients, as described elsewhere. In addition to the refraction measures, HOA up to 4th order in OSA standard were also obtained. The refraction and HOA measurements were obtained at a 4 mm circular pupil diameter. In order to allow direct comparison of all data, i.e. cyclo versus non-cyclo, we did not consider elliptical pupils. Statistical analyses were performed using SPSS (version 18.0, Chicago, USA). Paired t-test was used to determine significant differences between cycloplegic and non-cycloplegic conditions at each visual field angle. Additional pairwise comparisons were performed between on-axis, cycloplegic and non-cycloplegic autorefraction obtained from EM, Nidek autorefractor and subjective refraction. Assuming the within participant SD to be 0.30 D, a sample size of 24 subjects was needed to detect a paired difference of 0.18 ± 0.3 D between the test conditions with alpha (α) set to 0.05 and statistical power to 80%. The sample size did not accommodate for dropouts, as the study needed only one participant visit.

Results

The means of the selected power vectors and aberration coefficients are presented as a function of retinal eccentricity (−50° to 50°) for both the test conditions. The sign convention used throughout the manuscript for the field angles is positive for nasal and negative for temporal fields. In the interest of brevity, only statistically significant and/or relevant HOA with respect to horizontal meridian are represented in figures and tables. Fig. 2 showcases the differences in the measures of M and J180 as a function of horizontal visual field angle, for the two test conditions. At all field eccentricities, the mean spherical equivalent (M) measures obtained under the no-cycloplegia condition were significantly more myopic than under cycloplegia (p < 0.05). The differences were about 0.50 D centrally, increasing to about 1.00 D for peripheral field angles. Small (in the range of 0.25–0.50 D) yet statistically significant differences were found between the two test conditions, for J180 measures, particularly at peripheral field angles in the range of 30° to 40°. For the astigmatic component J45, there were no significant differences between the two conditions at any visual field angle.
Figure 2

Refractive components M (spherical equivalent in diopters) as a function of horizontal visual field angle (in degrees) for two test conditions: (a) non-cycloplegic state with +1.00 D fogging lens and (b) cycloplegic state. The shaded areas indicate 95% confidence intervals. The symbol ‘*’ indicate statistical significant differences (p < 0.05) at the indicated field angles. Sign convention used was postive for nasal and negative for temporal field angles.

Table 1 compares three Zernike coefficients: horizontal coma C(3, 1), horizontal trefoil C(3, 3) and primary spherical aberration C(4, 0), as a function of visual field eccentricity for the two test conditions. As noted, there were no statistically significant differences in the horizontal coma C(3, 1), which held true across most field eccentricities, except for nasal 10° and 20° field angles. In general, small and statistically insignificant changes were found in horizontal trefoil, with exceptions at 40° and 50° nasal field eccentricities (p < 0.05). Primary spherical aberration C(4, 0) was found to be significantly less positive in magnitude for the non-cycloplegic state than the corresponding measures obtained under cycloplegia. This relationship held true across the entire horizontal visual field, except at the 30° nasal and temporal visual field angles, for which the differences became insignificant (p > 0.05). There were no statistically significant differences in the horizontal secondary astigmatism and horizontal tetrafoil across the field eccentricities (p > 0.05), except for temporal 30° for the former and temporal 30° and central field angles for the latter. There were no significant differences found between the two test conditions (p > 0.05), for the following on-axis variables (central field measures): oblique astigmatic component (J45), vertical coma C(3, −1), vertical trefoil C(3, −3), vertical tetrafoil C(4, −4) and vertical secondary astigmatism C(4, −2).
Table 1

Aberration coefficients C(3, 1), C(3, 3) and C(4, 0) for 4 mm pupil diameter (in microns) as a function of horizontal visual field angle (in degrees) for cycloplegia and no-cycloplegia conditions. Sign convention for the field angles: positive for nasal and negative for temporal fields.

VariablesField angleNCycloplegia
No-cycloplegia
p-Value
MeanSDMeanSD
C(3, 1)Horizontal coma−50260.4070.1900.3810.1640.082
−40260.2500.1000.2460.0960.514
−30260.1540.0630.1530.0670.963
−20260.1050.0490.1170.1000.560
−10260.0460.0370.0430.0360.466
0260.0000.028−0.0040.0290.282
1026−0.0470.039−0.0560.0360.006
2026−0.0890.037−0.1020.0390.003
3026−0.1450.060−0.1580.0650.069
4026−0.2610.120−0.2540.1220.371
5026−0.6450.262−0.6010.2440.051



C(3, 3)Horizontal trefoil−50260.1180.1090.1020.0800.130
−40260.0430.0700.0410.0540.747
−30260.0250.0470.0270.0460.681
−20260.0130.0420.0420.1210.286
−10260.0070.0380.0110.0390.487
0260.0010.0360.0050.0360.334
1026−0.0080.044−0.0070.0450.755
2026−0.0110.054−0.0160.0560.392
30260.0030.0700.0030.0700.979
4026−0.0360.096−0.0110.0840.001
5026−0.2240.194−0.1630.1590.000



C(4, 0)Spherical aberration−5026−0.0460.045−0.0580.0480.083
−4026−0.0030.024−0.0140.0300.011
−3026−0.0060.021−0.0110.0200.208
−20260.0060.021−0.0140.0630.096
−10260.0160.0240.0080.0230.063
0260.0240.0220.0170.0240.002
10260.0240.0220.0180.0220.018
20260.0150.0210.0070.0280.013
3026−0.0030.033−0.0080.0370.244
4026−0.0290.047−0.0400.0460.095
5026−0.0720.080−0.1100.0650.025

All statistically significant p-values (i.e. p < 0.05) are italicized.

For the on-axis cycloplegic refraction, the Nidek autorefractor produced statistically significant lower values than the EM counterparts, for M (paired Δ = −0.13 ± 0.17 D, CI: −0.20 D to −0.06 D, p = 0.001) and J180 components (paired Δ = 0.05 ± 0.12 D, CI: 0.00–0.10 D, p = 0.047). However, all differences were clinically insignificant. The subjective refraction measures obtained under cycloplegia were also significantly lower than the corresponding EM measures, for M (paired Δ = −0.09 ± 0.21 D, CI: −0.17 D to 0.01 D, p = 0.044) and J180 components (paired Δ = 0.10 ± 0.16 D, CI: 0.04–0.17 D, p = 0.004), but again, were clinically insignificant. The oblique astigmatic component, J45, was not different when cycloplegic EM measures were compared with Nidek autorefractor and subjective refraction. When the non-cycloplegic measurements were considered, the EM reported clinically relevant higher levels of myopia than the Nidek autorefractor (paired ΔM = −0.41 ± 0.49 D, CI: −0.61 D to −0.22 D, p < 0.001) and subjective refraction (paired ΔM = −0.41 ± 0.52 D, CI: −0.62 D to −0.20 D, p < 0.001). No differences were observed for J180 and J45 components between the non-cycloplegic EM, Nidek and subjective measurements (Table 2).
Table 2

The comparison of the mean on-axis refraction decomposed into three Fourier components (M, J180 and J45 in Diopters) obtained in the following conditions: (i) EM and Nidek autorefractor under complete cycloplegia; (ii) EM and Nidek autorefractor, under no-cycloplegia; (iii) EM and subjective refraction under cycloplegia and (iv) EM and subjective refraction under no-cycloplegia. All units are in Diopters.

VariablesField angleNEM (cycloplegia)
Nidek (cycloplegia)
p-Value
MeanSDMeanSD
M026−1.581.51−1.451.500.001
J1800260.140.370.090.330.047
J45026−0.040.17−0.030.190.693

Discussion

The goal of the current study was to identify if the measurements obtained by the EM are influenced by its closed-view design. Overall, in our sample of young adults, the mean spherical equivalent (M) measures in the non-cycloplegic condition were found to be significantly more myopic than their cycloplegic counterparts (p < 0.05). The average differences between the cycloplegic and non-cycloplegic states are in the order of half a diopter centrally, reaching up to one diopter, across the horizontal visual field. The results suggests that the use of a +1.00 D fogging lens during non-cycloplegic measurements did not provide sufficient relaxation of accommodation and certainly was also not comparable to the measures obtained for the cycloplegic state. Querios et al.,32, 33 have demonstrated that successful relaxation of accommodation can be attained by the use of +2.00 D fogging spectacle lenses, while performing open-field central and peripheral refraction measurements. However, in contrast, in the current myopic participant sample, using the EM, we could not reproduce the claimed successful relaxation effects of accommodation. We offer two reasons for this lack of agreement. Firstly, the magnitude of the fogging lens used in the current study was exactly half of that proposed by Querios et al. As mentioned briefly in “Methods” section, the prototype EM can facilitate a six diopter range of defocus via its internal, telecentric, Badal system (−5.00 D to +1.00 D). This range of defocus limited our choice of fogging to a maximum of +1.00 D, without disturbing the telecentricity of the optical system. We plan to rectify this limitation by replacing the current Badal system with that of a much larger dioptric range for future studies. Secondly, considering the closed-view design of the EM (as opposed to the measures reported with open-field autorefractor), the very proximity of this instrument potentially created a weak stimulus for accommodation. Nevertheless, it is important to note that Querios and colleagues reported in their sub-group analysis that the paired differences between the cycloplegic state and the non-cycloplegic state with fogging lenses was significantly different only in the myopic group. Furthermore, contrary to our results, they reported lower mean spherical equivalent myopic refractive error with use of fogging lenses than those achieved with cycloplegia, indicating that either cycloplegia was ineffective or the measurement data were affected by the use of a fogging lens itself. Carkeet et al., investigated the effect of cycloplegia on ocular aberrations in young adults using Technolas Zywave aberrometer (Bausch and Lomb, USA) and found small yet statistically significant differences between higher order measures (3rd and 5th order) between cycloplegic and non-cycloplegic states. Hiraoka et al., confirmed similar results after assessing the effect of cycloplegia on ocular aberrations in children using the Topcon Hartmann-Shack aberrometer (KR-9000PW, Topcon, Tokyo, Japan) and 1% Atropine. Our results, reported from a custom-built instrument, are in agreement with these previous reports. However, it is worth noting that unlike in Hiraoka et al., we used 1% Tropicamide, for a shorter duration of cycloplegic effect. Further, it is important to note that the previous studies examined the effects of cycloplegia on the on-axis HOA while the current study extends this to horizontal peripheral visual field, encompassing −50° to 50°. Although small, it is interesting to note that there were significant differences between the Nidek autorefractor and EM measures even under complete cycloplegia. Such discrepancies in measurements could be plausibly explained by the differences in the measurement techniques. The EM is based on Hartmann–shack aberrometry principles in which a complete circular region of interest is considered for obtaining refraction measurements, while the Nidek autorefractor analyses projected infrared rings that only consider an annular region of the pupil. Further, it is worth noting that the EM measures were computed using Seidel defocus (i.e. paraxial curvature matched measures that includes defocus and primary spherical aberration) unlike the Nidek. It can be argued that the use of a single power fogging lens (+1 D) for testing various degrees of myopic participants in the study could be a limitation, as the fogging lens could have resulted in a slightly different level of accommodation control for different magnitudes of myopia. We chose to use one identical test case scenario across the experiment to reduce the possible effects of other confounding variables. Overall, our results suggest that the refraction values from the EM obtained for non-cycloplegic conditions could be slightly more myopic than their cycloplegic counterparts, despite the use of a +1.00 D fogging lens while performing the measurement. This shift can principally attributed to instrument-driven proximal accommodation. While performing subjective refraction, using the objective EyeMapper output as a start point, it is important to bear in mind that the start point could be geared towards more myopia or less hyperopia and appropriate subjective fogging techniques are needed for an optimal refraction endpoint. The EM used here was a prototype. The next generation capable of providing an increased dioptric range to facilitate greater levels of fogging is now under development. Further work may be required to assess whether the use of stronger fogging lenses would ameliorate accommodation effects due to the closed-view nature of the instrument, in particular, in studies that involve even younger participants for whose refraction measurements, accommodation control is much more crucial.

Conclusion

The use of +1.00 D fogging lens did not provide complete relaxation of accommodation which observation was clearly demonstrated when compared with the measurements produced under cycloplegia. Data from older individuals in the study confirm that proximal accommodation due to the closed view design of the EM plays an important role while performing refraction measurements. Accommodation affects the peripheral refraction more than central refraction for both, M and J180. A stronger fogging lens may be needed in vision studies involving younger populations where accommodation control can be crucial.

Conflicts of interest

Some of the authors are inventors on the patents that relate to the BHVI-EyeMapper: C. Fedtke (AU2011902736), K. Ehrmann (WO 2008/116270 A1, AU2011902736), D. Falk (AU2011902736), B. Holden (WO 2008/116270 A1).
  32 in total

Review 1.  Interventions to slow progression of myopia in children.

Authors:  Jeffrey J Walline; Kristina Lindsley; Satyanarayana S Vedula; Susan A Cotter; Donald O Mutti; J Daniel Twelker
Journal:  Cochrane Database Syst Rev       Date:  2011-12-07

2.  Peripheral refraction in high myopia with spherical soft contact lenses.

Authors:  Eva Kwok; Bhavna Patel; Simon Backhouse; John R Phillips
Journal:  Optom Vis Sci       Date:  2012-03       Impact factor: 1.973

3.  Measuring ocular aberrations in the peripheral visual field using Hartmann-Shack aberrometry.

Authors:  David A Atchison; Dion H Scott; W Neil Charman
Journal:  J Opt Soc Am A Opt Image Sci Vis       Date:  2007-09       Impact factor: 2.129

4.  Peripheral refraction in myopic children wearing orthokeratology and gas-permeable lenses.

Authors:  Pauline Kang; Helen Swarbrick
Journal:  Optom Vis Sci       Date:  2011-04       Impact factor: 1.973

5.  Peripheral refraction and ocular shape in children.

Authors:  D O Mutti; R I Sholtz; N E Friedman; K Zadnik
Journal:  Invest Ophthalmol Vis Sci       Date:  2000-04       Impact factor: 4.799

6.  Myopia and peripheral ocular aberrations.

Authors:  Ankit Mathur; David A Atchison; W Neil Charman
Journal:  J Vis       Date:  2009-09-21       Impact factor: 2.240

7.  Refractive error, axial length, and relative peripheral refractive error before and after the onset of myopia.

Authors:  Donald O Mutti; John R Hayes; G Lynn Mitchell; Lisa A Jones; Melvin L Moeschberger; Susan A Cotter; Robert N Kleinstein; Ruth E Manny; J Daniel Twelker; Karla Zadnik
Journal:  Invest Ophthalmol Vis Sci       Date:  2007-06       Impact factor: 4.799

8.  Influences of cycloplegia with topical atropine on ocular higher-order aberrations.

Authors:  Takahiro Hiraoka; Kazunori Miyata; Yasuko Nakamura; Takashi Miyai; Miyuki Ogata; Fumiki Okamoto; Tetsuro Oshika
Journal:  Ophthalmology       Date:  2012-10-16       Impact factor: 12.079

9.  Higher order ocular aberrations after cycloplegic and non-cycloplegic pupil dilation.

Authors:  Andrew Carkeet; Shajan Velaedan; Yong Kang Tan; Derrick Yu Jiang Lee; Donald T H Tan
Journal:  J Refract Surg       Date:  2003 May-Jun       Impact factor: 3.573

10.  The BHVI-EyeMapper: peripheral refraction and aberration profiles.

Authors:  Cathleen Fedtke; Klaus Ehrmann; Darrin Falk; Ravi C Bakaraju; Brien A Holden
Journal:  Optom Vis Sci       Date:  2014-10       Impact factor: 1.973

View more
  5 in total

1.  Instrument for fast whole-field peripheral refraction in the human eye.

Authors:  Enrique J Fernandez; Santiago Sager; Zhenghua Lin; Jiangdong Hao; Javier Roca; Pedro M Prieto; Zhikuang Yang; Weizhong Lan; Pablo Artal
Journal:  Biomed Opt Express       Date:  2022-04-21       Impact factor: 3.562

2.  Effect of cycloplegia on the measurement of refractive error in Chinese children.

Authors:  Tao Li; Xiaodong Zhou; Jie Zhu; Xiaojing Tang; Xiaoyan Gu
Journal:  Clin Exp Optom       Date:  2018-08-22       Impact factor: 2.742

3.  Agreement of wavefront-based refraction, dry and cycloplegic autorefraction with subjective refraction.

Authors:  Shahram Bamdad; Hamed Momeni-Moghaddam; Milad Abdolahian; David P Piñero
Journal:  J Optom       Date:  2020-09-28

Review 4.  Peripheral refraction and higher order aberrations.

Authors:  Dmitry Romashchenko; Robert Rosén; Linda Lundström
Journal:  Clin Exp Optom       Date:  2019-08-05       Impact factor: 2.742

5.  New Objective Refraction Metric Based on Sphere Fitting to the Wavefront.

Authors:  Mateusz Jaskulski; Andreí Martínez-Finkelshtein; Norberto López-Gil
Journal:  J Ophthalmol       Date:  2017-09-20       Impact factor: 1.909

  5 in total

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