Literature Search and Characteristics of Included Studies
A total of 238 studies were obtained from online databases according to the search strategy. After eliminating 117 duplicates, we reviewed abstracts of 121 studies, of which 93 studies were irrelevant. The final 28 studies were reviewed in full, and 16 studies were excluded because they did not meet the inclusion criteria. Finally. There are 12 studies included in this meta-analysis, with 9 RCTs and 3 cohort studies respectively. The screening process of eligible studies was shown in the flow diagram in Figure 1.
Table 1 shows the characteristics of 12 eligible studies. Five studies involved atropine monotherapy, including atropine at a low concentration of 0.01% and atropine at a medium concentration of 0.02% and 0.125%.[18][21] 8 articles involved atropine combined with ortho-k treatment, and the concentration of atropine included was 0.01%. [22][31] Zhao 2020 30 involved both atropine monotherapy and atropine combined with ortho-k therapy. Two studies (nozomi Kinoshita 2018 [24], nozomi Kinoshita 2020[27]) were conducted in Japan, and rest of studies were conducted in China. A total of 1218 myopic children were included in this study, with 671 in the atropine monotherapy group and 547 in AOK group. All children in the study had refractive error below -6.00D and were younger than 18 years old. 0.01% atropine was self-prepared in all studies. In the Yong Lyu study[21], the 0.02% atropine eye drops were made by diluting 1% atropine with saline. The 0.125% atropine used in Ren 2017[20]study was the lowest concentration of atropine available on the market.
Bias risk assessment
According to the Cochrane Collaboration's tool for assessing the risk of bias, , the quality of RCTs is high (Fig 2). Three studies [20][24][27]adopted randomization methods, but Only one study reported the allocation concealment[31]. Four Articles mentioned the blind method and explain applyment of the blind method in outcome evaluation[24][27][30][31]. Included studies completely report all outcome data In addition, other biases in six studies were low risk[24][26][27][29][31]. The Newcastle Ottawa scale was used to evaluate the cohort study, and results showed that moderate quality with four, seven and five scores in three cohort studies respectively( Table 2 ).
Atropine monotherapy group
Change in AL
Five studies[18][21][30] all reported and analyzed the changes in AL in the atropine monotherapy group and the ortho-k monotherapy group. There was statistical heterogeneity between the two groups (I2 =81%),which may be related to the different concentrations of atropine. Low concentration atropine and moderate concentration atropine were used as the first and second subgroups for analysis .Heterogeneity results showed that P < 0.0009, I2=79% > 50%. So random effect model was employed in this study.
Subgroup 1 (0.01% atropine) included three studies.[19][20][30]. Heterogeneity analysis showed that P = 0.83 , I2 = 0%. .The mean AL of atropine was 0.03 mm longer than that of ortho-k (z = 3.74, P = 0.0002), indicating that ortho-k control axial elongation more effectively than low concentration atropine monotherapy (Fig. 3).
Subgroup 2 ( 0.01% ~0.5% atropine) included two studies .[18][21] Heterogeneity analysis showed that P = 0.72, I2 = 0%. The mean AL of atropine was 0.09mm longer than that of ortho-k (z = 7.84, P < 0.00001), indicating that ortho-k control axial elongation more effectively than medium concentration atropine monotherapy(Fig. 3).
Atropine combined with ortho-k therapy group
Change in AL
There are 8 studies[24][27][29][31] on atropine combined with ortho-k and ortho-k monotherapy in slowing myopia progression. To investigate the effect of baseline diopter on axial elongation, we performed a subgroup analysis. We divided all included studies into two subgroups according to low myopia (- 1 ~ - 3D) and medium myopia (- 3 ~ - 6D).
There were 6 studies with baseline SER ranging from - 1 to - 3D ,146 patients in AOK group and 149 patients in ortho-k group.[24][27][29][31]1The combined results showed that the WMD in the low myopia subgroup was - 0.10 mm (95% CI, – 0.14, - 0.05), which represented that the mean axial elongation of AOK was 0.10 mm shorter than that of ortho-kwithin baseline diopter -3D.There were 3 studies[24][27][29] with baseline SER ranging from - 3 to - 6D,62 patients in AOK group and 62 patients in ortho-k group. The combined results showed that the WMD of the midium myopia subgroup was -0.04mm (95% CI, - 0.11, 0.03). There was no significant difference in the mean axial elongation between the AOK group and the ortho-k group (P = 0.30, I2 = 54%) (Fig. 4). However, due to the limited number of studies in midium myopia subgroup, extra data are needed to confirm the relationship between AOK and ortho-k in midum baseline myopia. It was suggested that the lower the baseline diopter, the better effect of AOK on the suppression of the axial elongation.
Publication Bias
A funnel plot demonstrated that all included studies have a certain degree of asymmetry, indicating that the research results may have potential publication bias. However, the results were as reference only because of the limited number of studies.(Fig. 5).
Sensitivity analysis
Subgroup 1 (low myopia “-1 to - 3D”): We performed sensitivity analysis (removal of a study one by one) to investigate the sources of heterogeneity. By omitting the Qi Tan 2019[26] study, heterogeneity was reduced from 88% to 44% and the overall effect was not inversed (Supplemental Table 3 ). The results are reliable. It may be related to that the follow-up time of this study is only one month, while other studies are more than one year.
Subgroup 2 (midium myopia“-3~ -6D”) showed that the heterogeneity decreased from 54% to 0% after excluding Tang 2020,[29]but the results remained unchanged. (Table 3 ). It suggested that Tang2020[29] took the axial length at 1 months after wearing OK as the baseline,Kinoshita2018 and 2020 [24][27]took the axial length at 3 months after wearing OK lensas the baseline. The central corneal thickness tends to stabilize after the first 1-2 months of OK therapy,[37][38] which may affect the experimental results.