Potential bias in observational PCSK9 inhibitor study
Joseph Edgar Blais, Angel Y S Wong

Abstract
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsLipoproteins and Cardiovascular Health · Computational Drug Discovery Methods
We read with interest the study by Huang et al.^1^ published in the European Heart Journal—Cardiovascular Pharmacotherapy.^1^ Their study attempts to provide insights into the effectiveness of proprotein convertase subtilisin/kexin type 9 (PCSK9) inhibitors against all-cause mortality and hospitalization-related outcomes. It is also an example of how researchers’ decisions about the design of their pharmacoepidemiological study can result in biased estimates. Several different biases can occur in both the design and analysis of observational studies that use real-world data. These biases often occur because of the growing availability of longitudinal electronic databases which permit investigators to look both forwards and backwards in calendar time.^2,3^ Looking into the future and making decisions about an individual’s eligibility to be included in a study is simply not possible in a randomized controlled trial or cohort study with primary data collection. Analysis decisions, such as identifying treatment groups one before the other or excluding people from the analysis on the basis of the number of visits, their diagnoses, or their treatment during follow-up can cause selection bias.^4,5^
In the study conducted by Huang et al.^1^ the authors defined the index date (time zero) as initiation of PCSK9 inhibitors and statins respectively. In the PCSK9 inhibitor treatment strategy, individuals were allowed to be prescribed a statin before the index date. Therefore, they were not incident (new) users of lipid lowering drugs, which means some individuals receiving the PCSK9 inhibitor treatment strategy would have a survival advantage over the statin treatment strategy. Moreover, as switchers were excluded from the analysis (i.e. for the PCSK9 inhibitor treatment strategy, individuals prescribed a statin ‘on or after the index date’) were excluded instead of censoring the follow-up time at treatment switching. This excludes the person-time for individuals in the statin treatment strategy that were subsequently prescribed a PCSK9 inhibitor and the person-time of individuals in the PCSK9 inhibitor treatment strategy that were subsequently prescribed a statin. In a randomized clinical trial, time zero is the date when participant eligibility is assessed, treatment strategies are assigned, and follow-up starts.^3^ As described by Hernán and colleagues, misalignment of these three decisions frequently occurs in observational studies and produces selection bias.^3^ The study design selected by Huang et al.^1^ suggests misalignment of time zero (specifically target trial emulation failure example #4 described by Hernán) and is an example of classical immortal time bias that could bias the results in favour of PCSK9 inhibitors.^5^
Mohyuddin and Prasad also describe how bias can be identified by inspecting the survival plots to understand how quickly an intervention exerts an effect on the outcome under study.^6^ Inspecting Figure 2 in Huang et al.^1^ shows a clear separation of the cumulative incidence of all-cause mortality immediately after 30 days of follow-up. The speed of this effect against all-cause mortality seems biologically improbable given that the landmark ODYSSEY OUTCOMES and FOURIER studies showed that the survival curves for all-cause mortality began to separate after 12 months of treatment.^7,8^ Lastly, defining the start of follow-up at 30 days after the index date causes selection bias in the hazard ratio because of depletion of susceptible individuals.^9^ The susceptibility of individuals to all-cause mortality does not remain constant over time as shown in Supplementary material online, Table S9 of the original paper. A higher risk of mortality might be expected, for example, within one month after a patient has a myocardial infarction who then goes on to initiate a statin at hospital discharge. This depletion of susceptibles may also contribute to the rapid separation of the survival curves at the start of follow-up.
How can these flaws in real-world data analysis be avoided? Clinical investigators who intend to use complex real-world data need to include team members with expertise in observational epidemiology and who are ideally experienced with the health system and data source used in the study. Using the target trial framework can help avoid common errors in study design.^3^ Furthermore, the questions listed in Lévesque et al.^5^ can also be used by researchers to avoid inducing immortal time bias.^5^ Study design diagrams and other visualisations have been promoted to explicitly describe the order of application of time windows for each temporal anchor and patient flow in real-world studies.^10^
An increase in the accessibility of electronic health records such as the TriNetX database should come with a commensurate increase in the responsibilities of authors, journal editors, and peer reviewers to ensure the application of appropriate methodologies to mitigate bias in real-world studies of cardiovascular pharmacotherapy. Researchers in pharmacoepidemiology are in danger of losing trust and credibility if we fail to appropriately respond to these methodological challenges when analyzing real-world data.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Huang CH, Wang SI, Fan FS, Lu HJ, Wei JC. Association of PCSK 9 inhibitors with mortality: insights from a retrospective cohort analysis. Eur Heart J Cardiovasc Pharmacother 2024;10:505–514.39054050 10.1093/ehjcvp/pvae 056 · doi ↗ · pubmed ↗
- 2Suissa S, Dell'Aniello S. Time-related biases in pharmacoepidemiology. Pharmacoepidemiol Drug Saf 2020;29:1101–1110.32783283 10.1002/pds.5083 · doi ↗ · pubmed ↗
- 3Hernán MA, Sauer BC, Hernández-Díaz S, Platt R, Shrier I. Specifying a target trial prevents immortal time bias and other self-inflicted injuries in observational analyses. J Clin Epidemiol 2016;79:70–75.27237061 10.1016/j.jclinepi.2016.04.014PMC 5124536 · doi ↗ · pubmed ↗
- 4Lund JL, Horváth-Puhó E, Komjáthiné Szépligeti S, Sørensen HT, Pedersen L, Ehrenstein V, Stürmer T. Conditioning on future exposure to define study cohorts can induce bias: the case of low-dose acetylsalicylic acid and risk of major bleeding. Clin Epidemiol 2017;9:611–626.29200891 10.2147/CLEP.S 147175 PMC 5703173 · doi ↗ · pubmed ↗
- 5Lévesque LE, Hanley JA, Kezouh A, Suissa S. Problem of immortal time bias in cohort studies: example using statins for preventing progression of diabetes. BMJ 2010;340:b 5087.20228141 10.1136/bmj.b 5087 · doi ↗ · pubmed ↗
- 6Mohyuddin GR, Prasad V. Detecting selection bias in observational studies-when interventions work too fast. JAMA Intern Med 2023;183:897–898.37306983 10.1001/jamainternmed.2023.2067 · doi ↗ · pubmed ↗
- 7Schwartz GG, Steg PG, Szarek M, Bhatt DL, Bittner VA, Diaz R, Edelberg JM, Goodman SG, Hanotin C, Harrington RA, Jukema JW, Lecorps G, Mahaffey KW, Moryusef A, Pordy R, Quintero K, Roe MT, Sasiela WJ, Tamby JF, Tricoci P, White HD, Zeiher AM; ODYSSEY OUTCOMES Committees and Investigators. Alirocumab and cardiovascular outcomes after acute coronary syndrome. N Engl J Med 2018;379:2097–2107.30403574 10.1056/NEJ Moa 1801174 · doi ↗ · pubmed ↗
- 8van Bruggen FH, Luijendijk HJ. Evolocumab’s long-term mortality risk unclear due to shortened follow-up of FOURIER. Am J Cardiovasc Drugs 2022;22:5–8.34008143 10.1007/s 40256-021-00480-y PMC 8748321 · doi ↗ · pubmed ↗
