Short-term efficacy of latanoprostene bunod for the treatment of open-angle glaucoma and ocular hypertension: a systematic literature review and a network meta-analysis.

BACKGROUND/AIMS: To assess the comparative efficacy of latanoprostene bunod (LBN), a novel prostaglandin analogue (PGA), to other medications for open-angle glaucoma and ocular hypertension on lowering intraocular pressure (IOP).
METHODS: A systematic literature review adapted from the Li et al (Ophthalmology, 2016) study was conducted. Medline, Embase and PubMed were searched for randomised controlled trials published between 1 January 2014 and 19 March 2020. Studies had to report IOP reduction after 3 months for at least two different treatments among placebo, PGAs (bimatoprost 0.01%, bimatoprost 0.03%, latanoprost, LBN, tafluprost, unoprostone) or apraclonidine, betaxolol, brimonidine, brinzolamide, carteolol, dorzolamide, levobunolol, timolol, travoprost. A Bayesian network meta-analysis was performed to provide the relative effect in terms of mean difference (95% credible interval) of IOP reduction and ranking probabilities. Surface under the cumulative ranking curve (SUCRA) was generated.
RESULTS: A total of 106 trials were included with data for 18 523 participants. LBN was significantly more effective than unoprostone (-3.45 (-4.77 to -2.12)). Although relative effect was not significative, compared with other PGAs, LBN numerically outperformed latanoprost (-0.70 (-1.83 to 0.43)) and tafluoprost (-0.41 (-1.87 to 1.07)), was similar to bimatoprost 0.01% (-0.02(-1.59 to 1.55)) and was slightly disadvantaged by bimatoprost 0.03% (-0.17 (-1.42 to 1.07)). LBN was significantly more efficient than the beta-blockers apraclonidine, betaxolol, brimonidine, brinzolamide, carteolol, dorzolamide and timolol. According to SUCRA, LBN was ranked second after bimatoprost 0.03%, followed by bimatoprost 0.01%.
CONCLUSION: LBN was significantly more effective than the PGA unoprostone and most of the beta-blockers. Compared with the most widely used PGAs, LBN numerically outperformed latanoprost and travoprost and was similar to bimatoprost 0.01%. © Author(s) (or their employer(s)) 2022. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.

Introduction

Glaucoma is a group of progressive optic neuropathies characterised by degeneration of retinal ganglion cells which may lead to vision loss and blindness.1 It is the number one cause of irreversible vision loss and the second leading cause of blindness worldwide.2 3 Primary open-angle glaucoma (POAG) is the most common form of the disease in North America with a prevalence of 3.3% (2.7 million people) in adults aged between 40 and 80 years in 2013.4

The goal of treatment is to reduce intraocular pressure (IOP), which is the only modifiable risk factor at this time.3 5 Initial treatment consists of topical therapies with several classes available, including prostaglandin analogues (PGAs), α-adrenergic agonist, beta-blockers and carbonic anhydrase inhibitors and parasympathomimetic agents.3 6 7 Among these, PGAs are the most effective medication because of their unmatched safety profile, IOP-lowering capabilities and their once-daily administration, with latanoprost, bimatoprost and travoprost being the most frequently used.8 Of note, in 2012, Lumigan (bimatoprost 0.03%) was discontinued and replaced by Lumigan RC (bimatoprost 0.01%) due to its favourable tolerability profile.9

In order to compare the different treatments a comprehensive assessment of their relative efficacy is crucial for clinicians and healthcare decision-makers,8 10 however, no head-to-head trials comparing all relevant competing therapies have been published. In the absence of direct evidence, the use of a network meta-analysis (NMA) may provide useful evidence.10 In 2016, Li et al published the results of a systematic review and an NMA which aimed to compare the effectiveness of first-line medications for patients with POAG or ocular hypertension (OH) and to provide relative ranking of these treatments. The authors conducted a systematic review in March 2014 in order to identify all randomised controlled trials (RCTs) comparing single active topical medication with no treatment/placebo or with another single topical treatment. Following a systematic review of 114 eligible trials, results of the NMA indicated that, compared with beta-blockers, α-adrenergic agonists and carbonic anhydrase inhibitors, PGAs were more efficacious in reducing IOP at 3 months. Authors also concluded that drugs within the PGA class, namely bimatoprost, latanoprost and travoprost were among the most efficacious, with intraclass difference found to be small and not clinically meaningful.11

Vyzulta (latanoprostene bunod (LBN) ophthalmic solution, 0.024% w/v), a novel nitric oxidedonating prostaglandin F2α analogue has received approbation for commercialisation in six different countries.12 The safety and efficacy of LBN has been well established through clinical studies (APOLLO and LUNAR studies),13 14 where LBN demonstrated enhanced efficacy compared with latanoprost and timolol.13–15 However, the effectiveness of LBN in comparison to other topical therapies other than latanoprost and timolol has not yet been evaluated. The objective of this study was to assess, through a systematic review and an NMA, the relative efficacy, as well as provide a relative ranking, of LBN compared with other topical medications, with a focus on PGAs, for the treatment of POAG and OH.

Material and methods

The pool of studies included in Li et al 11 previously described was considered and an exhaustive literature review was performed for studies published after 2013. An NMA was conducted according to a predefined protocol and was conformed to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses16 extension for NMA. The review question was established using the population, intervention, comparators, outcomes (PICO) framework. This systematic review was adapted from the work by Li et al previously described.11

Search strategy

MEDLINE, EMBASE and PubMed databases were searched on 19 March 2020 to identify RCTs published in English or French between 1 January 2014 and 19 March 2020. A manual search of reference lists was also performed to identify potentially relevant papers and systematic reviews. The same search strategy elaborated by Li et al was used but ‘latanoprostene bunod’ was added as a keyword.11 Detailed search strategies are presented in online supplemental appendix A.

Eligibility

Studies were selected if they reported relative efficacy between at least two different treatments (placebo, bimatoprost 0.01%, bimatoprost 0.03%, latanoprost, LBN, tafluprost, unoprostone, apraclonidine, betaxolol, brimonidine, brinzolamide, carteolol, dorzolamide, levobunolol, timolol or travoprost) in terms of IOP reduction after 3 months of usage. All eligibility criteria were defined a priori and were rigorously considered assuming the similarity assumption. Inclusion criteria included the following: RCTs with a parallel-group design (cross-over trials excluded); at least 60% of patients with a diagnostic of POAG and/or OH; trials that assess a monotherapy regimen (combination of medical treatments excluded); studies published in English and French between 1 January 2014 and 19 March 2020.Trials were excluded if they enrolled fewer than 10 participants in each group or if they evaluated a combination of medical treatments. Although no maximum or minimum duration of treatment was required, participants had to be followed for at least 28 days after randomisation.

Study selection and data extraction

Two reviewers independently screened the titles and abstracts of publications for potential eligibility. Using a predefined eligibility form (online supplemental appendix B), both reviewers screened the full text of all potentially eligible trials. Any disagreements were resolved by consensus or with the help of a third reviewer.

Data extraction was performed by two independent reviewers. Data extracted included: first author’s name, year of publication, trial design, location of trial, sample size, patients’ baseline characteristics, intervention characteristics and quantitative results with regard to treatment effect. For studies presenting multiple treatment durations, the duration closest to 3 months was used. If many IOP measures were available, the selection was made in this order: mean diurnal IOP, 24-hour mean IOP, peak IOP reduction and morning IOP. Any discrepancies were resolved by consensus or with the help of a third reviewer.

Quality assessment

As part of their systematic review, Li et al assessed the quality of included trials using the Cochrane Risk of Bias Tool, where the following seven methodological domains were graded as ‘low’, ‘high’ or ‘unclear’ risk of bias: sequence generation, allocation concealment, blinding of participants and personnel, blinding of outcome assessors, funding of the trial and financial relationship reported by the authors.11 17 Based on their work, new trials identified by the current systematic review were assessed using the same method.

Outcome definition

The primary outcome was defined as the mean reduction (MR) of IOP in continuous mmHg units after 3 months of treatment. The mean difference (MD) of the MR of IOP between two treatments with a 95% CI or credible interval (CrI) was calculated. An MD under 0 indicated that the treatment of reference performed a higher IOP reduction relative to its comparator and was therefore more effective.

Data synthesis and analysis

Using the ‘meta’ package in R, a pairwise meta-analysis (ie, direct comparisons) with a random-effect model was conducted for every treatment comparison with at least two trials. Statistical heterogeneity between studies was assessed using the I2 statistic, which describes the percentage of variability in effect estimates that is due to heterogeneity rather than sampling error (chance).18 Cochrane Handbook developed a rough guide for interpretation of I2: less than 40% might not be important, 30%–60% may represent moderate heterogeneity; 50%–90% may represent substantial heterogeneity and 75%–100% represented considerable heterogeneity.18 Pairwise comparisons with an I2 value greater than 65% were investigated to identify studies possibly causing heterogeneity.

An NMA, which combined direct and indirect comparisons, was conducted using a Bayesian random-effect model with Markov Chain Monte Carlo simulations executed with the ‘gemtc’ package in R.19 Using four parallel chains, 50 000 samples after 20 000-sample burn-in were obtained in each chain. Convergence of the model was assessed using the Brooks-Gelman-Rubin diagnostic in the ‘coda’ package in R. Consistency of the NMA, defined as a statistical discrepancy between direct and indirect comparison results, was evaluated using a node-splitting approach with the ‘gemtc’ package in R.19

The model ranked each treatment by their relative effect (probabilities of being more effective). Cumulative probability of being the most effective treatment was calculated. With that, the surface under the cumulative ranking curve (SUCRA) of each treatment is obtained.20 Specifically, SUCRA is a numeric presentation of the overall ranking and presents a single number, ranging from 0% to 100%, associated with each treatment, where 0% represents the least effective treatment and 100% represents the most effective treatment.21

Sensitivity analyses were conducted to assess the impact of heterogeneity between studies and inconsistency results by removing studies identified as possibly causing heterogeneity and including inconsistent combinations, respectively.

Supplementary analyses were conducted to evaluate the heterogeneity between baseline characteristics among trials included.

All concentrations of the same medication were combined in the same group except for bimatoprost 0.01% and bimatoprost 0.03%.

Results

Of the 2642 publications identified by the systematic review and the 114 studies used by Li et al, 106 RCTs met the a priori eligibility criteria and were included (figure 1; references of these RCTs are listed in online supplemental appendix C). Of these, 11 (10%) were published between 2014 and 2020. The total number of participants contributing to this network is 18 523 (complete characteristics of included studies are listed in online supplemental appendix D.

Figure 1

Organisational chart of the literature review. *Among the 128 excluded studies, seven were included in the Li et al publication. NMA, network meta-analysis; POAG, primary open-angle glaucoma.

Of the 106 trials, risk of selection bias (online supplemental appendix E) was rated as low for 54 (51%) and 33 (31%) studies when assessing sequence generation or allocation concealment, respectively, whereas the remaining trials were rated as having an ‘unclear risk’ except for one study with a ‘high risk’ in allocation concealment. Risk of performance bias, associated with blinding of participants, was rated as low (ie, reported blinding), high (ie, reported not blinding) or unclear (ie, not reported or unclear), for 42%, 37% and 21% of studies, respectively. Risk of detection bias, associated with blinding of the outcome assessor, was rated as low (ie, reported blinding), high (ie, reported not blinding) or unclear risk (ie, not reported) for 24%, 63% and 13% of studies, respectively. Of the 69 articles who reported funding for their research, 64 (93%) were funded by the industry. Of the 55 articles that reported financial relationship, 15 (27%) declared having no financial conflict of interest.

The 106 studies included compared 16 interventions (figure 2). A total of 138 direct comparisons were performed based on 93 two-arm trials, 11 three-arm trials and 2 four-arm trials. Results of the pairwise meta-analysis are presented in table 1. LBN was compared with timolol in two studies and latanoprost in one study. In both cases, LBN significantly lowers IOP more than the other treatments after 3 months (LBN vs timolol: MD (95% CI)=−1.42 (−1.84 to −1.01) and LBN vs latanoprost: −1.23 (−1.76 to −0.70)).

Figure 2

Network Graph. The nodes are weighted according to the number of participants randomised to that drug. The edges are weighted according to the number of direct comparison studies between drugs.

Table 1

Summary estimates for intraocular pressure at 3 months derived from the pairwise meta-analysis

Control	Experimental	Total no of studies	Meandifference*	95%CI		τ2†	I2‡
				Low	Up
Placebo	Bimatoprost 0.01%	1	−4.60	−5.60	−3.60	NA	NA
	Latanoprost	1	−3.10	−3.98	−2.22	NA	NA
	Unoprostone	1	−0.30	−1.50	0.90	NA	NA
	Betaxolol	2	−3.16	−4.17	−2.15	0.3	52%
	Brimonidine	1	−2.30	−3.99	−0.61	NA	NA
	Brinzolamide	1	−2.22	−3.48	−0.96	NA	NA
	Dorzolamide	4	−2.48	−3.84	−1.12	1.3	76%
	Levobunolol	2	−7.90	−8.94	−6.85	0.0	0%
	Timolol	4	−3.75	−4.75	−2.76	0.6	58%
Bimatoprost 0.01%	Latanoprost	2	1.02	0.68	1.37	0.0	0%
	Tafluprost	1	2.30	−0.91	5.51	NA	NA
	Travoprost	2	1.50	−1.98	4.97	5.2	80%
Bimatoprost 0.03%	Latanoprost	7	0.99	0.46	1.53	0.3	61%
	Travoprost	8	0.44	−0.52	1.40	1.4	86%
Latanoprost	Latanoprostene bunod	1	−1.23	−1.76	−0.70	NA	NA
	Tafluprost	3	−0.99	−1.92	−0.07	0.0	0%
	Unoprostone	6	2.90	2.16	3.63	0.3	37%
	Travoprost	7	−0.15	−1.30	1.00	1.9	87%
Apraclonidine	Timolol	2	−0.44	−3.91	3.03	5.6	89%
Betaxolol	Latanoprost	2	−1.84	−3.22	−0.47	0.0	0%
	Unoprostone	1	0.60	0.09	1.11	NA	NA
	Dorzolamide	2	−0.21	−0.82	0.40	0.0	0%
	Levobunolol	2	−4.65	−10.13	0.84	13.3	84%
	Timolol	6	−1.30	−2.46	−0.13	1.2	67%
Brimonidine	Latanoprost	5	−1.22	−2.13	−0.31	0.8	78%
	Betaxolol	1	2.00	0.90	3.10	NA	NA
	Brinzolamide	2	0.90	0.39	1.42	0.0	0%
	Timolol	4	0.42	0.04	0.81	0.0	0%
	Travoprost	1	−1.20	−3.77	1.37	NA	NA
Brinzolamide	Dorzolamide	2	−0.34	−0.84	0.16	0.0	0%
Carteolol	Levobunolol	1	−2.90	−4.59	−1.21	NA	NA
	Timolol	4	−0.27	−1.11	0.57	0.4	60%
Dorzolamide	Latanoprost	1	−2.90	−3.70	−2.10	0.0	NA
Levobunolol	Timolol	9	0.11	−0.40	0.62	0.1	15%
Timolol	Bimatoprost 0.03%	6	−2.06	−2.36	−1.75	0.0	0%
	Latanoprost	15	−1.18	−1.65	−0.70	0.6	76%
	Latanoprostene bunod	2	−1.42	−1.84	−1.01	0.0	0%
	Tafluprost	2	−0.50	−1.12	0.12	0.1	38%
	Unoprostone	2	0.94	−0.43	2.31	0.9	87%
	Brinzolamide	3	1.10	0.52	1.69	0.0	0%
	Dorzolamide	4	0.99	0.34	1.64	0.1	26%
	Travoprost	4	−0.89	−1.26	−0.52	0.0	0%
Travoprost	Tafluprost	1	−1.30	−2.93	0.33	NA	NA
Total	16 drugs	138§

, PGA

*Difference between the reduction in IOP during the study of the experimental drug and the control drug (mean difference under 0 favours the experimental drug). Results presented in bold are significant.

†τ2 describes the underlying between-study variability.

‡I2 is the percentage of variability in the treatment estimates which is attributable to heterogeneity.

§106 trials considered: 93 two-arm trials, 11 three-arm trials and 2 four-arm trials.

IOP, intraocular pressure; NA, not available; PGA, prostaglandin analogue.

Results of the NMA indicate that, when compared with placebo, all active drugs demonstrate an improved reduction of IOP at 3 months (table 2 and online supplemental appendix F). More specifically, the MDs in IOP reduction at 3 months for active drug in comparison to placebo range from −1.97 mm Hg for unoprostone to −5.59 mm Hg for bimatoprost 0.03% and are all statistically significative. Importantly, LBN shows the second greatest reduction in IOP vs placebo with an MD (95% CrI) of −5.42 mm Hg (−6.68 to –4.16). Furthermore, these results highlight the statistically significant superiority in efficacy of LBN compared with the PGA unoprostone (−3.45 (−4.77 to −2.12)) and the beta-blockers apraclonidine (−2.55 (−4.52 to −0.55)), betaxolol (−2.89 (−4.17 to −1.60)), brimonidine (−1.75 (−3.02 to −0.49)), brinzolamide (−2.88 (−-4.29 to −1.47)), carteolol (−2.17 (−3.65 to −0.69)), dorzolamide (−2.87 (−4.17 to −1.55)) and timolol (−1.69 (−2.80 to −0.58)). Although the relative effect was not significative, compared with other PGAs, LBN numerically outperformed latanoprost (−0.70 (−1.83 to 0.43)) and tafluoprost (−0.41 (−1.87 to 1.07)), was similar to bimatoprost 0.01% (−0.02 (−1.59 to 1.55)) and bimatoprost 0.03% demonstrated a slightly advantage over LBN (−0.17 (−1.42 to 1.07)) (table 2 and online supplemental appendix F). The model ranked each treatment by their relative effect (probabilities of being more effective) (table 3). According to these results, treatment with the higher probability of being ranked first is bimatoprost 0.03% with a probability of 37%, followed by LBN with a probability of 29%. LBN has a probability of 51% to be under the two best treatments and 70% to be under the three best treatments. Cumulative probability of being the most effective treatment was calculated and the cumulative ranking curve of each treatment (presented in online supplemental appendix G) was obtained to calculate the SUCRA. According to SUCRA results, LBN (SUCRA=88%) emerges as the second best treatment after bimatoprost 0.03% (94%) and followed in order by bimatoprost 0.01% (87%), tafluprost (78%), travoprost (73%), levobunolol (72%), latanoprost (68%), timolol (48%), brimonidine (47%), carteolol (38%), apraclonidine (30%), dorzolamide (23%), brinzolamide (22%), betaxolol (22%), unoprostone (11%) and placebo (0%).

Table 2

Mean difference in intraocular pressure at 3 months (95% credible interval) derived from the NMA

Placebo	−5.39	−5.59	−4.72	−5.42	−5.00	−1.97	−2.86	−2.53	−3.66	−2.53	−3.24	−2.55	−4.79	−3.73	−4.84
	(−6.60 to −4.21)	(−6.43 to−4.75)	(−5.39 to−4.04)	(−6.68 to−4.16)	(−6.15 to−3.88)	(−2.88 to−1.06)	(−4.65 to−1.12)	(−3.32 to−1.75)	(−4.48 to−2.85)	(−3.52 to−1.56)	(−4.40 to−2.09)	(−3.33 to−1.78)	(−5.64 to−3.95)	(−4.36 to−3.10)	(−5.69 to−3.98)
5.39	Bimatoprost 0.01%	−0.20	0.67	−0.02	0.39	3.42	2.53	2.86	1.73	2.85	2.15	2.84	0.60	1.66	0.55
(4.21 to 6.60)		(−1.41 to1.03)	(−0.44 to1.80)	(−1.59 to1.55)	(−1.06 to1.82)	(2.11 to 4.75)	(0.51 to 4.51)	(1.59 to 4.15)	(0.47 to 2.99)	(1.45 to 4.27)	(0.65 to 3.65)	(1.55 to 4.13)	(−0.70 to1.92)	(0.52 to 2.81)	(−0.60 to1.72)
5.59	0.20	Bimatoprost 0.03%	0.88	0.17	0.59	3.62	2.73	3.06	1.93	3.06	2.35	3.04	0.80	1.86	0.75
(4.75 to 6.43)	(−1.03 to1.41)		(0.26 to 1.48)	(−1.07 to1.42)	(−0.54 to1.69)	(2.69 to 4.55)	(0.94 to 4.47)	(2.18 to 3.95)	(1.09 to 2.76)	(2.01 to 4.10)	(1.18 to 3.50)	(2.12 to 3.94)	(−0.12 to1.71)	(1.24 to 2.48)	(0.08 to 1.43)
4.72	−0.67	−0.88	Latanoprost	−0.70	−0.29	2.75	1.85	2.18	1.05	2.18	1.47	2.17	−0.08	0.99	−0.12
(4.04 to 5.39)	(−1.80 to0.44)	(−1.48 to−0.26)		(−1.83 to0.43)	(−1.28 to0.69)	(2.01 to 3.49)	(0.13 to 3.53)	(1.46 to 2.92)	(0.41 to 1.69)	(1.27 to 3.09)	(0.41 to 2.53)	(1.41 to 2.92)	(−0.85 to0.70)	(0.58 to 1.39)	(−0.75 to0.52)
5.42	0.02	−0.17	0.70	Latanoprostene Bunod	0.41	3.45	2.55	2.89	1.75	2.88	2.17	2.87	0.62	1.69	0.58
(4.16 to 6.68)	(−1.55 to1.59)	(−1.42 to1.07)	(−0.43 to1.83)		(−1.07 to1.87)	(2.12 to 4.77)	(0.55 to 4.52)	(1.60 to 4.17)	(0.49 to 3.02)	(1.47 to 4.29)	(0.69 to 3.65)	(1.55 to 4.17)	(−0.68 to1.92)	(0.58 to 2.80)	(−0.69 to1.85)
5.00	−0.39	−0.59	0.29	−0.41	Tafluprost	3.04	2.14	2.47	1.34	2.47	1.76	2.46	0.21	1.28	0.17
(3.88 to 6.15)	(−1.82 to1.06)	(−1.69 to0.54)	(−0.69 to1.28)	(−1.87 to1.07)		(1.83 to 4.25)	(0.20 to 4.05)	(1.31 to 3.66)	(0.21 to 2.49)	(1.18 to 3.77)	(0.38 to 3.15)	(1.27 to 3.65)	(−0.97 to1.41)	(0.30 to 2.27)	(−0.92 to1.28)
1.97	−3.42	−3.62	−2.75	−3.45	−3.04	Unoprostone	−0.90	−0.56	−1.70	−0.57	−1.28	−0.58	−2.82	−1.76	−2.87
(1.06 to 2.88)	(−4.75 to−2.11)	(−4.55 to−2.69)	(−3.49 to−2.01)	(−4.77 to−2.12)	(−4.25 to−1.83)		(−2.73 to0.90)	(−1.50 to0.38)	(−2.63 to−0.75)	(−1.69 to0.56)	(−2.52 to−0.03)	(−1.57 to0.41)	(−3.83 to−1.81)	(−2.53 to−0.99)	(−3.82 to−1.92)
2.86	−2.53	−2.73	−1.85	−2.55	−2.14	0.90	Apraclonidine	0.33	−0.80	0.33	−0.38	0.31	−1.93	−0.86	−1.98
(1.12 to 4.65)	(−4.51 to−0.51)	(−4.47 to−0.94)	(−3.53 to−0.13)	(−4.52 to−0.55)	(−4.05 to−0.20)	(−0.90 to2.73)		(−1.43 to2.14)	(−2.54 to0.99)	(−1.51 to2.21)	(−2.28 to1.54)	(−1.46 to2.13)	(−3.70 to−0.13)	(−2.49 to0.81)	(−3.73 to−0.18)
2.53	−2.86	−3.06	−2.18	−2.89	−2.47	0.56	−0.33	Betaxolol	−1.13	0.00	−0.71	−0.02	−2.26	−1.20	−2.31
(1.75 to 3.32)	(−4.15 to−1.59)	(−3.95 to−2.18)	(−2.92 to−1.46)	(−4.17 to−1.60)	(−3.66 to−1.31)	(−0.38 to1.50)	(−2.14 to1.43)		(−1.99 to−0.29)	(−1.06 to1.03)	(−1.90 to0.46)	(−0.87 to0.81)	(−3.17 to−1.37)	(−1.88 to−0.52)	(−3.21 to−1.40)
3.66	−1.73	−1.93	−1.05	−1.75	−1.34	1.70	0.80	1.13	Brimonidine	1.13	0.42	1.11	−1.13	−0.06	−1.17
(2.85 to 4.48)	(−2.99 to−0.47)	(−2.76 to−1.09)	(−1.69 to−0.41)	(−3.02 to−0.49)	(−2.49 to−0.21)	(0.75 to 2.63)	(−0.99 to2.54)	(0.29 to 1.99)		(0.21 to 2.04)	(−0.75 to1.59)	(0.23 to 1.99)	(−2.05 to−0.21)	(−0.70 to0.57)	(−2.02 to−0.32)
2.53	−2.85	−3.06	−2.18	−2.88	−2.47	0.57	−0.33	0.00	−1.13	Brinzolamide	−0.71	−0.01	−2.26	−1.19	−2.30
(1.56 to 3.52)	(−4.27 to−1.45)	(−4.10 to−2.01)	(−3.09 to−1.27)	(−4.29 to−1.47)	(−3.77 to−1.18)	(−0.56 to1.69)	(−2.21 to1.51)	(−1.03 to1.06)	(−2.04 to−0.21)		(−2.01 to0.60)	(−0.99 to0.95)	(−3.34 to−1.17)	(−2.06 to−0.32)	(−3.36 to−1.24)
3.24	−2.15	−2.35	−1.47	−2.17	−1.76	1.28	0.38	0.71	−0.42	0.71	Carteolol	0.69	−1.55	−0.48	−1.59
(2.09 to 4.40)	(−3.65 to−0.65)	(−3.50 to−1.18)	(−2.53 to−0.41)	(−3.65 to−0.69)	(−3.15 to−0.38)	(0.03 to 2.52)	(−1.54 to2.28)	(−0.46 to1.90)	(−1.59 to0.75)	(−0.60 to2.01)		(−0.51 to1.89)	(−2.68 to−0.42)	(−1.46 to0.50)	(−2.76 to−0.42)
2.55	−2.84	−3.04	−2.17	−2.87	−2.46	0.58	−0.31	0.02	−1.11	0.01	−0.69	Dorzolamide	−2.24	−1.18	−2.29
(1.78 to 3.33)	(−4.13 to−1.55)	(−3.94 to−2.12)	(−2.92 to−1.41)	(−4.17 to−1.55)	(−3.65 to−1.27)	(−0.41 to1.57)	(−2.13 to1.46)	(−0.81 to0.87)	(−1.99 to−0.23)	(−0.95 to0.99)	(−1.89 to0.51)		(−3.18 to−1.29)	(−1.88 to−0.47)	(−3.21 to−1.35)
4.79	−0.60	−0.80	0.08	−0.62	−0.21	2.82	1.93	2.26	1.13	2.26	1.55	2.24	Levobunolol	1.06	−0.05
(3.95 to 5.64)	(−1.92 to0.70)	(−1.71 to0.12)	(−0.70 to0.85)	(−1.92 to0.68)	(−1.41 to0.97)	(1.81 to 3.83)	(0.13 to 3.70)	(1.37 to 3.17)	(0.21 to 2.05)	(1.17 to 3.34)	(0.42 to 2.68)	(1.29 to 3.18)		(0.38 to 1.75)	(−0.98 to0.89)
3.73	−1.66	−1.86	−0.99	−1.69	−1.28	1.76	0.86	1.20	0.06	1.19	0.48	1.18	−1.06	Timolol	−1.11
(3.10 to 4.36)	(−2.81 to−0.52)	(−2.48 to−1.24)	(−1.39 to−0.58)	(−2.80 to−0.58)	(−2.27 to−0.30)	(0.99 to 2.53)	(−0.81 to2.49)	(0.52 to 1.88)	(−0.57 to0.70)	(0.32 to 2.06)	(−0.50 to1.46)	(0.47 to 1.88)	(−1.75 to−0.38)		(−1.76 to−0.46)
4.84	−0.55	−0.75	0.12	−0.58	−0.17	2.87	1.98	2.31	1.17	2.30	1.59	2.29	0.05	1.11	Travoprost
(3.98 to 5.69)	(−1.72 to0.60)	(−1.43 to−0.08)	(−0.52 to0.75)	(−1.85 to0.69)	(−1.28 to0.92)	(1.92 to 3.82)	(0.18 to 3.73)	(1.40 to 3.21)	(0.32 to 2.02)	(1.24 to 3.36)	(0.42 to 2.76)	(1.35 to 3.21)	(−0.89 to0.98)	(0.46 to1.76)

, PGA

, others

Results presented in bold are significant.

Mean difference under 0 favours the drug in the column.

NAM, network meta-analysis; PGA, prostaglandin analogue.

Table 3

Ranking probabilities and SUCRA

Ranks	A	B	C	D	E	F	G	H	I	J	K	L	M	N	O	P
	Placebo	Bimatoprost 0.01%	Bimatoprost 0.03%	Latanoprost	Latanoprostene bunod	Tafluprost	Unoprostone	Apraclonidine	Betaxolol	Brimonidine	Brinzolamide	Carteolol	Dorzolamide	Levobunolol	Timolol	Travoprost
1	0.000	0.266	0.367	0.000	0.288	0.067	0.000	0.001	0.000	0.000	0.000	0.000	0.000	0.009	0.000	0.002
2	0.000	0.221	0.374	0.002	0.220	0.117	0.000	0.001	0.000	0.000	0.000	0.000	0.000	0.039	0.000	0.025
3	0.000	0.204	0.190	0.023	0.193	0.180	0.000	0.002	0.000	0.000	0.000	0.000	0.000	0.104	0.000	0.104
4	0.000	0.126	0.055	0.103	0.118	0.198	0.000	0.004	0.000	0.000	0.000	0.000	0.000	0.176	0.000	0.220
5	0.000	0.072	0.011	0.251	0.071	0.144	0.000	0.005	0.000	0.000	0.000	0.001	0.000	0.195	0.000	0.250
6	0.000	0.052	0.002	0.358	0.050	0.120	0.000	0.007	0.000	0.002	0.000	0.003	0.000	0.193	0.000	0.214
7	0.000	0.052	0.000	0.256	0.052	0.153	0.000	0.016	0.000	0.020	0.000	0.010	0.000	0.263	0.005	0.174
8	0.000	0.006	0.000	0.008	0.006	0.016	0.000	0.089	0.000	0.328	0.001	0.107	0.000	0.018	0.409	0.011
9	0.000	0.001	0.000	0.000	0.001	0.004	0.000	0.063	0.001	0.353	0.005	0.120	0.002	0.001	0.446	0.001
10	0.000	0.001	0.000	0.000	0.001	0.001	0.003	0.150	0.036	0.230	0.058	0.352	0.041	0.000	0.128	0.000
11	0.000	0.000	0.000	0.000	0.000	0.000	0.015	0.201	0.152	0.062	0.170	0.232	0.157	0.000	0.011	0.000
12	0.000	0.000	0.000	0.000	0.000	0.000	0.041	0.107	0.259	0.004	0.231	0.088	0.271	0.000	0.000	0.000
13	0.000	0.000	0.000	0.000	0.000	0.000	0.087	0.094	0.271	0.001	0.220	0.049	0.278	0.000	0.000	0.000
14	0.000	0.000	0.000	0.000	0.000	0.000	0.220	0.123	0.223	0.000	0.214	0.029	0.192	0.000	0.000	0.000
15	0.001	0.000	0.000	0.000	0.000	0.000	0.635	0.136	0.059	0.000	0.102	0.009	0.058	0.000	0.000	0.000
16	0.999	0.000	0.000	0.000	0.000	0.000	0.000	0.001	0.000	0.000	0.000	0.000	0.000	0.000	0.000	0.000
SUCRA*	0.0%	87.2%	93.5%	68.4%	87.6%	77.9%	10.6%	30.1%	22.2%	46.7%	22.3%	37.8%	22.7%	71.8%	48.5%	72.7%
Ranking according SUCRA*	16	3	1	7	2	4	15	11	14	9	13	10	12	6	8	5

, PGA

*Higher is the SUCRA, better is the rank of this treatment.

SUCRA, surface under the cumulative ranking curve.

Sensitivity analyses

A total of 10 direct comparisons were identified as possibly causing heterogeneity (online supplemental appendix H). When excluding these studies, the sensitivity analysis revealed no significant change in the NMA results (online supplemental appendix I). LBN was still significantly better than unoprostone and non-PGAs treatments, except for levobunolol and travoprost that was numerically superior. Although the relative effect was not significative, compared with other PGAs, LBN numerically outperformed latanoprost (−0.72 (−1.60 to 0.16)), tafluoprost (−0.60 (−1.80 to 0.61)) and bimatoprost 0.01% (-0.40 (−1.70 to 0.83)) and bimatoprost 0.03% demonstrated a slight advantage over LBN (0.13 (−0.88 to 1.10)).

The node-splitting approach allowed for the identification of two inconsistent nodes (levobunolol vs placebo and timolol vs levobunolol) (online supplemental appendix J). When excluding these nodes, the sensitivity analysis revealed no significant change in the NMA results. Compared with unoprostone (PGA) and other non-PGAs, results indicated that LBN was significantly better, excluding travoprost but including levobunolol (online supplemental appendix I). Compared with other PGAs, although the relative effect was not significative, LBN was still numerically superior to latanoprost (−0.66 (−1.60 to 0.31)), similar to bimatoprost 0.01% (0.09 (−1.30 to 1.50)) and disadvantaged by bimatoprost 0.03% (0.20 (−0.87 to 1.30)).

Four supplementary analyses were also conducted to evaluate the heterogeneity between baseline characteristics among trials by considering: (1) only studies published from 2000 onward, (2) studies with a washout period before randomisation, (3) studies that excluded prior glaucoma and cataract surgery, and (4) studies that excluded prior glaucoma laser. These analyses revealed that heterogeneity between baseline characteristics had no significant impact on the NMA results (online supplemental appendix K).

Also, the Brooks-Gelman-Rubin plot (online supplemental appendix L) illustrates that the NMA model converges.

Discussion

The objective of this study was to assess the relative efficacy of a new IOP-lowering medication, LBN, compared with other topical medications for the treatment of POAG and OH and to provide a relative ranking of these treatments. Findings from this NMA confirm that all drugs are more effective when compared with the placebo. Importantly, results also indicate that LBN is significantly more effective than unoprostone (PGA) and other non-PGAs drugs except levobunolol and travoprost for which LBN is numerically better although not significant. This demonstrates that LBN is more effective than timolol, which aligns with the conclusion drawn from the individual studies (APOLLO and LUNAR).13 14 Moreover, compared with other PGAs, LBN was numerically more effective than tafluoprost, similar to bimatoprost 0.01% and slightly disadvantaged by bimatoprost 0.03%.

This systematic review was adapted from the one conducted by Li et al that was previously published in a peer-reviewed journal.11 A clear research question was formed using the PICO framework and the analysis was conducted based on the predefined protocol. It should be noted that of the 114 trials eligible in the NMA published by Li et al, 19 were not included in this NMA. This is explained by the fact that our systematic review was limited to English or French publications, whereas Li et al did not impose any language restriction. Moreover, some full-text articles were not accessible via the databases exploited for this study. Nonetheless, results presented herein are consistent with the findings of Li et al. Indeed, when comparing PGAs in terms of IOP reduction at 3 months, the intraclass differences are relatively small and not significantly meaningful. In addition, this systematic review and NMA, which include the most recent PGA, namely LBN, provides new findings relevant to clinicians and decision-makers as it allows for the comparison of drugs that had not yet been evaluated in head-to-head trials.

It should be noted that there are some limitations associated with this NMA. First, although an NMA represents a powerful tool and may provide crucial information, an inherent limitation associated with NMA resides in the variability and the risk of biases of studies included. Due to possible variability between studies and between the comparisons made, a critical step when performing an NMA consists of validating the homogeneity and consistency assumptions. The sensitivity analyses conducted did not significantly alter the results, suggesting that the assumptions and conclusions made based on the statistical analysis are reliable and robust.

Second, this NMA focused on IOP reduction and did not include visual field outcomes. We acknowledge that the ultimate goal in the management of glaucoma consists of slowing or stopping structural damages leading to vision loss and that, consequently, visual field outcomes would be more clinically meaningful than IOP when comparing treatment response. However, due to the lag time between onset of optic neuropathy and clinically detectable visual field defects, the use of visual field outcomes to assess relative effectiveness of different interventions requires an extended time frame which poses challenges to the conduct of RCTs. Thus, although IOP does not measure structural of functional glaucomatous optic neuropathy, it remains the most commonly used surrogate endpoint of RCTs.10 20 22 Li et al reported that only 11% of trials included in their NMA reported any analysable visual field data. Moreover, the authors mentioned that since visual field data were reported in many different ways, the conduct of a pairwise meta-analysis or NMA would have been impossible.11 Finally, although our study provided a relative ranking of topical treatments for glaucoma based on IOP reduction at 3 months, the choice of treatment remains a multifactorial decision to take into consideration different factors, such as patient’s medical history and preference, risk factors and likelihood of compliance.

Conclusion

Results from the NMA showed that, LBN was significantly more effective than the PGA unoprostone and most of the beta-blockers. Although there was no significant relative effect, compared with the most widely used PGAs, LBN was numerically more efficient than latanoprost and tafluoprost, was similar to bimatoprost 0.01% and was slightly disadvantaged by bimatoprost 0.03%. LBN could potentially become a promising option for glaucoma patients.